M1.2 Matlab时域分析.pptVIP

控制系统的计算机辅助分析 主要内容： 一、概述 二、控制系统的数学描述与建模 三、模型的转换 四、模型连接和化简 五、控制系统的时域分析 控制系统的分析和设计流程 用Matlab绘制时间响应曲线 已知 MATLAB绘图命令 MATLAB提供了丰富的绘图功能 help graph2d可得到所有画二维图形的命令 help graph3d可得到所有画三维图形的命令 下面介绍常用的二维图形命令 6、文字标示 text(x,y,’字符串’) 在图形的指定坐标位置(x,y)处，标示单引号括起来的字符串。 gtext(‘字符串’) 利用鼠标在图形的某一位置标示字符串。 title(‘字符串’) 在所画图形的最上端显示说明该图形标题的字符串。 xlabel(‘字符串’)，ylabel(‘字符串’) 设置x，y坐标轴的名称。 输入特殊的文字需要用反斜杠（\）开头。 7、legend(‘字符串1’,‘字符串2’,…,‘字符串n’) 在屏幕上开启一个小视窗，然后依据绘图命令的先后次序，用对应的字符串区分图形上的线。 例exp2_5_.m 9、semilogx：绘制以x轴为对数坐标（以10为底），y轴为线性坐标的半对数坐标图形。 semilogy：绘制以y轴为对数坐标（以10为底），x轴为线性坐标的半对数坐标图形。 例exp2_17.m 用MATLAB解微分方程 MATLAB提供了ode23、ode45等微分方程的数值解法函数，不仅适用于线性定常系统，也适用于非线性及时变系统。 Lorenz chaotic attractor function xdot = lorenzeq(t,x) xdot = [ -8/3*x(1)+x(2)*x(3); -10*x(2)+10*x(3); -x(1)*x(2)+28*x(2)-x(3)]; %微分方程函数 %格式 % function xdot=filename(t,x) % xdot=[表达式1；表达式2；表达式3；...；表达式n-1] % 表达式1 对应 x1=x2 % 表达式2 对应 x2=x3 % 表达式3 对应 x3=x4 % ... % 表达式n-1 对应 xn-1=xn % 如例exp3_1.m % x(1)=Vo x(2)=iL % x(1)=x(2)/C % x(2)=(Vs-x(1)-R*x(2))/L t_final=100; x0=[0; 0; 1e-10]; [t,x] = ode45(lorenzeq,[0, t_final], x0); figure(1); plot(t,x); set(gcf,position,[7 305 501 380]) figure(2); plot3(x(:,1), x(:,2), x(:,3)); axis([10 40 -20 20 -20 20]) set(gcf,position,[516 305 501 380]) plot(t,x) plot3(x(:,1), x(:,2), x(:,3)); 例exp3_1.m  %微分方程函数 function xdot=elecsys(t,x) %状态导数 Vs=1; R=1.4; L=2; C=0.32; xdot=[x(2)/C;1/L*(Vs-x(1)-R*x(2))]; clear clc close t0=0; tfinal=15; x0=[0.5;0]; %初始化，电感电流为0，电容电压为0.5v %tol=0.001; %数值计算精度 [t,x]=ode45(elecsys,t0,tfinal,x0); %elecsys是系统微分方程的描述函数 figure(1) subplot(211) plot(t,x(:,1)) title(capacitor voltage) xlabel(time-sec) subplot(212) plot(t,x(:,2)) title(current of L) xlabel(time-sec) figure(2) vc=x(:,1); i=x(:,2); plot(vc,i); title(current versus capacitor voltage) xlabel(capacitor voltage) ylabel(current) Laplace 变换 inverse Laplace transform syms a s t w x g = 1/s^2 ilaplace(g) Returns t h =1/(s^2+4)^2 ilaplace(h) ans = 1/16*si