資源簡介
程序介紹小波分解的算法思想和小波能量譜的計算,在信號處理領域應用較廣

代碼片段和文件信息
%采樣頻率
fs=5000;
nfft=10240;
%信號
fid=fopen(‘a(chǎn).txt‘‘r‘);%故障
N=512;
xdata=fscanf(fid‘%f‘N);
fclose(fid);
xdata=(xdata-mean(xdata))/std(xdata1);
%功率譜
figure(1);
subplot(221);
plot(xdata);
xlabel(‘時間?t/s‘);
ylabel(‘幅值?A/W‘);
Y=abs(fft(xdatanfft));
subplot(222);
plot((0:nfft/2-1)/nfft*fsY(1:nfft/2));
xlabel(‘頻率?f/Hz‘);
ylabel(‘功率譜?P/W‘);
%--------------------------------------------------------------------------
%5層小波包分解;
T=wpdec(xdata3‘db1‘);
plot(T);%最后一句為顯示小波包四層分解樹結(jié)構(gòu)
%-----------------------------------------------------------------------
figure(2);
subplot(211);
%重構(gòu)第一層各系數(shù)
y10=wprcoef(T[10]);
y11=wprcoef(T[11]);
figure(2);
subplot(221);
plot(y10);
xlabel(‘第一層低頻系數(shù)時間?t/s‘);
ylabel(‘第一層低頻系數(shù)幅值?A/W‘);
Y10=abs(fft(y10nfft));
subplot(222);
plot((0:nfft/2-1)/nfft*fsY10(1:nfft/2));
xlabel(‘第一層低頻系數(shù)頻率?f/Hz‘);
ylabel(‘第一層低頻系數(shù)功率譜?P/W‘);
subplot(223);
plot(y11);
xlabel(‘第一層高頻系數(shù)時間?t/s‘);
ylabel(‘第一層高頻系數(shù)幅值?A/W‘);
Y11=abs(fft(y11nfft));
subplot(224);
plot((0:nfft/2-1)/nfft*fsY11(1:nfft/2));
xlabel(‘第一層高頻系數(shù)頻率?f/Hz‘);
ylabel(‘第一層高頻系數(shù)功率譜?P/W‘);
%------------------------------------------------------------------
%-----------------------------------------------------------------
%重構(gòu)第二層各系數(shù)
y20=wprcoef(T[20]);
y21=wprcoef(T[21]);
y22=wprcoef(T[22]);
y23=wprcoef(T[23]);
figure(3);
subplot(421);
plot(y20);
xlabel(‘Y10分解的低頻系數(shù)時間t/s‘);
ylabel(‘幅值A/W‘);
Y20=abs(fft(y20nfft));
subplot(422);
plot((0:nfft/2-1)/nfft*fsY20(1:nfft/2));
xlabel(‘Y10分解的低頻系數(shù)頻率f/Hz‘);
ylabel(‘功率譜P/W‘);
subplot(423);
plot(y21);
xlabel(‘Y10分解的高頻系數(shù)時間t/s‘);
ylabel(‘幅值A/W‘);
Y21=abs(fft(y21nfft));
subplot(424);
plot((0:nfft/2-1)/nfft*fsY21(1:nfft/2));
xlabel(‘Y10分解的高頻系數(shù)頻率f/Hz‘);
ylabel(‘功率譜P/W‘);
subplot(425);
plot(y22);
xlabel(‘Y11分解的低頻系數(shù)時間t/s‘);
ylabel(‘幅值A/W‘);
Y22=abs(fft(y22nfft));
subplot(426);
plot((0:nfft/2-1)/nfft*fsY22(1:nfft/2));
xlabel(‘Y11分解的低頻系數(shù)頻率f/Hz‘);
ylabel(‘功率譜P/W‘);
subplot(427);
plot(y23);
xlabel(‘Y11分解的高頻系數(shù)時間t/s‘);
ylabel(‘幅值A/W‘);
Y23=abs(fft(y23nfft));
subplot(428);
plot((0:nfft/2-1)/nfft*fsY23(1:nfft/2));
xlabel(‘Y11分解的高頻系數(shù)頻率f/Hz‘);
ylabel(‘功率譜P/W‘);
%------------------------------------------------------------------
%-----------------------------------------------------------------
%重構(gòu)第三層各系數(shù)
y30=wprcoef(T[30]);
y31=wprcoef(T[31]);
y32=wprcoef(T[32]);
y33=wprcoef(T[33]);
y34=wprcoef(T[34]);
y35=wprcoef(T[35]);
y36=wprcoef(T[36]);
y37=wprcoef(T[37]);
figure(4);
subplot(921);
plot(y30);
xlabel(‘Y20分解的低頻系數(shù)時間t/s‘);
ylabel(‘幅值A/W‘);
Y30=abs(fft(y30nfft));
subplot(922);
plot((0:nfft/2-1)/nfft*fsY30(1:nfft/2));
xlabel(‘Y20分解的低頻系數(shù)頻率f/Hz‘);
ylabel(‘功率譜P/W‘);
subplot(923);
plot(y31);
xlabel(‘Y20分解的高頻系數(shù)時間t/s‘);
ylabel(‘幅值A/W‘);
Y31=abs(fft(y31nfft));
subplot(924);
plot((0:nfft/2-1)/nfft*fsY31(1:nfft/2));
xlabel(‘Y20分解的高頻系數(shù)頻率f/Hz‘);
ylabel(‘功率譜P/W‘);
subplot(925);
plot(y32);
xlabel(‘Y21
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????4855??2006-05-23?22:44??xiaobobao.m
-----------?---------??----------?-----??----
?????????????????4855????????????????????1
- 上一篇:小波能譜熵
- 下一篇:springboot入門視頻代碼教材資源.txt
評論
共有 條評論