资源描述：
《信号与系统matlab实验3连续时间lti分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验三连续时间LTI系统分析姓名学号班级通信一班一、实验目的（一）掌握使用Matlab进行连续系统时域分析的方法1、学会使用符号法求解连续系统的零输入响应和零状态响应2、学会使用数值法求解连续系统的零状态响应3、学会求解连续系统的冲激响应和阶跃响应（二）掌握使用Matlab进行连续时间LTI系统的频率特性及频域分析方法1、学会运用MATLAB分析连续系统地频率特性2、学会运用MATLAB进行连续系统的频域分析（三）掌握使用Matlab进行连续时间LTI系统s域分析的方法1、学会运用MATLAB求拉普拉斯变换（LT）2、学会运用MATLAB求拉普拉斯反变换（ILT）3、学会在MATLAB环
2、境下进行连续时间LTI系统s域分析二、实验原理及实例分析（一）连续系统时域分析（详细请参见实验指导第二部分的第5章相关部分）（二）连续时间LTI系统的频率特性及频域分析（详细请参见实验指导第二部分的第8章相关部分）（三）拉普拉斯变换及连续时间系统的s域分析（详细请参见实验指导第二部分的第10、11章相关部分）12三、实验过程（一）熟悉三部分相关内容原理（二）完成作业已知某系统的微分方程如下：其中，为激励，为响应。1、用MATLAB命令求出并画出时系统的零状态响应和零输入响应（零状态响应分别使用符号法和数值法求解，零输入响应只使用符号法求解）；>>eq='D2y+3*Dy+2*y=0';>
3、>cond='y(0)=1,Dy(0)=2';>>yzi=dsolve(eq,cond);yzi=simplify(yzi);>>eq1='D2y+3*Dy+2*y=Dx+3*x';eq2='x=exp(-3*t)*Heaviside(t)';cond='y(-0.01)=0,Dy(-0.001)=0';yzs=dsolve(eq1,eq2,cond);yzs=simplify(yzs.y)yzs=heaviside(t)*(-exp(-2*t)+exp(-t))>>yt=simplify(yzi+yzs)yt=-3*exp(-2*t)+4*exp(-t)-exp(-2*t)*heavi
4、side(t)+exp(-t)*heaviside(t)12>>subplot(3,1,1);>>ezplot(yzi,[0,8]);gridon;>>title('rzi');>>subplot(3,1,2);>>ezplot(yzs,[0,8]);>>gridon;>>title('rzs');>>subplot(3,1,3);>>ezplot(yt,[0,8]);gridon;>>title('完全响应')sys=tf([1,3],[1,3,2]);t=ts:dt:te;f=exp(-3*t).*uCT(t);y=lsim(sys,f,t);plot(t,y),gridon;axi
5、s([0,8,-0.02,0.27]);12xlable('Time(sec)'),ylable('y(t)');title('零状态响应')2、使用MATLAB命令求出并画出系统的冲激响应和阶跃响应（数值法）；用卷积积分法求系统的零状态响应并与（1）中结果进行比较；t=0:0.001:4;sys=tf([1,3],[1,3,2]);h=impulse(sys,t);g=step(sys,t);subplot(2,1,1);plot(t,h),gridon;xlable('Time(sec)'),ylable('h(t)');title('冲激响应');subplot(2,1,2);12
6、plot(t,g),gridon;xlable('Time(sec)'),ylable('g(t)');title('阶跃响应')_dt=0.01;t1=0:dt:8;f1=exp(-3*t1);t2=t1;sys=tf([1,3],[1,3,2]);f2=impulse(sys,t2);[t,f]=ctsconv(f1,f2,t1,t2,dt)function[f,t]=ctsconv(f1,f2,t1,t2,dt)f=conv(f1,f2);f=f*dt;ts=min(t1)+min(t2);te=max(t1)+max(t2);t=ts:dt:te;subplot(1,1,1)pl
7、ot(t,f);gridon;axis([min(t),max(t),min(f)-abs(min(f)*0.2),max(f)+abs(max(f)*0.2)]);12title('卷积结果')3、使用MATLAB命令求出并画出此系统的幅频特性和相频特性；使用频域分析法求解系统的零状态响应并与（1）中结果进行比较；>>w=-3*pi:0.01:3*pi;b=[1,3];a=[1,3,2];H=freqs(b,a,w);subplo
