《通信系统的计算机模拟第六讲精品文稿.ppt》由会员分享，可在线阅读，更多相关《通信系统的计算机模拟第六讲精品文稿.ppt（57页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、通信系通信系统的的计算机模算机模拟第六第六讲1第1页，本讲稿共57页第五讲第五讲 回顾回顾带通信号低通仿真模型带通信号低通仿真模型l复包络：频域，同相和正交复包络：频域，同相和正交分量的谱分量的谱l能量与功率能量与功率l随机带通信号的正交模型随机带通信号的正交模型l随机复包络信号的若干结论随机复包络信号的若干结论l信噪比信噪比与系统的低通仿真模型与系统的低通仿真模型l线性时不变系统线性时不变系统l多载波信号多载波信号l非线性与时变系统非线性与时变系统2第2页，本讲稿共57页第五章第五章 滤波器模型与仿真方法滤波器模型与仿真方法l滤波器是组成通信系统的许多子系统中的重要部分滤波器是组成通信系统的
2、许多子系统中的重要部分l这些滤波器中的许多是这些滤波器中的许多是模拟模拟的，为了便于仿真，必须将它们的，为了便于仿真，必须将它们映射映射为合适为合适的等价数字滤波器的等价数字滤波器近似并且都会在仿真结果中引入误差近似并且都会在仿真结果中引入误差l根据定义，滤波器具有频率选择性，并且具有有限时间冲激响应或根据定义，滤波器具有频率选择性，并且具有有限时间冲激响应或者无限时间冲激响应。者无限时间冲激响应。l由于滤波器具有频率选择性，它们引入了记忆由于滤波器具有频率选择性，它们引入了记忆?，而由于这种记忆，需要滤，而由于这种记忆，需要滤波器过去的输入和波器过去的输入和/或输出来计算给定时刻的滤波器输出
3、。或输出来计算给定时刻的滤波器输出。l因此，滤波器需要存储，而存储和恢复采样值会大大地增加仿真程因此，滤波器需要存储，而存储和恢复采样值会大大地增加仿真程序的计算量，反过来，这又增加了仿真程序的运行时间，因此，我序的计算量，反过来，这又增加了仿真程序的运行时间，因此，我们寻求们寻求较好的仿真结构较好的仿真结构以获得减少计算量的算法。以获得减少计算量的算法。l不是滤波器设计不是滤波器设计3第3页，本讲稿共57页5.1 概述概述l数字滤波器的经典综合方法通常都是基于模拟滤波器原型之上的数字滤波器的经典综合方法通常都是基于模拟滤波器原型之上的l仿真模型所要求的数字滤波器的设计，仿真模型所要求的数字滤
4、波器的设计，从模拟滤波器的从模拟滤波器的S（拉普拉斯）域传（拉普拉斯）域传递函数着手，来寻求递函数着手，来寻求等价的数字滤波器等价的数字滤波器等价的数字滤波器等价的数字滤波器，等价于模拟原型的数字滤波器。等价于模拟原型的数字滤波器。l基本方法是基本方法是基于时域判据或频域判据基于时域判据或频域判据基于时域判据或频域判据基于时域判据或频域判据来衡量其等价性。来衡量其等价性。l通过要求通过要求数字滤波器的输出数字滤波器的输出数字滤波器的输出数字滤波器的输出与作为设计基础的与作为设计基础的模拟原型的采样输出相匹模拟原型的采样输出相匹模拟原型的采样输出相匹模拟原型的采样输出相匹配配配配，来，来运用时域
5、判据进行滤波器的设计运用时域判据进行滤波器的设计运用时域判据进行滤波器的设计运用时域判据进行滤波器的设计。l“冲激不变法冲激不变法”是使是使数字滤波器的冲激响应数字滤波器的冲激响应与与模拟原型的冲激响应的采模拟原型的冲激响应的采样相匹配样相匹配的一种设计方法的一种设计方法l对于阶跃不变数字滤波器，对于阶跃不变数字滤波器，数字滤波器的阶跃响应数字滤波器的阶跃响应则与对应的模拟则与对应的模拟原型滤波原型滤波器阶跃响应的采样相匹配器阶跃响应的采样相匹配。l模拟原型和得到的数字滤波器在模拟原型和得到的数字滤波器在时域上是等价的时域上是等价的，它们在频域上也它们在频域上也近似近似近似近似等价等价。（采样
6、频率。（采样频率混叠）混叠）4第4页，本讲稿共57页双线性双线性z变换变换l最常用方法：最常用方法：双线性双线性z变换：变换：l模拟原型模拟原型（双线性双线性z变换）变换）数字等价模型数字等价模型l它允许对频率响应上的点进行匹配，使得模拟原型和它允许对频率响应上的点进行匹配，使得模拟原型和所得数字滤波器在某些指定频率值上具有相同的频率所得数字滤波器在某些指定频率值上具有相同的频率响应（包括幅度和相位）。响应（包括幅度和相位）。l双线性变换法以引入双线性变换法以引入非线性频率畸变非线性频率畸变（Warping）为）为代价消除了混叠误差。代价消除了混叠误差。5第5页，本讲稿共57页IIR and
7、FIRl基于模拟原型的综合方法给出无限时间冲激响应（基于模拟原型的综合方法给出无限时间冲激响应（infinite-duration impulse response,IIR）数字滤波器。）数字滤波器。l原因：模拟滤波器的冲激响应会渐近趋于零，严格地讲，这种冲激响应在持续原因：模拟滤波器的冲激响应会渐近趋于零，严格地讲，这种冲激响应在持续时间上是无限长的。时间上是无限长的。l为获得用于仿真程序的有限时间冲激响应（为获得用于仿真程序的有限时间冲激响应（finite-duration impulse response,FIR）数字滤波器，由标准）数字滤波器，由标准IIR设计方法所产生的冲激响应可能设
8、计方法所产生的冲激响应可能需要需要截短（截短（truncate）。l通过在通过在FIR滤波器的冲激响应中滤波器的冲激响应中包含足够多的项数包含足够多的项数，由截短引起的误差，由截短引起的误差可以降低到一个允许的水平下。可以降低到一个允许的水平下。6第6页，本讲稿共57页Pure数字滤波器数字滤波器l数字滤波器的一个重要特性是，在没有对应的模拟滤波器时也可设计数字滤波器的一个重要特性是，在没有对应的模拟滤波器时也可设计数字滤波器。数字滤波器。l允许近似给定的幅度响应，但要保持一个完全线性的相位响应，这种滤波允许近似给定的幅度响应，但要保持一个完全线性的相位响应，这种滤波器为器为FIR滤波器，可通
9、过横向延时线（滤波器，可通过横向延时线（Transversal Delay Line,TDL）结构实现。结构实现。l这种滤波器有多种设计方法，这种滤波器有多种设计方法，最基本的方法是将最基本的方法是将以采样频率为周期以采样频率为周期的的期望幅度响应期望幅度响应展开展开成为傅里叶级数，所得傅里叶级数的系数就确定了数成为傅里叶级数，所得傅里叶级数的系数就确定了数字滤波器的冲激响应。字滤波器的冲激响应。l这是频率采样的一个例子这是频率采样的一个例子因为在各频率点上对期望频率响应进行因为在各频率点上对期望频率响应进行“采样采样”，对这些频率采样点进行傅里叶反变换（，对这些频率采样点进行傅里叶反变换（I
10、FFT）之后，就获得）之后，就获得了滤波器的冲激响应。了滤波器的冲激响应。l将滤波器的输入信号与冲激响应进行卷积运算，实现了仿真滤波器将滤波器的输入信号与冲激响应进行卷积运算，实现了仿真滤波器模型并且获得了滤波器的输出信号。模型并且获得了滤波器的输出信号。7第7页，本讲稿共57页数字滤波器分类8第8页，本讲稿共57页 5.2 IIR与与FIR滤波器滤波器 l数字滤波器通常根据冲激响应的持续时间数字滤波器通常根据冲激响应的持续时间分为分为IIR滤波器和滤波器和FIR滤波器。滤波器。l和冲激响应分类紧密联系的是滤波器的实和冲激响应分类紧密联系的是滤波器的实现或结构。现或结构。9第9页，本讲稿共57
11、页5.2.1 IIR滤波器滤波器l数字滤波器：计算数字滤波器：计算过去的过去的N个输出个输出采样点采样点yn-k(1 k N）、当前的输入当前的输入采样点采样点xn，以及，以及过去的过去的N个输入个输入采样点采样点xn-k(1 k N)的加权和作为当前的输出采样点。的加权和作为当前的输出采样点。任何权值任何权值ak或或bk(k 0)具有非零值，处理器就会有记忆，因而具有频具有非零值，处理器就会有记忆，因而具有频率选择性，我们称之为率选择性，我们称之为“滤波器滤波器”。在时变系统中，其中的一个或多个权值还会是序数在时变系统中，其中的一个或多个权值还会是序数n的函数。对（的函数。对（5-1）式两边
12、进行）式两边进行Z变换，可以得到传递函数变换，可以得到传递函数H（z）。）。记住记住z变换是一种线性运算，和的变换是一种线性运算，和的z变换等于变换等于z变换的和，变换的和，k个采样周个采样周期的延时等价于乘以因子期的延时等价于乘以因子z-k。（5-1）10第10页，本讲稿共57页l由上式得传递函数为由上式得传递函数为此即线性时不变滤波器的传递函数的一般形式（5-2）（5-3）11第11页，本讲稿共57页滤波器的冲激响应以及频率响应滤波器的冲激响应以及频率响应l冲激响应（记为冲激响应（记为hn）是传递函数）是传递函数H（z）的逆变换，用的逆变换，用exp(j2fT)fT)替代传递函数中的替代传
13、递函数中的变量，就可得频率响应。变量，就可得频率响应。数字滤波器的冲激响应记作数字滤波器的冲激响应记作hn，只要令，只要令x（n)=(n)就可得到它。就可得到它。（5-4）（5-5）12第12页，本讲稿共57页Why IIRl由于式（由于式（5-1）的递归性）的递归性(yn是是yn-1的函数），冲激响应通常的函数），冲激响应通常是无限长的，是无限长的，IIR滤波器。冲激响应滤波器。冲激响应hn是序数是序数n的离散函数，的离散函数，而频率响应是连续变量而频率响应是连续变量f的连续函数。的连续函数。13第13页，本讲稿共57页 5.2.2 FIR滤波器滤波器l对于对于k1有有 ak=0,则得到则得
14、到FIR滤波器。滤波器。FIR数字滤波器的传递数字滤波器的传递函数如下：函数如下：其对应的冲激响应下如：仅在取值范围仅在取值范围0nN内，冲激响应具有非零值，最多有个内，冲激响应具有非零值，最多有个N+1非非零项，它是有限长的。零项，它是有限长的。（5-6）（5-7）（5-8）14第14页，本讲稿共57页5.2.3 综合与仿真综合与仿真 l在仿真程序中使用滤波器涉及两种截然不同的运算操作在仿真程序中使用滤波器涉及两种截然不同的运算操作l综合：给定了滤波要求，并确定了满足这些要求的滤波器的传递函综合：给定了滤波要求，并确定了满足这些要求的滤波器的传递函数数H（z），这就建立了仿真模型。），这就建
15、立了仿真模型。l综合运算所得结果通常表示成两个向量：一个包含分母系数，另一个则包含分综合运算所得结果通常表示成两个向量：一个包含分母系数，另一个则包含分子系数。这两个向量定义了传递函数式，以及给定滤波器的输入产生滤波器输子系数。这两个向量定义了传递函数式，以及给定滤波器的输入产生滤波器输出的算法。出的算法。l即使采用复数算法，综合操作的计算量通常也不大，因为即使采用复数算法，综合操作的计算量通常也不大，因为滤波器的综合运算只进滤波器的综合运算只进行一次，因而发生在仿真程序主循环的外部。行一次，因而发生在仿真程序主循环的外部。l第二种运算操作涉及第二种运算操作涉及在每一个仿真时间步，计算滤波器的
16、输出在每一个仿真时间步，计算滤波器的输出。这。这种运算在蒙特卡罗仿真程序中要重复进行数百万次，甚至数十亿次，如果想种运算在蒙特卡罗仿真程序中要重复进行数百万次，甚至数十亿次，如果想在合理的仿真运行时间内完成，要尽可能地减小这种运算的计算量。在合理的仿真运行时间内完成，要尽可能地减小这种运算的计算量。15第15页，本讲稿共57页5.3 IRR与与FIR滤波器的实现滤波器的实现 lTarget:最小化计算量最小化计算量仿真时间短仿真时间短5.3.1 直接直接II型与转置直接型与转置直接II型的实现型的实现lIIR数字滤波器的一种有效方法是采用转置直接数字滤波器的一种有效方法是采用转置直接型结构。型
17、结构。l转置直接转置直接型滤波器结构的信号流图，及推导出这种结构的直接型滤波器结构的信号流图，及推导出这种结构的直接型型结构的信号流图。结构的信号流图。16第16页，本讲稿共57页17第17页，本讲稿共57页产生转置直接产生转置直接型滤波器的规则型滤波器的规则l1重画原有（直接重画原有（直接型）信号流图，保持结构（所有链路的相对位置型）信号流图，保持结构（所有链路的相对位置不变）。不变）。l2将每次链路中的信号流反向。将每次链路中的信号流反向。l3对新链路赋上与原有链路相同的运算（乘以一个常数、时延等）。对新链路赋上与原有链路相同的运算（乘以一个常数、时延等）。l4如有需要，翻转（从左到右）新
18、构造的信号流图，使得输入如有需要，翻转（从左到右）新构造的信号流图，使得输入/输输出信号流的方向与原始信号流图一致（注：信号流的方向通常是从左出信号流的方向与原始信号流图一致（注：信号流的方向通常是从左向右）。向右）。l新获得的信号流图叫做转置直接新获得的信号流图叫做转置直接型（型（DF）结构，其传递函数与原始信）结构，其传递函数与原始信号流图相同。号流图相同。18第18页，本讲稿共57页转置直接转置直接型的优势型的优势l5-2所示的四阶转置直接所示的四阶转置直接型滤波器。（型滤波器。（可以很方便地扩展到高阶滤波器）可以很方便地扩展到高阶滤波器）。l给定输入采样，要计算输出信号给定输入采样，要
19、计算输出信号l第一步，计算状态变量。对于四阶滤波器，在我们的公式中有五个第一步，计算状态变量。对于四阶滤波器，在我们的公式中有五个状态变量。这五个状态变量的计算式如下：状态变量。这五个状态变量的计算式如下：（5-9）（5-10）（5-11）（5-12）（5-13）19第19页，本讲稿共57页直接对应着直接对应着Matlab程序程序20第20页，本讲稿共57页lx和和y分别代表滤波器的当前输入和当前输出。分别代表滤波器的当前输入和当前输出。l式（式（5-9）到式（）到式（5-13）采样步长的延时，）采样步长的延时，不需要存储和恢不需要存储和恢复复l这是这是MATLAB子程序子程序filter要基
20、于转置直接要基于转置直接型结构的原型结构的原因。因。l在首次进入仿真循环之前必须进行初始化，输出一个瞬态在首次进入仿真循环之前必须进行初始化，输出一个瞬态响应。这段时间通常称为响应。这段时间通常称为“稳定时间稳定时间”，其大小为滤波器，其大小为滤波器带宽倒数的若干倍。带宽倒数的若干倍。21第21页，本讲稿共57页矩阵形式矩阵形式 l数字滤波器的状态方程通常使用矩阵形式来表达数字滤波器的状态方程通常使用矩阵形式来表达 Wn 和和Wn-1 是的是的k1的列向量，分别表示当前及过去的状态变量。的列向量，分别表示当前及过去的状态变量。Fc和和Fd是是kk的系数矩阵，的系数矩阵，B是是k1的列向量，用来
21、将输入的列向量，用来将输入xn耦合耦合到状态变量。单输出信号到状态变量。单输出信号yn的输出方程为：的输出方程为：这里这里C是一个的是一个的1k行向量行向量（5-14）（5-15）22第22页，本讲稿共57页(5-16)23第23页，本讲稿共57页可算图可算图l在主对角线及其上方元素全部为零，具有这种属在主对角线及其上方元素全部为零，具有这种属性的信号流图称为性的信号流图称为“可算图可算图”：状态变量可以依：状态变量可以依次计算出来次计算出来l可算的充要条件：每一个闭合回路中的路径上至可算的充要条件：每一个闭合回路中的路径上至少含有一个延时单元，即少含有一个延时单元，即z1。24第24页，本讲
22、稿共57页Maltab生成生成l对于大量不同的模拟原型，可以轻易算出滤波器系数对于大量不同的模拟原型，可以轻易算出滤波器系数ak和和bk。l尽管生成的这些滤波器系数通常会用尽管生成的这些滤波器系数通常会用MATLAB子程序子程序filter，但是使用下列三行简单的但是使用下列三行简单的MATLAB程序，就可以逐个采样点程序，就可以逐个采样点地仿真地仿真IIR滤波器。滤波器。25第25页，本讲稿共57页Maltab生成生成l其中其中a和和b在仿真循环外确定在仿真循环外确定l参数参数sreg代表长度为代表长度为order+1的移位寄存器，其中的移位寄存器，其中order是滤波器阶数。是滤波器阶数。
23、l可使用可使用MATLAB滤波器综合子程序，如滤波器综合子程序，如butter（巴特沃思滤波器）、（巴特沃思滤波器）、cheby1（切比雪夫（切比雪夫1型滤波器）及型滤波器）及elliptic（椭圆滤波器），来计算滤波器系（椭圆滤波器），来计算滤波器系数向量数向量a和和b。l必须小心以确保存放分子系数和分母系数的向量保持相同的长度，必须小心以确保存放分子系数和分母系数的向量保持相同的长度，如果这两个向量长度不相等，较短的向量应该进行补零，使之与较如果这两个向量长度不相等，较短的向量应该进行补零，使之与较长向量的长度相等。长向量的长度相等。26第26页，本讲稿共57页例例5-1 来确定一个四阶巴
24、特沃思滤来确定一个四阶巴特沃思滤波器的冲激响应。波器的冲激响应。ln=40;%number of sampleslorder=4;%filter orderlb,a=butter(order,0.1?);%prototypel%The following program segment is the block processing implementation.lin1=1,zeros(1,n-1);%input vectorlout1=filter(b,a,in1);%output vectorl%The following program segment is the sample-by
25、-sample implementation.lsreg=zeros(1,order+1);%initialize shift registerlfor k=1:nl if k=1l in=1;%impulse inputl elsel in=0;lend27第27页，本讲稿共57页lout=b(1)*in+sreg(1,1);%determine outputl sreg=in*b-out*a+sreg;%update registerl sreg=sreg(1,2:(order+1),0;%shiftl out2(k)=out;%create output vectorlendlsubpl
26、ot(2,1,1)lindex=0:n-1;lstem(index,out1)lxlabel(Sample Index)lylabel(Block Processing)lsubplot(2,1,2)lstem(index,out2)lxlabel(Sample Index)lylabel(Serial Processing)28第28页，本讲稿共57页Result29第29页，本讲稿共57页 5.3.2 FIR滤波器的实现滤波器的实现 l令令ak0，可从图，可从图5-2直接得出直接得出FIR数字滤波器的实现策略。抽头延时线或数字滤波器的实现策略。抽头延时线或横向延时线（横向延时线（TDL）结
27、构，它们用带抽的延时线（级联式单元）结构，它们用带抽的延时线（级联式单元z1）来实）来实现，其中的抽头乘上了权值。现，其中的抽头乘上了权值。30第30页，本讲稿共57页 5.4 IIR滤波器：综合方法与滤波器特性滤波器：综合方法与滤波器特性 l引言：引言：l无限时间冲激响应（无限时间冲激响应（IIR）滤波器通常根据模拟原）滤波器通常根据模拟原型（巴特沃思、切比雪夫、椭圆等）来进行设计，型（巴特沃思、切比雪夫、椭圆等）来进行设计，并采用递归结构实现。并采用递归结构实现。l除射频端以外，除射频端以外，DSP物理硬件的滤波器已经是数物理硬件的滤波器已经是数字形式的了，模拟原型为起点的数字滤波器综字形
28、式的了，模拟原型为起点的数字滤波器综合的经典方法，就变得不太重要。合的经典方法，就变得不太重要。l根据模拟原型设计仍是基本和简单方法根据模拟原型设计仍是基本和简单方法31第31页，本讲稿共57页 5.4.1 冲激不变滤波器冲激不变滤波器 l假设数字滤波器的单位冲激响应，除了幅度比例有差别外，等于假设数字滤波器的单位冲激响应，除了幅度比例有差别外，等于模拟滤波器冲激响应的采样模拟滤波器冲激响应的采样该滤波器称为该滤波器称为模拟滤波器的冲激模拟滤波器的冲激不变实现。滤波器的传递函数不变实现。滤波器的传递函数表示L1拉普拉斯反变换，Z表示z变换，Ha(s)表示模拟原型的传递函数，T表示采样周期，同时
29、还是幅度比例系数。（5-17）32第32页，本讲稿共57页l一阶低通滤波器，其直流增益一阶低通滤波器，其直流增益(f=0)为为1。l例例5-2 由式（由式（5-17），首先对进行拉普拉斯反变换，由此），首先对进行拉普拉斯反变换，由此得到如下冲激响应：得到如下冲激响应：（5-17）（5-18）u(t)为单位阶跃函数，表示所得的滤波器是因果的。进行z变换之前，对冲激响应在采样时刻进行采样。例子33第33页，本讲稿共57页lz变换为变换为进行求和并乘上采样周期得：的冲激不变等价(5-22)34第34页，本讲稿共57页5.4.2 阶跃不变滤波器阶跃不变滤波器 l阶跃不变滤波器常用于仿真系统中阶跃不变滤
30、波器常用于仿真系统中l阶跃不变滤波器阶跃不变滤波器:数字滤波器的单位阶跃响应等数字滤波器的单位阶跃响应等价于模拟滤波器阶跃响应的采样。价于模拟滤波器阶跃响应的采样。l模拟滤波器阶跃响应的频域表示为模拟滤波器阶跃响应的频域表示为Ha(s)/s，其，其阶跃不变滤波器的传递函数如下：阶跃不变滤波器的传递函数如下：(5-23)35第35页，本讲稿共57页阶跃不变滤波器阶跃不变滤波器 l左边是传递函数为左边是传递函数为H(z)的数字滤波器对采样后的单位阶跃信号的响应，的数字滤波器对采样后的单位阶跃信号的响应，右边表示原型模拟滤波器单位阶跃响应采样的右边表示原型模拟滤波器单位阶跃响应采样的z变换。变换。l
31、再次采用（再次采用（5-18）式作为模拟原型。这个模拟原型的阶跃响应如下：）式作为模拟原型。这个模拟原型的阶跃响应如下：使用部分分式展开可得(5-23)(5-24)(5-25)36第36页，本讲稿共57页l取拉普拉斯反变换，并在取拉普拉斯反变换，并在t=nT处采样得处采样得再进行z变换得阶跃不变数字滤波器的传递函数为(5-26)(5-27)(5-28)OR(5-29)37第37页，本讲稿共57页Result将上式写成式（5-3）定义的标准形式有(5-30)38第38页，本讲稿共57页 5.4.3 双线性双线性z变换滤波器变换滤波器 l引言引言l双线性双线性z变换滤波器可能是仿真程序中使用最变换
32、滤波器可能是仿真程序中使用最普遍的滤波器了普遍的滤波器了l采用双线性变换方法的数字滤波器综合很直采用双线性变换方法的数字滤波器综合很直观并且完全是代数形式的观并且完全是代数形式的l双线性变换不会表现出混叠现象。双线性变换不会表现出混叠现象。?39第39页，本讲稿共57页 双线性双线性z变换滤波器变换滤波器综合方法综合方法l简单的代数变换，将模拟原型简单的代数变换，将模拟原型Ha（s)映射为一个数字滤波器，映射为一个数字滤波器，数字滤波器的传递函数定义为数字滤波器的传递函数定义为lC为常数，因为Ha（s)通常定义为的多项式之比，（5-31）定义的代数映射刚好给出的合适形式，即的z-1多项式之比。
33、l具有理想的属性，why？lexp(sdT)来替换双线性z变换中的变量z，sd是数字滤波器的复频域变量：（5-31）（5-32）40第40页，本讲稿共57页 双线性双线性z变换滤波器变换滤波器综合方法综合方法cont1l用用j2f fa a替换替换s,j2f fd d替换替换sd（5-33）得到数字滤波器和模拟原型的关系调节C值使得对任意采样频率都有fafd的方法称为预畸变（5-35）（5-34）或者41第41页，本讲稿共57页 双线性双线性z变换滤波器变换滤波器综合方法综合方法cont2fa，fd显然是非线性关系，余切函数期望：模拟数字 传递函数的形状保持不变尽可能线性，幅度相位不变x0,近
34、似线性变换要求 fs是采样频率.满足的近似线性关系所需C值。（5-35）线性关系要求（5-36）如要求fa fd（5-37）42第42页，本讲稿共57页例例5-4l为了说明双线性变换综合方法，我们再次使用（为了说明双线性变换综合方法，我们再次使用（5-18）定义的简单一阶模拟原型。）定义的简单一阶模拟原型。标准形式（5-18）（5-40）（5-39）（5-41）43第43页，本讲稿共57页例例5-5l这个例子将对前面三个例子中设计的三种滤波器的响应进行这个例子将对前面三个例子中设计的三种滤波器的响应进行比较。假设采样频率为比较。假设采样频率为100Hz（Ts0.01s），参数为），参数为a=2
35、1010，从而使模拟滤波器的，从而使模拟滤波器的3dB频率为频率为10Hz。冲激不变滤波器的传递函数为 阶跃不变滤波器的传递函数为（5-45）（5-43）（5-42）（5-44）44第44页，本讲稿共57页How about 双线性双线性z变换变换l第一种第一种l如果不采用频率预畸变，那么模拟、数字在低频段非常吻合如果不采用频率预畸变，那么模拟、数字在低频段非常吻合将C、a的值代入（5-48）（5-47）（5-46）45第45页，本讲稿共57页How about 双线性双线性z变换变换l第二种第二种l设选择设选择C值，使得频率响应在模拟滤波器的值，使得频率响应在模拟滤波器的3dB频率处相匹配频
36、率处相匹配（5-51）（5-50）（5-49）（5-48）46第46页，本讲稿共57页三个数字滤波器比较三个数字滤波器比较l%T=0.01;lf=0:0.1:50;lz=exp(-i*2*pi*f*T);%see(5.4)la0=0.239057;a1=0.239057;b1=0.521886;%bilinear invariantlnum=a0+a1*z;lden=1-b1*z;lampx=abs(num./den);la0=0.628319;b1=0.533488;%impulse invariantlnum=a0;lden=1-b1*z;lampy=abs(num./den);la0=1
37、.0;a1=0.533488;b1=0.533488;%step invriantlnum=(a0-a1)*z;lden=1-b1*z;lampz=abs(num./den);lplot(f,ampx,f,ampy,f,ampz)lxlabel(Frequency-Hz)lylabel(Amplitude Response)47第47页，本讲稿共57页Result单位增益零点48第48页，本讲稿共57页例例56l前面的例子为非常简单的一阶原型，前面的例子为非常简单的一阶原型，ak,bk可可以不用计算机确定以不用计算机确定l考虑一个五阶椭圆滤波器，具有考虑一个五阶椭圆滤波器，具有1dB通带波通带
38、波纹纹（Passband Ripple）和一个）和一个最小最小20dB的的阻带衰减。阻带衰减。l也可以用前面的方法做，分析方法比较复杂也可以用前面的方法做，分析方法比较复杂l也可以采用也可以采用Matlab的综合工具的综合工具49第49页，本讲稿共57页lfs=100;%set sampling frequencylfc=20;%set cuttoff frequencylf=0:0.1:50;%define frequency vectorlb,a=ellip(5,1,20,2*pi*fc,s);%synthesize elliptic filterlh=freqs(b,a,2*pi*f);
39、%amp.resp.of analog filterlbz1,az1=impinvar(b,a,fs);%impulse invariant digital filterlh1=freqz(bz1,az1,f,fs);%amplitude response of abovelbz2,az2=bilinear(b,a,fs);%bilinear z filter(not prewarped)lh2=freqz(bz2,az2,f,fs);%amplitude response of abovelbz3,az3=bilinear(b,a,fs,fc);%bilinear z filter(prew
40、arped)lh3=freqz(bz3,az3,f,fs);%amplitude response of abovelsubplot(211)%subplot 1lplot(f,abs(h),f,abs(h1)%plotlxlabel(Frequency-Hz)%label x axislylabel(Amplitude Response)%label y axislsubplot(212)%subplot 2lplot(f,abs(h2),f,abs(h3)%plotlxlabel(Frequency-Hz)%label x axislylabel(Amplitude Response)%l
41、abel y axis50第50页，本讲稿共57页51第51页，本讲稿共57页特殊情况：梯形积分法特殊情况：梯形积分法l离散时间积分由如下的表达式定义：离散时间积分由如下的表达式定义：两边进行z变换，可得梯形积分器的传递函数（5-52）（5-53）（5-55）（5-54）52第52页，本讲稿共57页特殊情况：梯形积分法特殊情况：梯形积分法如果基于双线性Z变换综合算法（5-55）（5-56）或者，积分运算等价于除以s，53第53页，本讲稿共57页5.4.4 IIR数字滤波器的计算机辅助设计数字滤波器的计算机辅助设计 l由于由于IIR方法能给出较小型的结构（较少的延时单元），它方法能给出较小型的结
42、构（较少的延时单元），它有时也能派上用场有时也能派上用场l计算机辅助设计方法的一个优点在于：无须等价模拟原型就可计算机辅助设计方法的一个优点在于：无须等价模拟原型就可以开发数字滤波器以开发数字滤波器l也便于快速而系统地研究实现中的折中，比如在复杂度以及精度也便于快速而系统地研究实现中的折中，比如在复杂度以及精度（实际与目标幅度响应之间的误差）。（实际与目标幅度响应之间的误差）。54第54页，本讲稿共57页例例5-7order=20,%degree of polynomialsf=0 0.1 0.12 0.23 0.25 0.35 0.37 1;%frequency pointsamp=1 1 0 0 1 1 0 0;%amplitude responseb,a=yulewalk(order,f,amp);%synthesize filterfreqz(b,a)%display results黑色区域是过渡带55第55页，本讲稿共57页56第56页，本讲稿共57页5.4.5 IIR滤波器的误差源滤波器的误差源57第57页，本讲稿共57页