双线性变换法IIR数字滤波器设计
用双线性变换法设计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数字滤波器实验四双线性变换法设计IIR数字滤波器（实验学时：2学时）⼀、实验⽬的1、掌握IIR数字滤波器设计的⼀般⽅法和步骤；2、学会利⽤给定的技术指标，设计满⾜设计要求的IIR模拟原型滤波器；3、掌握⽤双线性变换法设计IIR数字滤波器的原理与⽅法；4、掌握数字滤波器的计算机仿真⽅法；5、学会⽤所设计的滤波器对实际信号进⾏滤波处理，获得对数字滤波的感性认识。
⼆、实验基本原理与⽅法⽬前，IIR数字滤波器设计的最通⽤的⽅法是借助于模拟滤波器的设计⽅法。
模拟滤波器的设计已经有了⼀套相当成熟的⽅法，它不但有完整的设计公式，⽽且还有较为完整的图表供查询，因此，充分利⽤这些已有的资源将会给数字滤波器的设计带来很⼤的⽅便。
IIR数字滤波器的设计步骤如下：1、按⼀定规则将给出的数字滤波器的技术指标转换成模拟低通滤波器的技术指标；2、根据转换后的技术指标设计模拟低通滤波器()G s。
常⽤的模拟低通滤波器主要有巴特沃什模拟滤波器、切⽐雪夫模拟滤波器、椭圆滤波器等。
各种模拟滤波器的特点及设计⽅法参看教材第三章内容，在Matlab中，各种模拟滤波器的设计都有相应的函数来实现。
3、将设计出的模拟滤波器按⼀定规则转换成数字滤波器()H z。
常⽤的转换⽅法有两种：冲激响应不变法和双线性映射法。
这两种⽅法的特点及转换公式可参考教材第三章的内容，在Matlab中，⼆者分别由impinvar函数和bilinear函数实现。
若设计的数字滤波器是低通的，那么上述设计⼯作可以结束，若所设计的是⾼通、带通或带阻滤波器，那么还有步骤（4）；（4）将⾼通、带通或带阻数字滤波器的技术指标先转化为模拟低通滤波器的技术指标，然后按上述步骤（2）设计出模拟低通数字滤波器，再将其转换为所需的数字滤波器。
三、实验内容及步骤本实验主要包括以下三个⽅⾯的内容：1、设计⼀巴特沃什模拟低通滤波器，技术指标参数如下：在通频带内频率低于10Hz 时，最⼤衰减⼩于1dB ；在阻带内频率⼤于15Hz 时，最⼩衰减⼤于15dB 。
双线性变换法设计IIR数字滤波器通信专业
双线性变换法设计IIR 数字滤波器一、 实验目的1) 熟悉用双线性变换法设计IIR 数字滤波器的原理和方法。
2) 掌握数字滤波器设计的计算机仿真方法。
3) 观察对心电 图信号的实际滤波作用，获得数字滤波的感性认识。
二、 实验原理在滤波器的设计过程中，逼近是一个重要的环节，所谓逼近就是根据性能指标的要求，对理想特性进行逼近，以求得一个因果、稳定且可实现的传递函数。
目前模拟滤波器的频域设计理论已经发展的相当成熟，它不仅具有简单而严格的设计公式，而且设计参数已经表格化了。
借助模拟滤波器的逼近方法，用模拟域到数字域的某种变换,完成数字滤波器的逼进,这是一类简单而行之有效的方法。
双线性变换法采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-/T~/T 之间，再用z=e 转换Z 平面上。
也就是说，第一步先将整个S 平面压缩映射到S1平面的-/T~/T 一条横带里；第二步再通过标准变换关系z=e 将此横带变换到整个Z 平面上去。
这样就使S 平面与Z 平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，影射关系如图7-1所示。
图 7-1 双线性变换的映射关系为了将S 平面的整个虚轴j Ω压缩到S1平面j Ω1轴上的-π/T 到π/T 段上，可以通过以下的正切实现) 2tan(21T T Ω=Ω （7-1）式中T 仍是采样间隔。
当Ω1由-π/T 经过0变化到-π/T 时，Ω由-∞+∞-变化到经过0,也即映射了整个j Ω轴。
将式（7-1）写成Tj T j T j T j eee e T j 11112/2/2Ω-ΩΩ-Ω+-=Ω (7-2)将此关系到解析延拓到整个S 平面和S1平面，令j Ω=s ，j Ω1=s1，则得T s Ts ee T Ts Teee e T s T s T s T s T s 111111112)2tanh(2212/2/2/2/--+-⋅==+-⋅=-- (7-3) 再将s1平面通过以下标准变换关系映射到Z 平面： Ts e1z =(7-4)从而得到S 平面和Z 平面的单值映射关系为11112--+-⋅=z z T s (7-5) sTsT s T s T z -+=-+=222121 (7-6)式（7-5）与式（7-6）是s 平面与z 平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双母性变换。
实验十二 用双线性变换法设计IIR数字滤波器
用双线性变换法设计IIR数字滤波器[实验目的]熟悉模拟巴特沃茨滤波器设计和用双线性变换法设计IIR数字滤波器的方法。
[实验原理]利用双线性变换设计IIR数字滤波器，首先要设计出满足指标要求的模拟滤波器的传递函数Ha(S),然后由Ha(S)通过双线性可得所要设计的IIR数字滤波器的系统函数H(Z)。
[实验内容及要求]用双线性变换法设计IIR滤波器是IIR滤波器设计的经典方法，首先根据模拟滤波器的指标设计出来相应的模拟滤波器，然后将设计好的模拟滤波器转换为满足给定技术指标的数字滤波器。
在MATLAB的数字信号处理工具箱中提供了相应的设计函数，常用的有：1、butterworth数字和模拟2、零极点增益模型到传递函数模型的转换函数3、低通到低通模拟滤波器的转换函数4、双线性变换函数5、Butterworth数字和模拟滤波器设计函数6、计算butterworth低通原型模拟滤波器的极点和增益函数7、计算离散时间系统的频率响应函数用双线性变换法设计BUTTERWORTH低通IIR数字滤波器，要求使用buttord，butter，biliner。
滤波器技术指标：取样频率1HZ，通带内临界频率0.2HZ，通带内衰减小于1db，阻带临界频率0.3HZ，阻带内衰减大于25DB。
解：wp=0.2*2*pi;ws=0.3*2*pi;rp=1;rs=25;fs=1;ts=1/fs;wp2=2*fs*tan(wp/2*ts);ws2=2*fs*tan(ws/2*ts);[n,wn]=buttord(wp2,ws2,rp,rs,'s')[z,p,k]=buttap(n)[bap,aap]=zp2tf(z,p,k)[b,a]=lp2lp(bap,aap,wn)[bz,az]=bilinear(b,a,fs);[h,w]=freqz(bz,az);subplot(2,1,1);plot(w/pi,abs(h));grid onxlabel('频率');ylabel('幅度');subplot(2,1,2);plot(w/pi,20*log10(abs(h)));grid onxlabel('频率');ylabel('幅度');n = 6wn = 1.7043z =[]p =-0.2588 + 0.9659i-0.2588 - 0.9659i-0.7071 + 0.7071i-0.7071 - 0.7071i-0.9659 + 0.2588i-0.9659 - 0.2588ik =1bap =0 0 0 0 0 0 1 aap =1.0000 3.8637 7.4641 9.1416 7.4641 3.8637 1.0000 b =24.5076a =1.0000 6.5850 21.6809 45.2558 62.976555.5591 24.5076频率幅度频率幅度。
实验：用双线性变换法设计IIR数字滤波器
（设计性 4学时）一.实验目的：(1）熟悉用双线性变换法设计IIR 数字滤波器的原理与方法。
(2）掌握数字滤波器的计算机仿真方法。
(3）通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。
二.实验内容及步骤:(1) 用双线性变换法设计一个巴特沃斯低通IIR 数字滤波器，设计指标参数为：在通带内频率低于时,最大衰减小于1dB;在阻带内[ , pi] 频率区间上，,最小衰减大于15dB ；(2) 以 为采样间隔，打印出数字滤波器在频率区间[ 0, ]上的幅频响应特性曲线；(3) 用所设计的滤波器对实际心电图信号采样序列(在本实验后面给出)进行仿真滤波处理，并分别打印出滤波前后的心电图波形图，观察总结滤波作用与效果。
（4）采用不同阶数的Butterworth 低通滤波器，比较滤波效果。
三.实验步骤：（1）复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计IIR 数字滤波器的内容，按照教材例6.4.2，用双线性变换法设计数字滤波器系统函数H （z ）。
方法一：教材例6.4.2种已求出满足本实验要求的数字滤波器系统函数：方法二：根据设计指标，调用MATLAB 信号处理工具箱函数buttord 和butter ，也可得到H （z ）。
（2）编写滤波器仿真程序，计算H(z)对心电图信号采样序列x(n)的相应序列y(n)。
（3）在通过计算机上运行仿真滤波程序，并调用通用绘图子程序，完成实验内容（2）和（3）。
本实验要用的MATLAB 绘图函数参阅教材。
四.,思考题:用双线性变换法设计数字滤波器过程中,变换公式： s=11z1z 1T 2--+-中T 的取值,对设计结果有无影响? 为什么? 五.实验报告要求（1）简述实验目的及原理；（2）由所打印的特性曲线及设计过程简述双线性变换法的特点；（3）对比滤波前后的心电图信号波形,说明数字滤波器的滤波过程与滤波作用；(4) 简要回答思考题.六:心电图信号采样序列 x(n)：人体心电图信号在测量过程中往往受到工业高频干扰，所以必须经过低通滤波处理后，才能作为判断心脏功能的有用信息。
用双线性变换法设计IIR
实验三 用双线性变换法设计IIR 滤波器一、实验目的1、了解两种工程上最常用的变换方法：脉冲响应不变法和双线性变换法。
2、掌握双线性变换法设计IIR 滤波器的原理及具体设计方法，熟悉双线性法设计低通、带通和高通IIR 滤波器的计算机程序。
3、观察用双线性变换法设计的滤波器的品与特性，并与脉冲响应不变法相比较，了解双线性变换法的特点。
4、熟悉用双线性变换法涉及数字Butterworth 和Chebyshev 滤波器的全过程。
5、了解多项式呈几何多项式乘方运算的计算机变成方法。
二、实验原理与方法从模拟滤波器设计IIR 数字滤波器具有四种方法：微分—差分变换法、脉冲响应不变法、双线性变换法、匹配Z 变换法：在工程上常用的是其中两种：脉冲响应不变法、双线性变换法。
脉冲响应不变法需要经历如下步骤：由已知系统传输函数)(s H 计算系统冲激响应)(t h ：对)(t h 进行等间隔取样得到)()(nT h t h =,由)(t h 获得数字滤波器的系统响应)(z H 。
这种方法非常直观，其算法宗旨是保证所设计的IIR 滤波器的脉冲响应和响应模拟滤波器的冲激响应在采样个点上完全一致。
而双线性变换法的设计准则是使数字滤波器的频率响应与参考模拟滤波器的频率响应相似。
脉冲响应不变法一个重要的特点是频率坐标的变换式现行的，其缺点是有频谱的周期延拓效应，存在频谱混淆的现象。
为了克服脉冲响应不变法可能产生的频谱混淆，提出了双线性变换法，它依靠双线性变换式：ss z z z s -+=+-=--11,1111 其中，jwre z j s =Ω+=,σ 建立起s 平面和z 平面的单值映射关系，数字频域和模拟频域之间的关系：Ω==Ωarctg w w tg 2),2( （3—1）由上面的关系式可知,当∞>-Ω时，w 中止在折叠频率π=w 处，整个Ωj 周单值的对应于单位圆的一周。
因此双线性变换法不同于脉冲响应不变法，不存在频谱混淆的问题。
实验六 用双线性变换法设计IIR数字滤波器
实验六用双线性变换法设计IIR数字滤波器一、实验目的学会运用MATLAB设计数字低通、带通、高通、带阻滤波器的设计方法。
二、实验涉及的matlab子函数bilinear功能：双线性变换——将s域映射到z域。
调用格式：[numd,dend]= bilinear (num,den,Fs)，将模拟域系统函数转换为数字域的系统函数，Fs为采样频率。
三、实验原理下面举例说明用双线性变换法设计各种数字滤波器的过程。
例1、采用双线性变换法设计一个巴特沃斯数字低通滤波器，要求：wp=0.25*pi，rp=1db，ws=0.4*pi，as=15db，滤波器采样频率Fs=100hz。
MATLAB源程序为：%数字滤波器指标wpd=0.25*pi; %滤波器的通带截止频率wsd=0.4*pi; %滤波器的阻带截止频率Rp=1;As=15; %输入滤波器的通阻带衰减指标%转换为模拟原型滤波器指标Fs=100;T=1/Fs;wp=(2/T)*tan(wpd/2);ws=(2/T)*tan(wsd/2);%模拟原型滤波器计算[n,wc]=buttord(wp,ws,Rp,As,'s') %计算阶数n和截止频率[z0,p0,k0]=buttap(n); %归一化切比雪夫1型原型设计ba=k0*poly(z0); %求原型滤波器系数baa=poly(p0); %求原型滤波器系数a[ba1,aa1]=lp2lp(ba,aa,wc);%变换为模拟低通滤波器%用双线性变换法计算数字滤波器系数[bd,ad]=bilinear(ba1,aa1,Fs) %双线性变换%求数字系统的频率特性[H,w]=freqz(bd,ad);dbH=20*log10(abs(H)/max(abs(H))); %化为分贝值subplot(2,2,1),plot(w,abs(H));ylabel('|H|');title('幅度响应');axis([0,pi,0,1.1]);grid subplot(2,2,2),plot(w,angle(H));ylabel('\phi');title('相位响应');axis([0,pi,-4,4]);grid subplot(2,2,3),plot(w,dbH);title('幅度响应(dB)'); ylabel('dB');xlabel('频率');axis([0,pi,-40,5]);grid subplot(2,2,4),zplane(bd,ad); axis([-1.1,1.1,-1.1,1.1]);title('零极图');运行结果为： n = 5wc = 103.2016bd = 0.0072 0.0362 0.0725 0.0725 0.0362 0.0072 ad = 1.0000 -1.9434 1.9680 -1.0702 0.3166 -0.0392 则所求滤波器的系统函数为54321-5432-10.0392z 166z 3.01.0702z z 68.911.9434z 1z 072.00z 362.00z 725.00z 725.000.0362z 0.0072H(z)---------+-+-+++++=|H |幅度响应φ相位响应幅度响应(dB)d B频率-1-0.500.51Real PartI m a g i n a r y P a r t零极图例2、采用双线性变换法设计一个椭圆数字高通滤波器，要求通带250hz ，1db ，阻带150hz ，20db ，滤波器采样频率为Fs=1000hz 。
实验五：双线性变换法的设计IIR-数字滤波器
实验五、双线性变换法设计IIR 数字滤波器一、实验目的：1、熟悉用双线性变换法设计IIR 数字滤波器的原理与方法。
2、掌握数字滤波器的计算机仿真方法。
3、熟悉Batterworth 滤波器设计方法及特点 二、实验原理（一）、IIR 数字滤波器的设计步骤：① 按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标； ② 根据模拟滤波器技术指标设计为响应的模拟低通滤波器；③ 跟据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器；④ 如果要设计的滤波器是高通、带通或带阻滤波器，则首先把它们的技术指标转化为模拟低通滤波器的技术指标，设计为数字低通滤波器，最后通过频率转换的方法来得到所要的滤波器。
在MATLAB 中，经典法设计IIR 数字滤波器主要采用以下步骤：IIR 数字滤波器设计步骤（二）、用模拟滤波器设计数字滤波器的方法 1、冲激响应不变法：冲激响应不变法是从时域出发，要求数字滤波器的冲激响应h (n ) 对应于模拟滤波器h (t ) 的等间隔抽样。
优点：时域逼近良好；保持线性关系。
缺点：频域响应混叠。
只适用于限带低通滤波器和带通滤波器2、双线性变换法优点：克服了频域混叠模拟滤波器原型 buttap,cheb1ap频率变换 模拟离散化 bilinear,impin varIIR 数字滤波器/Tπ/T π-3/Tπ3/Tπ-j ΩσjIm (z)Re(z)1S 平面Z 平面1S ~S T Tππ-将整个平面压缩变换到平面一个的带状区域缺点：高频时会引起畸变1）冲激响应不变法impinvar格式：[BZ,AZ]= impinvar （B,A,Fs ）功能：把具有[B,A]模拟滤波器传递函数模型转换为采样频率为Fs 的数字滤波器的传递函数模型[BZ,AZ]，Fs 默认值为1。
例：一个4阶的Butterworth 模拟低通滤波器的系统函数如下：12251)(234++++=s s s s s H a试用冲激响应不变法求出Butterworth 模拟低通数字滤波器的系统函数。
IIR数字滤波器设计及实现
实验三IIR数字滤波器设计及实现一、实验目的(1)熟悉用脉冲响应不变法和双线性变换法设计IIR数字滤波器的原理与方法；(2)学会调用MATLAB信号处理工具箱中滤波器设计函数设计IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。
二、实验原理设计IIR数字滤波器一般采用脉冲响应不变法和双线性变换法。
脉冲响应不变法：根据设计指标求出滤波器确定最小阶数N和截止频率Wc；计算相应的模拟滤波器系统函数；将模拟滤波器系统函数:'转换成数字滤波器系统函数双线性变换法：根据数字低通技术指标得到滤波器的阶数N；取合适的T值，几遍校正计算相应模低通的技术指标--；根据阶数N查表的到归一化低通原型系统函数。
,将"' Q 代入。
‘去归一化得到实际的,/ ：' ；用双线性变换法将：’转换成数字滤波器三、实验内容及步骤1、用脉冲响应不变法设计(1)根据设计指标求出滤波器确定最小阶数N和截止频率Wcclear;close all;clc; % 开始准备fp=3400;fs=5000;Fs=22050;Rp=2;Rs=20;T=1/Fs; % T=1s 的模拟滤波器设计指标W1p=fp/Fs*2; W1s=fs/Fs*2; % 求归一化频率[N, Wn] = buttord(W1p, W1s, Rp, Rs, 's'; % 确定 butterworth 的最小阶数 N 和频率参数Wn 得到结果为:N 二7Wn 二 0.3266 即：该设计指标下的模拟滤波器最小阶数为N=7,其截至频率为Wn =0.3266；(2)计算相应的模拟滤波器系统函数打:, clear;close all;clc; % 开始准备fp=3400;fs=5000;Fs=22050;Rp=2;Rs=20;T=1/Fs; % T=1s 的模拟滤波器设计指标W1p=fp/Fs*2; W1s=fs/Fs*2; % 求归一化频率[N, Wn] = buttord(W1p, W1s, Rp, Rs, 's'; % 确定 butterworth 的最小阶数 N 和频率参数 Wn[B,A]=butter(N,1,'s' %计算相应的模拟滤波器系统函数得到结果为： B = 1.0e-003 * 0 00 0 0 0 0 0.3966 A =1.0000 1.4678 1.0773 0.5084 0.1661 0.0375 0.0055 0.0004 >>(3)将模拟滤波器系统函数转换成数字滤波器系统函数 clear;close all;clc; % 开始准备fp=3400;fs=5000;Fs=22050;Rp=2;Rs=20;T=1/Fs; % T=1s 的模拟滤波器设计指标W1p=fp/Fs*2; W1s=fs/Fs*2; % 求归一化频率[N, Wn] = buttord(W1p, W1s, Rp, Rs, 's'; % 确定 butterworth 的最小阶数 N 和频率参数Wn[B,A]=butter(N,1,'s' ； %计算相应的模拟滤波器系统函数 [Bz,Az]=impinvar(B,A %用脉冲相应不变法将模拟滤波器转换成数字滤波器 sys=tf(Bz,Az,T； %得到传输函数‘‘‘‘‘ Bz =1.0e-004 *-0.0000 0.0045 0.2045 0.8747 0.7094 0.1090 0.0016 0Az =1.0000 -5.5415 13.2850 -17.8428 14.4878 -7.1069 1.9491 -0.2304>>>>即：由Bz和Az可以写出数字滤波器系统函数为：Transfer function:-9.992e-015 z~7 + 4.454e-007 z~6 + 2.045e-005 z~5 + 8.747e-005 z~4 + 7.094e-005 z"3 + 1.09e-005 z~2+ 1.561e-007 z z 7 - 5.541 z 6 + 13.28 z 5 - 17.84 z 4 + 14.49 z 3 - 7.107 z 2 + 1.949 z - 0.2304Sampling time: 4.5351e-005>>(4)绘图clear;close all;clc; % 开始准备fp=3400;fs=5000;Fs=22050;Rp=2;Rs=20;T=1/Fs; % T=1s 的模拟滤波器设计指标W1p=fp/Fs*2; W1s=fs/Fs*2; % 求归一化频率[N, Wn] = buttord(W1p, W1s, Rp, Rs, 's'; % 确定butterworth 的最小阶数N 和频率参数Wn[B,A]=butter(N,Wn,'s'; %计算相应的模拟滤波器系统函数[Bz,Az]=impinvar(B,A; %用脉冲响应不变法将模拟滤波器转换成数字滤波器sys=tf(Bz,Az,T；%得到传输函数‘ [H,W]=freqz(Bz,Az,512,Fs; % 生成频率响应参数plot(W,20*log10(abs(H; % 绘制幅频响应grid on; %加坐标网格得到结果为：观察实验结果图可看到：在频率为3402Hz处频率为衰减2.015db,在频率为5017Hz处幅度衰减21.36db。
双线性变换法IIR数字滤波器设计
双线性变换法IIR 数字滤波器设计一、实验目的1、掌握用双线性变换法设计低通IIR 数字滤波器的基本原理和算法。
2、掌握用双线性变换法设计高通和带通IIR 数字滤波器的基本原理和算法。
3、进一步了解数字滤波器和模拟滤波器的频率响应特性。
二、实验原理与计算方法１、双线性变换法设计IIR 低通数字滤波器的基本原理和算法双线性变换法设计数字滤波器，采用了二次映射的方法，就是先将整个s 平面压缩到s 1平面的一个Tj Tj ππ~-的横形条带范围内，然后再将这个条带映射到z 平面上，就能建立s 平面到z 平面的一一对应关系。
对于低通数字滤波器，映射关系为zzT z z T s ++-=+-=--11211211 (1) 其中T 为抽样周期。
用双线性变换法设计低通IIR 数字滤波器的基本步骤，首先根据设计要求确定相应的模拟滤波器的传递函数)(s H a ，再应用(1)式得数字滤波器的传递函数)(z Hzz T s a s H z H ++-==112)()( (2)通常可以给定的参数为：低通数字滤波器通带边界频率p p f πΩ21=、阻带边界频率s s f πΩ21=和对应的通带衰减函数p α、阻带衰减函数s α。
s 1平面中的模拟角频率1Ω与数字角频率ω的关系为线性关系T 1Ωω=，在计算模拟滤波器的阶数N 、极点s i 和传递函数)(s H a 之前，应作预畸变处理22tan 22tan 21T f TT TπΩΩ== (3)模拟滤波器的阶数N 、极点s i 和传递函数)(s H a 的计算方法与冲激响应不变法相同，可以采用Butterworth 逼近或Chebyshev 逼近。
２、双线性变换法设计IIR 高通、带通、带阻数字滤波器的基本原理和算法由于双线性变换法获得的数字滤波器频率响应特性中不会出现混叠现象，因此可以适用于高通、带通和带阻滤波器的设计。
IIR 数字滤波器的设计通常要借助于模拟低通滤波器的设计，由原型低通滤波器到其他形式（高通、带通、带阻）IIR 数字滤波器的频带变换有模拟频带变换法和数字频带变换法。
试验三用双线性变换法设计IIR滤波器
实验三 用双线性变换法设计I ＩＲ滤波器一、实验目的１、了解两种工程上最常用的变换方法：脉冲响应不变法和双线性变换法。
2、掌握双线性变换法设计IIR 滤波器的原理及具体设计方法，熟悉双线性法设计低通、带通和高通II Ｒ滤波器的计算机程序。
3、观察用双线性变换法设计的滤波器的品与特性，并与脉冲响应不变法相比较,了解双线性变换法的特点。
4、熟悉用双线性变换法涉及数字Butterworth 和Chebyshev 滤波器的全过程。
5、了解多项式呈几何多项式乘方运算的计算机变成方法。
二、实验原理与方法从模拟滤波器设计IIR 数字滤波器具有四种方法:微分—差分变换法、脉冲响应不变法、双线性变换法、匹配Z 变换法：在工程上常用的是其中两种：脉冲响应不变法、双线性变换法。
脉冲响应不变法需要经历如下步骤：由已知系统传输函数)(s H 计算系统冲激响应)(t h :对)(t h 进行等间隔取样得到)()(nT h t h =，由)(t h 获得数字滤波器的系统响应)(z H 。
这种方法非常直观,其算法宗旨是保证所设计的IIR 滤波器的脉冲响应和响应模拟滤波器的冲激响应在采样个点上完全一致。
而双线性变换法的设计准则是使数字滤波器的频率响应与参考模拟滤波器的频率响应相似。
脉冲响应不变法一个重要的特点是频率坐标的变换式现行的，其缺点是有频谱的周期延拓效应,存在频谱混淆的现象。
为了克服脉冲响应不变法可能产生的频谱混淆，提出了双线性变换法，它依靠双线性变换式：ss z z z s -+=+-=--11,1111 其中,jw re z j s =Ω+=,σ 建立起s 平面和z 平面的单值映射关系，数字频域和模拟频域之间的关系：Ω==Ωarctg w w tg 2),2( (3—1）由上面的关系式可知,当∞>-Ω时,w 中止在折叠频率π=w 处，整个Ωj 周单值的对应于单位圆的一周。
因此双线性变换法不同于脉冲响应不变法，不存在频谱混淆的问题。
双线性变换法设计IIR数字滤波器
双线性变换法设计IIR数字滤波器一、实验目的1. 掌握利用双线性变换法设计IIR数字滤波器的原理和具体方法。
2. 深理解利用双线性变换法设计IIR数字滤波器的优缺点及使用范围。
二、实验内容1. 利用巴特沃思模拟滤波器,通过双线性变换法设计巴特沃思数字滤波器,数字滤波器的技术指标为通带临界频率fp=0.2πrad ，通带内衰减小于rp=1dB；阻带临界频率fs=0.3πrad，阻带内衰减大于rs=10dB；采样频率为FS=π/64.程序代码T=64/pi; %设置采样周期为π/64fs=1/T; %采样频率为周期倒数wpz=0.2;wsz=0.3;%设置归一化通带和阻带截止频率Wp=2*tan(wpz*pi/2);Ws=2*tan(wsz*pi/2);Ap=1;As=10; %设置通带最大和最小衰减[N,Wc]=buttord(Wp,Ws,Ap,As,'s'); %调用butter函数确定巴特沃斯滤波器阶数[B,A]=butter(N,Wc,'s'); %调用butter函数设计巴特沃斯滤波器W=linspace(0,pi,400*pi); %指定一段频率值[D,C]=bilinear (B,A,fs); %调用双线性变换法[Hz,W]=freqz(D,C,W);%返回频率响应sys=tf(D,C,T) % 给出传输函数H（z）subplot(3,1,1);plot(W,abs(Hz)/abs(Hz(1))); %绘出巴特沃斯数字低通滤波器的幅频特性曲线grid on;title('巴特沃斯数字滤波器1');xlabel('Frequency/Hz');ylabel('Magnitude');subplot(3,1,2);plot(W,20*log10(abs(Hz)));grid on;title('巴特沃斯数字滤波器1');xlabel('Frequency/Hz');ylabel('dB');[Nd,wdc]=buttord(wpz,wsz,Ap,As);[Bdz,Adz]=butter(Nd,wdc);[Hzz,W]=freqz(Bdz,Adz,W);subplot(3,1,3);plot(W,abs(Hzz)/abs(Hzz(1)));grid on;title('巴特沃斯数字滤波器2');xlabel('Frequency/Hz');ylabel('Magnitude');总结与体会在课设之前，我对滤波器设计中的函数基本上处于一种模糊状态。
实验二 IIR数字滤波器设计
实验二 IIR 数字滤波器设计一．实验目的1.掌握双线性变换法设计IIR 数字滤波器的原理及具体设计方法，熟悉用双线性变换法设计低通、带通和高通IIR 数字滤波器的计算机编程。
2.观察用双线性变换法设计的数字滤波器的频域特性，了解双线性变换法的特点。
3.熟悉用双线性变换法设计数字Butterworth 和Chebyshev 滤波器的全过程。
4. 通过观察对实际心电图信号的滤波作用，获得数字滤波工程应用的认识。
二.实验原理与方法1. IIR 数字滤波器可以借助于模拟滤波器设计，即先设计一个适于技术要求的原型模拟滤波器，再按一定的准则用映射的方法将模拟原型的传递函数Ha(s)变换为数字滤波器的系统函数H(z)，从而完成数字滤波器的设计任务。
这是一类简单而有效的方法，因为模拟滤波器理论已经相当成熟，有大量公式图表可以利用。
2. 双线性变换法的设计准则是使数字滤波器的频率响应与参考模拟滤波器的频率响应相似。
由双线性变换式 1111z s z ---=+ 建立s 平面与z 平面的单值映射关系，频率变换关系为()2tg ωΩ=。
s 平面的频率轴j Ω单值对应于z 平面上的单位圆j z e ω=，因此不存在频率混叠问题。
由于Ω与ω间的非线性关系，使各个临界频率位置发生非线性畸变，可以通过预畸变校正。
用双线性变换法设计数字滤波器时，先将数字滤波器的各临界频率经过频率预畸变求得模拟原型滤波器的各临界频率，设计模拟原型传递函数，通过双线性变换，正好将这些频率点映射到所需位置上。
双线性变换法设计数字低通滤波器步骤如下：(1)确定数字滤波器的性能指标，包括：通带、阻带临界频率，通带内最大衰减，阻带内最小衰减，采样周期 T 。
(2)确定相应的数字频率。
(3)计算预畸的模拟低通原型临界频率。
(4)计算低通原型阶数N 和3dB 频率ΩC ，求得传递函数Ha(s)。
(5)用低通变换公式1111z s z---=+代入Ha(s)，求得数字滤波器系统函数H(z)。
(完整word版)用双线性变换法设计IIR数字低通滤波器课程设计
V=课程设计报告书姓名:班级：学号:时间:设计题目用双线性变换法设计IIR数字低通滤波器设计要求1. 通过实验加深对双线性变换法设计IIR滤波器基本方法的了解.2. 了解MATLAB有关双线性变换法的子函数。
3．掌握用双线性变换法设计数字低通滤波器的方法。
本次课程设计是采用双线性变换法基于MATLAB设计一个IIR数字低通滤波器, 其中要求通带截止频率为ωp=0.25π;通带最大衰减Rp=1dB;阻带最小衰减As=15dB;阻带截止频率ωs=0.4π;滤波器采样频率Fs=100Hz.设计过程摘要: 根据IIR滤波器的特点, 在MATLAB坏境下用双线性变换法设计IIR数字滤波器。
利用MATLAB设计滤波器, 可以随时对比设计要求和滤波器特性调整参数, 直观简便, 极大的减轻了工作量, 有利于滤波器设计的最优化。
1.关键词：双线性变换法 , 数字滤波器 , MATLAB , IIR2.设计原理与步骤1.1设计原理滤波器的种类很多, 从功能上可分为低通、高通、带通和带阻滤波器, 每一种又有模拟滤波器和数字滤波器两种形式。
如果滤波器的输人和输出都是离散时间信号, 则该滤波器的冲击响应也必然是离散的, 这种滤波器称之为数字滤波器。
数字滤波器是一种用来过滤时间离散信号的数字系统, 通过对抽样数据进行数学处理来达到频域滤波的目的。
数字滤波器也是具有一定传输选择特性的数字信号处理装置, 其输入、输出均为数字信号, 实质上是一个由有限精度算法实现的线性时不变离散系统。
IIR数字滤波器采用递归型结构, 即结构上带有反馈环路。
IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成, 可以组合成直接型、正准型、级联型、并联型四种结构形式, 都具有反馈回路。
数字滤波器根据其冲激响应函数的时域特性, 可分为两种, 即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。
IIR 数字滤波器的特征是, 具有无限持续时间冲激响应, 需要用递归模型来实现, 其差分方程为：（1-1）（1-2）设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H(z), 使其频率响应H(z)满足所希望得到的频域指标, 即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。
数字信号处理实验三 ——用双线性变换法设计IIR数字滤波器
实验内容及步骤
5、设计一个CB型带通DF，满足：通带边界频率为100Hz～500Hz， 通带衰耗小于3dB，过渡带宽20Hz，阻带衰耗大于15dB，取样频率 fs=2000Hz。写出其模拟滤波器方程及数字滤波器方程。
6、人体心电图信号在测量过程中往往受到工业高频干扰，所以必须 经过低通滤波处理后，才能作为判断心脏功能的有用信息。下面程序 中给出一实际心电图信号采样序列样本x(n)，其中存在高频干扰。在 实验中，以x(n)作为输入序列，滤除其中的干扰成分。我们用双线性 变换法设计一个巴特沃斯低通IIR数字滤波器对上述信号进行处理。 直接运行程序，结果输出滤波器幅频特性曲线图，有噪声的心电图采 集信号波形图和经过三级二阶滤波器后的心电图信号波形图，可以看 出低通滤波器除信号中高频噪声的滤波效果。
1 T j T
z
2
1 T
j
2 T
1
2
2
1 j T
z
2 1
1 j T
2
1 T j T
z
2
1 T
j
2 T
1
2
2
j
j Im[z]
-1
1 Re[z]
0
0
s平面
z平面
（2）整个左半平面映射成单 位圆的内部，因此这是一个 稳定的变换。 （3）虚轴以一一对应的方式 映射成单位圆，因此在频域 中不会出现混叠。
[bz,az]=bilinear(bt,at,FS); %双线性变换为数字滤波器
[H,W] = freqz(bz,az); %求解数字滤波[器z,的p,频k率]=响b应uttap(n); plot(W*FS/(2*pi),abs(H));grid; [b,a] = zp2tf(z,p,k); % 零极点增益模型到传递 xlabel('频率/Hz');ylabel('幅值'); 函数模型的转换
实验五双线性变换法设计IIR数字滤波器
实验五IIR数字滤波器设计一、实验目的（1）熟悉用双线性变换法设计IIR数字滤波器的原理与方法；（2）学会调用MATLAB信号处理工具箱中滤波器设计函数设计各种IIR数字滤波器，学会根据滤波需求确定滤波器指标参数。
（3）掌握IIR数字滤波器的MATLAB实现方法。
（4）通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。
二、实验原理设计IIR数字滤波器一般采用间接法（脉冲响应不变法和双线性变换法），应用最广泛的是双线性变换法。
基本设计过程是：①先将给定的数字滤波器的指标转换成过渡模拟滤波器的指标；②设计过渡模拟滤波器；③将过渡模拟滤波器系统函数转换成数字滤波器的系统函数。
MATLAB信号处理工具箱中的各种IIR 数字滤波器设计函数都是采用双线性变换法。
第六章介绍的滤波器设计函数butter、cheby1 、cheby2 和ellip可以分别被调用来直接设计巴特沃斯、切比雪夫1、切比雪夫2和椭圆模拟和数字滤波器。
本实验要求读者调用如上函数直接设计IIR数字滤波器。
本实验的数字滤波器的MATLAB实现是指调用MATLAB信号处理工具箱函数filter对给定的输入信号x(n)进行滤波，得到滤波后的输出信号y(n）。
三、实验内容及步骤（1）调用信号产生函数mstg产生由三路抑制载波调幅信号相加构成的复合信号st，该函数还会自动绘图显示st的时域波形和幅频特性曲线，如图5.1所示。
容易看出，图5.1中三路调幅信号的载波频率分别为250Hz 、500Hz 、1000Hz 。
抑制载波单频调幅信号的数学表示式为0001()cos(2)cos(2)[cos(2())cos(2())]2c c c s t f t f t f f t f f t ππππ==-++ 其中，cos(2)c f t π称为载波，f c 为载波频率，0cos(2)f t π称为单频调制信号，f 0为调制正弦波信号频率，且满足0c f f >。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
双线性变换法IIR 数字滤波器设计一、实验目的1、掌握用双线性变换法设计低通IIR 数字滤波器的基本原理和算法。
2、掌握用双线性变换法设计高通和带通IIR 数字滤波器的基本原理和算法。
3、进一步了解数字滤波器和模拟滤波器的频率响应特性。
二、实验原理与计算方法１、双线性变换法设计IIR 低通数字滤波器的基本原理和算法双线性变换法设计数字滤波器，采用了二次映射的方法，就是先将整个s 平面压缩到s 1平面的一个Tj Tj ππ~-的横形条带范围内，然后再将这个条带映射到z 平面上，就能建立s 平面到z 平面的一一对应关系。
对于低通数字滤波器，映射关系为zzT z z T s ++-=+-=--11211211 (1) 其中T 为抽样周期。
用双线性变换法设计低通IIR 数字滤波器的基本步骤，首先根据设计要求确定相应的模拟滤波器的传递函数)(s H a ，再应用(1)式得数字滤波器的传递函数)(z Hzz T s a s H z H ++-==112)()( (2)通常可以给定的参数为：低通数字滤波器通带边界频率p p f πΩ21=、阻带边界频率s s f πΩ21=和对应的通带衰减函数p α、阻带衰减函数s α。
s 1平面中的模拟角频率1Ω与数字角频率ω的关系为线性关系T 1Ωω=，在计算模拟滤波器的阶数N 、极点s i 和传递函数)(s H a 之前，应作预畸变处理22tan 22tan 21T f TT TπΩΩ== (3)模拟滤波器的阶数N 、极点s i 和传递函数)(s H a 的计算方法与冲激响应不变法相同，可以采用Butterworth 逼近或Chebyshev 逼近。
２、双线性变换法设计IIR 高通、带通、带阻数字滤波器的基本原理和算法由于双线性变换法获得的数字滤波器频率响应特性中不会出现混叠现象，因此可以适用于高通、带通和带阻滤波器的设计。
IIR 数字滤波器的设计通常要借助于模拟低通滤波器的设计，由原型低通滤波器到其他形式（高通、带通、带阻）IIR 数字滤波器的频带变换有模拟频带变换法和数字频带变换法。
（1）模拟频带变换法首先将给定的对数字滤波器(DF)的技术要求转换为一个低通模拟滤波器(AF)的技术要求，根据这种要求用某种逼近设计出原型的低通模拟滤波器(LP AF)，计算出模拟滤波器的阶数N 、极点s i 和传递函数)(s H a ，再按照双线性变换的变换关系，将模拟滤波器的传递函数)(s H a 转换为数字滤波器的传递函数)(z H 。
表8-1中列出了将给定的对数字滤波器(DF)的技术要求直接转换为对一个低通模拟滤波器(AF)的技术要求的频率预畸变校正关系和转换公式。
例：数字高通滤波器的设计首先将给定的数字高通滤波器的技术指标根据公式转换为模拟低通滤波器的技术指标，利用cheb1ord(Wp,Ws,ap,as,'s')函数求出chebyshev模拟低通滤波器的阶数N，再利用cheb1ap(N,ap)函数求出模拟低通滤波器系统函数)(sHa的零极点，zp2tf(z,p,k)函数将零极点转换为系统函数系数；然后利用lp2hp由模拟低通滤波器的系统函数得到模拟带通滤波器的系统函数，bilinear函数则用于实现双线性变换法由模拟带通滤波器系统函数)(sHa 计算数字数字带通滤波器系统函数)(zH的系数。
（2）数字频带变换法首先将给定的对数字滤波器(DF)的技术要求转换为一个低通模拟滤波器(AF)的技术要求，用双线性变换法将原型的低通模拟滤波器(LP AF)映射为低通数字滤波器，再将数字低通滤波器根据相应的变换公式经频带变换到各型数字滤波器。
例：数字高通滤波器的设计函数[bhp,ahp]=zmapping(blp,alp,Nz,Dz)用来实现从数字低通滤波器得到数字高通滤波器的有理函数。
%数字滤波器技术指标>>wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;%对应的模拟滤波器技术指标>>T=1;Fs=1/T;Wp=(2/T)*tan(wp/2); =(2/T)*tan(ws/2);>>[cs,cd]=afd_chb1(Wp,Ws,Rp,As); %Chebyshev模拟滤波器>> [blp,alp]=bilinear(cs,cd,Fs) %双线性变换>>wphp=0.6*pi; %数字高通滤波器截止频率%低通－高通频带变换>>alpha=-(cos((wplp+wphp)/2))/(cos((wplp+wphp)/2))>>Nz=-[alpha,1];Dz=[1,alpha];>> [bhp,ahp]=zmapping(blp,alp,Nz,Dz) %数字高通滤波器的系统函数系数（3）IIR数字滤波器的设计可利用MATLAB提供的函数直接设计相应的数字滤波器。
函数buttord 和cheb1ord 用来根据给定的技术指标求出滤波器的阶数N 和边界频率wn ，butter 和cheby1则根据阶数和边界频率设计相应的数字滤波器。
输入的参数不同则所设计的滤波器类型不同。
[N,wn]=buttord(wp,ws,Rp,As); [N,wn]=cheb1ord(wp,ws,Rp,As); [b,a]=butter(N,wn); [b,a]=cheby1(N,Rp,wp);三、实验内容（1）Chebyshev IIR 数字带通滤波器满足如下技术指标低阻带边界频率Hz 1001=s f ，高阻带边界频率Hz 6001=s f ，阻带衰减函数dB 18=s α 低通带边界频率Hz 2001=p f ，高通带边界频率Hz 4002=p f ，通带波动dB 2=δ 抽样频率Hz f 2000=,记录所得的模拟滤波器的阶数N ，画出模拟滤波器和数字滤波器的频率响应的幅频和相频特性曲线。
A. 实验代码：fs1=100; fs2=600; fp1=200; fp2=400; fsa=2000; As=18; Rp=2;T=1./fsa %对应的模拟滤波器技术指 w1=2.*pi.*(fp1./fsa) %Chebyshev 模拟滤波器 w2=2.*pi.*(fp2./fsa)wp1=2*pi*fp1*T;wp2=2*pi*fp2*T ws2=2.*pi.*(fs2./fsa)cosw0=(sin(w1+w2))./(sin(w2)+sin(w1)) w0=acos(cosw0);bw=wp2-wp1 Wp=(cosw0-cos(w2))./sin(w2) Ws=(cosw0-cos(ws2))./sin(ws2)[N,omgn]=cheb1ord(Wp,Ws,Rp,As,'s') %返回模拟低通滤波器阶数N 和边界频率 [z,p,k]=cheb1ap(N,Rp) %得系统函数零极点 [blp,alp]=zp2tf(z,p,k); %由零极点得系数[bhp,ahp]=lp2bp(blp,alp,w0,bw) %模拟低通到模拟带通[bdf,adf]=bilinear(bhp,ahp,1) %双线性变换将模拟带通滤波器转换成数字带通滤波器[BPA,wa]=freqs(bhp,alp,fsa); [BPD,wd]=freqz(bdf,adf,fsa); subplot(2,2,1); plot(abs(BPA));title('模拟带通滤波器幅频特性') subplot(2,2,2); plot(angle(BPA));title('模拟带通滤波器相频特性') subplot(2,2,3); plot(abs(BPD));title('数字带通滤波器幅频特性') subplot(2,2,4); plot(angle(BPD));title('数字带通滤波器相频特性')实验截图：（2）Chebyshev IIR 数字带阻滤波器满足如下技术指标低阻带边界频率Hz 10001=s f ，高阻带边界频率Hz 20001=s f ，阻带衰减函数dB 20=s α 低通带边界频率Hz 5001=p f ，高通带边界频率Hz 30002=p f ，通带波动dB 3=δ 抽样频率kHz f 10=,记录所得的模拟滤波器的阶数N ，画出模拟滤波器和数字滤波器的频率响应的幅频和相频特性曲线。
实验代码：fs1=1000; fs2=2000; fp1=500; fp2=3000; fsa=10000;As=18;Rp=2;T=1./fsa %对应的模拟滤波器技w1=2.*pi.*(fp1./fsa);%Chebyshev模拟滤波器w2=2.*pi.*(fp2./fsa)wp1=2*pi*fp1*T;wp2=2*pi*fp2*Tws2=2.*pi.*(fs2./fsa)cosw0=(sin(w1+w2))./(sin(w2)+sin(w1))w0=acos(cosw0);bw=wp2-wp1Wp=(cosw0-cos(w2))./sin(w2)Ws=(cosw0-cos(ws2))./sin(ws2)[N,omgn]=cheb1ord(Wp,Ws,Rp,As,'s') %返回模拟低通滤波器阶数N和边界频率?n [z,p,k]=cheb1ap(N,Rp) %得系统函数零极点[blp,alp]=zp2tf(z,p,k) %由零极点得系数[bhp,ahp]=lp2bs(blp,alp,w0,bw) %模拟低通到模拟带通[bdf,adf]=bilinear(bhp,ahp,1) %双线性变换将模拟带通滤波器转换成数字带通滤波器[BPA,wa]=freqs(bhp,alp,fsa)[BPD,wd]=freqz(bdf,adf,fsa)subplot(2,2,1);plot(abs(BPA));title('模拟带阻滤波器幅频特性')subplot(2,2,2);plot(angle(BPA));title('模拟带阻滤波器相频特性')subplot(2,2,3);plot(abs([BPD]));title('数字带阻滤波器幅频特性')subplot(2,2,4);plot(angle([BPD]));title('数字带阻滤波器相频特性')实验截图：Welcome !!! 欢迎您的下载，资料仅供参考！。