求高手解答怎么求包络线
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('原始信号');
%运用小波方法滤波
=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('功率谱');
这个程序对么？画出来的对吗？ 这样求是不对的{:{40}:}
页:
[1]