|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
rt,这是我写的小波程序,信号比较简单,目的是想弄清楚一些基本概念!
clc
N=512
n=1:N
fs=4
x=sin(1.2*pi*n/fs)+sin(0.15*pi*n/fs)
[c,l]=wavedec(x,3,'db4')
cA3=appcoef(c,l,'db4',3)
[cd1,cd2,cd3]=detcoef(c,l,[1,2,3])
subplot(511)
plot(x)
subplot(512);plot(cA3)
subplot(513);plot(cd3)
subplot(514);plot(cd2);
subplot(515);plot(cd1)
如果大家运行了一下,可以发现其中的L似乎不对,我得出的l =70 70 133 259 512,为什么不是理论上的64 64 128 256 512?以小弟的理解频带是以fs/2开始往下划分的,也即v0空间对应的最高频是fs/2,程序中信号频率是0.6hz和0.075hz,而fs/2=2hz,那么在d1和a3中应该无频率成分,理论上只有在d2和d3中有频率成分,请问为什么会出现图中这种情况?谢谢大家··· |
|