bishe.asv

来自「以上是我大学毕业设计」· ASV 代码 · 共 24 行

load a.txt;
xdata= a;
figure(1);
plot(xdata);
fbh=cwt(xdata,1:2:64,'cgau4','plot');
%采样频率
s=real(fbh)
plot
%5层小波包分解
T=wpdec(s,5,'db2');
%重构低频信号
y1=wprcoef(T,[5,2]);
%y1的波形
figure(2);
subplot(2,2,1);
plot(y1);
xlabel('时间 t/n');
ylabel('电压 V/v');
%y1的功率谱
Y1=abs(fft(y1,nfft));
subplot(2,2,2);
plot((0:nfft/2-1)/nfft*fs,Y1(1:nfft/2));
xlabel('频率 f/Hz');
ylabel('功率谱 P/W');