资源描述：
《离散系统的matlab的实现》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、离散系统的MATLAB实现一、设计目的通过该设计，理解系统的单位抽样响应，频率响应及零极点增益的概念及意义。二、课程设计环境计算机、MATLAB软件三、设计内容及其主要MATLAB函数如图所示的一个离散系统该离散系统对应的输入输出查分方程为：y(n)-0.4y(n-1)-0.5y(n-2)=0.2x(n)+0.1x(n-l)求该系统的单位抽样响应，频率响应及零极点增益。1、单位抽样响应系统的单位抽样响应是当输入信号为输入信号为单位抽样信号时系统的输岀响应。MATLAB中有两个函数可以计算系统的单位抽样响应:filte
2、r函数和impz函数。1)filtei•函数filter函数是利用递归滤波器或非递归滤波器对数据进行滤波。因为一个离散系统可以看做一个滤波器。系统的输出就是输入经过滤波器滤波的结果。filtei•函数的格式为：y=filter(b,a,x)此函数是对由a和b组成的系统对输入进行滤波，如果输入为单位抽样信号5("),那么输出就是系统的单位抽样响应。2)impz函数impz函数的格式为：impz(b,a)可直接给出系统的单位抽样响应。2、频率响应H(八由向量a和b组成的系统的频率响应为MH@")=f1+工水)以伙MATL
3、AB中的freqz函数使用基于FFT的算法来计算由向量a和b组成的系统的频率响应。其一般用法为：【h,f】=fwqz(b,a,n,Fs)其中向量a和b由离散系统决定，Fs为采样频率，在【0,空】2频率范围内选取n个频率点，并记录在f中。由于freqz函数是从用基2的FFT算法，n常取2的幕次方，这样可以提高运算速度。3、零极点增益利用MATLAB中roots的函数可求得系统的零、极点，从而得到系统的零极点增益表示。其一般用法为：zr=roots(b)zp=roots(a)1、按题目要求进行设计，写出设计报告，给出源程
4、序。2、分别利用filter函数和impz函数画出给定的离散系统的单位抽样响应。3、画出给定离散系统的频率响应。4、画出系统的零、极点图。五、设计总体方案1・、离散系统单位抽样响应(1)利用filtei•函数计算系统的单位抽样响应。其一般格式为:y=filter(b,a,x)然后利用plot函数可以画出给定的离散系统的单位抽样响应图。其一般格式为：plot(x,y)还可以对图形进行修饰，例如plot(x,y,'s‘)表1线性、点型和颜色标识符颜色标识符线型和点型y黄占八、、m品红0圆圈c青XX号r红++号g绿—实线b
5、蓝*星号w白■■虚线k黑■点划线一长划线表2其他修饰函数图例标注标题X轴y轴legend^文text(x,y,'文title('文xlabel('文ylabel('文字')字')字，)字')字')(2)利用impz函数可直接画出系统的单位抽样响应，其一般格式为impz(b,a)2、离散系统的频率响应利用freqz函数计算由向量3和匕组成的系统的频率响应，其一般格式为：【h,f]=freqz(b,a,n,Fs)再利用plot函数画出离散系统的幅频特性和相频特性图，其中h为y轴，f为x轴。3、零、极点增益利用roots函
6、数可计算系统的零、极点。其一般格式为：zr=roots(b)zp=roots(a)利用polar函数画出zr的幅频特性曲线和相频特性曲线。其中polar函数是极坐标绘图函数。在画图时需要用到打开图形保持holdon和关闭图形保持holdoffo再利用zplane函数得出相应离散体统的零极点增益图。与roots函数得出来的图相对照，看是否一致。zplane函数的一般格式为：zplane(b,a)六、设计参数由图中给出的离散系统对应的输入差分方程：y(n)-0.4y(n-1)-0.5y(n-2)=0.2x(n)+0.1x
7、(n-l)得出：b=[0.2,0.1];a=[lr0.4,-0.5]在用filter函数中x的选值只要有一个1就行，其他的数值必须为零。在用freqz函数中，【h,f】=freqz(b,a,n,Fs)。Fs为采样频率，n的取值在【0,竺】频率范围内，并记录在f中。由于freqz函数是2从用基2的FFT算法，n常取2的幕次方，这样可以提高运算速度。所以可取：Fs=20；n=2A3・七、源程序和实验图Vl＞用filter函数»a=[l,-0.4r0.5];»b=[0.2,0.1J；»x=[O,O,l,O,O];»y=fi
8、lter(b,a,x)00.20000.18000.1720»plot(x,y；o*)离散系统的单位抽样响应v2＞用impz函数»impz(b,a)离散系统的单位抽样响应n(samples)6oo.2864218o.「「「「O..Ooooooapa一-duJY420ooo.o.v3＞求频率响应»a=[l,-0.4,-0.5];»b=[0・2,0
