系统时域分析的仿真实验.doc

实验名称：系统时域分析的仿真实验 一、实验目的 1、掌握Matlab求解LTI系统响应的方法。 2、会用simulink建立系统模型。 二、上级实验 1、连续系统 程序1： num=1; den=[1,2,2,1]; G=tf(num,den); %多项式模型 Str_title=阶跃响应冲激响应零输入响应零状态响应; %字符串数组，调用为标题 for k=1:4 switch k case 1 subplot(221); %阶跃响应 step(G) title(str_title(1:4)); case 2 %冲激响应 subplot(222); impulse(G) title(str_title(5:8)); case 3 %零输入响应 subplot(223); sys=ss(G); %转为状态空间模型 x0=[1,1,1]; %初始状态 initial(sys,x0) title(str_title(9:13)); case 4 %零状态响应 subplot(224); [u,t]=gensig(square,5,20,0.01); % 方波square 正弦波sin lsim(G,u,t) title(str_title(14:18)); end end xlabel(时间); ylabel(幅值); 程序2： syms t; %定义函数变量 x(t)=sin(0.1*t)+sin(0.5*t)+sin(10*t);%输入信号的函数 Lap=laplace(x(t)); %拉氏变换 syms s; %定义函数变量 Hs=1/(s^3+2*s^2+2*s+1); %系统函数 Y=Hs.*Lap; %频域下的输出 y=ilaplace(Y,s,t); %时域下的输出 subplot(211); fplot(y,[0,50]); %函数图像显示 title(拉氏变换求系统响应); xlabel(时间); ylabel(幅值); num=1; den=[1,2,2,1]; G=tf(num,den); %多项式模型 t=0:0.01:50; u=sin(0.1*t)+sin(0.5*t)+sin(10*t); subplot(212); lsim(G,u,t); axis([0 50 -2.5 2.5]); title(系统响应); xlabel(时间); ylabel(幅值); 图形1： 图形2： 2、离散系统 程序1： num=1； den=[2 -1 -1]; G=tf(num,den,0.1,Variable,z); N=50; %作图点数 n=0:N-1; hn=impz(num,den,n); %冲激函数 subplot(221); stem(n,hn,filled); title(冲激响应) grid on; gn=stepz(num,den,n); %阶跃函数 subplot(222); stem(n,gn,filled); title(阶跃响应) grid on; n=0:50; nl=length(n); y01=[0 1]; %y的初始状态 x01=[0 0]; %x的初始状态 x1=zeros(1,nl);