|
楼主 |
发表于 2008-12-4 13:55
|
显示全部楼层
回复 板凳 mochen 的帖子
不好意思我错了 我做出来了 是程序错了
下面我重新给出程序 采样频率是400Hz 取2048个点!
n=1:2048;
T=1/400;
t=T*n;
s=sin(10*pi*t)+sin(70*pi*t)+sin(150*pi*t)+sin(300*pi*t);
[c,l]=wavedec(s,3,'db40');
[cd1,cd2,cd3]=detcoef(c,l,[1,2,3]);
ca3=appcoef(c,l,'db40',3);
figure(1)
subplot(511);plot(1:length(s),s);title('s');
subplot(512);plot(1:l(1),ca3);title('ca3');
subplot(513);plot(1:l(2),cd3);title('cd3');
subplot(514);plot(1:l(3),cd2);title('cd2');
subplot(515);plot(1:l(4),cd1);title('cd1');
figure(2)
y=fft(s,2048) %我做了s的频谱来验证一下杨建国老师的结论!
k=1:2048;
plot(400*k/2048,abs(y)); |
|