自动控制原理实验指导书

1、.自动控制原理实验指导书王娜 编写电气工程与自动化学院自动化系2017年11月实验一 控制系统的时域分析实验目的 1、熟悉并掌握Matlab操作环境和基本方法，如数据表示、绘图等命令； 2、掌握控制信号的拉氏变换与反变换laplace和ilaplace，控制系统生成模型的常用函数命令sys=tf(num,den)，会绘制单位阶跃、脉冲响应曲线；3、会构造控制系统的传递函数、会利用matlab函数求取系统闭环特征根； 4、会分析控制系统中 对系统阶跃、脉冲响应的影响。 实验内容及步骤 1、矩阵运算a) 构建矩阵：A=1 2;3 4; B=5 5;7 8; 解： A=1 2;3 4A = 1 2

2、3 4 B=5 5;7 8B = 5 5 7 8b) 已知A=1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4 ，求矩阵A的特征值、特征多项式和特征向量.解： A=1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4; V,D=eig(A)V = 0.4181 -0.4579 - 0.3096i -0.4579 + 0.3096i -0.6044 0.6211 -0.1757 + 0.2740i -0.1757 - 0.2740i 0.0504 0.5524 0.7474 0.7474 -0.2826 0.3665 -0.1592 - 0.067

3、5i -0.1592 + 0.0675i 0.7432 D = 13.0527 0 0 0 0 -4.1671 + 1.9663i 0 0 0 0 -4.1671 - 1.9663i 0 0 0 0 2.1815 p=poly(A)p =-6.9000 -77.2600 -86.1300 604.55002. 基本绘图命令 a) 绘制余弦曲线y=cos(x)，x0，2 解： x=linspace(0,2*pi); y=cos(x); plot(x,y)b) 线形为点划线、颜色为红色、数据点标记为加号； 解：hold on;plot(x,y,r-.+)c) 加网格线 解：grid ond) 标注

4、控制：x 、y坐标轴名称和标题“y=cos(t)”； 解：xlabel(x); ylabel(y); title(y=cos(x)3. 常用拉氏变换和反变换的命令F=laplace(f)：f(t)的拉氏变换，结果为F(s)，默认变量为s ；f=ilaplace(F)：F(s)的拉氏反变换，结果为f(t)，变量为t；例1-1试求函数的拉氏变换式，并用拉氏反变换观察变换结果。解：MATLAB程序如下： clear; %清除所有变量syms t A w b s %定义符号变量t, A, w, b, sft= A*sin(w*t+b); %定义f(t)的符号函数ft的表达式Fs =laplace(ft

5、) %求ft的拉氏变换式Fs, 即F(s)运行结果：Fs =A*(cos(b)*w/(s2+w2)+sin(b)*s/(s2+w2)可利用拉氏反变换对上述结果进行检验：ft=ilaplace(Fs) %求Fs的拉氏反变换式ft运行结果：ft =sin(t)即 f(t)=L-1F(s)=L-11/(s2+1)=sin(t)4.求系统的单位阶跃响应说明：step(num,den),其中num：传递函数分子表达式，den：传递函数分母表达式，幂次由高到低排列。例1-1：若已知单位负反馈前向通道的传递函数为，试作出其单位阶跃响应曲线，准确读出其动态性能指标，并记录数据。解：1）作单位阶跃响应曲线mat

6、lab参考程序graph.m如下：sys=tf(100,1 5 0);sysc=feedback(sys,1);step(sysc);grid on;2)运行程序得到系统的单位阶跃曲线如下：3)在曲线图中空白区域，单击鼠标右键，在快捷菜单中选择“Characteristics”命令，可以显示动态性能指标“Peak Response”(峰值Cp)，“Setting Time”（调节时间ts）、“Rise Time”（上升时间tr）和稳态值“Steady State”，如图：4）单击鼠标右键，在出现的快捷菜单中选择“Properties”命令，显示属性编辑对话框，如图：5）在“Option”选项卡

