《(改进)上机实验6连续系统的复频域分析.doc》由会员分享，可在线阅读，更多相关《(改进)上机实验6连续系统的复频域分析.doc（9页珍藏版）》请在金锄头文库上搜索。
1、上机实验六 连续系统的复频域分析一、 实验目的（1） 了解连续系统复频域分析的基本实现方法；（2） 掌握相关调用格式及作用。二、 实验原理复频域分析法主要有两种，即留数法和直接的拉普拉斯变换法，利用MATLAB进行这两种分析的基本原理如下。1. 基于留数函数的拉普拉斯变换设LTI系统的传递函数为若的零极点分别为，则可以表示为利用MATLAB的residue函数可以求解。2. 直接的拉普拉斯变换法经典的拉普拉斯变换分析法，即先从时域变换到复频域，在复频域经国处理后，在利用拉普拉斯反变换从复频域变换到时域，完成对时域问题的求解，涉及的函数有laplace函数和ilaplace函数等。三、 涉及的M
2、ATLAB函数1. residue函数功能：按留数法，求部分公式展开系数。调用格式：r,p,k,=residue(num,den):其中num,den分别是多项式系数按降序排列的行向量。2.Laplace函数功能：用符号推理求解拉氏变换。调用格式：L=laplace（F）：F为函数，默认为变量t的函数，返回L为s的函数。在调用函数时，要用syms命令定义符号变量t；3llaplace函数功能：用符号推理求解拉氏变换。调用格式：L=llaplace（F）4.Ezplot函数功能：用符号型函数的绘图函数。调用格式：Ezplot(f):f为符号型函数Ezplot(f,min,max):可指定横轴范围
3、Ezplot(f,xmin,xmax,ymin,ymax):可指定横轴范围和纵轴范围Ezplot(x,y):绘制参数方程的图像，默认x=x(t),y=y(t),0t clear;syms t1;syms s;a=1 2;b=1 4 3;sys=tf(a,b);t=0:0.1:5;y1=impulse(sys,t);y2=step(sys,t);y3=laplace(cos(20*t1);y4=laplace(exp(-t1);Hs=sym(s+2)/(s2+4*s+3);y3=y3*Hs;y4=y4*Hs;y3=ilaplace(y3);y4=ilaplace(y4);y3=vpa(y3,2)
4、;y4=vpa(y4,2);subplot(1,3,1);plot(t,y1,t,y2);text(2,0.6,阶跃响应);text(1,0.2,冲激响应);subplot(1,3,2);ezplot(y3,0,2);hold on;ezplot(cos(20*t),0,2);axis(0 1 -1.5 1.5);text(0.2,1.2,Input Voltage);text(0.2,-0.07,Output Voltage);subplot(1,3,3);ezplot(y4,0,10);hold on;ezplot(exp(-t),0,10);axis(0 10 0 1.5);text(1
5、.6,0.5,Input Voltage);text(2,0.25,Output Voltage);2. 若系统的传递函数为 试确定其零点，画出零极点分布图，确定其阶跃响应。clear;a=2 0 1 -3 1 4;b=5 2 -6 -3 5 2 -4 2 -1;subplot(2,2,1);zplane(a,b);legend(零点,极点);t=0:0.1:5;sys=tf(a,b);y=step(sys,t);subplot(2,2,2);plot(t,y);title(阶跃响应);xlabel(时间(t);ylabel(y(t);3.若某系统的传递函数为试确定其零极点，画出零极点分布图，
6、确定其冲击相应clear;a=1.65 -0.331 -576 90.6 19080;b=1 0.996 463 97.8 12131 8.11 0;subplot(2,2,3);zplane(a,b);legend(零点,极点);t=0:0.1:10;sys=tf(a,b);y=impulse(sys,t);subplot(2,2,4);plot(t,y);title(冲激响应);xlabel(时间(t);ylabel(y(t);六、 实验要求1在计算机中输入程序，验证并记录实验结果，并将实验结果存入指定存储区域。2对于设计性实验，应自行编制完整的实验程序，重复验证实验的过程，并在实验报告中给出完整的自编程序。实现对信号的模拟，并得出试验结果。3在实验报告中写出完整的MATLAB程序，并给出实验结果。
