数字信号处理实验报告——离散系统的频率响应分析和零极点分布.docxVIP

数字信号处理实验报告——离散系统的频率响应分析和零极点分布 实验 3 离散系统的频率响应分析与零、极点分布 实验目的 :加深对离散系统的频率响应分析与零、极点分布的概念理解。 实验原理 : 离散系统的时域方程为 N M d k y ( n k ) p k x ( n k ) k 0 k 0 其变换域分析方法如下 : y[n] x[ n] h[n] x[m]h[n m] Y(e j ) X (e j )H (e j ) 频域 m H ( ej ) p(e j ) p0 p1e j 系统的频率响应为 D(e j ) d0 d1e j  ... pM e jM ... d Ne jN y[n] x[n] h[n] x[m]h[n m] Y( z) X ( z) H (z) Z 域 m H ( z) p(z) p0 p1 z 1 ... pM z M D(z) d0 d1z 1 ... dN z N 系统的转移函数为 M p z i M z 1 ) (1 i i k i 1 H (z) 0 K N N i z 1 ) 分解因式 dk z i (1 , 其中 i 与 i 称为零、极 i 0 i 1 点。 MATLAB中 , 可以用函数 [z,p,K]=tf2zp(num,den) 求得有理分式形式的系统转移函数的零、极点 , 用函数 zplane(z,p) 绘出零、极点分布图 ; 也可以用函数 zplane(num,den) 直接绘出有理分式形式的系统转移函数的零、极点分布图。 另外 , 在 MATLAB中 , 可以用函数 [r,p,k]=residuez(num,den) 完成部分分式 展开计算 ; 可以用函数 sos=zp2sos(z,p,K) 完成将高阶系统分解为 2 阶系统的串 联。 实验内容 : 求系统 0.0528 0.797z 1 0.1295z 2 0.1295z 3 0.797z 4 0.0528z 5 H (z) 1.8107z 1 2.4947z 2 1.8801z 3 0.9537z 4 0.2336z 5 1 的零、极点与幅度频率响应。 实验要求 : 编程实现系统参数输入 , 绘出幅度频率响应曲线与零、极点分布图。 实验过程 : 数字信号处理实验报告——离散系统的频率响应分析和零极点分布 画出零极点图 : num=[0 、 0528,0 、 797,0 、 1295,0 、 1295,0 、 797,0 、 0528]; den=[1,-1 、 8107,2 、 4947,-1 、 8801,0 、 9537,-0 、2336]; [z,p,K]=tf2zp(num,den); zplane(b,a); title(零极点图 ); 绘制幅度频率响应曲线 : num=[0 、 0528,0 、 797,0 、 1295,0 、 1295,0 、 797,0 、 0528]; den=[1,-1 、 8107,2 、 4947,-1 、 8801,0 、 9537,-0 、2336]; [H,W]=freqz(num,den); M=abs(H); subplot(2,1,1); xlabel(Mormalized Frequency(x Πrad/sample)) ; ylabel(Magritude(dB)); title(幅度频率响应曲线 ); grid; subplot(2,1,2); plot(W); xlabel(Mormalized Frequency(x Πrad/sample)); ylabel(Phase(degree)); grid; 数字信号处理实验报告——离散系统的频率响应分析和零极点分布