7、的“Show setting time within”文本框中，设置时间误差带2%或5%。6）读图中数据可得到系统稳态值为1，动态性能指标为：上升时间tr=0.127s,超调量Mp=44%，峰值时间tp=0.321s，调节时间ts=1.41s。7）已知二价震荡环节的传递函数G(s)= ，其中，从0变化到1.25，求此系统的单位阶跃响应曲线、脉冲响应曲线和斜坡响应曲线。解：参考函数如下：(1)系统单位阶跃响应曲线的程序代码:syms sfor zeta=0,0.25,0.5,0.75,1.0,1.25; wn=0.4; wn=sym(num2str(wn); zet=sym(num2str(ze

8、ta); if zeta=0 figure(1) ezplot(ilaplace(wn2/s/(s2+wn2),0 80 ); grid on title(xi=0) hold on elseif zeta=1 ezplot(ilaplace(wn2/s/(s+wn)2),0 80 ); hold on; else ezplot(ilaplace(wn2/s/(s2+2*zet*wn*s+wn2),0 80 ); hold on; end end grid on title(xi:0,0.4,0.7,0.9,1.0,1.5, ) axis(0 80 0 1.8) gtext(wn=0.4)绘图

9、结果显示：系统脉冲响应曲线的程序代码:syms sfor zeta=0,0.25,0.5,0.75,1.0,1.25; wn=0.4; wn=sym(num2str(wn); zet=sym(num2str(zeta); if zeta=0 figure(1) ezplot(ilaplace(wn2/(s2+wn2),0 80); grid on title(xi=0)hold on; elseif zeta=1 ezplot(ilaplace(wn2/(s+wn)2),0 80); hold on; else ezplot(ilaplace(wn2/(s2+2*zet*wn*s+wn2),0

10、 80); hold on; end end grid on title(xi:0,0.25,0.5,0.75,1.0,1.25, ) axis(0 80 -0.4 0.4) gtext(0) gtext(0.25) gtext(0.5) gtext(0.75) gtext(1.0) gtext(1.25)绘图结果显示：系统斜坡响应曲线的程序代码:syms sfor zeta=0,0.25,0.5,0.75,1.0,1.25; wn=0.4; wn=sym(num2str(wn); zet=sym(num2str(zeta); if zeta=0 figure(1) ezplot(ilapla

11、ce(wn2/s2/(s2+wn2),0 80); grid on title(xi=0) hold on elseif zeta=1 ezplot(ilaplace(wn2/s2/(s2+2*wn*s+wn2),0 80); hold on; else ezplot(ilaplace(wn2/s2/(s2+2*zeta*wn*s+wn2),0 80); hold on; end end grid on title(xi:0, 0.25,0.5,0.75,1.0, 1.25) axis(0 80 0 80) gtext(0) gtext(0.25) gtext(0.5) gtext(0.75)

12、 gtext(1.0) gtext(1.25)绘图结果显示：由上至下为zeta=0，0.25，0.5，0.75，1.0，1.25【分析】：可见，当wn一定时，系统随着阻尼比zeta的增大，闭环极点的实部在s左半平面的位置更加远离原点，虚部减小到0，超调量减小，调节时间缩短，稳定性更好。实验要求 1. 利用所学知识，编写实验内容中1.到4.的相应程序，并将所用到的命令行、中间变量和最终结果及产生的图形写在实验报告上。2. 在例1-1的7）中，令zeta=0.25保持不变，编写程序并绘制出wn=10,30和50时，对应系统的单位阶跃、脉冲和斜坡响应曲线。实验二 控制系统的频域分析一、实验目的1.以

13、二阶系统为例，掌握控制系统频率特性的基本原理；2.掌握控制系统Nyquist图的绘制方法，并加深理解控制系统奈奎斯特稳定性判据的实际应用；3.掌握利用matlab函数求取控制系统频域指标相角裕度、幅值裕度的方法；4.掌握控制系统伯德图的绘制方法，并会利用伯德图分析控制系统稳定性。二、实验原理1.对数频率特性曲线又称频率特性的对数坐标图或伯德图，又两张图组成，一张为对数幅频特性，其纵坐标为20log|G(jw)|,单位为分贝，用符号dB表示。20log|G(jw)|常用L（w）表示。另一张是相频特性图，其纵坐标为（），两张图的纵坐标均按线性分度，横坐标是角频率w，采用lg（w）分度（为了能在一张

14、图上同时展示出频率特性的低频和高频部分）。故坐标点w不得为零。1到10的距离等于10到100的距离。这个距离表示十倍频程，用dec表示。2.对数稳定判据对数频率特性曲线是奈氏判据移植于对数频率坐标的结果。若G（jw）H（jw）包围（-1，j0）点，即G（jw）H（jw）在点（-1，j0）左边有交点，在Bode图中表现为L（w）0分贝所在的频段范围内， 与-180线有交点。对数频率稳定性判据的内容为：闭环系统稳定的充分必要条件为：当w从零变化为+时，在开环系统对数幅频特性曲线L（w）0分贝的频段内，相频特性穿越（2k+1）（k=0，1、2，.）的次数N为P/2,其中N=N+-N-，N+为正穿越次

15、数，N-为负穿越次数。P为开环传函的正实部极点数。3. 稳定裕度1） 相角裕度 。当开环相频特性曲线（奈氏曲线）的幅值为1时，其相位角与-180（即负实轴）的相角差，称为相角裕度，即： ，式中，为奈氏曲线与单位圆相交处的频率，称为幅值穿越频率或剪切频率。当w=wc时，有| G（jw）H（jw）|=1. 相角裕度的含义是，对于闭环稳定系统，如果开环相频特性再滞后度，则系统将变为临界稳定。当 0时，相位裕度为正，闭环系统稳定。当=0时，表示奈氏曲线恰好通过（-1，j0）点，系统处于临界稳定状态。当1,即20lgKg（dB）0时，闭环系统稳定。当Kg=1时，系统处于临界稳定状态。当Kg1,即20lg

16、Kg（dB）0时，闭环系统不稳定。三、实验内容及步骤（1）产生频率向量 频率向量可由logspace( )函数来构成。此函数的调用格式为=logspace(m,n,npts)，此命令可生成一个以10为底的指数向量(10m10n )，点数由npts任意选定。 （2）绘制系统奈氏曲线 奈氏图由nyquist函数绘制，调用格式如下：格式1：nyquist(num,den):作nyquist图，角频率向量的范围自动设定，默认w（-，+）。 注意：在自动控制理论中，幅频特性L(w)为w的偶函数，相频特性 为w的奇函数，则w从零变化到+与从零变化到-的曲线是关于实轴对称的。格式2：nyquist（num,

17、den,w）:作开环系统的奈氏曲线，角频率向量w的范围可以人工给定，生成格式为w=logspace(d1,d2,n),表示将变量w作对数等分，d1和d2表示变量范围10d110d2，n为等分点数。格式3：re,im,w=nyquist(num,den):返回变量不做曲线，其中re为频率响应的实部，im为频率响应的虚部，w是频率点。例2-1 已知 ，绘制nyquist图，判定系统的稳定性。解：num=0.5;den=1 2 1 0.5;figure(1);nyquist(num,den)由于横坐标角频率的范围不够，从图中很难看出w从-变化到+时的相角，通过重新设置坐标范围显示全部范围的曲线，在当

18、前图形Figure1窗口中选择“Eidt” 菜单选项下的命令“Axes Properties”选项，在图形下方会增加一个坐标设置对话框，如图所示：根据实际需要更改该对话框的参数，使图形完全显示w从-变化到+时，系统奈氏曲线的形状。为应用奈氏曲线稳定性判据对闭环系统判稳，必须知道G(s)H(s)不稳定根的个数p是否为0，可以通过求其特征方程的根函数roots求得:p=1 2 1 0.5;roots(p)结果显示，系统有三个特征根：-1.5652，-0.2174+j0.5217，-0.2174-j0.5217.而且特征根的实部全为负数，都在s平面的左半平面，是稳定根，故p=0.由此，系统奈氏曲线没

19、有包围且远离（-1，j0）点，且p=0，因此系统闭环稳定。（3）绘制连续系统伯德图格式1：bode（num,den,w）:使用给定角频率w绘制系统的bode图。其中w为对数等分，生成格式为w=logspace(d1,d2,n),表示将变量w作对数等分，d1和d2表示变量范围10d110d2，n为等分点数。格式2：bode(num,den):在当前图形窗口中直接绘制系统的bode图，角频率w的范围自动设定。格式3：mag,phase,w=bode(num,den):返回变量格式，不做图，计算系统bode图的输出数据，其中输出变量mag是系统bode图的幅值向量mag=|G(jw)|,注意此幅值不

20、是分贝值，须用mag（dB）=20*log(mag)转换；phase为bode图的幅角向量，phase=G(jw),单位为（）；w是系统bode图的频率向量，单位是rad/s。例2-2 已知控制系统开环传递函数 ，绘制其bode图。解：参考程序如下：num=10;den=1 2 10;bode(num,den) %显示系统的伯德图运行结果如下：num=10;den=1 2 10;bode(num,den)例2-3 在上述系统bode图中，确定谐振峰值的大小Mr和谐振频率wr。解：m,p,w=bode(num,den);%返回变量格式，得到(m,p,w)向量mr=max(m) %由最大值函数得到

21、m的最大值wr=spline(m,w,mr)%由插值函数spline求得谐振频率。运行结果为：谐振峰值Mr=1.6667，谐振频率wr=2.8284rad/s其中spline（）为3次样条函数插值，函数格式为：spline(x,y,xi),式中，y=f(x),xi为等分值。（4）计算系统的稳定裕度，包括增益裕度Gm和相位裕度Pm。格式1：margin(num,den)给定开环系统的数学模型，作bode图，并在图上标注增益裕度Gm和对应频率wg，相位裕度Pm和对应频率wc。格式2：Gm,Pm,wg,wc=margin(num,den)返回变量格式，不做图。格式3：Gm,Pm,wg,wc=marg

22、in(m,p,w)给定频率特性的参数向量：幅值m、相位p和频率w，由插值法计算Gm及wg、Pm及wc。例2-4 已知单位负反馈系统的开环传函，求系统的稳定裕度，并分别用格式2和格式3计算，比较误差。解：k=2;z=;p=0 -1 -2;num,den=zp2tf(z,p,k);Margin(num,den);Gm1,Pm1,wg1,wc1=margin(num,den) %格式2求出系统稳定裕度m,p,w=bode(num,den);Gm2,Pm2,wg2,wc2=margin(m,p,w) %格式3求出系统稳定裕度程序运行后显示系统的bode图如图所示，并在图的上方标出了稳定裕度，同时格式2

23、、3求出稳定裕度显示在命令窗口中。分析：比较以下两组数据：由格式2计算出的数据：Gm1=3；wg1=1.4142rad/s，Pm1=32.6133，wc1=0.7494rad/s由格式3计算出的数据：Gm2=3；wg2=1.4141rad/s，PM=32.6138，wc2=0.7492rad/s由于格式2计算时的w的点数是自动给定的，因此直接求取的稳定裕度与采用格式3经插值后得到的稳定裕度比较，后者的精确度高一些，这两种格式计算出的Gm不是分贝值，可以转化成dB，20lgGm=20lg3dB=9.54dB.（5）绘制尼柯尔斯图(Nichols图) 尼柯尔斯图是将线性非时变系统在不同频率下的增益

24、分贝值及相位绘在一直角坐标系的图上，尼柯尔斯图将二种波德图（波德增益图及波德相位图）结合成一张图，而频率只是曲线中的参数，不直接在图中显示。尼柯尔斯图，其中可看到增益裕度及相位裕度。尼柯尔斯图可以用来分析系统的稳定性，以及增益裕度、相位裕度等有关系统相对稳定性。在尼柯尔斯图上可以看到相位-180度，增益0dB的点。找出尼柯尔斯图对应相位-180度的点1：若此点在增益0dB的点上方，表示其增益大于0dB，对应的单位回授系统不稳定。若此点在增益0dB的点下方，表示其增益小于0dB，对应的单位回授系统稳定，而两者的距离即为增益裕度。而根据尼柯尔斯图对应增益0dB度的点也可以判断是否稳定，及相位裕度：若此点在相位-180度点左方，表示其相位小于-180度，对应的单位回授系统不稳定。若此点在相位-180度点右方，表示其相位大于-180度，对应的单位回授系统稳定，而两者的距离即为相位裕度。函数调用格