[已答复] 傅立叶滤波

[复制链接]
fan040908 发表于 2009-11-28 01:18:59
clc
n=0:255;
xn=0.1*n.*exp(-0.03*n);
noise=(1+0.5*rand(1,length(n))).*cos(2*pi*5000*n);
y=xn+noise;
subplot(221)
plot(xn)
subplot(222)
plot(y)
Y=fft(y,256);
H=[Y(1:49),zeros(1,157),Y(208:256)];
x=ifft(H,256);
subplot(2,2,[3,4])
plot(abs(x))
这个是书本上的一个例子,说的是傅立叶滤波的,可是出来的结果不是很一致也,课本上的滤波效果很好,只是尾巴出现一点点的波动,可是我出来的,就愁人了啊。是不是rand这个随机函数的频谱是在整个频率范围内存在的呢?noise=(1+0.5*rand(1,length(n))).*cos(2*pi*5000*n);虽然我已经把它搬到了高频,可是出来的效果真的是...愁人哦,各位师兄有空看看下哈,随便给小弟解答下下哈

4 条回复


songzy41 发表于 2009-11-28 15:46:44
原帖由 fan040908 于 2009-11-28 01:18 发表
傅立叶滤波
clc
n=0:255;
xn=0.1*n.*exp(-0.03*n);
noise=(1+0.5*rand(1,length(n))).*cos(2*pi*5000*n);
y=xn+noise;
subplot(221)
plot(xn)
subplot(222)
plot(y)
Y=fft(y,256);
H=[Y(1:49),zeros(1,157),Y(208:256)];
x=ifft(H,256);
subplot(2,2,[3,4])
plot(abs(x))

错误是在下语句:
H=[Y(1:49),zeros(1,157),Y(208:256)];
其一H长不是256;其二H的实部和虚部不构成共轭对称。现改为如下,LZ试一下:
H=[Y(1:49),zeros(1,159),Y(209:256)];

fan040908 发表于 2009-11-29 12:20:54
???那个实部和虚部不对称的话有什么后果呢?

urion 发表于 2009-11-29 19:31:36
请问2楼 这个滤波效果还是不怎么好

songzy41 发表于 2009-11-29 20:37:50
原帖由 fan040908 于 2009-11-29 12:20 发表
???那个实部和虚部不对称的话有什么后果呢?

这就会出现一些你不想要的频率分量。
原帖由 urion  于 2009-11-29 19:31 发表
请问2楼 这个滤波效果还是不怎么好

是的,因为噪声是宽带信号,用了带通滤波器同样带入了不少噪声分量。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
热门教程
站长推荐
快速回复 返回顶部 返回列表