信号与系统matlab课后作业_北京交通大学分析.docx

信号与系统MATLAB平时作业学院：电子信息工程学院班级：姓名：学号：教师：钱满义MATLAB习题M3-1一个连续时间LTI系统满足的微分方程为y’’(t)+3y’(t)+2y(t)=2x’(t)+x(t)(1)已知x(t)=e-3tu(t)，试求该系统的零状态响应yzs(t)；(2)用lism求出该系统的零状态响应的数值解。利用(1)所求得的结果，比较不同的抽样间隔对数值解精度的影响。解：由于则将带入原方程式化简得所以又因为则该系统的零状态响应(2)程序代码1、ts=0;te=5;dt=0.1;sys=tf([2 1],[1 3 2]);t=ts:dt:te;x=exp(-3*t).*(t=0);y=lsim(sys,x,t)2、ts=0;te=5;dt=1;sys=tf([2 1],[1 3 2]);t=ts:dt:te;x=exp(-3*t).*(t=0);y1=-0.5*exp(-3*t).*(exp(2*t)-6*exp(t)+5).*[t=0];y2=lsim(sys,x,t)plot(t,y1,r-,t,y2,b--)xlabel(Time(sec))legend(实际值,数值解)用lism求出的该系统的零状态响应的数值解在不同的抽样间隔时与（1）中求出的实际值进行比较表格 1 抽样间隔为1序号12345实际值0.0976-0.0189-0.0178-0.0082-0.0032数值解0.442090.000570.003970.000990.00016表格 2 抽样间隔为0.1间际值0.15170.22960.25960.25980.24250.21590.1854数值解0.15290.23130.26160.26180.24440.21760.1868891011121314150.15420.12460.09760.07380.05320.0360.02160.010.15540.12550.09830.07430.05360.03620.02180.010116171819202122230.0008-0.0065-0.012-0.016-0.0189-0.0208-0.022-0.0230.0008-0.0065-0.0121-0.0162-0.0191-0.021-0.0221-0.0232425262728293031-0.0225-0.0222-0.0216-0.0208-0.0199-0.0188-0.0178-0.017-0.0227-0.0224-0.0218-0.021-0.02-0.019-0.0179-0.0173233343536373839-0.0156-0.0145-0.0134-0.0124-0.0115-0.0106-0.0097-0.009-0.0157-0.0146-0.0135-0.0125-0.0116-0.0106-0.0098-0.0094041424344454647-0.0082-0.0075-0.0068-0.0062-0.0057-0.0052-0.0047-0.004-0.0082-0.0075-0.0069-0.0063-0.0057-0.0052-0.0048-0.004484950-0.0039-0.0036-0.0032-0.0039-0.0036-0.0033将两种结果画在同一幅图中有图表 1 抽样间隔为1图表 2 抽样间隔为0.1图表 3 抽样间隔为0.01当抽样间隔dt减小时，数值解的精度越来越高，从图像上也可以看出数值解曲线越来越逼近实际值曲线，直至几乎重合。M3-4conv函数可计算起点为k=0的两个序列的卷积。利用此函数和卷积的性质，编写计算非零起点任意两个序列的卷积的程序，并利用该程序计算下面两个序列的卷积。,解：程序代码m=-2; %序列x最左边的位置n=-1; %序列h最左边的位置x=[0.85 0.53 0.21 0.67 0.84 0.12];h=[0.68 0.37 0.83 0.52 0.71];z=conv(x,h)N=length(z);stem(m+n:N+m+n-1,z)grid图表 4 树根图得出序列的卷积M3-8某离散时间LTI系统满足差分方程试利用impz函数求其单位脉冲响应，并画出前30点的图。解：程序代码k=0:30;a=[1 0.7 -0.45 -0.6];b=[0.8 -0.44 0.36 0.02];h=impz(b,a,k);stem(k,h);grid解得单位脉冲响应为图表 5 单位脉冲响应k012345h[k]0.8-11.42-0.9440.6998-0.06276789101112-0.20760.537-0.50690.4719-0.23