数字信号处理实验三用双线性变换法设计IIR数字滤波器Word文档下载推荐.docx
《数字信号处理实验三用双线性变换法设计IIR数字滤波器Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《数字信号处理实验三用双线性变换法设计IIR数字滤波器Word文档下载推荐.docx（6页珍藏版）》请在冰豆网上搜索。
H（z）
16
0.0007378（1z）
A=0.09036
由（3.1）式和（3.2）式可见，滤波器H（z）由三个二阶滤波器Hi（z）、H2（z）和
H3（z）级联组成，如图3.1所示。
图3.1滤波器H（z）的组成
⑵编写滤波器仿真程序，计算H（z）对心电图信号抽样序列x（n）的响应序列y（n）
设yk（n）为第k级二阶滤波器Hk（z）的输出序列，yk-1（n）为输入序列，如图3.1所示。
由（3.2）式可得到差分方程：
yk（n）=Ayk-1（n）+2Ayk-1（n-1）+Ayk-1（n-2）+Bkyk（n-1）+Ckyk（n-2）——（3.3）
当k=1时，yk-1（n）=x（n）。
所以H（z）对x（n）的总响应序列y（n）可以用顺序迭代算法得到。
即依次对k=1,2,3，求解差分方程（3.3），最后得到y3（n）=y（n）。
仿真程序就是实现上述求解差分方程和顺序迭代算法的通用程序。
也可以直接调
用MATLABfilter函数实现仿真。
（3）在通用计算机上运行仿真滤波程序，并调用通用绘图子程序，完成“2.实验
内容⑵和（3）”。
本实验要用的MATLAB绘图函数见附录。
3.实验内容
（1）用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。
设计指标参数为：
在通带内频率低于0.2n时，最大衰减小于1dB;
在阻带内［0.3n,n］频率区间上，最小衰减大于15dB。
以x=[ones（1,10）,zeros（1,40）]信号验证。
⑵以0.02n为抽样间隔，打印出数字滤波器在频率区间［0,冗/2］上的幅频响应特性曲线。
数濾波器旳颉率响应
程序代码:
ws=0.3*pi;
Ap=1;
As=15;
T=1;
WP=（2/T）*tan（wp/2）;
WS=（2/T）*tan（ws/2）;
[N,Wc]=buttord（WP,WS,Ap,As,'
s'
）;
[bt,at]=butter（N,Wc,'
）
[bd,ad]=bilinear（bt,at,1/T）
%数字滤波器的频率响应
[hd,wd]=freqz（bd,ad,60）;
%画图
figure
（1）;
plot（wd/pi,abs（hd））,grid;
axis（[0,1,0,1]）
xlabel（'
w/pi'
ylabel（'
abs（h）'
）%定义X轴与Y轴的标称含义
title（'
数滤波器的频率响应'
）%定义图形名称
n=[1:
56];
%输入系列的长度
x=[-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];
（3）用所设计的滤波器对实际心电图信号抽样序列（在本实验后面给出）进行仿
真滤波处理，并分别打印出滤波前后的心电图信号波形图，观察总结滤波作
用与效果
时域
議isH旗
程序代码：
wp=0.2*pi;
WP=（2/T）*tan（wp/2）;
WS=（2/T）*tan（ws/2）;
[hd,wd]=freqz（bd,ad,n）;
x=[-4,
-2,0,-4,
-6,
-4,-2,-4,-6,
-6,-4,
-4,-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,
-2,-2,-2,0];
y=filter（bd,ad,x）;
figure
（1）
subplot（2,1,1）;
plot（n,x）;
inputsequence'
n'
x（n）'
subplot（2,1,2）;
plot（n,y）;
outputsequence'
y（n）'
k=[0:
2*pi/63:
2*pi];
figure
（2）;
X1=fft（x,64）;
X1=abs（X1）;
stem（k,X1）;
Y=fft（y,64）;
stem（k,abs（Y））;
4.思考题
用双线性变换法设计数字滤波器过程中，变换公式
中T的取值，对设计结果有无影响？
为什么?
答：
无影响，因为数字滤波器的传输函数H（ejw）以2为周期，最高频率在处，因此，s<
，按照线性关系ss/T,那么一定满足s/T,因此T可
以任选。
在数字滤波器的设计过程中，设计问题起始于数字滤波器上的传输函数
X（ejwk）
Hc（s）再通过公式变换得到数字滤波器的传输函数
X（e），这样取值T对设计
器的传输函数Hc（s）的技术要求并设计出模拟滤波，然后模拟滤波器的传输函数
的影响就被抵消，因此T的取值对设计结果没有影响。
121212~
（11.2686Z0.7051z）（11.0106z0.3583z）（10.9044z0.2155z）
升级会员 