如何设计截止频率为a的低通滤波器

serendipity_w

谢谢！

serendipity_w

谁能给发个m文件

HolySaint

学会搜索啊

http://forum.vibunion.com/forum/thread-49432-1-1.html

serendipity_w

clear all;
close all;
fs=500;
t=1/fs:1/fs:10;
x=sin(2*pi*1*t)+sin(2*pi*2*t)+sin(2*pi*4*t)+sin(2*pi*6*t);%频率成分由1，2，4，6HZ组成
s1=sin(2*pi*1*t)+sin(2*pi*2*t)+sin(2*pi*4*t);
N=length(t);
df=fs/N;
f1=[0:length(t)-1]*fs/N;
F=fft(x);
low_pass=rectpuls(f1,4.1);%滤除Hx的成分
res1=conv(x, real(ifft(low_pass)));%时域卷积的结果
res2=real(ifft(F.*low_pass));%频域乘积的反变换
res1应该和res2，s1一样才对啊，画出图来为什么不一样啊？
figure;
subplot(311);plot(t,res1(1:length(t)));grid on;%
subplot(312);plot(t,res2);grid on;
subplot(313);plot(t,s1);grid on;

大家看看折算不算是一个低通滤波器啊，还有为什么最后换出来的三个图都不一样啊？
untitled.fig (129.96 KB, 下载次数: 16)

2007-8-16 16:16 上传

点击文件名下载附件

[ 本帖最后由 serendipity_w 于 2007-8-16 16:16 编辑 ]

破凰

首先取得滤波器的带宽不对，其次不能这样来设计时域滤波器的系数。你设计的滤波器的频率响应如下图（幅值太小）

频域滤波参考此帖
FFT滤波的物理意义

serendipity_w

看了你的回帖有些所得，不过还是不太明白，
低通滤波器到底要怎么设计那？

破凰

可以找一本数字信号处理的书来看