华东理工大学MATLAB资料.ppt

b [1 1/0.9];a [1 0.9]; [Fh,w] freqz b,a ; [Gd,w] grpdelay b,a ; subplot 311 plot w/pi,abs Fh ;ylabel |H w | ; grid on; axis [0 max w/pi 0 1.5] subplot 312 plot w/pi,angle Fh ; ylabel ang[H w ] ; grid on; subplot 313 plot w/pi,Gd ;ylabel grd[H w ] ; grid on; b [1 0 -1/0.81];a [1 0 -0.81]; [Fh,w] freqz b,a ; [Gd,w] grpdelay b,a ; subplot 311 plot w/pi,abs Fh ;ylabel |H w | ; grid on; axis [0 max w/pi 0 1.5] subplot 312 plot w/pi,angle Fh ; ylabel ang[H w ] ; grid on; subplot 313 plot w/pi,Gd ;ylabel grd[H w ] ; grid on; 6、最小相位系统 1 逆系统的概念： 对于一个系统函数为H z 的LSI系统，其逆系统定义为：系统函数Hi z 的逆系统与H z 级联后，总的系统函数是1。 2 最小相位系统概念的引出 一个因果稳定的系统H z ，如果要求它的逆系统Hi z 也是因果稳定的，那么要求该系统H z 的全部零点和极点都在单位圆内。 H z 的全部零点和极点都在单位圆内的系统称为最小相位系统。 ? 3 频率响应的补偿 —— 非最小相位系统分解的由来 G z 失真系统 Hd z 补偿系统 Hc z x n y n 信号处理的目标：G z 1 或 x n y n 信号要求： G z Hd z ·Hc z ∴ Hd z 与Hc z 互为逆系统。 对于一般的Hd z ，它是因果稳定的，但不一定是最小相位系统，可能有在单位圆外的零点。 非最小相位系统的分解 G z 失真系统 Hd z 补偿系统 Hc z x n y n 4 最小相位系统名称的含义及其性质 最小相位滞后性质 任何非最小相位系统H z 都能表示为： 其相位响应有如下关系： 解释：前面已经说明：全通系统的相位响应是非正（群延迟为正），意味着全通系统部分总是使相位滞后。因此，Hmin z 则具有最小相位滞后的特性，故称最小相位系统。 最小群延迟性质 解释：全通系统的群延迟总是为正，因此，最小相位系统Hmin z 具有最小群延迟的特性。 最小能量延迟性质 从时域分析，最小相位系统时域响应的波形延迟和能量延迟均最小。 例：一个数字滤波器具有以下性质 1 属于高通滤波器，具有一个零点，一个极点。 2 极点距z平面原点的距离为r 0.9。 3 直流信号不能通过该系统。 求： 1 画出该滤波器的零极点图。 3 画出该滤波器的幅频特性曲线。 2 为归一化频率特性，使|H ej? | 1，求滤波器的 系统函数。 4 写出系统的差分方程。 解： 1 系统零极点图： Re[z] jIm[z] 0 1 因为是高通，极点为-0.9 直流分量不通过，零点为1 2 系统函数： 3 滤波器的幅频特性 4 系统差分方程 |H ej? | ? ? 2? 0 1 例：将下面的系统函数所表示的高通滤波器转换为一个 陷波器。陷波器阻止频率为?0 ?/4，及其谐波通过 求： 1 写出系统差分方程。 2 画出系统零极点图。 3 画出两个滤波器的幅度特性曲线。 解： 1 先分析原滤波器： Re[z] jIm[z] 0 1 阻止频率为?0 0的高通滤波器 要求?0 ?/4及其谐波分量不通过，则应增加?0 ?/4及其谐波分量处的零极点。修改后的系统函数及差分方程： 2 画出系统零极点图。 Re[z] jIm[z] 0 1 3 画出两个滤波器的幅度特性曲线。 傅立叶和他的贡献 傅立叶分析是信号处理中最基本、最核心的内容之一 ，这种分析方法得名于法国数学家B.J.Fourier 1768 -1830 。实际上，傅立叶分析不仅应用于信号处理领域， 作为数学的一个主要分支，傅立叶分析在很多工程领域都 有广泛的应用。 在傅立叶之前，人们很早就已经认识到可以用三角函数来表示复杂的周期信号，甚至在古巴比伦和古埃及时期，人们就开始用这种方法来预测天文事件。 在近代科学史上，18世纪伟大的数学家欧拉重新使这种方法走进