控制系统仿真及CAD试题（研2012）.docVIP

控制系统仿真及CAD试题（研2012） 1、4-2某反馈控制系统的开环传递函数为 试绘制其根轨迹。 解：在MATLAB命令窗口中输入下列命令： ng=1.0; dg=poly([0,-4,-2+4j,-2-4j]); rlocus(ng,dg) 运行结果为： 4-3已知某系统传递函数为 试绘制其伯德图。 解：分子分母同乘100*200得到 在Matlab窗口中输入下列命令： k=80*200; num=[1 100]; den=conv([2.5 100], [(1/200) 2*0.3 200]); w=logspace(-1,1,100); [m,p]=bode(k*num,den,w); subplot(2,1,1); semilogx(w,20*log10(m)); grid; xlabel(Frequency(rad/sec)); ylabel(Gain(dB)); subplot(2,1,2); semilogx(w,p); grid; xlabel(Frequency(rad/s)); ylabel(Phase(deg)); 可绘制该系统的伯德图如下所示。 4-4设控制系统具有如下的开环传递函数 试求取当K=10时的相角裕度和幅值裕度，并画出其伯德图。 解：在MATLAB命令窗口中输入下列命令： k=10; num=1; den=poly([0,-1,-5]); [m,p,w]=bode(k*num,den); subplot(2,1,1); semilogx(w,20*log10(m)); grid on; ylabel(Gain(dB)); subplot(2,1,2); semilogx(w,p); grid on; xlabel(Frequency(rad/sec)); ylabel(Phase(deg)); [gm,pm,wcg,wcp]=margin(m,p,w) 这里gm,wcg为幅值裕度值与相应的频率 pm,wcp为相角裕度值与相应的频率，运行结果为：gm =3.0000，pm = 25.4489，wcg =2.2361，wcp = 1.2241。因此，系统的幅值裕度和相角裕度分别为3dB和25.4489°。系统的伯德图如下所示。 4-13对于高阶系统的设计问题，往往要进行降阶近似处理，并要验证近似效果。已知某高阶系统模型为 经简化处理后，模型等效为 试比较两个模型在单位阶跃信号作用下的响应情况，并分析近似效果。 解：在Matlab命令窗口中输入下列命令： num1=[35,10861,13285,82402,278376,511812,422964,194480]; den1=[1,33,437,3017,11870,27470,37492,28880,9600]; sys1=tf(num1,den1); num2=[35,284.98]; den2=[1,10.114,12.31]; sys2=tf(num2,den2); step(sys1,-,sys2,:r,15) legend(原模型, 简化后模型) ; grid on 程序运行结果如下： 从曲线中可以看出，降阶后系统响应无超调，调整时间缩短，系统响应变快，但是简化前后终值有差异。 2-5用四阶龙格-库塔法求解题2-3数值解，并与前两题结果相比较。 解：四阶龙格-库塔法表达式，其截断误差为同阶无穷小，当h步距取得较小时，误差是很小的. 程序如下： h=0.1; disp(四阶龙格-库塔方法求解函数数值解为); disp(y=); y=1; for t=0:h:1 disp(y); k1=-y; k2=-(y+k1*h/2); k3=-(y+k2*h/2); k4=-(y+k3*h); y=y+(k1+2*k2+2*k3+k4)*h/6; end 得到结果 四阶龙格-库塔方法求解函数数值解为 y= 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.4066 0.3679 (2)比较这几种方法： 对于四阶龙格-库塔方法 真值 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.4066 0.3679 龙库 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493