matlab实验大作业,matlab实验的基本操作大作业

matlab实验的基本操作大作业

MATLAB仿真 期末大作业 姓 名： 龚文娟 班 级： 0308412 学 号： 030841203 指导教师： 李绍武 2011春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由串联，其中： A表示自己学号最后一位数(可以是零)，K为开环增益。要求： (1) 设K=1时，建立控制系统模型，并绘制阶跃响应曲线(用红色虚线，并标注坐标和标题)；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； (2) 在第(1)问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace模块将命令窗口的阶跃响应数据导入Simulink模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink模型窗口绘制阶跃响应曲线，用out1或者to workspace模块将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 (3)用编程法或者rltool法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，位置稳态误差小于等于0.1rad，开环系统截止频率，相角裕度大于等于45度，幅值裕度大于等于10dB。 MATLAB分析： (1) K=1时，在命令窗口建立系统模型，用step函数绘制阶跃响应曲线。 命令窗口的程序： num=[4]; den=[1 1 0] G=tf(num,den) sys=feedback(G,1) t=[0:0.01:10]; step(sys, r: ); %用红色虚线 xlabel( t ); ylabel( y ); title( 阶跃响应曲线 ) % 标注坐标和标题 [y,t]=step(sys) 程序运行结果如下： 用游离鼠标法求取上升时间、超调量、调节时间、峰值时间 上升时间= 0.865-0.232=0.633s 超调量=(1.44-1)/1*100%=44% 调节时间=5.38s 峰值时间=1.66s (2) 用from workspace模块将命令窗口的阶跃响应数据导入Simulink模型窗口，用示波器显示阶跃响应曲线。 (3) 用编程法设计串联校正环节，首先计算超前校正参数的子函数的程序代码如下： function Gc=cqjz(G,kc,yPm) %超前校正 G=tf(G); [mag,phase,w]=bode(G*kc); Mag=20*log10(mag); [Gm,Pm.Wcg,Wcp]=margin(G); fam=(yPm-getfield(Pm, Wcg ))*pi/180; alf=(1+sin(fam))/(1-sin(fam)) Mn=-10*log10(alf); Wcgn=spline(Mag,w,Mn); T=1/(Wcgn*sqrt(alf)) Tz=alf*T; Gc=tf([Tz 1],[T 1]); 其次，主函数MATLAB代码如下： num=1; den=conv([1,0],[1,1]); G=tf(num,den) kc=10; yPm=45+12 Gc=cqjz(G,kc,yPm) G=G*kc; GGc=G*Gc Gyclose=feedback(G,1) Gxclose=feedback(GGc,1) figure(1); step(Gxclose, b: ); hold on step(Gyclose, r ); grid gtext( 校正前的 ); gtext( 校正后的 ); figure(2); bode(G, b: ); hold on bode(GGc, r ); grid gtext( 校正前的 );gtext( 校正后的 );gtext( 校正前的 );gtext( 校正后的 ); figure(3); nyquist(Gxclose, b: ); hold on nyquist(Gyclose, r ); grid gtext( 校正前的 );gtext( 校正后的 ); 运行结果如下： %超前校正环节传递函数 %校正后系统闭环传递函数 系统的校正前后的单位阶跃响应曲线和校正前后的开环Bode图： 系统校正前后的闭环Nyquist图如下所示： 最后，验证已校正系统的相角裕度和幅值裕度。 γ’=39.2979°+ 17.9642°=57.2621°＞45° 幅值裕度仍为＋∞dB。从而满足此系统的性能要求。