第四次过程性考核

1、控制系统CAD 作业四控制系统CAD课程的平时作业共布置4次记分作业，每次作业满分2.5分，共10分。在完成记分作业的前提下，辅导教师还可根据学生学习的实际情况，再布置一些有针对性的作业，作为对记分作业的补充，供学生更好地掌握学习内容。一、填空题1称为系统的幅频特性，它描述系统对_不同频率特性_输入信号的稳态响应_幅值_衰减（或放大）的特性。称为系统的_相频_特性，它描述系统对不同频率输入信号的稳态响应及相位滞后或超前的特性。2若sys为系统传递函数，则在函数mag，phase，w=bode（sys）中函数的返回值mag和phase分别表示_频率特性的幅频值_和_相位_。3若sys为系统传递函

2、数，则函数margin（sys）的功能是_绘制系统的Bode图_，它和函数bode（sys）的区别在于所绘制的Bode图上带有稳定裕量标记_。4若sys为系统传递函数，则函数nyguist（sys）的功能是_绘制系统的Nyquist曲线_。5常用的串联校正有_相位超前校正、相位滞后校正、 和相位滞后-超前校正_三种方式。6某串联校正装置的传递函数为，该串联校正装置提供的是_超前（微分）_校正，它能改善系统的_稳定性和快速性_。7MATLAB提供的单输入单输出系统可视化设计工具SISO Design Tool是一种基于_根轨迹分析_或_频域性能指标_的系统设计工具。8欲用SISO Design

3、Tool对系统进行校正设计，需在MATLAB命令窗口输入_rltool(sys)_命令启动SISO Design Tool。二、简答题1已知单位负反馈系统的开环传递函数为，利用bode（）函数可分别绘制出K取不同值时系统的Bode图，如题图4-1所示，请根据Bode图分析系统开环增益K对系统稳定性的影响。系统的幅频特性曲线随着K的增大，位置提高，因此曲线2所对应的K相对较大，从图中还可知，曲线2所对应的幅值裕量约为-10（dB），相位裕量约为-23度，对应的闭环系统为不稳定。2根据上题的Bode图，简述用鼠标移动的方法求取稳定裕量的步骤。在Bode图曲线上任意一点用鼠标左键单击，曲线上便会标注

4、出“”，同时系统会显示单击点的对数幅频（或相频）值以及对应的频率。此时若按住鼠标左键并拖动图中的“”，则显示的数据随着“”位置的改变而变化。 在开环对数频率特性上，幅值裕量的定义为： g是相角 ()=-180°； 相位裕量的定义为：。（c是的频率），据此即可求出系统的幅值裕量和相位裕量。3在使用SISO Design Tool对系统进行校正设计时，需要调整校正装置提供的零、极点位置，系统提供了哪两种方法？简述各自的操作步骤。在使用SISO Design Tool对系统进行校正设计时，若要调整校正装置提供的零、极点位置，除了可以在显示系统根轨迹的“SISO Design for SIS

5、O Design Task”窗口用鼠标拖移零极点外，还可以再在“Control and Estimation Tools Manager”窗口中选择“Compensator Editor”选项卡，在Location文本框中能准确设定所添加的开环零点和极点的位置。三、编程题1若系统的单位阶跃响应为，试求系统的频率特性表达式提示：利用MATLAB符号运算功能，分别对输入、输出信号进行拉氏变换，并对结果进行简化，MATLAB中的拉氏变换函数为laplace（）。syms t s w;>> Y=1-1.8*exp(-4*t)+0.8*exp(-9*t);>> R=1;>&

6、gt; YS=laplace(Y,t,s);>> RS=laplace(R,t,s);>> GS=YS/RS;>> GSS=simple(GS);GSF=factor(GS);2已知系统的开环传递函数为，试分别计算当和时开环频率特性的幅频值和幅角值。clc;clf;clear>> w=0.5,2;>> num=10;den=conv(conv(1 0,2 1),1 0.5 1);>> G=tf(num,den);>> for j=1:length(w)H=freqresp(G,w(j);mag=abs(H);p

7、hase=angle(H)*180/pi;disp('w=',num2str(w(j)disp('幅频=',num2str(mag)3已知系统的开环传递函数为，试：（1）绘制对数幅频特性曲线以及相频特性曲线；（2）确定系统相位裕量、幅值裕量以及对应的截止频率和相位穿越频率。clc;clf;clear>> num=10;den=(conv(conv(1 0,0.1 1),0.04 1);Gk=tf(num,den)bode(Gk);grid>> Gm,Pm,Wg,Wc=margin(Gk)Gmdb=20*log10(Gm)4设负反馈系统的开

8、环传递函数为，试确定系统闭环稳定时的临界值。clc;clear den=conv(1 0,1 -10); t=0:0.01:100; for i=1:length(t) num=10*t(i) 1; Gk=tf(num,den); Gm,Pm,Wg,Wc=margin(Gk); if Gm=1|Pm=0 t(i) break end end 5已知单位反馈系统的开环传递函数为，试绘制出系统的闭环幅频特性曲线，并求出系统的谐振峰值和谐振频率。clc;clf;clearnum=5;den=conv(conv(1 0,0.1 1),1 1);Gk=tf(num,den)Gb=feedback(Gk,

9、1,-1)w=0:0.01:100;>> mag,phase=bode(Gb,w);>> M,n=max(mag);>> Mr=20*log10(M)wr=w(n)bode(Gb,w);grid6已知某单位反馈系统的开环传递函数为，校正装置的传递函数为，试绘制出系统校正前和校正后的对数幅频特性曲线以及校正装置的对数幅频特性曲线，并求出校正后的相位裕量。clc;clf;clearnum0=10;den0=conv(1 0,2 1);numc=conv(10 1,2 1);denc=conv(100 1,0.2 1);G0=tf(num0,den0);Gc=tf

10、(numc,denc);G=G0*Gc;Gm,Pm,Wg,Wc=margin(G)bode(G0,'b',G,'r',Gc,'g');gridgtext('校正前');gtext('校正后');gtext('校正装置')gtext('校正前');gtext('校正后');gtext('校正装置')7单位负反馈系统的开环传递函数为，若要求其相位裕量，试确定串联校正装置的传递函数（提示：串联超前校正）。clc;clf;clearnum0=1000;den0

11、=conv(1 0,conv(0.001 1,0.1 1);G0=tf(num0,den0);Gm0,Pm0,Wg0,Wc0=margin(G0);r=45;ep=10;r0=Pm0;p=(r-r0+ep)*pi/180;a=(1+sin(p)/(1-sin(p);mag0,phase0,w=bode(G0);mag20=20*log10(mag0);Mcm=-10*log10(a);Wn=spline(mag20,w,Mcm);T=1/(Wn*sqrt(a);numc=a*T 1;denc=T 1;Gc=tf(numc,denc);G=G0*Gc;Gm,Pm,Wg,Wc=margin(G)bode(G0,'b',G,'r',Gc,'g');grid;gtext('校正前');gtext('校正后');gtext('校正装置')gtext('校正前');gte