声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1014|回复: 0

[滤波] IIR椭圆带通滤波 频带设置问题 附完整简洁的程序和问题描述

[复制链接]
发表于 2010-2-4 19:27 | 显示全部楼层 |阅读模式

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

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

x
%需求描述:设计一个0.5-30HZ的带通滤波器
%问题描述:原应出现0.5-30HZ内的5条谱线,实际上一条也没出现,如图所示。

clc;
clear all;
fs=512;% 采样频率(HZ)


%%%%%%%%%%%%%%%%%%%%%%%%%%
%我认为以下4行代码是问题关键,但不知如何设置才满足0.5-30HZ的滤波需求。

wp1=0.5;%允许最低频率(角频率)
wp2=30;%允许最高频率(角频率)
ws1=0.1; %最低截止频率(角频率)
ws2=31;%最高截止频率(角频率)

%%%%%%%%%%%%%%%%%%%%%%%%%



wp=[wp1 wp2];
ws=[ws1 ws2];
rp=2 ;%波纹系数
rs=60 ;%波纹系数
[n,wn]=ellipord(wp/(fs/2),ws/(fs/2),rp,rs) ;%获取固有频率和阶
[b,a]=ellip(n,rp,rs,wn) ;%返回传递函数参数
%freqz(b,a,512,1000);%1000是频率点数
[h w]=freqz(b,a);
figure(1);
plot(w*fs/(2*pi),abs(h));grid;%频率响应
xlabel('频率响应')
t=[0:1/512:100];

%这是待滤波信号,全由简单的正弦谐波构成。
y1=sin(t)+sin(10*t)+sin(20*t)+sin(30*t)+sin(100*t)+sin(200*t);


y=filter(b,a,y1) ;%用椭圆滤波器滤波,y为滤波后的结果
z=fft(y);
N=length(t);
Fs=512;
%f0=1/(2*pi);
F=([0:N-1])*Fs/N;
out=abs(z);
figure(2)
plot(F(1:N/4),out(1:N/4));
xlabel('HZ')

[ 本帖最后由 sdf4940 于 2010-2-4 19:28 编辑 ]
频谱.jpg
频率响应.jpg
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 17:53 , Processed in 0.063043 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表