上机实验6连续系统的复频域分析111.docVIP

上机实验6 连续系统的复频域分析 一、实验目的 （1）了解连续系统复频域分析的基本实现方法； （2）掌握相关函数的调用格式及作用。 二、实验原理 复频域分析方法主要有两种，即留数法和直接的拉普拉斯变换法，利用matlab进行两种分的基本原理如下。 1.基于留数函数的拉普拉斯变换法 设LTI系统的传递函数为 H(s)=B(s)/A(s) 若H(s)的零点分别为，…，，，…，，则可以表示为 H(s) H(s)= 利用matlab的residue函数可以求解。 2.直接的拉普拉斯变换法 经典的拉普拉斯变换分析法，即先从时域变换到复频域，在复频域经过处理之后，再利用拉普拉斯反变换从复频域换到时域，完成对时间问题的求解，涉及的函数有 Laplace和ilaplace函数等。 二、涉及的matlab函数 1.residue函数 功能：按留数法，求部分分式展开系数。 调用格式： [r,p,k]=residue(num,den);其中num，den分别是B(s)，A(s)多项式系数按降序排列的行向量。 2.laplace函数 功能：用符号推理求解拉普拉斯变换。 调用格式： L=laplace(F):F为函数，默认为变量t的函数，返回L为s的函数。在调用时，要用sym命令定义符号变量t； 3.ilaplace函数 功能：符号推理求解拉氏反变换。 调用格式：L=ilaplace(F) 4.ezplot函数 功能：用符号型函数的绘图函数。 ezplot(f):f为符号型函数。 ezplot(f,[min max]):可指定横轴范围 ezplot(f,[xmin xmax,ymin ,ymax]):可指定横轴范围和纵轴范围 ezplot(x,y):绘制参数方程的图像，默认x=x(t),y=y(t),0t2*pi; 5.roots函数 功能：求多项式的根； 调用格式： r=roots(c):其中c是多项式的系数向量(自高次到低次)，r为根向量，注意，matlab默认根为列向量。 四、实验内容与方法 1.验证性实验 （1）系统零极点求解. 已知H(s) ==(s^2-1)/(s^3+2s^2+3s+2),画出H(s)的零极点图。 Matlab程序1： b=[1 0 -1]; a=[1 2 3 2]; zs=roots(b); ps=roots(a); plot(real(zs),imag(zs),go,real(ps),imag(ps),mx,markersize,12); grid on legend(零点,极点); 运行结果： 直接调用zplane函数绘制： Matlab程序2： b=[ 1 0 -1]; a=[1 2 3 2]; zplane(b,a); grid on legend(零点,极点); 运行结果： （2）一个线性非时变电路的转移函数为： H(s)==10^4(s+6000)/(s^2+875s+88*10^6) 若=12.5cos(8000t)V,求的稳态响应。 ?稳态滤波法求解。 MATLAB程序： w=8000; s=j*w; num=[0,1e4,6e7];den=[1,875,88e6]; H=polyval(num,s)/polyval(den,s); mag=abs(H) mag = 4 phase=angle(H)/pi*180 phase = 36.8699 t=2:1e-6:2.002; vg=12.5*cos(w*t); vo=12.5*mag*cos(w*t+phase*pi/180); plot(t,vg,t,vo);grid; text(0.25,0.85,output voltage,sc); text(0.07,0.35,input voltage,sc); title(稳态滤波输出); ylabel(电压(v)),xlabel(时间（s)); 系统的稳态响应如右图： （3）拉氏变换法求解。 MATLAB程序： Syms t m; Hs=sym((10^4*(s+6000))/(s^2+875*s+88*10^6)); Vs=laplace(12.5*cos(8000*t));Vos=Hs*Vs; Vo=ilaplace(Vos); V0=vpa(Vo,4); ezplot(Vo,[1,1+5e-3]);hold on; ezplot(12.5*cos(8000*t),[1,1+5e-3]);axis([1,1+2e-3,-50,50]); grid; 得到稳态