求高手解答怎么求包络线

可乐比海盗

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('功率谱');

这个程序对么？画出来的对吗？

开心一族

这样求是不对的