论坛 › 信号处理方法 › 小波变换的问题

qqvirile 发表于 2008-12-3 20:31

小波变换的问题

我想要处理信号x=sin(10*pi*t)+sin(70*pi*t)+sin(150*pi*t)+sin(300*pi*t);但是我小波变化之后的结果看不出任何信息！请各位高手指点
n=400;
t=1:n;
x=sin(10*pi*t)+sin(70*pi*t)+sin(150*pi*t)+sin(300*pi*t);
=wavedec(x,3,'sym4');
=detcoef(c,l,);
ca3=appcoef(c,l,'sym4',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');

qqvirile 发表于 2008-12-3 20:32

cd1的图

这是cd1的图

mochen 发表于 2008-12-4 09:31

应该是你的显示那几句程序有问题

qqvirile 发表于 2008-12-4 12:41

回复 板凳 mochen 的帖子

我还是觉得程序本身有问题 我的绘图程序不应该有错误,我换个小波函数 但是结果是一样的 还是会有无数的频率分量！我很差异这个问题 我的信号很简单的 哪来的那么多的频率分量啊！！请继续指点 期待大家的指点！
我对小波的理论已经了解很多了 但是具体到作图的时候还是会有很多的困惑 希望大家一起来关注！

qqvirile 发表于 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);
=wavedec(s,3,'db40');
=detcoef(c,l,);
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));

graduate 发表于 2008-12-5 16:33

三级小波分解得到的 bookkeeping vector L什么作用？？

mingmingtree 发表于 2012-3-20 15:38

楼主现在还在学小波么？多想能够交流一下啊！！

查看完整版本: 小波变换的问题