自动控制lab2MATLAB用于时域分析.docxVIP

MATLAB用于时域分析王少丹通信1301201308030104实验目的在系统的时域分析中，可利用MATLAB完成系统的输出响应分析、稳定性分析、求动态性能指标以及误差分析等工作。实验原理在MATLAB中，可以通过单输入单输出系统的传递函数，进行系统的脉冲响应、阶跃响应、一般输入相应等时域分析。如果给定系统的传递函数为G(s)=，则其时域响应可以由以下函数得到。单位阶跃响应 step(num,den)或y=step(num,den,t)单位脉冲响应 impulse(num,den) 或 y=impulse(num,den,t)一般输入响应 y=lsim(num,den,u,t)注意，时间t是事先定义的矢量；u为输入信号。此外，还可以求出系统的超调量、调节时间以及稳态误差。实验内容一阶系统响应二阶系统响应若ω确定，系统的瞬态响应和ξ的取值有关。稳定性分析既可以通过对系统输出响应曲线的分析来判别系统的稳定性，也可以通过对系统闭环特征方程的求解来进行稳定性判别。求动态性能指标0型，稳态误差分析实例分析SIMULINK用于系统仿真实验代码%单位阶跃响应曲线sys1=tf([100],[1 0]);sys2=tf([0.1],[1]);sys=feedback(sys1,sys2);figurestep(sys)%二阶系统不同zeta值时的单位阶跃响应t=0:0.1:12;num=[1];zeta1=0;den1=[1 2*zeta1 1];zeta3=0.3;den3=[1 2*zeta3 1];zeta5=0.5;den5=[1 2*zeta5 1];zeta7=0.7;den7=[1 2*zeta7 1];zeta9=1.0;den9=[1 2*zeta9 1];[y1,x,t]=step(num,den1,t);[y3,x,t]=step(num,den3,t);[y5,x,t]=step(num,den5,t);[y7,x,t]=step(num,den7,t);[y9,x,t]=step(num,den9,t);figureplot(t,y1,t,y3,t,y5,t,y7,t,y9)grid on;%系统稳定性分析den=[1 2 8 12 20 16 16];roots(den)t=0:0.01:2;num=[1000]den=[1 34.5 1000];[y,x,t]=step(num,den,t);plot(t,y);%求超调量maxy=max(y);yss=y(length(t));pos=100*(maxy-yss)/yss%求峰值时间for i=1:1:201if y(i)==maxy,n=i;endendtp=(n-1)*0.01%求调节时间for i=1:1:201if(y(i)1.05y(i)0.95),m=i;break;endendts=(m-1)*0.01%求单位脉冲响应及稳态误差t=0:0.1:15;[num1,den1]=cloop([1],[1 1]);[num2,den2]=cloop([1],[1 1 0]);[num3,den3]=cloop([4 1],[1 1 0 0]);y1=impulse(num1,den1,t);y2=impulse(num2,den2,t);y3=impulse(num3,den3,t);figuresubplot(311);plot(t,y1);subplot(312);plot(t,y2);subplot(313);plot(t,y3);er1=0-y1(length(t))er2=0-y2(length(t))er3=0-y3(length(t))%求单位阶跃响应及稳态误差t=0:0.1:20;[num1,den1]=cloop([1],[1 1]);[num2,den2]=cloop([1],[1 1 0]);[num3,den3]=cloop([4 1],[1 1 0 0]);y1=step(num1,den1,t);y2=step(num2,den2,t);y3=step(num3,den3,t);figuresubplot(311);plot(t,y1);subplot(312);plot(t,y2);subplot(313);plot(t,y3);er1=0-y1(length(t))er2=0-y2(length(t))er3=0-y3(length(t))%求单位斜坡响应及稳态误差t=0:0.1:20;t1=0:0.1:100;[num1,den1]=cloop([1],[1 1]);[num2,den2]=cloop([1],[1 1 0]);[num3,den3]=cloop([4 1],[1 1 0 0]);y1=step(n