没有合适的资源?快使用搜索试试~ 我知道了~
IIR数字滤波器和FIR数字滤波器设计
4星 · 超过85%的资源 需积分: 17 23 下载量 159 浏览量
2011-04-13
09:09:38
上传
评论
收藏 262KB DOC 举报
温馨提示
试读
14页
有关滤波器的设计,用双线性变换法设计IIR数字滤波器,用窗函数法设计 FIR数字滤波器。
资源推荐
资源详情
资源评论
实验三:用双线性变换法设计 IIR 数字滤波器
1. 实验目的
(1) 熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法。
(2) 掌握数字滤波器的计算机仿真方法。
(3) 通过观察对实际心电图信号的滤波作用,获得数字滤波的感性知识。
2. 实验内容
(1) 用双线性变换法设计一个巴特沃斯低通 IIR 数字滤波器。设计指标参数为:在通带内频
率低于 0.2π 时,最大衰减小于 1dB;在阻带内[0.3π, π] 频率区间上,最小衰减大于 15dB。
(2) 以 0.02π 为抽样间隔,打印出数字滤波器在频率区间[0, π/2]上的幅频响应特性曲线。
(3) 用所设计的滤波器对实际心电图信号抽样序列(在本实验后面给出)进行仿真滤波处理 ,
并分别打印出滤波前后的心电图信号波形图,观察总结滤波作用与效果。
3. 实验步骤
(1) 复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计 IIR 数字滤波器的内容,按照
“2. 实验内容(1)”的要求设计满足指标的数字滤波器函数 H(z)。
现给出满足本实验要求的数字滤波器系统函数:
H(z) =
= —— (3.1)
式中 H
k
(z) = , k = 1 , 2 , 3 —— (3.2)
A = 0.09036
B
1
= 1.2686 , C
1
=-0.7051
B
2
= 1.0106 , C
2
=-0.3583
B
3
= 0.9044 , C
3
=-0.2155
由 (3.1) 式和(3.2)式可见,滤波器 H(z)由三个二阶滤波器 H
1
(z)、H
2
(z)和 H
3
(z)级联组成,
如图 3.1 所示。
图 3.1 滤波器 H(z)的组成
(2) 编写滤波器仿真程序,计算 H(z)对心电图信号抽样序列 x(n)的响应序列 y(n)。
设 y
k
(n)为第 k 级二阶滤波器 H
k
(z)的输出序列,y
k-1
(n)为输入序列,如图 3.1 所示。由
(3.2)式可得到差分方程:
y
k
(n) = A y
k-1
(n) + 2A y
k-1
(n-1) + A y
k-1
(n-2) + B
k
y
k
(n-1) + C
k
y
k
(n-2) —— (3.3)
当 k = 1 时,y
k-1
(n) = x(n)。所以 H(z)对 x(n)的总响应序列 y(n)可以用顺序迭代算法得到。即
依次对 k = 1,2,3,求解差分方程 (3.3),最后得到 y
3
(n) = y(n)。仿真程序就是实现上述求
H
1
(z)
H
2
(z)
H
3
(z)
x(n)
y
1
(n) y
2
(n)
y
3
(n) = y(n)
H(z)
解差分方程和顺序迭代算法的通用程序。也可以直接调用 MATLAB filter 函数实现仿真。
(3) 在通用计算机上运行仿真滤波程序,并调用通用绘图子程序,完成“ 2.实验内容(2)和(3)
”。
本实验要用的 MATLAB 绘图函数见附录。
4. 思考题
用双线性变换法设计数字滤波器过程中,变换公式
s =
中 T 的取值,对设计结果有无影响?为什么?
5. 心电图信号抽样序列 x(n)
人体心电图信号在测量过程中往往受到工业高频干扰,所以必须经过低通滤波处理后,
才能作为判断心脏功能的有用信息。下面给出一实际心电图信号抽样序列样本 x(n),其中
存在高频干扰。在实验中,以 x(n)作为输入序列,滤除其中的干扰成分。
{x(n)}={-4, -2, 0, -4, -6, -4, -2, -4, -6, -6,
-4, -4, -6, -6, -2, 6, 12, 8, 0, -16,
-38,-60, -84, -90, -66, -32, -4, -2, -4, 8,
12, 12, 10, 6, 6, 6, 4, 0, 0, 0,
0, 0, -2, -4, 0, 0, 0, -2, -2, 0,
0, -2, -2, -2, -2, 0,}
程序 1:
b=1.2;a=[1,-0.9,0.3]; %y(n)-0.9y(n-1)+0.3y(n-2)=1.2x(n)
x=[ones(1,10),zeros(1,40)];%利用函数 ones( )和 zeros(),产生一个 50 点序列,前 10 个点值
为 1,后 40 个点值为 0.
y=filter(b,a,x);%信号 x 通过滤波器后的输出
n=0:49;
%创建图形窗口 1,绘制 x(n)和 y(n)的时域序列图,图中要求定义 X 轴与 Y 轴的标称含义,
图形的名称
figure(1)
subplot(2,1,1);stem(n,x,'.');axis([0 50 0 1]);title('x(n)的时域序列图');
subplot(2,1,2);stem(n,y,'.');axis([0 50 0 4]);title('y(n)的时域序列图');
%对 x(n)和 y(n)做 64 点 FFT
N=64;p=0:63;
X1=fft(x,N); X2=fft(y,N);Y1=abs(X1);Y2=abs(X2);
%创建图形窗口 2,绘制 X(k)和 Y(k)的幅频特性,横轴显示为(0-2);并观察结果,分析滤
波器特性
figure(2)
k=[0:63];
Wk=(2*k)/N;
Wk=Wk/pi;
subplot(2,1,1);stem(p,Y1);title('Y1 的幅频特性')%定义图形名称
subplot(2,1,2);stem(p,Y2);title('Y2 的幅频特性')%定义图形名称
%创建图形窗口 3,绘制滤波器幅频特性
[h,w]=freqz(b,a);%计算频率响应
figure(3);
plot(w/pi,abs(h));grid;
xlabel('w/pi');ylabel('abs(h)')%定义 X 轴与 Y 轴的标称含义
title('滤波器的频率响应')%定义图形名称
结果:
剩余13页未读,继续阅读
资源评论
- 想飞的星星2013-03-19比较实用。可以使用。
飞雪3363
- 粉丝: 0
- 资源: 3
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功