完整word版控制系统CAD课程形考作业四参考答案

1、控制系统CAD形成性考核册-12 -控制系统CAD作业4姓名：学号：得分：教师签名：、填空题1. A（豹）=|G（jG）|称为系统的幅频特性，它描述系统对Xo输入信号的稳态响应衰减（或放大）的特性。蚁=NG（j称为系统的特性，它描述系统对不同频率输入信号的稳态响应，相位滞后（* 0）的特性。（不同频率、幅值、相频特性）2.若sys为系统传递函数，在函数mag,phase,w=bode（sys）中，函数的返回值 mag和phase分别表示。（频率特性的幅频值、相位）3.若sys为系统传递函数，函数margin （sys）的功能是:，它和函数bode（sys）的区别在于:。（绘制系统的 Bode图

2、、所绘制的Bode图上带有稳定裕量标记）4.若sys为系统传递函数，函数 nyguist （sys）的功能是:（绘制系统的Nyquist曲线）二、简答题1. 已知单位负反馈系统的开环传递函数为:G（s）二gs+DG+S，利用bodeO函数可分别绘制出K取不同值时系统的 Bode图如题图7-1所示，请根据Bode图，分析系统开环增益K对系统稳定性的影响。Bode Diagram501Frequency (rad/sec)P题图7-1 不同K取值下的系统Bode图（系统的幅频特性曲线随着K的增大，位置提高，因此曲线 2所对应的K相对较大，从图中还可知，曲线 2所对应的幅值裕量约为-10（ dB），

3、相位裕量约为-23度，对应的闭环系统为不稳定。）2. 在上题的Bode图中，简述用鼠标移动的方法求取稳定裕量的步骤。在Bode图曲线上任意一点用鼠标左键单击，曲线上便会标注出“ 同时系统会显示单击点的对数幅频（或相频）值以及对应的频率。此时若按住鼠标左键并拖动图中的“；则显示的数据随着“ ”位置的改变而变化。在开环对数频率特性上，幅值裕量的定义为:h=-20lg|Gk(j%)|，(Q 是相角 M3)=-180o 的频率)相位裕量的定义为：Y=Wc）（兀）。（3c是20lg|Gk（jeo）|=0的频率），据此即可求出系统的幅值裕量和相位裕量二、编程题1.若系统的单位阶跃响应为:y（t） =1 1

