马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
t=0:0.005:1*pi;
fs=20000;
s=4*sin(2*200*pi*t).*(sin(2*4500*pi*t))+25*(sin(2*4500*pi*t));
subplot(411);plot(t,s);title('原始信号');
%运用小波方法滤波
[c,l]=wavedec(s,1,'db10');
d1=wrcoef('d',c,l,'db10');
a1=0;
subplot(412);plot(d1);title('重构高频信号');
%hilbert变换,求包络线
y=hilbert(d1);
y1=abs(y);
subplot(413);plot(y1);title('包络线');
%FFT求包络谱
nfft=1024;
p=abs(fft(y1,nfft));
figure(2);
plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));
title('包络谱');xlabel('频率');ylabel('功率谱');
这个程序对么?画出来的对吗? |