matlab由频率响应计算差分方程,MATLAB在数字信号处理中的应用3.ppt

MATLAB在数字信号处理中的应用3.ppt

2.2.2 ROC and pole-zero 例10： (1) 用因式分解形式表示上面的z 变换 (2) 画出零点和极点 (3) 确定收敛域 常用z 变换得到的是有理多项式 常常需要进行因式分解，求出极零点和收敛域 2.2.2 ROC and pole-zero 方法一：符号运算 MATLAB中求因式分解：factor, 将因式分解展开成多项式：expand syms z % construct symbolic object 'z' NUM = 2*z^4 + 16*z^3 + 44*z^2 + 56*z + 32; DEN = 3*z^4 + 3*z^3 - 15*z^2 + 18*z - 12; % prime factors NUMFAC = factor(NUM) DENFAC = factor(DEN) % verify by expanding iNUM = expand(NUMFAC) iDEN = expand(DENFAC) 例10： (1) 用因式分解形式表示上面的z 变换 2.2.2 ROC and pole-zero 例10： (1) 用因式分解形式表示上面的z 变换 方法一：数值运算(求极点和零点) MATLAB中求极零点：roots, 根据极零点求多项式系数：poly num = [2,16,44,56,32]; den = [3,3,-15,18,-12]; % determine zeros and poles zeros = roots(num) poles = roots(den) % polynomial coefficients inum = poly(zeros); iden = poly(poles); gain = num(1)*iden(1)/inum(1)/den(1) 2.2.2 ROC and pole-zero MATLAB中画零极点图：zplane, zplane(zeros,poles), zplane(num,den) 例10： (2) 画出零点和极点 >> zplane(zeros,poles) >> zplane(num,den) (3) 确定收敛域 略 2.2.3 inverse z-Transform 例11：求的z 反变换。 MATLAB中函数：[R,p,C] = residuez(b,a) R：留数；p：极点；C：直接项 思考：在MATLAB中如何求z反变换 方法一：iztrans 方法二：部分分式展开法 b = 1; a = poly([0.9,0.9,-0.9]); [R,p,C] = residuez(b,a) 2.2.3 inverse z-Transform 例11：求的z 反变换。 求幂级数系数可以通过求的单位脉冲响应 在MATLAB中可以通过filter函数实现 方法二：长除法(幂级数法) [delta,n] = impseq(0,0,7); ix = filter(b,a,delta) x = 0.75*0.9.^n+0.5*n.*0.9.^n+0.25*(-0.9).^n 长除法可以用于验证z变换和z反变换的正确性 2.2.3 inverse z-Transform 例11：求的z 反变换。 R = 0.2500 0.2500 + 0.0000i 0.5000 - 0.0000i p = -0.9000 0.9000 + 0.0000i 0.9000 - 0.0000i C = [] 2.2.4 system function 系统函数是单位脉冲响应的z变换：ztrans, 单位脉冲响应是系统函数的z反变换：iztrans, residuez, filter MATLAB中由系统函数计算单位脉冲响应：[x,n] = impz(b,a,N) 作用等同于filter(长除法)，输出为列向量 例11：求的z 反变换。 [ix,n] = impz(b,a,7); n = n'; ix = ix' x = 0.75*0.9.^n+0.5*n.*0.9.^n+0.25*(-0.9).^n 2.2.4 system function 系统函数的分子系数就是差分方程中 x 的系数——b 系统函数的分母系数就是差分方程中 y 的系数——a MATLAB中对于差分方程的计算可以用于系统函数的计算 思考：如何由差分方程计算系统函数 2.2.4 system function M