《精》连续时间系统的LTI系统的时域仿真.docVIP

信号与系统课程设计报告书 课题名称 连续时间系统的LTI系统的时域仿真 ------零输入响应与零状态响应 姓 名 王莹 学 号 院、系、部 电气工程系 专 业 电子信息工程 指导教师 孙秀婷 康朝红 2011年 12 月25日 零输入响应与零状态响应 王莹 一、设计目的 进一步理解拉普拉斯变换及其性质和运算，学会用Matlab，并通过这一次的设计，进一步提高自己的实践动手能力。掌握信号经过LTI系统的时域分析方法。?巩固已经学过的知识，加深对知识的理解和应用，加强学科间的横向联系，学会应用MATLAB对实际问题进行仿真。 二、设计要求（1）根据实际问题建立系统的数学模型； （2）用MATLAB描述此系统； （3）仿真实现并绘制输出信号的波形。要求用两种方法仿真实现完全响应。对仿真结果进行比较，并与理论值比较。 题目要求 1、已知离散时间系统的差分方程为：2y(n) - y(n-1) - 3y(n-2)=2x(n) - x(n-1)，x(n)= u(n) ，y(-1)=1，y(-2)=3 ， 求系统的零输入响应、零状态响应和全响应 2、 初始条件为,； 四、题目求解 1、解：将差分方程Z变换得： …………………………………….(1) 依题意有：x(-1)=0,x(-2)=0,y(-1)=1,y(-2)=3 ,X(z)= 将上式变形如下： ………..(2) …………………………….(3) 易得系统函数为H(z)= ① 零输入时 零输入时，x(n)=0,差分方程右边为0，z变换后应为 = = 将Y(z)进行Z反变换，得到其零输入响应为： y(n)= ② 零状态时 零状态时，将y(-1)=0,y(-2)=0代入上面的式（2）中，得 Y(z)= X(z)= = = 将其Z反变换，得到零状态响应为： y(n)= ③ 全响应 与上面同理，y(-1)=1,y(-2)=3 将上面式(3)变形得： Y(z)= = Z反变换得全响应为 Y(n)= 五、由所解结果可得编写程序： ⒈fen=[2 -1 0]; %系统函数分子的系数 mu=[2 -1 -3]; %系统函数分母的系数 n=0:50; nl=length(n); %求零输入响应 y01=[1 3]; %y的初始状态 x01=[0 0]; %x 的初始状态 x1=zeros(1,nl); %求1*n1的全0矩阵 zi1=filtic(fen,mu,y01,x01); %为filter函数准备初始值 y1=filter(fen,mu,x1,zi1); %求零输入响应 subplot(311); %将当前窗口分裂成三部分，在1处绘制图形 stem(n,y1,r.); title(零输入响应); grid on; %求零状态响应 y02=[0 0]; x02=[0 0]; x2=0.5.^n; zi2=filtic(fen,mu,y02,x02); y2=filter(fen,mu,x2,zi2); subplot(312); %在2处绘制图形 stem(n,y2,r.); title(零状态响应); grid on; %求全响应 y03=[1 3]; x03=[0 0]; x3=0.5.^n; zi3=filtic(fen,mu,y03,x03); y3=filter(fen,mu,x1,zi3); subplot(313); %在3处绘制图形 stem(n,y3,r.); title(全响应); grid on; ⒉a=[6,4,5]; n=length(a)-1; Y0=[1,5]; b=[1,3,4]; dt=0.01;te=6; t=0:dt:te; u=sin(3*t)+cos(2*t); p=roots(a);V=rot90(vander(p));c=V\Y0; y1=zeros(1,length(t)); for k=1:n y1=y1+c(k)*exp(p(k)*t);end te=t(end); dt=te/(length(t)-1); [r,p,k]=residue(b,a); h=r(1)*exp(p(1)*t)