(第八次实验报告.docVIP

信号与系统实验报告 姓名：___杨___超___ 班级：___0101110___ 学号：_2011210309__ 指导老师：___郑丹玲____ 实验八 LTI离散系统的频域分析 实验目的 利用MatLab绘制LTI离散系统的零极图； 根据离散系统的零极点分布，分析系统单位响应h(k)的时域特性； 利用Matlab求解LTI离散系统的幅频特性和相频特性。 实验原理 1、LTI离散系统可采用线性常系数差分方程来描述，其中y(k)为系统输出信号，f(k)为系统输入信号。 2、离散系统的稳定性 离散系统稳定的充要条件：系统函数H(z)的所有极点均位于z平面的单位圆内。对于三阶以下的低阶系统，利用求根公式可方便的求出离我三系统的极点位置，判断系统的稳定性。对于高阶系统，手工求解极点位置则非常困难，这是可利用Matlab来实现。 3、离散系统的频率响应 4、分析离散系统特性常用的函数 （1）roots()函数 求系统函数的零点和极点位置，调用命令格式如下： p=roots(A) A为待求根的关于z的多项式的系数构成的行向量，返回向量p则是包含该多项式所有根位置的列向量。 （2）tf2zp()函数 调用格式为:[q,p,c]=tf2zp(b,a) b和a分别表示H（z）的分子和分母多项式的系数向量，q和p分别为求得的零点和极点，c为H（z）中的常数。 （3）zplane()函数 应用在zplance函数可以得到系统函数H(z)的零极点分布图，其调用格式为： zplance(b,a) （4）impz()函数 根据H（z）零极点分布，绘制系统的单位响应h(k)的序列波形，其调用格式如下： impz(b,a) （5）freqz()函数 计算离散系统频率响应的数值，再利用函数abs(),angle()及plot()函数，可绘制出系统的幅频特性曲线和相频特性曲线。Freqz()的调用格式： 【H，w】=freqz(b,a,N) 三、实验内容 figure(1) b=[1 0]; a=[1 -0.8]; subplot(2,1,1) zplane(b,a); subplot(2,1,2) impz(b,a); figure(2) b=[1 0]; a=[1 0.8]; subplot(2,1,1) zplane(b,a); subplot(2,1,2) impz(b,a); figure(3) b=[1 0]; a=[1 -1.2 0.72]; subplot(2,1,1) zplane(b,a); subplot(2,1,2) impz(b,a); figure(4) b=[1 0]; a=[1 -1]; subplot(2,1,1) zplane(b,a); subplot(2,1,2) impz(b,a); figure(5) b=[1 0]; a=[1 -1.6 1]; subplot(2,1,1) zplane(b,a); subplot(2,1,2) impz(b,a); figure(6) b=[1 0]; a=[1 -1.2]; subplot(2,1,1) zplane(b,a); subplot(2,1,2) impz(b,a); figure(7) b=[1 0]; a=[1 -2 1.36]; subplot(2,1,1) zplane(b,a); subplot(2,1,2) impz(b,a); a=[1 0 0 0 0 0 0 0 -0.9]; b=[1 0 0 0 0 0 0 0 -1]; h=freqz(b,a,400,whole); h1=abs(h); h2=angle(h); subplot(2,1,1) plot(h1); xlabel(角频率); ylabel(幅度); title(幅频特性曲线); subplot(2,1,2) plot(h2); xlabel(角频率); ylabel(); title(幅频特性曲线); 四、思考题 1、从实验内容1的结果进行分析，说明系统函数的零极点分布于系统时域特性的关系，那个系统不稳定？ 答：系统函数的所有极点均位于z平面的单位圆内，则系统的稳定性高。第四、五、六、七个系统不稳定。 2、实验内容2中的离散系统是一个说明类型的滤波器，具有什么功能？ 答：该系统是低通滤波器，它具有通低频阻高频的功能。 五、实验心得 通过这次试验，我了解了如何用MatLab进行信号过滤，如何画出LTI离散系统的零极图，还能根据离散系统的零极点分布，分析系统单位响应h(k)的时域特性，最后，还能利用Matlab求解L