控制系统CAD课程形考作业三参考答案

《控制系统CAD》形成性考核册-1-控制系统CAD作业3第5章一、填空题1.在控制系统时域分析时，常用的时域动态性能指标有：上升时间tr、、和。（峰值时间tm、调节时间ts、超调量σ%）2.MATLAB中step()函数是函数。它可以用来或求取。（单位阶跃响应、绘制单位阶跃响应波形图，输出量的数值向量）3.impulse函数不带输出变量时，可在当前窗口中直接绘制出系统的。（单位冲激响应曲线）4.在MATLAB中，欲求取系统的时域性能，可以利用鼠标和菜单的操作，在单位阶跃响应曲线上求取，也可以求取。（编制简单程序）5.要判别系统的稳定性，可使用的MATLAB函数有、、和四个函数。若要直接求出系统闭环特征方程的根，应使用函数。（roots()、zpkdata()、tf2zp()、pzmap()、roots()）6.若在Simulink仿真集成环境中分析系统在斜坡信号和阶跃信号作用下系统的稳态误差，在Simulink建模时，可分别选用和模块作为系统的输入信号模块。（Step、Ramp）二、写出下列程序（命令）的功能或执行结果姓名：学号：得分：教师签名：《控制系统CAD》形成性考核册-2-1.G=zpk([],[-1+3*I,-1-3*i],3)step(G)绘制出3()(13)(13)Gssjsj的单位阶跃响应曲线2.num=7*[1,5];den=conv([1,0,0],conv([1,10],[1,1]));[gm,pm,wg,wc]=margin(num,den)求系统)1)(10()5(7)(2sssssG的幅值裕度和相角裕度3.num=5*[1,5,6];den=[1,6,10,8];impulse(num,den);gridon;求三阶系统的8106)65(5)(232ssssssG单位阶跃响应和单位冲激响应。三、编程题1.典型二阶系统传递函数为：222()2nnnGsss，阻尼比ζ的值为0.4，试编制程序绘制出自然振荡角频率ωn从1变化到10（每次递增2）的单位阶跃响应曲线，并对系统的动态性能的变化作定性分析。clc;clf;clearwn=[1,3,5,7,9];%自然振荡角频率dr=0.4;%阻尼比t=linspace(0,20,200);%将时间t在0-20间200等分holdon;%打开图形叠加功能（同一个坐标上画多条曲线)forj=1:5num=wn(j)^2;den=[1,2*dr*wn(j),0];%开环传递函数分母多项式系数Gk=tf(num,den);%求得开环传递函数Gb=feedback(Gk,1);%求得闭环传递函数《控制系统CAD》形成性考核册-3-step(Gb);grid;%画图gtext(strcat('wn=',num2str(j)))%添加鼠标文字注释，strcat()为字符串连接函数end随着自然振荡角频率ωn的增大，系统的响应速度加快。2.设单位反馈系统的开环传递函数：249()10.5kGsss，试求系统的单位阶跃响应的性能指标σ％、tS（2％误差带）。clc;clearnum=49;den=[110.50];Gk=tf(num,den)%求得开环传递函数Gb=feedback(Gk,1,-1)%求得单位反馈的闭环传递函数[y,t]=step(Gb);%求取单位阶跃响应%计算调节时间Cend=dcgain(Gb);%求取系统的终值n=length(t);%求得仿真时间t序列的长度while(y(n)0.98*Cend)&(y(n)1.02*Cend)%通过循环求取单位阶跃响应首次进入±2%h(∞)内的所需时间序列号n=n-1;endTs=t(n)%求得调节时间%计算超调量[Y,i]=max(y);%求得系统响应的峰值Os=100*(Y-Cend)/Cend%求得超调量Ts=0.8074（s）Os=2.8335（%）3.单位负反馈系统的开环传递函数为：27(1)()(3)(45)ksGsssss，试求在单位斜坡输入信号作用下的稳态误差，试写出MATLAB命令序列（需首先判断系统的稳定性）。num=[77];den=[conv(conv([10],[13]),[145])];Gk=tf(num,den);Gb=feedback(Gk,1,-1)Transferfunction:7s+7《控制系统CAD》形成性考核册-4--------------------------------s^4+7s^3+17s^2+22s+7r=roots([1717227])r=-4.0287-1.2632+1.5198i-1.2632-1.5198i-0.4449闭环特征根的实部均为负，系统稳定。rss=tf([10],[01]);Kv=dcgain(rss*Gk)essv=1/Kvessv=2.1429第6章一、填空题1.所谓根轨迹，是指控制系统开环传递函数的某一参数（如开环增益）从变化到时，在s平面上移动的轨迹。（零、无穷大、闭环特征根）2.根轨迹法的实质就是依据反馈系统中开、闭环间的关系，由的零、极点直接寻求的总体规律。（传递函数、开环传递函数、闭环根轨迹）3.根据给定的根轨迹增益向量k，绘制系统sys的部分根轨迹的MATLAB函数是。（rlocus(sys,k)）4.计算与根轨迹上任一闭环极点相对应的根轨迹增益的MATLAB函数是。（rlocfind()）二、编程分析题1.已知单位反馈系统的开环传递函数，试用MALAB编程绘出系统根轨迹。⑴()(0.21)(0.51)KGssss(2)(1)()(21)KsGsss《控制系统CAD》形成性考核册-5-(1)num=1;den=[conv(conv([10],[0.21]),[0.51])]sys=tf(num,den);rlocus(sys);-20-15-10-505-15-10-5051015RootLocusRealAxisImaginaryAxis(2)num=[12];den=[conv([10],[21])]sys=tf(num,den);rlocus(sys);-7-6-5-4-3-2-101-2-1.5-1-0.500.511.52RootLocusRealAxisImaginaryAxis2.已知控制系统的开环传递函数为32324)37kssGssss（，试用MALAB分析该系统的根轨迹的起点和终点。num=[1104];den=[1370];G=tf(num,den);rlocus(G)；p=roots(den)《控制系统CAD》形成性考核册-6-z=roots(num)p=0-1.5000+2.1794i-1.5000-2.1794iz=-2.00000.5000+1.3229i0.5000-1.3229i-2.5-2-1.5-1-0.500.51-2.5-2-1.5-1-0.500.511.522.5RootLocusRealAxisImaginaryAxis图中显示了该系统的根轨迹。可以看到，该系统有3个开环极点和3个开环零点，因此根轨迹有3个分支，它们的起点是开环极点0,-1.5+2.18j和-1.5-2.18j，终点是开环零点-2,0.5+1.32和0.5-1.32j，根轨迹的一个分支从极点0开始，终于零点-2；另两条分支分别从极点-1.5+2.18j和-1.5-2.18j开始，从圆弧变化，最后分别终于零点0.5+1.32和0.5-1.32j。3.已知单位负反馈控制系统的开环函数为()(6)(3)kKGssss，请利用MATLAB分析该系统临界稳定时的K值与系统的闭环极点，并求根轨迹与虚轴交点处的根轨迹增益。num=1den=[conv(conv([10],[16]),[13])]G=tf(num,den);rlocus(G);k=0:0.1:100;[kp]=rlocfind(G)Selectapointinthegraphicswindow《控制系统CAD》形成性考核册-7-selected_point=0.0533+4.2391ik=163.6757p=-9.01690.0084+4.2605i0.0084-4.2605i-20-15-10-505-15-10-5051015RootLocusRealAxisImaginaryAxis