4、.8ed +0肌耳（t AO），试求系统的频率特性表MA TLAB达式(提示：利用MATLAB符号运算功能，分别对输入、输出信号进行拉氏变换， 并对结果化简, 中的拉氏变换函数为laplace()。syms t s w;Y=1-1.8*ex p(-4*t)+0.8*ex p(-9*t);R=1;YS=la place(Y,t,s);RS=la place(R,t,s);GS=YS/RS;GSS=sim ple(GS);GSF=factor(GS);GW=subs(GSF,s,i*w) dis p(w=,num2str(w(j)36即：G(J)=：(9 + )(4 +血)102.已知系统的开环传

5、递函数为G(s)H(S)二心 4、/ 2 4、，试分别计算当s(2s+1)(s +0.5s+1)GW =36/(i*w + 9)*(i*w + 4)3=0.5和3=2时开环频率特性的幅频值A( 3和幅角值 y 3。clc;clearw=0.5,2;3 =0.5、2num=10;den=conv(conv(1 0,2 1),1 0.5 1);G=tf(num,den);求得传递函数for j=1:length(w)H=freqres p(G,w(j);求得频率响应mag=abs(H);求复数H的模phase=angle(H)*180/pi;求复数H的幅角并把单位转化为度dis p(幅频=,num

6、2str(mag)dis p(幅角=,num2str( phase)endw=0.5幅频=17.8885幅角=-153.4349w=2幅频=0.38348幅角=32.47123.已知系统的开环传递函数为:10Gk(s)=k s(0.1s+1)(0.04s+1) o试：（1）绘制对数幅频特性曲线以及相频特性曲线;（2）确定系统相位裕量 丫和幅值裕量h以及对应的截止频率 oc和相位穿越频率COgoclc;clf;clearnum=10;den=conv(conv(1 0,0.1 1),0.04 1);开环传递函数分母多项式系数Gk=tf(num,den)求得开环传递函数bode(Gk);gridG

7、m, Pm,Wg,Wc=margin(Gk)求得系统的稳定裕量以及对应的频率Gmdb=20*log10(Gm)将幅值裕量转化为分贝值Gm = 3.5000Pm = 35.7873Wg = 15.8114Wc = 7.6120Gmdb = 10.88144.设负反馈系统的开环传递函数为:Gk(s)10(Ts+1)s（s10），试确定系统闭环稳定时T的临界值。clc;clearden=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=marg

8、in(Gk);endif Gm=1| Pm=0endans=1t(i)break5.已知单位反馈系统的开环传递函数为:Gk(S2s(0.1s+1)(s+1)，试绘制出系统的闭环幅频特性曲线，并求出系统的谐振峰值和谐振频率。clc;clf;clear求取开环传递函数num=5;den=conv(conv(1 0,0.1 1),1 1);Gk=tf(num,den)Gb=feedback(Gk,1,-1)求取闭环传递函数w=0:0.01:100;设定分析的频率起止范围mag ,p hase=bode(Gb,w);求得频率响应M,n=max(mag);求取谐振峰值幅值Mr=20*log10(M)wr

9、=w(n)求取谐振频率bode(Gb,w);grid绘制闭环对数频率特性Mr =12.6957wr = 2.1500、填空题1.常用的串联校正有三种方式O (相位超前校正、相位滞后校正、相位滞后-超前校正)0 4s +12.某串联校正装置的传递函数为：Gc(s)=帛，该串联校正装置提供的是.校正。它能改善系统的O 超前(微分)、稳定性和快速性3. MATLAB提供的单输入/单输出系统可视化设计工具SISO Design Tool 是一种基于的系统设计工具。(根轨迹分析、频域性能指标)4. 欲使用SISO Design Tool对系统进行校正设计，需在MA TLAB命令窗口输入二、简答题命令启动

10、 SISO Design Tool 。 rltool(sys)在使用SISO Design Tool对系统进行校正设计时，若要调整校正装置提供的零、极点位置,系统提供了哪两种方法？请简述各自的操作步骤。(在使用SISO Design Tool对系统进行校正设计时，若要调整校正装置提供的零、极点位置，除了可以在显示系统根轨迹的“ SISO Design for SISO Design Task窗口用鼠标拖移零极点外，还可以再在“Control andEstimation Tools Manager窗口中选择Compensator Editor选项卡，在Location文本框中能准确设定所添加的开

11、环零点和极点的位置。)三、编程题1.已知某单位反馈系统开环传递函数为:10Gk(s)=聞，校正装置的传递函数为:C , 、(10s+1)(2s+1)Gc(s)=(205+1)(0 2s+1)，试绘制出系统校正前和校正后的对数幅频特性曲线以及校正装置的对数幅频特性曲线，并求出校正后的相位裕量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(numc,denc)Gk=G0*GcGm, Pm,Wg,Wc=margin(Gk);bode(G0,

12、b,Gk,r,Gc,g);gridPm = 73.6339 (度)2.单位负反馈系统的开环传递函数为:1000Gk(s)=，若要求其相位裕量s(0.001s+1)(0.1s+1)丫45o，试确定串联校正装置的传递函数（提示：串联超前校正）。clc;clf;clearnum0=1000;den0=conv(1 0,conv(0.001 1,0.1 1);G0=tf(num0,den0)Gm0, Pm 0,Wg0,Wc0=margin(G0);求得校正前系统的稳定裕量以及对应的频率r=45;e p=10;r0=Pm0;p=(r-r0+e p)* pi/180;计算需附加的相角a=(1+sin( p)”(1-sin( p);计算amag0, phase0,w=bode(G0);求得校正前系统的幅频mag20=20*log10(mag0);Mcm=-10*log10(a);Wm=s pline(mag20,w,Mcm);求得3mT=1/(Wm*sqrt(a);numc=a*T 1;denc