华东理工大学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)。实际上，傅立叶分析不仅应用于信号处理领域， 作为数学的一个主要分支，傅立叶分析在很多工程领域都 有广泛的应用。 在傅立叶之前，人们很早就