系统仿真matlab.doc

复习1 MATLAB命令窗口的使用, 输入并分析以下例子求导：syms a x;f=2*a*x2+5*x+6;diff(f);diff(f,x);diff(f,a)积分：syms a x;f=2*a*x+5；int(f);int(f,x);int(f,a);解方程: Syms x;f=x2+2*X+2;solve(f,x);x1,x2=solve(x1+2*x2=4,2*x1+3*x2=6)解微分方程：1）dsolve(D2y=1-y)2）dsolve(D2y=x+Dy, y(0)=1, Dy(0)=0, x)拉氏变换： syms ty=laplace(2*exp(-t)+1)反拉氏变换：syms sy=ilaplace(s+2)/(s2+4*s+3)2 在同一坐标系中分别绘制两个函数曲线。（至少有三种方法） 在不同窗口分别绘制两个函数曲线 （1）y1=sin2t+cost （2）y2=sin2t 3 输入以下传递函数，并求传函零极点。4 系统结构图如下：各环节传递函数为：求系统的开环传递函数和闭环传递函数。5 已知某单位反馈系统的开环传递函数为要求：绘制系统带栅格线的闭环根轨迹，并确定使系统产生重实根和纯虚根的开环增益k。6. 已知某系统的开环传递函数为：要求（1）绘制系统的奈奎斯特曲线，判断闭环系统的稳定性，求出系统的单位阶跃响应。（2）如系统稳定，求其稳定裕度、谐振振幅，谐振频率。1、syms a x;f=2*a*x2+5*x+6;diff(f)diff(f,x)diff(f,a)f1=2*a*x+5;int(f1)int(f1,x)int(f1,a)f2=x2+2*x+2;solve(f,x)dsolve(D2y=1-y)dsolve(D2y=x+Dy, y(0)=1, Dy(0)=0, x)syms t;y=laplace(2*exp(-t)+1)syms s;y=ilaplace(s+2)/(s2+4*s+3)2、t=0:pi/50:2*pi;y1=sin(2*t)+cos(t);y2=sin(2*t);figure(1)plot(t,y1)hold onfigure(2)plot(t,y2)3、num=3 0 9;den=conv(1 0 5,conv(1 2 1,conv(1 2,conv(1 2,1 2);z,p=tf2zp(num,den) 4、 G=tf(12 24 12 20,2 4 6 2 2);num1=5 3;den1=1 0;Gc=tf(num1,den1);den2=1 1000;H=tf(1000,den2);Pk=series(Gc,G);Gk=series(Pk,H)Gb=Pk/(1+Gk)5、num=1;den=conv(1 0,conv(0.01 1,0.02 1);rlocus(num,den)k poles=rlocfind(num,den)6、num=26;den=conv(1 6,1 -1);w=logspace(-1,2);G=tf(num,den);figure(1)nyquist(G)hold onfigu