小波包分解并求能量谱

程序介绍小波分解的算法思想和小波能量谱的计算，在信号处理领域应用较广

%采样频率
fs=5000;
nfft=10240;
%信号
fid=fopen（‘a.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）;%最后一句为显示小波包四层分解树结构

%-----------------------------------------------------------------------
figure（2）;
subplot（211）;
%重构第一层各系数
y10=wprcoef（T[10]）;
y11=wprcoef（T[11]）;
figure（2）;
subplot（221）;
plot（y10）;
xlabel（‘第一层低频系数时间 t/s‘）;
ylabel（‘第一层低频系数幅值 A/W‘）;
Y10=abs（fft（y10nfft））;
subplot（222）;
plot（（0:nfft/2-1）/nfft*fsY10（1:nfft/2））;
xlabel（‘第一层低频系数频率 f/Hz‘）;
ylabel（‘第一层低频系数功率谱 P/W‘）;
subplot（223）;
plot（y11）;
xlabel（‘第一层高频系数时间 t/s‘）;
ylabel（‘第一层高频系数幅值 A/W‘）;
Y11=abs（fft（y11nfft））;
subplot（224）;
plot（（0:nfft/2-1）/nfft*fsY11（1:nfft/2））;
xlabel（‘第一层高频系数频率 f/Hz‘）;
ylabel（‘第一层高频系数功率谱 P/W‘）;
%------------------------------------------------------------------

%-----------------------------------------------------------------
%重构第二层各系数
y20=wprcoef（T[20]）;
y21=wprcoef（T[21]）;
y22=wprcoef（T[22]）;
y23=wprcoef（T[23]）;
figure（3）;
subplot（421）;
plot（y20）;
xlabel（‘Y10分解的低频系数时间t/s‘）;
ylabel（‘幅值A/W‘）;
Y20=abs（fft（y20nfft））;
subplot（422）;
plot（（0:nfft/2-1）/nfft*fsY20（1:nfft/2））;
xlabel（‘Y10分解的低频系数频率f/Hz‘）;
ylabel（‘功率谱P/W‘）;
subplot（423）;
plot（y21）;
xlabel（‘Y10分解的高频系数时间t/s‘）;
ylabel（‘幅值A/W‘）;
Y21=abs（fft（y21nfft））;
subplot（424）;
plot（（0:nfft/2-1）/nfft*fsY21（1:nfft/2））;
xlabel（‘Y10分解的高频系数频率f/Hz‘）;
ylabel（‘功率谱P/W‘）;
subplot（425）;
plot（y22）;
xlabel（‘Y11分解的低频系数时间t/s‘）;
ylabel（‘幅值A/W‘）;
Y22=abs（fft（y22nfft））;
subplot（426）;
plot（（0:nfft/2-1）/nfft*fsY22（1:nfft/2））;
xlabel（‘Y11分解的低频系数频率f/Hz‘）;
ylabel（‘功率谱P/W‘）;
subplot（427）;
plot（y23）;
xlabel（‘Y11分解的高频系数时间t/s‘）;
ylabel（‘幅值A/W‘）;
Y23=abs（fft（y23nfft））;
subplot（428）;
plot（（0:nfft/2-1）/nfft*fsY23（1:nfft/2））;
xlabel（‘Y11分解的高频系数频率f/Hz‘）;
ylabel（‘功率谱P/W‘）;
%------------------------------------------------------------------

%-----------------------------------------------------------------
%重构第三层各系数
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分解的低频系数时间t/s‘）;
ylabel（‘幅值A/W‘）;
Y30=abs（fft（y30nfft））;
subplot（922）;
plot（（0:nfft/2-1）/nfft*fsY30（1:nfft/2））;
xlabel（‘Y20分解的低频系数频率f/Hz‘）;
ylabel（‘功率谱P/W‘）;
subplot（923）;
plot（y31）;
xlabel（‘Y20分解的高频系数时间t/s‘）;
ylabel（‘幅值A/W‘）;
Y31=abs（fft（y31nfft））;
subplot（924）;
plot（（0:nfft/2-1）/nfft*fsY31（1:nfft/2））;
xlabel（‘Y20分解的高频系数频率f/Hz‘）;
ylabel（‘功率谱P/W‘）;
subplot（925）;
plot（y32）;
xlabel（‘Y21

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4855  2006-05-23 22:44  xiaobobao.m

----------- ---------  ---------- -----  ----

                 4855                    1