新·MATLAB命令绘出当激励信号为x(n)=u(n)时,该系统的零状态响应及系统的零状态响应的11个数值解.pptVIP

姓名： 班级： 学号： 指导老师： 一.设计题目： 37.已知某LTI系统的差分方程为y(n)-0.9y(n-1)=x(n).试用MATLAB命令绘出当激励信号为x(n)=u(n)时，该系统的零状态响应及系统的零状态响应的11个数值解，与理论结果对照。 二．设计目的： 掌握零状态响应的解法。 1. 差分方程的变换域解法： 设LTI离散系统的激励为f(k),响应为y(k),描述n阶系统的后向差分方程为 式中系数 均为实数，设f(k)是在k=0时接入（即若k0,f(k)=0）,且系统的初始状态为y(-1),y(-2),······,y(-n+1)，y(-n),令Z[y(k)]=Y(Z),Z[f(k)]=F(Z). 根据单边z变换的移位特性可得： Z[y(k-i)]= + (1)Z[f(k-j)]= + = (2) 对差分方程两边取z变换,并将式（1），式（2）代入，得： 由上式可得： Y(z)= + F(z) 式中： A(z)= , B(z)= 是仅与差分方程系数有关的z的多项式； 系统的零状态响应为： yf(k)= Z-1[Yf(z)]= Z-1 [ F(z)] 2.计算过程： 由于A(z) = 1- B(z)=1 (z)= X(z) = 由部分分式法展开求解： = = + 由留数定理求出系数A,B A=Res[ (z),0.9] = (z-0.9)︱z=0.9= =-9 B=Res[ (z),1]= (z-1)︱z=1= =10 所以：(z)= +取其逆变换得零状态响应： (n)=[(-9) +10]u(n)=10(1- )u(n)  代入数值求得： (0)=10(1- )u(0)=1 (1)=10(1- )u(1)=1.9 (2)=10(1- )u(2)=2.71 同理可得： (3)=3.439 (4)=4.0951 (5)=4.68559 (6)=5.217031 (7)=5.6953 (8)=6.1258 (9)=6.5132 (10)6.8619 四.MATLAB源程序代码： num=[1]; den=[1,-0.9]; n=0:10; x=(n=0); y=filter(num,den,x) subplot(2,1,1);stem(x);title(LTI系统的输入序列); subplot(2,1,2);stem(y);title(LTI系统的零状态响应); 数值解： t =0 1 2 3 4 5 6 7 8 9 10 y =1.0000 1.9000 2.7100 3.4390 4.0951 4.6856 5.2170 5.6953 6.1258 6.5132 6.8619 图形： 1.开始的时候对离散域的Z域解法不是很清楚，造成了计算后结果的差异。真正学懂后解决就很简单了。 2.在MATLAB上面编程也出现了很多问题，主要是自己对这门软件的了解不够深，并且对这门软件在信号处理方面的这一大块的函数及调用都了解的很少，通过看书查阅资料才慢慢了解了MATLAB在信号处理这方面强大的功能。 3.对函数的格式运用缺乏理解。 七.自我感受： 通过本次课程设计，我从中学到了很多，认识到自己的不足，要想真正的用MATLAB这门软件将信号处理的完美无缺就需要我们专门去学习这门软件，将它各个方面的知识都要掌握，这也让我认识到了自己的不足。但通过设计也督促我去查阅更多的资料，让我从各个不同的方面去考虑解决问题的方法，也了解了课堂中没有注意到问题及自己在实践与理论相结