离散系统的MATLAB实现.docxVIP

离散系统的MATLAB实现 、实验目的 学习利用matlab求解系统频率响应的方法。 学习利用matlab求解系统输出响应的方法。 加深对离散系统频率响应概念的理解。 二、设计内容及主要MATLAB函数 离散系统对应的输入输出差分方程为： y(n)-0.4y(n-1)-0.5y(n-2)=0.2x(n)+0.1x(n-1) 求该系统的单位抽样响应，频率响应及零极点增益。 单位抽样响应 系统的单位抽样响应是当输入信号为单位抽样信号时系 统的输出响应。 MATLAB中有两个函数可以计算系统的单位抽样响应： filter函数和impz函数。 1) filter 函数 filter函数是利用递归滤波器或非递归滤波器对数据 进行滤波。因为一个离散系统可以看做一个滤波器。系统 的输出就是输入经过滤波器滤波的结果。filter函数的格式 为： y=filter(b,a,x) 此函数是对由a和b组成的系统对输入进行滤波，如 果输入为单位抽样信号8(n),那么输出就是系统的单位 抽样响应。 2) impz函数 impz函数的格式为： impz(b,a)可直接给出系统的单位抽样响应。 频率响应H (ejw)) 由向量a和b组成的系统的频率响应为： M Zb(r)e-jwr r=0 H(ejw)= N 1+ Za(k)e-jkw k=1 MATLAB中的freqz函数使用基于FFT的算法来计算由 向量a和b组成的系统的频率响应。其一般用法为： [h,f]=freqz(b,a,n,Fs) 其中向量a和b由离散系统决定，Fs为采样频率，在 [0,Fs]频率范围内选取n个频率点，并记录在f中。由于freqz 函数是采用基2的FFT算法，n常取2的幂次方，这样可以 提高运算速度。 零极点增益 利用MATLAB中roots的函数可求得系统的零、极点，从 而得到系统的零极点增益表示。其一般用法为: zr=roots(b) zp=roots(a) 三、源程序 b=[0.20.1]; a=[-0.4-0.5]; x=[1,zeros(1,60)]; y=filter(b,a,x) stem(y) fs=1000; [h,f]=freqz(b,a,64,fs) abs(h) plot(f,abs(h)) p=angle(h) ph=p*180/p plot(f) roots(b) roots(a) zplane(b,a) 四、仿真结果 单位抽样响应图: 系统的频率响应图:  零点图:  极点图:二兑卯1怠晏盗T3| 口田|口回X|Pile Edit View Insert Tools Desktc-p 7fi ndow Help-1.5 -1■0.5 0Real Part0.5Note new toolbar butt 极点图: 二兑卯1怠晏盗T3| 口田|口回 X| Pile Edit View Insert Tools Desktc-p 7fi ndow Help -1.5 -1 ■0.5 0 Real Part 0.5 Note new toolbar butt oils : data bruihi ng linked plots Fl ：±7 vi d e 0 AJrtJE^E- 五、设计总结 在编程调试的过程中，不可避免地遇到了很多的困难。比如， 函数使用方法不正确，数据输入方式不正确等等。但在不断的修 改与调试下，成功的编好了程序，并且调试成功，做出了仿真结 果。学到了很多基本的MATLAB应用知识，也了解了 MATLAB 在信号与系统中的一些应用，更好地理解了系统的单位抽样响 应，频率响应及零极点增益的概念及意义。同时，我也更加熟练 地掌握了各相关函数的使用方法，增强了 MATLAB的实际应用 知识。