%小波包分解后RBF預(yù)測
%clc
clear
close?all

%---------------------------------------------------
%?產(chǎn)生訓(xùn)練樣本與測試樣本
clear?all;
clc;
n=4;
wpname=‘db8‘;

m=34;
b1=[];
a1=zeros（116）;
E1=zeros（116）;
for?i1=1:m
%?????mid=fopen（‘E:\研究生\機(jī)車軸承數(shù)據(jù)\株洲滾動軸承故障數(shù)據(jù)\株洲測試200509二進(jìn)制全部\050927-174712_12.BIN‘‘r‘）;
%?????s=fread（mid‘single‘）;?
%?fs=200;
%?sj?=?load（‘G:\數(shù)據(jù)\ASCE數(shù)據(jù)整理\二期數(shù)據(jù)\data\Shaker\Random\shm03s‘）;
%?data1=sj.dasy.DA02（1:3000:）‘;
%?t=（1:length（data1））/（fs/2）;
%?figure（1）
%?plot（tdata1‘b‘）
%?xlabel（‘時間?t/s‘）;
%?ylabel（‘幅值?A‘）;
%?%?title（‘原始數(shù)據(jù)‘）;
%?f=calculatef（data1fs）;
%?%?t=（1:length（data1））/fs;
%?sf=f（1:）;
%?data=sf;
%?figure（2）
%?plot（data‘b‘）
%?xlabel（‘?dāng)?shù)據(jù)點(diǎn)‘）;
%?ylabel（‘頻率?f/Hz‘）;
%?????fclose（mid）;
????year=‘2004‘;
????month=‘03‘;
????day=‘04‘;
????hour=‘09‘;
????minite=‘27‘;
????second=‘46‘;
????pathx=‘D:\論文\論文參考資料\nasa\bearing_IMS\3rd_test\4th_test\txt\‘;
????filename0=[year‘.‘month‘.‘day‘.‘hour‘.‘minite‘.‘second];
????filename=[pathxfilename0];
????x0=load（filename）;
????w1=load（filename）;??
????y11=w1（::）;
%?????w1=load（‘C:\Users\Administrator\Desktop\幾種算法\ccec_data\12k_Fan_End_Bearing_Fault_data\B007_0.mat‘）;
%?????y11=w1.X118_DE_time;
%?????yy=fft1（y11）;
%?????plot（yy）
%?????y21=w1.X118_FE_time;
%?????y31=w1.X118_BA_time;
????y11=y11（（600*i1-599）:（600*i1））;
????Data1=y11-mean（y11）;
????Data1=detrend（Data1）;
%?????T?=?wpdec（Data4‘sym4‘）;
%?????E?=?wenergy（T）;
????wpt1=wpdec（Data1nwpname）;?%對數(shù)據(jù)進(jìn)行小波包分解
????for?j1=1:2^n?%wpcoef（wpt1[ni-1]）是求第n層第i個節(jié)點(diǎn)的系數(shù)
????E1（j1）=norm（wpcoef（wpt1[nj1-1]）2）;%求第i個節(jié)點(diǎn)的范數(shù)平方，其實(shí)也就是平方和
????end
????%?disp（‘每個節(jié)點(diǎn)的能量E（i）‘）;
????%?E
????%?disp（‘小波包分解總能量E_total‘）;
????E_total1=sum（E1）;?%求總能量
????for?k1=1:2^n
????a1（k1）=?E1（k1）/E_total1;%求每個節(jié)點(diǎn)的概率
????end
%?????a=[1a];
????b1=[b1;a1];??
end
%?[x1y1]=size（b1）;
save?b1?b1
B=reshape（b1.‘1544）;
%?B=reshape（b1.‘1?15000）;
input1=[];
output1=[];
for