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

《数字信号处理实验报告——离散系统的频率响应分析和零极点分布.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、实验3离散系统的频率响应分析和零、极点分布实验目的：加深对离散系统的频率响应分析和零、极点分布的概念理解。实验原理：离散系统的时域方程为其变换域分析方法如下：频域系统的频率响应为Z域系统的转移函数为分解因式，其中和称为零、极点。在MATLAB中，可以用函数[z，p，K]=tf2zp（num，den）求得有理分式形式的系统转移函数的零、极点，用函数zplane（z，p）绘出零、极点分布图；也可以用函数zplane（num，den）直接绘出有理分式形式的系统转移函数的零、极点分布图。另外，在MATLAB中，可以用函数[r，p，k]=residuez（

2、num，den）完成部分分式展开计算；可以用函数sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。实验内容：求系统的零、极点和幅度频率响应。实验要求：编程实现系统参数输入，绘出幅度频率响应曲线和零、极点分布图。实验过程：画出零极点图：>>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('零极点图'

3、);>>绘制幅度频率响应曲线：>>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('MormalizedFrequency(xΠrad/sample)');>>ylabel('Magritude(dB)');>>title('幅度频率响应曲线');>>grid;>>subplot(2,1,2

4、);>>plot(W);>>xlabel('MormalizedFrequency(xΠrad/sample)');>>ylabel('Phase(degree)');>>grid;>>