声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1001|回复: 0

[小波] 小波分解程序有一处错误,请大家帮忙看看,谢谢!

[复制链接]
发表于 2009-5-28 10:00 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
load wd8;
data1=data(:,4);
data1=data1-mean(data1);         %隔直流
s1=data1;
N=1024;
fs=1280;       %采样频率
Y=fft(data1(1:N));
pyy=Y.*conj(Y)/N;
w='db3';
[c,l]=wavedec(s1,6,w);
for i=1:6
    A(i,:)=wrcoef('a',c,l,w,i);
    D(i,:)=wrcoef('d',c,l,w,i);
end
figure(1);
f=fs*(0:N/2-1)/N;
plot(f,pyy(1:N/2));
title('wd8信号的傅里叶变换系数','fontsize',18);
figure(2);
t=1:1000;                %为什么此处改为t=(1:N)/fs,就不能出来小波分解的系数图了???
subplot(7,2,1);plot(t,s1(t),'r');
title('Orig.signal and approx. 1 to6.');
subplot(7,2,2);plot(t,s1(t),'r');
title('Orig.signal and details 1 to6.');
for i=1:6
    subplot(7,2,2*i+1);plot(t,A(6-i+1,t),'b');
    subplot(7,2,2*i+2);plot(t,D(6-i+1,t),'g');
end


,在此程序中,fft没有问题,频率坐标也对应,上面t定义为1:1000时,可以正常画图,但是坐标和时间不对应啊,那样只代表两个采样点,其时间并不是1秒,请大家帮忙解释一下,急求,谢谢各位了!
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-13 20:32 , Processed in 0.106800 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表