MATLAB实验教案.docxVIP

实验一 离散系统的时域分析和复频域分析1.实验目的（1）掌握在时域求系统响应的方法。（2）掌握时域离散系统的时域特性。（3）通过实验判断系统稳定性（4）掌握利用Z变换对系统进行复频域分析。（5）掌握系统零、极点的绘制方法。（6）通过复频域分析系统稳定性、频率特性。（7）熟悉Z变换的应用2.实验设备计算机MATLAB R2012a仿真软件3.实验原理（1）离散系统的时域分析在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，利用filter函数或conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。系统的时域特性是指系统的线性移不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应，或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号、输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的。系统的稳态输出是指当n→∞时系统的输出。如果系统稳定，信号加入系统后，系统输出的开始一段称为暂态效应，随着n的加大，幅度趋于稳定，达到稳态输出。注意在以下实验中均假设系统的初始状态为零。（2）离散系统的复频域分析离散系统的时域方程为 其变换域分析如下 频域系统频率响应为 Z域系统的转移函数为分解因式=,其中，和称为零、极点。在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阶系统的串联。判断Z域因果LTI稳定性：当且仅当系统函数的全部极点位于单位圆内时，因果LTI系统是稳定的。4.实验内容及步骤（1）给定矩形脉冲 是脉冲响应为的LTI系统的输入，求输出 。利用filter函数对差分方程进行数值求解，画出输入序列和脉冲响应。（2）已知差分方程： ① 画出在 的脉冲响应。② 画出在 的单位阶跃响应。③ 判断由 表征的这个系统的稳定性。（3） 一个线性时不变系统由差分方程描述如下： ① 画出 内的系统的脉冲响应，并确定稳定性。② 如果这个系统输入为 ,求在范围内的响应 。（4）求下列序列的Z变换并用“zplane”画出零、极点分布图。 ，该序列的Z变换为： （5）求下列直接型系统函数的零、极点，并将它转换成二阶节形式 （6）给定因果系统 要求:① 求系统函数并画出零、极点示意图。② 画出系统的幅频特性 和相频特性 。③求脉冲响应 并画出序列图。5.程序清单如下(1)内容1:调用filter解差分方程b=[1];a=[1,-0.9];n=-5:50;x=stepseq(0,-5,50)-stepseq(10,-5,50);h=0.9.^n;h(1:5)=0;y=filter(b,a,x);subplot(3,1,1);stem(n,x);ylabel(x(n));axis([-5 50 min(x) max(x)]);subplot(3,1,2);stem(n,h);ylabel(h(n));axis([-5 50 min(h) max(h)]);subplot(3,1,3);stem(n,y);xlabel(n);ylabel(y(n));axis([-5,50,-0.5,8]);（2）涉及的功能函数stepseqfunction[x,n]=stepseq(n0,n1,n2)n=[n1:n2];x=[(n-n0)=0];(3)内容2：稳定性分析b=[1];a=[1,-1,0.9];x=impseq(0,-20,120);n=[-20:120];h=filter(b,a,x);subplot(2,1,1);stem(n,h);title(脉冲响应);xlabel(n);ylabel(h(n))x=stepseq(0,-20,120);s=filter(b,a,x);subplot(2,1,2);stem(n,s);title(阶跃响应);xlabel(n);ylabel(u(n))（4）涉及的功能函数impseqfunction[x,n]=im