自动控制原理的MATLAB仿真与实践第3章 线性系统的时域分析.pptVIP

3）绘制不同参数时的单位阶跃响应uc(t)。 clear?; E=1; wn=3; t=0:0.02:6 ; %设置参数E和wn, 时间采样点t et=[0.1,0.2,0.5,0.95] ; %设置4个不同的阻尼向量et for i=1:4 ; %对et循环计算 uc=E*wn^2*(1/wn^2+1/(4*et(i)^2*wn^2-4*wn^2)^(1/2)*(1/(-et(i)*wn+1/2*(4*et(i)^2*wn^2-4*wn^2)^(1/2))*exp((-et(i)*wn+1/2*(4*et(i)^2*wn^2-4*wn^2)^(1/2)).*t)-1/(-et(i)*wn-1/2*(4*et(i)^2*wn^2-4*wn^2)^(1/2)).*exp((-et(i)*wn-1/2*(4*et(i)^2*wn^2-4*wn^2)^(1/2)).*t))) ; plot(t,uc,k-); hold on ; %绘制响应曲线并保持图形 end; * title(uc(t)的二阶欠阻尼响应曲线) gtext(et=0.1) %在图中对应的响应曲线的合适位置 %点击鼠标左键添加阻尼参数标示et=0.1 gtext(et=0.2) gtext(et=0.5) gtext(et=0.95) xla