信号与系统实验报告2.docxVIP

实验二连续线性时不变系统分析姓名罗治民班级电子12--BF学号 14122502289日期2014.6.17成绩一、实验目的：深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义。掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。实验条件：计算机一台，MATLAB软件。三、实验原理：1.连续系统的时域响应（1）单位冲激响应h(t)impulse(sys)计算系统的冲激响应。sys由函数tf(b,a)。h=impulse(sys, t)；微分方程：该系统的单位冲激响应h(t)。a=[1,5,6]; b=[2,8];sys=tf(b,a);t=0:0.1:10;h=impulse(sys,t);plot(h);xlabel(t); title(h(t));（2）单位阶跃响应g(t)step(sys)计算系统的阶跃响应。sys由函数tf(b,a)，g=step(sys, t)；微分方程：该系统的阶跃响应响应h(t)。a=[1,5,6]; b=[2,8];sys=tf(b,a);t=0:0.1:10;g=step(sys,t);plot(g);xlabel(t); title(g(t));(3) 零输入响应yzi(t):dsolve（’equ1’，’equ2’，…）计算零输入响应。微分方程：输入为y(0-)=1，y’(0-)=1,时系统的零输入响应。MATLAB程序：yi=dsolve(D2y+3*Dy+2*y=0,y(0)=1,Dy(0)=1);t=0:0.1:10；ezplot(yi,t); title(yi); ylabel(yi);（4）零状态响应yzs(t)lsim(sys, x, t)系统的零状态响应。参数: sys可由函数tf(b,a)获得x为输入信号，t为定义的时间向量。微分方程：，输入为时系统的零状态响应。a=[1,5,6]; b=[2,8];sys=tf(b,a);t=0:10/300:10;x=exp(-t);y=lsim(sys,x,t); plot(t,y);（5）连续系统的全响应y(t)的计算 全响应：MATLAB提供了专门用于求解连续系统这些响应的函数dsolve()，其调用格式为：dsolve（’equ1’，’equ2’，…）2．连续系统的系统函数零极点分析：MATLAB中提供了roots函数计算系统的零极点，提供了pzmap函数绘制连续系统的零极点分布图。系统函数为：计算其零极点，画出分布图。b=[2,3,1];a=[1,2,2,1];z=roots(b)；p=roots(a)；sys=tf(b,a);pzmap(sys)；3．连续系统的频率响应MATLAB中freqs函数可以分析连续系统的频响，格式如下：H=freqs(b,a,w)；计算系统在指定频率点向量w上的频响H；w为频率点向量。[H,w]=freqs(b,a)；%自动选取200个频率点计算频率响应。已知某连续系统的系统函数为：分析系统的幅频率特性。b=[1]; a=conv([1,1],[1,1,1]);[H,w]=freqs(b,a);plot(w,abs(H));xlabel(Frequency(rad/s));ylabel(Amplitude);title(Magnitude response);四、实验步骤及结果测试：（包括程序、结果图）1．描述某线性时不变系统的微分方程为：且f(t)=t2，y(0-)=1，y’(0-)=1；试求系统的单位冲激响应、单位阶跃响应、全响应、零状态响应、零输入响应、自由响应和强迫响应。编写相应MATLAB程序，画出各波形图。解：冲激响应：t=0:0.1:10; b=[1,2]; a=[1,3,2];sys=tf(b,a); h=impulse(sys,t);plot(t,h);xlabel(t);title(h(t));阶跃响应：t=0:0.1:10; b=[1,2]; a=[1,3,2]; sys=tf(b,a); g=step(sys,t);plot(t,g);xlabel(t);title(g(t));全响应：t=0:0.1:10;yt=dsolve(D2y+3*Dy+2*y=2*t+2*t^2,y(0)=1,Dy(0)=1);ezplot(yt,t);xlabel(t);ylabel(y) ;title(y(t));零状态响应：a=[1,3,2];b=[1,2];sys=tf(b,a);t=0:10/300:10;x=square(t);ys=lsim(sys,x,t);plot(t,ys);xlabel(t);title(ys(t))零输入响应：t=0:0.1:10;yit=dsolve(D2y+3*Dy+2*y=0,y(0