《系统的相位结构数学信号处理课件.ppt》由会员分享，可在线阅读，更多相关《系统的相位结构数学信号处理课件.ppt（79页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第5章 离散时间系统的相位、结构与逆系统5.1离散时间系统的相频响应；5.2FIR系统的线性相位；5.3具有线性相位系统的零点分布；5.4全通系统和最小相位系统；5.5谱分解；5.6FIR系统的结构；5.7离散时间系统的Lattice结构；5.8逆系统5.9MATLAB文件5.1 离散时间系统的相频响应幅频响应相频响应如果：我们称其为线性相位。若：也称线性相位对输入，有假定：所以：输出是输入的简单移位，移位的大小正比于因此不会发生失真。例：令则：没有发生相位失真具有线性相位例：令若：则：发生了相位失真如果令：再令：则：则：由于：定义：如果系统的相频响应不是线性的，那么系统的输出将不再是输入信号
2、作线性移位后的组合，因此，输出将发生失真。定义：为系统的群延迟(GroupDelay,GD)为系统的相位延迟(PhaseDelay,PD)显然，若系统具有线性相位，则其GD为常数。若：则：即：相位延迟反映了载波信号的延迟，而群延迟反映了输出包络的延迟。思考：如何实现对信号的零相位滤波？若要保证系统是因果的，又如何实现？5.2 FIR 系统的线性相位在绝大部分信号处理的场合，人们都期盼系统具有线性相位，但是，如何实现线性相位？对FIR系统，如果保证：则该系统具有线性相位。上述对称有四种情况：第一类FIR系统偶对称奇对称第二类FIR系统1.为奇数令：并利用的对称性，有第一类FIR系统令：令：实数最
3、后有：相位增益所以，只要保证滤波器的系数偶对称，该滤波器必然具有线性相位。2.为偶数令：则：第二类FIR系统：3.为奇数4.为偶数请掌握四种情况下线性相位表达式的推导方法。的线性组合，在时，易取得最大值，因此这一类滤波器易体现低通特性，且是偶函数。通过频率移位，又可体现高通、带通、带阻特性。所以，经典的低通、高通、带通和带阻滤波器的都是偶对称的。说明：第一类FIR系统是的线性组合，在时，的值为零，且是奇函数。这一类滤波器都是作为特殊形式的滤波器，如Hilbert变换器、差分器等。第二类FIR系统是最好取为奇数，以便以中心点为对称。思考：四类滤波器的对称点在何处5.3 具有线性相位系统的零点分布
4、所以，的零点也是的零点,反之亦然令：则：的零点分布:零点分布可能有四种情况：1.不在实轴也不在圆上，应是一对共轭零点，模1;2.不在实轴，但在圆上，也是一对共轭零点；模1；3.在实轴但不在圆上，无共轭，角度0，模1;4.在实轴，但在圆上，无共轭，角度0，模1;四个零点同时存在,构成四阶系统.在单位圆内把该式展开，其系数也是对称的，是具有线性相位的子系统。无共轭零点,有镜象零点无镜象对称零点,有共轭零点.一个具有线性相位的FIR数字滤波器的转移函数可表示为上述四类FIR子系统的级联，即：很容易证明，每一个子系统的系数都是对称的，因此它们都具有线性相位。无镜象零点,也无共轭零点.5.4 全通系统和
5、最小相位系统如果一个系统的幅频响应对所有的频率都等于1(或一个常数),即则称系统为全通系统。最简单的全通系统，纯延迟全通系统一阶全通系统：镜像对称二阶全通系统：一对位于单位圆内的共轭极点，一对共轭零点和极点以单位圆为镜像对称。高阶全通系统：高阶全通系统的另一种表示形式：即：对该全通系统，请自己证明：1.是IIR系统(不考虑纯延迟形式）；2.极点数和零点数相等；3.极点和零点是以单位圆镜像对称的；4.极点都在单位圆内,零点都在单位圆外；5.全通系统的群延迟始终为正值。全通系统的特点：IIR系统的无限长，无法对称，即无法作到线性相位。在实际中，可以用一个全通系统和IIR系统相级联，在不改变幅频响应
6、的情况下对相频响应做矫正，使其接近线性相位。全通系统的应用：全通系统还广泛应用在系统分析及一些特殊滤波器的设计方面（如功率互补IIR滤波器组）一阶全通系统极零图幅频相频抽样响应三阶全通系统一个离散系统，其极点必须在单位圆内，但对零点没有限制，如果：1.所有的零点都在单位圆内：最小相位系统；2.所有的零点都在单位圆外：最大相位系统；3.单位圆内、外都有零点：混合相位系统。最小相位系统1.在具有相同幅频响应的因果的稳定的滤波器2.集合中,最小相位滤波器具有最小的相位偏移；最小相位系统的性质：例：作为作业，请证明如下两个系统具有相同的幅频响应：那一个是最小相位系统幅频相频2.在所有具有相同幅频响应的
7、离散系统中,最3.小相位系统的具有最小的延迟；令：累计能量有：所以，最小相位系统的单位抽样响应又称最小延迟序列。思考：具有线性相位的FIR系统是否是最小相位系统？例.三个系统：它们具有相同的幅频响应，试判断，那一个是最小相位系统？最大相位系统？混合相位系统？请注意：为保证系统具有相同的幅频响应（相同的定标），的表达式。3.设为最小相位系统令：构成一对Hilbert变换则：和复倒谱：Cepstrum4.对于稳定因果系统，当且仅当其是最小相位5.系统时,该系统才有逆系统（InverseSystem)。令：记：的逆系统Deconvolution(反卷积）Systemidentification(系统
8、辨识）5.任一非最小相位的因果系统的转移函数均可由一个最小相位系统和一个全通系统级联而成,即：由于最小相位系统有着以上特殊的性质，因此有着广泛的应用，特别是在信号的建模与系统辨识方面。要理解，具有相同幅频响应的系统，它们所对应的转移函数可以是不相同的，区别就在于相位（或零点的位置）。那么，如何由一个最小相位系统得到具有相同幅频响应的最大相位、混合相位系统？5.5 谱分解（Spectral factorization)令：显然，具有线性相位。将一个转移函数的极零点重新分配，得到两个转移函数，这一过程（或方法）就称为“谱分解”。最常用的是将具有线性相位系统的转移函数作分解，并且往往是分解成两个具有
9、相同幅频响应的子系统。=1.0000，4.0500，8.1000，14.9956，27.7248，43.2996，51.1831，43.2996，27.7248，14.9956，8.1000，4.0500，1.0000例.令显然，该系统具有线性相位，共有12个零点：下图是对作谱分解的结果，可以看出，分解后的两个系统具有相同的幅频响应。谱分解的目的是想得到因果的、符合某种要求的系统，这在信号建模、特殊滤波器的设计中经常要用到。分解的一般方法是：令一个系统是最小相位系统；则另一个系统必然是最大相位系统。这样，两个系统都有着相同的幅频响应。另外一种分解方法是得到两个混合系统，目的是保证它们都具有线性
10、相位。5.6 FIR 系统的结构直接实现:一、直接实现和级联实现级联实现:乘法量减少一半二、具有线性相位的FIR系统的结构FIR系统该系统实际上是一个N点平均器。IIR系统三、FIR系统的递归实现及梳状滤波器该系统可由一FIR系统和一个一阶IIR系统级联而成，极零点抵消后，仍是一FIR系统。令令IIR实现梳状滤波器N点平均器思路：用DFT系数表示系统函数四、频率抽样实现令：梳状滤波器N个一阶IIR系统则：可按上述级联方式得到系统的信号流图：该结构一方面反映了Z变换、DTFT、DFT之间的关系，另一方面，给出了FIR滤波器设计的一种有效方法。5.7 离散时间系统的 Lattice 结构Latti
11、ce结构又称“格形”结构，是一种非常新颖、有特色的结构，在基于模型的功率谱估计、语音信号处理、自适应滤波方面有着重要的应用。对一个FIR系统，其Lattice结构是：反射系数Lattice结构的基本单元1.全零点系统(FIR)的Lattice结构如何实现滤波器系数和的相互转换定义：Lattice结构中的基本关系：是Lattice结构中第m个上、下结点相对输入端的转移函数。得到由低阶倒高阶，或由高到低的递推关系。得到时域递推关系：低到高阶高到低阶MATLAB中有相应的m文件。例：看作是FIR系统的逆形式。2.全极点系统(IIR)的Lattice结构的求解方式同FIR系统Lattice结构的计算方
12、法,只是将多项式的系数换成.系数系数及及注意：在递推求解的过程中，反射系数有关反射系数的更多讨论见第12章信号建模。3.极零系统的Lattice结构上半部对应全极系统上半部对应全极系统下半部对应全零系统下半部对应全零系统两组Lattice系数求出同全极系统；递推求解5.8 离散系统的状态变量描述描述：差分方程、转移函数、线性卷积1.LSI系统的状态变量与状态方程转移函数、差分方程、中间变量的关系1.“状态”指系统内一组变量,它包含了系统全部2.过去的信息,由这一组变量和现在与将来的3.输入，可求出现系统现在和将来的全部输出；4.2.可用于分析多输入、多输出系统；如何选择状态变量？有着不同的方法
13、。方法之一是选择作为系统的状态。定义一组新的变量相互关系状态方程输出方程上述内容讨论了如何由差分方程转换为状态方程。当然，反过来也可以。两边取两边取Z变换：变换：2.由状态方程求系统的转移函数状态方程输出方程3.由状态方程求输出及单位抽样响应抽样响应为：零输入解零状态解例对系统，当时，即是系统的单位抽样响应，显然，该序列称为Fibonacci序列。试利用状态方程求。解：1fiftfilt.m本文件实现零相位滤波。其调用格式是：y=filtfilt(B,A,x)。式中B是的分子多项式，A是分母多项式，x是待滤波信号，y是滤波后的信号。2grpdelay.m求系统的群延迟。调用格式gdw=grpd
14、elay(B,A,N),或gdF=grpdelay(B,A,N,FS)式中B和A仍是的分子、分母多项式，gd是群延迟，w、F是频率分点，二者的维数均为N；FS为抽样频率，单位为Hz。5.9与本章内容有关的MATLAB文件3tf2latc.m和latc2tf.m：实现转移函数和Lattice系数之间的相互转换。tf2latc的调用格式是：(1)k=tf2latc(b),(2)k=tf2latc(1,a),(3)k,c=tf2latc(b,a),其中（1）对应全零系统，（2）对应全极系统，（3）对应极零系统。latc2tf的调用格式和tf2latc正好相反。需要说明的是，tf2latc求出的Lat
15、tice系数k和本书求出的k差一个负号，这是由于我们在图中用的是k。4.latcfilt.m用来实现Lattice结构下的信号滤波。调用格式是：(1)y,g=latcfilt(k,x)：对应全零系统(2)y,g=latcfilt(k,1,x)：对应全极系统(3)y,g=latcfilt(k,c,x)：对应极零系统x是待滤波的信号，y是用Lattice结构作正向滤波的输出，g是作反向滤波的输出。若输入x是则输出y是的系数；g是的系数。5.tf2ss.m和ss2tf.m实现转移函数和相应状态变量之间的转换。二者的调用格式分别是：A,B,C,D=tf2ss(b,a),b,a=ss2tf(A,B,C,D)。式中b,a分别是分子、分母多项式的系数向量，A,B,C及D的定义见书。6.sos2ss.m实现由转移函数的二阶级联形式转换为状态变量表示。调用格式：A,B,C,D=sos2ss(sos,g),A,B,C,D的定义 见 书。有 关 sos和 g的 说 明 见 2.8节 关 于tf2sos.m的说明。