机械工程控制基础上机实习指导书.docxVIP

第一篇 编程练习第一章 系统的时间响应分析利用MATLAB分析时间响应用MATLAB分析时间响应在MATLAB中，可以用impluse函数、step函数和lsim函数对线性连续系统的时间响应进行仿真计算。其中impluse函数用于生成单位脉冲响应；step函数用于生成单位阶跃响应；lsim函数用于生成任意输入的时间响应。有关的说明分别以图1.1.1、图1.1.2、图1.1.3所示。图1.1.1 impluse函数图1.1.2 step函数图1.1.3 lism函数【例】系统的传递函数为求系统在时间常数不同取值时的单位脉冲响应、单位阶跃响应、任意输入的时间响应。令，应用impluse函数可以生成单位脉冲响应；应用step函数可以生成单位阶跃响应；应用lsim函数可以生成任意输入的时间响应。【源代码1】t=[0:0.01:0.8]; %仿真时间区段%%nG=[50];tao=0;dG=[0.05 1+50*tao 50];G1=tf(nG,dG);tao=0.0125;dG=[0.05 1+50*tao 50];G2=tf(nG,dG); %三种tao值下的系统模型%tao=0.025;dG=[0.05 1+50*tao 50];G3=tf(nG,dG);%[y1,T]=impulse(G1,t);[y1a,T]=step(G1,t);[y2,T]=impulse(G2,t);[y2a,T]=step(G2,t); %系统响应%[y3,T]=impulse(G3,t);[y3a,T]=step(G3,t);%subplot(121),plot(T,y1,--,T,y2,-.,T,y3,-);legend(tao=0,tao=0.0125,tao=0.025)xlabel(t(sec)),ylabel(x(t)),grid on; %生成图形%subplot(122),plot(T,y1a,--,T,y2a,-.,T,y3a,-);legend(tao=0,tao=0.0125,tao=0.025)grid on;xlabel(t(sec)),ylabel(x(t));【生成的图形1】图1.1.4 单位脉冲响应曲线（左）和单位阶跃响应曲线（右）【源代码2】%t=[0:0.01:1];u=sin(2*pi.*t) %仿真时间区段和输入%%tao=0.025;nG=[50];dG=[0.05 1+50*tao 50];G=tf(nG,dG); %求系统模型%%y=lsim(G,u,t); %求系统响应%%plot(t,u,--,t,y,-,t,u-y,-.,linewidth,1)legend(u(t),xo(t),e(t)) %产生图形%grid on;xlabel(t(sec)),ylabel(x(t));【生成的图形2】图1.1.5 任意输入的时间响应及误差曲线图形利用MATLAB求系统的瞬态性能指标在求出系统的单位阶跃响应后，根据系统瞬态性能指标的定义，可以得到系统的上升时间、峰值时间、最大超调量和调整时间等性能指标。【源代码3】t=0:0.001:1;yss=1;dta=0.02; %设定仿真时间段和误差限%%nG=[50];tao=0;dG=[0.05 1+50*tao 50];G1=tf(nG,dG);tao=0.0125;dG=[0.05 1+50*tao 50];G2=tf(nG,dG);tao=0.025;dG=[0.05 1+50*tao 50];G3=tf(nG,dG); %计算三种时间常数下，系统的单位阶跃响应%y1=step(G1,t);y2=step(G2,t);y3=step(G3,t);%r=1;while y1(r)yss;r=r+1;end %上升时间%tr1=(r-1)*0.001;%[ymax,tp]=max(y1);tp1=(tp-1)*0.001; %峰值时间%%mp1=(ymax-yss)/yss; %最大超调量%%s=1001;while y1(s)1-dta y1(s)1+dta;s=s-1;end %调整时间%ts1=(s-1)*0.001;%r=1;while y2(r)yss;r=r+1;endtr2=(r-1)*0.001;[ymax,tp