实验三 用双线性变换法设计IIR数字滤波器
实验三IIR数字滤波器实施方案及软件实现
实验三IIR数字滤波器设计及软件实现1．实验目地（1）熟悉用双线性变换法设计IIR数字滤波器地原理与方法；（2）学会调用MATLAB信号处理工具箱中滤波器设计函数（或滤波器设计分析工具fdatool）设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数.b5E2R。
（3）掌握IIR数字滤波器地MATLAB实现方法.（3）通过观察滤波器输入输出信号地时域波形及其频谱，建立数字滤波地概念.2．实验原理设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛地是双线性变换法.基本设计过程是：①先将给定地数字滤波器地指标转换成过渡模拟滤波器地指标；②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器地系统函数.MATLAB信号处理工具箱中地各种IIR数字滤波器设计函数都是采用双线性变换法.第六章介绍地滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2和椭圆模拟和数字滤波器.本实验要求读者调用如上函数直接设计IIR数字滤波器.p1Ean。
本实验地数字滤波器地MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定地输入信号x(n)进行滤波，得到滤波后地输出信号y(n）.DXDiT。
3. 实验内容及步骤（1）调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成地复合信号st，该函数还会自动绘图显示st地时域波形和幅频特性曲线，如图10.4.1所示.由图可见，三路信号时域混叠无法在时域分离.但频域是分离地，所以可以通过滤波地方法在频域分离，这就是本实验地目地.RTCrp。
图10.4.1 三路调幅信号st 地时域波形和幅频特性曲线（2）要求将st 中三路调幅信号分离，通过观察st 地幅频特性曲线，分别确定可以分离st 中三路抑制载波单频调幅信号地三个滤波器（低通滤波器、带通滤波器、高通滤波器）地通带截止频率和阻带截止频率.要求滤波器地通带最大衰减为0.1dB,阻带最小衰减为60dB.5PCzV 。
用双线性变换法设计IIR数字滤波器
用双线性变换法设计IIR数字滤波器一.实验目的：1）熟悉用双线性变换法设计IIR数字滤波器的原理和方法。
2）掌握数字滤波器设计的计算机仿真方法。
3）观察对心电图信号的实际滤波作用，获得数字滤波的感性认识。
二.实验内容：1）用双线性变换法设计巴特沃斯IIR数字滤波器。
设计指标为：通带内频率低于0.3πrad时，最大衰减小于2dB；在阻带[0.6π,π]频率区间上,最大衰减大于20dB。
2）绘出数字滤波器的幅频响应特性曲线、相频特性曲线和零极图。
三.实验程序：wp=0.3*pi;ws=0.6*pi;Rp=2;As=20;Fs=100;T=1/Fs;Omgp=(2/T)*tan(wp/2);%原型通带频率预修正Omgs=(2/T)*tan(ws/2);%原型阻带频率预修正[n,Omgc]=buttord(Omgp,Omgs,Rp,As,'s');%计算阶数n和截至频率[ba1,aa1]=butter(n,Omgc,'s');%求模拟滤波器系数[bd,ad]=bilinear(ba1,aa1,Fs)%用双线性变换法求数字滤波器系数b,a[sos,g]=tf2sos(bd,ad);%由直接型变换为级联型[H,w]=freqz(bd,ad);%求数字系统的频率特性subplot(2,2,1);plot(w/pi,abs(H));title('幅度响应');axis([0,1,0,1.1]);subplot(2,2,2);plot(w/pi,angle(H)/pi);title('相位响应');axis([0,1,-1,1]);subplot(2,2,3);zplane(bd,ad);title('零极图');axis([-1.1,1.1,-1.1,1.1]);四.实验结果：程序运行结果：bd =0.0780 0.2339 0.2339 0.0780ad =1.0000 -0.7958 0.5020 -0.0824五．总结：通过这个实验，对设计数字滤波器的整个过程有了比较好的掌握。
实验三 用双线性变换法设计IIR数字滤波器
实验_三_题目_用双线性变换法设计IIR 数字滤波器 第16周星期_3_第6,7节一. 实验目的(1)熟悉用双线性变换法设计IIR 数字滤波器的原理与方法。
(2)掌握数字滤波器的计算机仿真方法。
(3)通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。
二、实验容、方法、设计程序及实验结果(1)复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计IIR 数字滤波器的容，用双线性变换法设计数字滤波器系统函数()z H 。
其中满足本实验要求的数字滤波器系统函数为：()()()()()212121612155.09044.013583.00106.117051.02686.1110007378.0-------+-+-+-+=zz z z z zz z H ()z H k k ∏==31（3.1）式中： ()()3211212121，，，k zC z B z z A z H k k k =--++=---- （3.2）2155.09044.03583.00106.17051.02686.109036.0332211-==-==-===C B C B C B A ，，，根据设计指标，调用MATLAB 信号处理工具箱buttord 和butter ，也可以得到()z H 。
由公式（3.1）和（3.2）可见，滤波器()z H 由三个二阶滤波器()z H 1、()z H 2和()z H 3级联而成，如图3-1所示。
()n y图3-1 滤波器z H 的组成此参数下的程序如下：%实验三，用双线性变换法设计 IIR数字滤波器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];k=1; %控制滤波循环变量close all; %关闭全部绘图窗口figure(1); %创建绘图窗口subplot(2,2,1); %定位子图 1n=0:55; %横坐标stem(n,x,'.'); %画出枝干图axis([0,56,-100,50]); %调整坐标xlabel('n'); %标注横坐标ylabel('x(n)'); %标注纵坐标title('心电图信号采集序列x(n)');%命名该子图B=[0.09036,2*0.09036,0.09036];%H1 滤波器的分子系数矩阵A=[1.2686,-0.7051]; %H1滤波器的分母系数矩阵A1=[1.0106,-0.3583]; %H2滤波器的分母系数矩阵A2=[0.9044,-0.2155]; %H3滤波器的分母系数矩阵while(k<=3)y=filter(B,A,x); %进行滤波x=y; %重新赋值X 进行下一次滤波k=k+1; %控制循环变量if k==2A=A1;else A=A2;endendsubplot(2,2,3); %定位子图3stem(n,y,'.');axis([0,56,-100,50]);xlabel('n');ylabel('y(n)');title('三级滤波后的心电图信号(原坐标)');subplot(2,2,2)stem(n,y,'.');axis([0,56,-15,5]);xlabel('n');ylabel('y(n)');title('调整坐标后的心电图信号');%求数字滤波器的幅频特性A=[0.09036,2*0.09036,0.09036];%滤波器的分子系数矩阵B1=[1,-1.2686,0.7051]; %H1滤波器的分母系数矩阵B2=[1,-1.0106,0.3583]; %H2滤波器的分母系数矩阵B3=[1,-0.9044,0.2155]; %H3滤波器的分母系数矩阵[H1,w]=freqz(A,B1,100); %进行滤波器幅频特性分析[H2,w]=freqz(A,B2,100);[H3,w]=freqz(A,B3,100);H4=H1.*(H2); %点积H=H4.*(H3);db=20*log10(abs(H)+eps);subplot(2,2,4)plot(w/pi,db);axis([0,0.5,-50,10]);xlabel('w');ylabel('|H(e^j^w)|');grid on; %显示方格title('滤波器的幅频响应曲线');（2）用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。
用双线性变换法设计IIR数字低通滤波器课程设计
V=课程设计报告书姓名：班级：学号：时间：和j Ω3π / T π / T－3π / T－π / Too σ－11jIm[z ]Re[z ]Z 平面S 平面图2脉冲响应不变法的映射关系由（1-9）式，数字滤波器的频率响应和模拟滤波器的频率响应间的关系为 (1-10)这就是说，数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。
正如采样定理所讨论的，只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以内时，即 (1-11)才能使数字滤波器的频率响应在折叠频率以内重现模拟滤波器的频率响应，而不产生混叠失真，即|ω|<π (1-12)但是，任何一个实际的模拟滤波器频率响应都不是严格限带的，变换后就会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真，如图2所示。
这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。
当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。
这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。
－3π－2π……)j (a ΩH Ωoo －π2π3ππω＝Ω T)(e j ωH T π2T πT πT π2－图3脉冲响应不变法中的频响混叠现象对某一模拟滤波器的单位冲激响应h a (t )进行采样，采样频率为f s ，若使f s 增加，即令采样时间间隔（T =1/f s ）减小，则系统频率响应各周期延拓分量之间相距更远，因而可减小频率响应的混叠效应。
⎪⎭⎫ ⎝⎛-=∑∞-∞=T k j H T e H k a j πωω21)(2||sT Ω=≥Ωπ0)(=Ωj H a ⎪⎭⎫ ⎝⎛=T j H T e H a j ωω1)(《信息处理综合课程设计》课程设计行政班级学号__ 姓名_选课班级__ 任课教师成绩________。
实验三 IIR数字滤波器的设计
《数字信号处理》（2010-2011学年第1学期）论文成绩：实验三 IIR数字滤波器的设计学生姓名：闫春遐所在院系：电子信息工程学院自动化系年级专业：2008级自动化系学号：00824049指导教师：王亮完成日期：2010年10月26日实验三 IIR数字滤波器的设计一、实验目的（1）掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR 数字滤波器的计算机编程。
（2）观察双线性变换及脉冲响应不变法设计滤波器的频域特性，了解双线性变换法及脉冲响应不变法的特点。
（3）熟悉巴特沃斯滤波器、切比雪夫和椭圆滤波器的频率特性。
二、实验内容实验中有关变量的定义：f通带边界频率cr f阻带边界频率δ通带波动At最小阻带衰减f采样频率sT采样周期上机实验内容：（1）c f=0.3 kHz，δ=0.8 dB，r f=0.2 kHz，At=20 dB，T=1 ms；设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。
解答：>> wc=2*1000*tan(2*pi*300/(2*1000));>> wr=2*1000*tan(2*pi*200/(2*1000));>> [N wn]=cheb1ord(wc,wr,0.8,20,'s');>> [B A]=cheby1(N,0.8,wn,'high','s');>> [num den]=bilinear(B,A,1000);>> [h w]=freqz(num,den);>> f=w/pi*500;>> plot(f,20*log10(abs(h)));grid;axis([0,500,-80,10]);xlabel('频率/Hz');ylabel('幅度/dB');（2）c f=0.2 kHz， =1 dB，r f=0.3 kHz，At=25 dB，T=1 ms；分别用脉冲响应不变法及双线性响应变换法设计一巴特沃斯数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。
【精品】用双线性变换法设计IIR数字低通滤波器课程设计
V=课程设计报告书姓名：班级:学号：时间:设计过程摘要：根据IIR滤波器的特点,在MATLAB坏境下用双线性变换法设计IIR 数字滤波器.利用MATLAB设计滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化.关键词：双线性变换法,数字滤波器，MATLAB，IIR1.设计原理与步骤1。
1设计原理滤波器的种类很多，从功能上可分为低通、高通、带通和带阻滤波器，每一种又有模拟滤波器和数字滤波器两种形式。
如果滤波器的输人和输出都是离散时间信号，则该滤波器的冲击响应也必然是离散的,这种滤波器称之为数字滤波器。
数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。
数字滤波器也是具有一定传输选择特性的数字信号处理装置，其输入、输出均为数字信号,实质上是一个由有限精度算法实现的线性时不变离散系统。
IIR数字滤波器采用递归型结构，即结构上带有反馈环路。
IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成,可以组合成直接型、正准型、级联型、并联型四种结构形式,都具有反馈回路。
数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR)数字滤波器和有限长冲激响应（FIR）数字滤波器。
IIR数字滤波器的特征是，具有无限持续时间冲激响应,需要用递归模型来实现，其差分方程为:（1-1）(1—2）设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H（z）,使其频率响应H(z)满足所希望得到的频域指标,即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数.本次课程设计采用双线性变换法。
1。
2设计步骤：（1）将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标；(2)根据转换后的技术指标设计模拟低通滤波器H （s ）；（3）在按一定规则将H （s)转换为H （z ）；若所设计的数字滤波器是低通的，那么上述设计工作可以结束，若所设计的是高通、带通或者带阻滤波器，那么还有步骤：(4）将高通、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标，然后按上述步骤(2）设计出模拟低通滤波器H(s)，再将H （s)转换为所需的H(z ）.2。
实验三、IIR数字滤波器设计和应用
丫L«Z 构4090459 R 昭L 静IVERSITY 学院物电学院专业、班级10电子实验课程名称数字信号处理实验 ______________ 教师及职称杨卫平教授_________________ 开课学期2013 至2013 学年下学期填报时间2013 年5 月10 日云南师范大学教务处编印Normalized Freque ncy ( X ： rad/sample)(2).设计一个数字带通滤波器，要求在100~200Hz 通带内纹波不大于 3dB ，通带两边各50Hz外是阻带，衰耗不小于40dB 。
抽样频率为1000Hz 。
Wp=[100 200]/500;Ws=[100-50 200+50]/500; Rp=3;Rs=40;[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs) [b,a]=cheby1(N,Rp,Wn); freqz(b,a,512,1000);title('Chebyshev Type I Ban dpass Filter'); axis([0,500,-80,0]);.se^rngHT esaB —11 1-500 00.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.90 0.1 0.2 0.30.4 0.5 0.6 0.7 0.8 0.91Normalized Freque ncy (心;rad/sample))BoceanrrnaaM-100 -200 -300 -400-203 •实验设备及材料个人计算机一台及 MATLAB 仿真软件。
_______________________________________________________4.实验方法步骤及注意事项 实验方法步骤：先打开电脑，然后再打开MATLAB 仿真软件，在Blank M-File 中输入程序，然后再编译运行程序，直到程序能编译运行为止。
实验三：用双线性变换设计IIR数字滤波器
实验三：用双线性变换设计IIR数字滤波器pb03023058 甘小华一、实验目的1、熟悉用双线性变换法设计IIR数字滤波器的原理与方法。
2、掌握数字滤波器的计算机方针方法。
3、通过观察对实际心电图信号的滤波作用，获得数字滤波器的感性知识。
二、实验原理冲激响应不变法设计IIR数字滤波器会引起频域的混叠，因此有了双线性变换设计IIR 数字滤波器。
双线性变换所采用的由s平面到z平面的映射为s=(z-1) / (z+1)即z=(1+s) / (1-s)三、实验内容1、用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。
设计参数：在通带内频率低于0.2*pi时，最大衰减小于1dB;在阻带内[0.3*pi,pi]频率区间上，最小衰减大于15dB。
2、以0.02*pi为采样间隔，打印出数字滤波器在频率区间[0,pi/2]上的幅频响应特性曲线。
3、用所设计的滤波器对实际心电图信号序列进行仿真滤波处理,可采用filter函数。
4、求得满足本实验要求的数字滤波器系统函数H(z)=H1(z)+H2(z)+H3(z)Hi(z)=A(1+2/z+1/zz) / (1-Bi/z-Ci/zz)式中A=0.09036B1=1.2686 C1=-0.7051B2=1.0106 C2=-0.3583B3=0.9044 C3=-0.21555、实验所用的matlab程序为：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]k=1;close all;figure(1);subplot(2,2,1);n=1:56;stem(n,x,'.');axis([0 56 -100 50]);hold on;n=0:60;m=zeros(61);plot(n,m);xlabel('n');ylabel('x(n)');title('心电图信号采样序列x(n)');B=[0.09036 2*0.09036 0.09036];A=[1.2686 -0.7051];A1=[1.0106 -0.3583];A2=[0.9044 -0.2155];while(k<=3)y=filter(B,A,x);x=y;if k==2A=A1;endif k==3A=A2;endk=k+1;endsubplot(2,2,3);n=1:56;stem(n,y,'.');axis([0 56 -15 5]);hold on;n=0:60m=zeros(61);plot(n,m);xlabel('n');ylabel('y(n)');title('三级滤波后的心电图信号');A=[0.09036,0.1872,0.09036];B1=[1,-1.2688,0.7051];B2=[1,-1.0106,0.3583];B3=[1,-0.9044,0.2155];[H1,w]=freqz(A,B1,100);[H2,w]=freqz(A,B2,100);[H3,w]=freqz(A,B3,100);H4=H1.*(H2);H=H4.*(H3);mag=abs(H);db=20*log10((mag+eps)/max(mag)); subplot(2,2,2);plot(w/pi,db);axis([0,0.6,-50,10]);xlabel('w/pi');ylabel('20log10|H(jw)|');title('滤波器的幅频响应曲线');6、 运行该程序，得到下面三幅图n x (n )(b)ny (n )(c)w/pi20l o g 10|H (j w )|(a)图(a) 滤波器的幅频衰减曲线(b) 心电图信号采样序列 (c) 三级滤波后的心电图信号7、 实验总结：（1） 双线性变换的特点：双线性变换采用的是非线性频率映射，消除了冲激响应不变法中的频域混叠现象，并且幅频响应的衰减比较快，这正是我们所希望的。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验＿三_题目_用双线性变换法设计IIR数字滤波器第１6周星期_３_第６,7节一.实验目得（1)熟悉用双线性变换法设计IIR数字滤波器得原理与方法。
(２）掌握数字滤波器得计算机仿真方法。
(３)通过观察对实际心电图信号得滤波作用,获得数字滤波得感性知识。
二、实验内容、方法、设计程序及实验结果(1)复习有关巴特沃斯模拟滤波器设计与用双线性变换法设计ＩIR数字滤波器得内容,用双线性变换法设计数字滤波器系统函数。
其中满足本实验要求得数字滤波器系统函数为:()()()()()212121612155.09044.013583.00106.117051.02686.1110007378.0-------+-+-+-+=zzzzzzzzHﻩﻩﻩﻩﻩﻩﻩ(3、1) 式中: ﻩ(3、２)根据设计指标,调用MATＬＡＢ信号处理工具箱butｔｏrd与ｂutｔer,也可以得到。
由公式(3、1)与（３、２)可见,滤波器由三个二阶滤波器、与级联而成,如图3-1所示。
%x=、、、k=1；%控制滤波循环变量clｏｓe alｌ; %关闭全部绘图窗口fｉgure(1);%创建绘图窗口subplot（2,2,1)；%定位子图１n=0:55; %横坐标sｔem(n,x,'、'）; ％画出枝干图ａxiｓ([0,56，-100,5０]）; %调整坐标ｘlａbeｌ（'n'); %标注横坐标ylaｂｅl(＇ｘ(ｎ)');%标注纵坐标titｌｅ('心电图信号采集序列x（n)'）;%命名该子图B=[0、090３６,2*0、09036,0、０9０36]；%H1 滤波器得分子系数矩阵Ａ=[1、２686,-0、7０51]; %H１滤波器得分母系数矩阵A1＝[１、0１0６，-0、３58３]；%H2滤波器得分母系数矩阵A2=［０、9044,-0、21５5］; ％H３滤波器得分母系数矩阵ｗhiｌe(k<=3)y=fiｌteｒ（B,A,x); %进行滤波图3-1 滤波器得组成x=ｙ; %重新赋值X进行下一次滤波k=k+１; %控制循环变量if k＝＝２A＝A1；eｌｓe Ａ=A2;endeｎdsubplot(2，２,３); ％定位子图3stem(n，ｙ,'、');axiｓ([0,56，－10０,5０]);ｘlａbel（＇n');yｌabel('ｙ（n）');ｔiｔｌe（'三级滤波后得心电图信号(原坐标)');suｂploｔ(2,2,２)stｅｍ(n，y,'、');axis(［0,５6,－１5，5]);ｘlabel（'n');ｙｌabel(＇y(n)'）；tｉtle（'调整坐标后得心电图信号');％求数字滤波器得幅频特性A=[０、０９036,2*0、09036，0、0９036］;%滤波器得分子系数矩阵B1＝［１,-1、２6８6,０、7051］；％H1滤波器得分母系数矩阵B2=［1,-１、0106,0、３5８3］; %H2滤波器得分母系数矩阵B３=［1,－0、9０44，０、２155］;%H3滤波器得分母系数矩阵[H1，w]=freqz（A,B1，1０0)；%进行滤波器幅频特性分析[H2,w］=freqｚ(Ａ，B2，100);［H3，w]=ｆreqz(A,B3,100);Ｈ４=Ｈ1、*(Ｈ２）；%点积H=H４、*(H３）;ｄb=2０*log1０(aｂs(H）＋eps）；subplot（2,2,4)plot(w/ｐi,ｄb);axiｓ(［0，0、5,-５0,1０]);xlａbel（'w');ylabel('|H(e^j^ｗ)｜＇);grｉd oｎ; %显示方格tｉtlｅ(＇滤波器得幅频响应曲线')；（2）用双线性变换法设计一个巴特沃斯低通IIＲ数字滤波器。
设计指标参数为:在通带内频率低于0、2π时,最大衰减小于１dB;在阻带内[0、3π,π］频率区间上,最小衰减大与15dB。
参数如下截图:程序设计如下:x=[-４,－２,0,-４,-6，-４,-2，-4,－６,-6,-4,-４，-6,-６,-2,6,12,８,0,-１6,、、、－38，－６０，-84，-9０，-66,-32，－4,-2,-4，8,12,12,10,6,6,６,4,、、、0，0,0,0,0,-2,-４,０，0，0，-２，-2,0,0，-2,－2，-2,-２，０]；k＝1; %控制滤波循环变量％关闭全部绘图窗口figurｅ（2); %创建绘图窗口subplot(２,2，1); ％定位子图 1n=0:55; %横坐标stｅm(n,x,'、＇）; ％画出枝干图ａxｉｓ([０,５６,-10０,5０］）；%调整坐标xlａbｅl('n')； %标注横坐标ylabｅl('x（n)')； %标注纵坐标title('心电图信号采集序列x(n)')；％命名该子图B=[0、0007,0、０0４4,0、０11１,0、０１48,０、０111,0、0044，0、０0０7]; %H1 滤波器得分子系数矩阵Ａ＝［１、0００0,-3、1８38,4、6228,-３、７80０,1、8１39，－0、4８01,０、0５4５］； %H1滤波器得分母系数矩阵ｙ=filtｅr(B,Ａ,ｘ)；ｓuｂpｌot(2,２,3); ％定位子图3sｔem(n,y,'、');ａxｉｓ([0，56,－100,5０］）;xlabel（'n');ylabel('y(n)');tｉｔlｅ('滤波后得心电图信号(原坐标）');suｂploｔ(2,２,2)stｅm(n,ｙ,'、');axis(［0,5６,－100，30])；xlaｂｅl('n'）;ylabel('y(n）＇);titｌe('调整坐标后得心电图信号')；%求数字滤波器得幅频特性A=[0、０007,0、0０４４,0、0111,0、0148,０、0１11,０、００４４，0、00０7]；%滤波器得分子系数矩阵B=[1、0000，-3、1８38,4、62２８,-3、7８00,1、8139,－０、4801,0、0545]; %Ｈ1滤波器得分母系数矩阵[Ｈ,w］=frｅqz（A,B,１00)；％进行滤波器幅频特性分析dｂ＝2０＊loｇ10（ａｂs(H)+ｅｐｓ）；subpｌot(2,２,4)pｌｏt(w/pi,db)；axis([0,0、5,－5０,1０]);xlａbel('w'）;ｙlaｂeｌ（'|H（ｅ^j＾w)｜');ｇrid on; %显示方格ｔｉtle（'滤波器得幅频响应曲线');三、(1)、双线性变换法得特点:1．模拟滤波器经过双线性变换后,不存在频率特性得混叠失真,因而对模拟滤波器得频率响应函数Ha（s）无限带要求,而且能够直接用于设计低通、高通、带通、带阻等各种类型得数字滤波器。
2．与冲激响应不变法中模拟频率与数字频率之间得线性关系ω=ΩTs不同得就是,双线性变换法中模拟滤波器得频率与所转换成得数字滤波器得频率之间就是非线性关系，但就是,如果事先进行频率预畸变,这种非线性关系不会使所设计得数字滤波器得幅频特性受到影响。
由于频率得非线性关系会产生相频特性失真，所以若对数字滤波器得相位特性要求较严,则不宜采用这种变换方法。
3．双线性变换方法比较容易,不需要将模拟系统函数进行部分分式分解。
（２）、数字滤波得过程:顺序迭代、设yｋ（n)为第k级二阶滤波器Hｋ(z）得输出序列， y ｋ－1 (ｎ）为输入序列。
由式3－１可得到差分方程:ｙk(n)=Ayｋ-1 (n)+2Ａｙk-1（ｎ-1)+Aｙk-1（n-2)＋B k y k(ｎ-1）+C kｙk（n-2）当k=1 时,yk－1（ｎ)＝x(n)。
所以H(z)对ｘ（n)得总响应序列y(n)可以用顺序迭代算法得到。
即依次对k=1，２,3，求解差分方程(3、3),最后得到y3(n)=ｙ(n)。
数字滤波得作用:通过对有噪声得心电图采集信号波形图与经过三级二阶滤波器滤波后得心电图信号波形图得对比分析,可以瞧出低通滤波器滤除信号中高频噪声得滤波效果,信号变得更平缓。
四、思考题答：Ｔ取值对结果没有影响。
设数字滤波器ｗ1处衰减为δ１,ｗ２处衰减为δ2双线性变换公式:由于就是低通滤波器:得:模拟巴特沃思低通滤波器得幅度平方函数为：令:α,N就是一个与T无关得量。
查表得归一化模拟低通滤波器为:则去归一化模拟低通滤波器为：数字系统函数为:由以上各式得:α就是一个与T无关得量。
H(z)就是一个与T无关得函数，所以Ｔ取值对结果没有影响。