|
我用matlab验证过了,代码如下:
t = 1:256;
y = sin(2*pi*t*1.25/200) + randn(1,256);
[c,l] = wavedec(y,5,'db4');
a5 = wrcoef('a',c,l,'db4',5);
yr = sin(2*pi*t*1.25/200);
plot(t,y,'-k',t,a5,'-r',t,yr,'-b')
验证的结果表明,采用db4小波滤波的效果还是蛮不错的,附图中,黑色的表示原始曲线;红色的滤波后的曲线,蓝色的表示理论曲线。
但是我想说明的是:你的边界延拓一般要采用matlab的方法;同时单支重构后采用FFT变换的第3个点近似,同直接计算FFT变换取第三个点近似结果是严格一致的(用模拟数据计算过),所以你不能用FFT计算的第三个点表示你计算得到的近似的1.25Khz的信号。
[ 本帖最后由 songjianjiang 于 2010-5-18 11:54 编辑 ] |
评分
-
1
查看全部评分
-
|