如何设计截止频率为a的低通滤波器
谢谢！ 谁能给发个m文件 学会搜索啊http://forum.vibunion.com/forum/thread-49432-1-1.html 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=*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;
大家看看折算不算是一个低通滤波器啊，还有为什么最后换出来的三个图都不一样啊？
[ 本帖最后由 serendipity_w 于 2007-8-16 16:16 编辑 ] 首先取得滤波器的带宽不对，其次不能这样来设计时域滤波器的系数。你设计的滤波器的频率响应如下图（幅值太小）
频域滤波参考此帖
FFT滤波的物理意义
回复 #5 破凰 的帖子
看了你的回帖有些所得，不过还是不太明白，低通滤波器到底要怎么设计那？ 可以找一本数字信号处理的书来看
页:
[1]