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

serendipity_w · 发表于 2007-8-16 13:21

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

谢谢！

serendipity_w · 发表于 2007-8-16 13:21

谁能给发个m文件

HolySaint · 发表于 2007-8-16 14:19

学会搜索啊

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

serendipity_w · 发表于 2007-8-16 16:15

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)

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

破凰 · 发表于 2007-8-16 17:58

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

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

serendipity_w · 发表于 2007-8-16 20:15

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

破凰 · 发表于 2007-8-16 20:43

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

账号		自动登录	找回密码
密码			我要加入

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

本帖被以下淘专辑推荐:

评分

回复 #5 破凰的帖子

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

本帖被以下淘专辑推荐:

评分

回复 #5 破凰 的帖子

回复 #5 破凰的帖子