自动控制原理MATLAB仿真实验报告材料
MATLAB与控制系统仿真实验报告
MATLAB与控制系统仿真实验报告第一篇：MATLAB与控制系统仿真实验报告《MATLAB与控制系统仿真》实验报告2013-2014学年第 1 学期专业：班级：学号：姓名：实验三 MATLAB图形系统一、实验目的：1.掌握绘制二维图形的常用函数。
2.掌握绘制三维图形的常用函数。
3.熟悉利用图形对象进行绘图操作的方法。
4.掌握绘制图形的辅助操作。
二、实验原理：1，二维数据曲线图（1）绘制单根二维曲线plot(x,y);（2）绘制多根二维曲线plot(x,y)当x是向量，y是有一维与x同维的矩阵时，则绘制多根不同颜色的曲线。
当x，y是同维矩阵时，则以x，y对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。
（3）含有多个输入参数的plot函数plot(x1,y1,x2,y2,…,xn,yn)（4）具有两个纵坐标标度的图形plotyy(x1,y1,x2,y2)2,图形标注与坐标控制1）title(图形名称)；2）xlabel（x轴说明）3）ylabel（y轴说明）4）text（x，y图形说明）5）legend（图例1，图例2，…）6）axis（[xmin xmax ymin ymax zmin zmax]）3, 图形窗口的分割 subplot（m,n,p）4,三维曲线plot3（x1,y1,z1,选项1，x2,y2，选项2，…,xn,yn,zn,选项n）5，三维曲面mesh(x,y,z,c)与surf(x,y,z,c)。
一般情况下，x，y，z是维数相同的矩阵。
X，y是网格坐标矩阵，z是网格点上的高度矩阵，c用于指定在不同高度下的颜色范围。
6，图像处理1）imread和imwrite函数这两个函数分别用于将图象文件读入matlab工作空间，以及将图象数据和色图数据一起写入一定格式的图象文件。
2）image和imagesc函数这两个函数用于图象显示。
为了保证图象的显示效果，一般还应使用colormap函数设置图象色图。
(最新版)自动控制原理MATLAB仿真实验报告
实验一 MATLAB及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性；二、预习要点1、系统的典型响应有哪些？2、如何判断系统稳定性？3、系统的动态性能指标有哪些？三、实验方法（一）四种典型响应1、阶跃响应：阶跃响应常用格式：1、；其中可以为连续系统，也可为离散系统。
2、；表示时间范围0---Tn。
3、；表示时间范围向量T指定。
4、；可详细了解某段时间的输入、输出情况。
2、脉冲响应：脉冲函数在数学上的精确定义：其拉氏变换为：所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式：①；②③（二）分析系统稳定性有以下三种方法：1、利用pzmap绘制连续系统的零极点图；2、 利用tf2zp 求出系统零极点；3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.四、实验内容 (一) 稳定性1． 系统传函为()27243645232345234+++++++++=s s s s s s s s s s G ，试判断其稳定性2． 用Matlab 求出的极点。
%Matlab 计算程序num=[3 2 5 4 6];den=[1 3 4 2 7 2];G=tf(num,den);pzmap(G);p=roots(den) 运行结果： p =-1.7680 + 1.2673i -1.7680 - 1.2673i 0.4176 + 1.1130i 0.4176 - 1.1130i -0.2991P ole-Zero MapReal AxisI m a g i n a r y A x i s-2-1.5-1-0.500.5-1.5-1-0.50.511.5图1-1 零极点分布图由计算结果可知，该系统的2个极点具有正实部，故系统不稳定。
自动控制原理 matlab实验报告
自动控制原理实验（二）一、实验名称：基于MATLAB的控制系统频域及根轨迹分析二、实验目的：（1）、了解频率特性的测试原理及方法；（2）、理解如何用MATLAB对根轨迹和频率特性进行仿真和分析；（3）、掌握控制系统的根轨迹和频率特性两大分析和设计方法。
三、实验要求：（1）、观察给定传递函数的根轨迹图和频率特性曲线；（2）、分析同一传递函数形式，当K值不同时，系统闭环极点和单位阶跃响应的变化情况；（3）、K值的大小对系统的稳定性和稳态误差的影响；（4）、分析增加系统开环零点或极点对系统的根轨迹和性能的影响。
四、实验内容及步骤（1）、实验指导书：实验四（1）、“rlocus”命令来计算及绘制根轨迹。
会出根轨迹后，可以交互地使用“rlocfind”命令来确定点击鼠标所选择的根轨迹上任意点所对应的K值，K值所对应的所有闭环极点值也可以使用形如“[K, PCL] = rlocfind(G1)”命令来显示。
（2）、波特图：bode(G1, omga)另外，bode图还可以通过下列指令得出相位和裕角：[mag,phase,w] = bode(sys)（3）、奈奎斯特图：nuquist(G, omega)（2）课本：例4-1、4-2、4-7五实验报告要求（1）、实验指导书：实验四思考题请绘制下述传递函数的bode图和nyquist图。
1. 根据实验所测数据分别作出相应的幅频和相频特性曲线；2. 将思考题的解题过程(含源程序)写在实验报告中。
幅频特性曲线相频特性曲线Gs = zpk([10], [-5; -16; 9], 200)subplot(1, 2, 1)bode(Gs)gridsubplot(1, 2, 2)nyquist(Gs)grid（2）课本：例4-1、4-2、4-7图像结果：程序：Gs = zpk([-1], [0; -2; -3],1) rlocus(Gs)图像结果：程序：Gs = zpk([-2], [-1-j; -1+j],1) rlocus(Gs)程序：K=[0.5 1 2]for i=1:1:3num=[1,1,0,0]; den=[1,1,K(i)]; sys=tf(num,den); rlocus(sys); hold ongrid onend图像结果：目标：改变增益K和转折频率依次调节源程序：k1=[4.44,10,20];num=[1,2];den=conv([1,1],[1,2,4]);%一阶转折频率 1/T（wn1=2,wn2=1）二阶转折频率 wn3=wn'=2，伊布西塔=1/2 num1=[1,1];den1=conv([1,2],[1,2,4]);%一阶转折频率 1/T（wn1=1,wn2=2）二阶转折频率 wn3=wn'=2，伊布西塔=1/2 t=[0:0.1:7]; %for i=1:3g0=tf(k1(i)*num,den);g=feedback(g0,1);[y,x]=step(g,t);c(:,i)=y;g1=tf(k1(i)*num1,den1);g(1)=feedback(g1,1);[y1,x]=step(g(1),t);c1(:,i)=y1;endplot(t,c(:,1),'-',t,c(:,2),'-',t,c(:,3),'-',t,c1(:,1),'-',t,c1(:,2), '-',t,c1(:,3),'-');gridxlabel('Time/sec'),ylabel('out')结果分析：在本题中（1）改变k值：k值越大，超调量越大，调节时间越长，峰值时间越短，稳态误差越小（2）改变转折频率：超调量，调节时间，峰值时间，稳态误差同样有相应的变化。
自动控制原理MATLAB仿真实验报告
实验一 MATLAB 及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点1、 系统的典型响应有哪些2、 如何判断系统稳定性3、 系统的动态性能指标有哪些 三、实验方法（一） 四种典型响应1、 阶跃响应：阶跃响应常用格式：1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。
2、),(Tn sys step ；表示时间范围0---Tn 。
3、),(T sys step ；表示时间范围向量T 指定。
4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。
2、 脉冲响应：脉冲函数在数学上的精确定义：0,0)(1)(0〉==⎰∞t x f dx x f其拉氏变换为：)()()()(1)(s G s f s G s Y s f ===所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式： ① )(sys impulse ； ②);,();,(T sys impulse Tn sys impulse③ ),(T sys impulse Y =（二） 分析系统稳定性 有以下三种方法：1、 利用pzmap 绘制连续系统的零极点图；2、 利用tf2zp 求出系统零极点；3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.四、实验内容 (一) 稳定性1． 系统传函为()27243645232345234+++++++++=s s s s s s s s s s G ，试判断其稳定性2． 用Matlab 求出253722)(2342++++++=s s s s s s s G 的极点。
%Matlab 计算程序num=[3 2 5 4 6];den=[1 3 4 2 7 2];G=tf(num,den);pzmap(G);p=roots(den)运行结果： p =+ - + -P ole-Zero MapReal AxisI m a g i n a r y A x i s-2-1.5-1-0.500.5-1.5-1-0.50.511.5图1-1 零极点分布图由计算结果可知，该系统的2个极点具有正实部，故系统不稳定。
自动控制原理MATLAB仿真实验(于海春)
自动控制原理MATLAB仿真实验（于海春）实验一典型环节的MATLAB仿真一、实验目的1．熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法。
2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。
3．定性了解各参数变化对典型环节动态特性的影响。
二、SIMULINK 的使用MATLAB中SIMULINK是一个用来对动态系统进行建模、仿真和分析的软件包。
利用SIMULINK功能模块可以快速的建立控制系统的模型，进行仿真和调试。
1．运行MATLAB软件，在命令窗口栏“>>”提示符下键入imulink命令，按Enter键或在工具栏单击按钮，即可进入如图1-1所示的SIMULINK仿真环境下。
2．选择File菜单下New下的Model命令，新建一个imulink仿真环境常规模板。
图1-1SIMULINK仿真界面图1-2系统方框图3．在imulink仿真环境下，创建所需要的系统。
以图1-2所示的系统为例，说明基本设计步骤如下：1）进入线性系统模块库，构建传递函数。
点击imulink下的“Continuou”，再将右边窗口中“TranferFen”的图标用左键拖至新建的“untitled”窗口。
2）改变模块参数。
在imulink仿真环境“untitled”窗口中双击该图标，即可改变传递函数。
其中方括号内的数字分别为传递函数的分子、分母各次幂由高到低的系数，数字之间用空格隔开；设置完成后，选择OK，即完成该模块的设置。
3）建立其它传递函数模块。
按照上述方法，在不同的imulink的模块库中，建立系统所需的传递函数模块。
例：比例环节用“Math”右边窗口“Gain”的图标。
4）选取阶跃信号输入函数。
用鼠标点击imulink下的“Source”，将右边窗口中“Step”图标用左键拖至新建的“untitled”窗口，形成一个阶跃函数输入模块。
5）选择输出方式。
控制工程基础实验——Matlab仿真实验报告
实验一：Mat lab 仿真实验1.1直流电机的阶跃响应给直流电机一个阶跃，直流电机的传递函数如下:画出阶跃响应如下:Step Resp onse零极点分布:POle-ZeroMap0.8 0.60.4-0.4 -0.6 -0.8g m-0.2 -1-10000-9000 -8000 -7000 -6000 ReaWi@ -4000 -3000 -2000 -1000s A 0.2G(s)=50(0.1s 1)(1 10*s 1)分析：直流电机的传递函数方框图如下:所以传递函数可以写成:n (s) 1/C EU a (S )FaS 2 T m S T式中，T mJ^,T a =L分别为电动机的机电时间常数与电磁时间常数。
一般 C M C ER相差不大。
而试验中的传递函数中，二者相差太大，以至于低频时:（低频时）0.1s 1所以对阶跃的响应近似为:x °(t) = 50(1 - e 处)G(s)二 _______ 50(0.1s 1)(1 10，s 1)直流电机传递函数方块图1.2直流电机的速度闭环控制如图1-2,用测速发电机检测直流电机转速，用控制器Gc(s)控制加到电机电枢上 的电压。
图1-2氏流电机速度闭环控1.2.1假设Gc(s)=100,用matlab 画出控制系统开环Bode 图，计算增益剪切频率、相位 裕量、相位剪切频率、增益裕量Bode Diagram5d^MTno0050O■ ■90sa&—80 ^17021、10幅值裕量Gm1 2 310 10 10Freque ncy (rad/sec)=11.1214410510610相位裕量Pm = 48.1370幅值裕量对应的频率值(相位剪切) wcg =3.1797e+003相位裕量对应的频率值(幅值剪切)wcp =784.3434从理论上，若G c(s) =100，那么开环传递函数为：G(s)=100-4(0.1s 1)(0.001s 1)(1 10 s 1)1001 (0.1 )2 J (0.001 )2 1 (1 10A )21 1 1.G(j H 一[tan (0.1 ) tan (0.001 ■) tan (0.0001 ■)]令G(j%)|=1，假设J1+(0.佃J托0.1灼，J+(1汉10鼻国)2屯1 得：c=786.15 继而，.G(j c)二-[tan d(0.1 c) tan'(0.001 c) tan'(0.0001 c)] = 48.06 1.2.2 通过分析bode图，选择合适的K p作为G c(s)，使得闭环超调量最小。
自动控制原理MATLAB分析与设计仿真实验报告(最终版)
兰州理工大学《自动控制原理》MATLAB分析与设计仿真实验报告学院：电气工程与信息工程学院专业班级： 13级自动化3班姓名：学号：时间： 2015年12月Step ResponseTime (seconds)A m p l i t u d e1234567891000.511.5System: sys1Rise time (seconds): 1.17System: sys1P eak amplitude: 1.41Overshoot (%): 40.6At time (seconds): 2.86System: sys1Final value: 1第三章 线性系统的时域分析法一、教材第三章习题3.5设单位反馈系统的开环传递函数为G(s)=0.41(0.6)s s s ++(1)试求系统在单位阶跃输入下的动态性能。
(2)忽略闭环零点的系统在单位阶跃输入下的动态性能。
(3)对(1) 和(2)的动态性能进行比较并分析仿真结果。
(1)A :程序如下。
B :系统响应曲线如下图。
Step Response Time (seconds)A m p l i t u d e01234567891000.20.40.60.811.21.4System: sys1Final value: 1System: sys1Settling time (seconds): 8.08System: sys1P eak amplitude: 1.16Overshoot (%): 16.3At time (seconds): 3.63System: sys1Rise time (seconds): 1.64（2）A :程序如下。
B :系统响应曲线如下图。
(3) A :程序如下。
B 响应曲线如下图。
阶跃响应t (sec)c (t )0123456789100.20.40.60.811.21.4System: sysRise Time (sec): 1.46System: sys1Rise Time (sec): 1.64System: sys1P eak amplitude: 1.16Overshoot (%): 16.3At time (sec): 3.63System: sys P eak amplitude: 1.18Overshoot (%): 18At time (sec): 3.16System: sys1Final Value: 1System: sys1Settling Time (sec): 8.08System: sysSettling Time (sec): 7.74120,0.1ττ==120.1,0ττ==分析：忽略闭环零点时，系统的峰值时间，调节时间，上升时间均为增大的，而超调量减小。
自动控制原理MATLAB仿真实验报告
实验一 典型环节的MATLAB 仿真 一、实验目的1．熟悉MATLAB 桌面和命令窗口，初步了解SIMULINK 功能模块的使用方法。
2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。
3．定性了解各参数变化对典型环节动态特性的影响。
二、实验内容① 比例环节1)(1=s G 和2)(1=s G ；Simulink 图形实现：示波器显示结果：② 惯性环节11)(1+=s s G 和15.01)(2+=s s GSimulink 图形实现：示波器显示结果：③ 积分环节s s G 1)(1Simulink 图形实现：示波器显示结果：④ 微分环节s s G )(1Simulink 图形实现：波器显示结果：⑤ 比例+微分环节（PD ）2)(1+=s s G 和1)(2+=s s G1）、G1（s ）=s+2Simulink 图形实现：示波器显示结果：2）、G2（s）=s+1 Simulink图形实现：示波器显示结果：⑥ 比例+积分环节（PI ）s s G 11)(1+=和s s G 211)(2+=1）、G1（1）=1+1/sSimulink 图形实现：示波器显示结果：2）G2（s）=1+1/2s Simulink图形实现：示波器显示结果：三、心得体会通过这次实验我学到了很多，对课本内容加深了理解，熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法，加深对各典型环节响应曲线的理解，这为对课程的学习打下了一定基础。
实验二线性系统时域响应分析一、实验目的1．熟练掌握step( )函数和impulse( )函数的使用方法，研究线性系统在单位阶跃、单位脉冲及单位斜坡函数作用下的响应。
2．通过响应曲线观测特征参量ζ和nω对二阶系统性能的影响。
3．熟练掌握系统的稳定性的判断方法。
二、实验内容1．观察函数step( )的调用格式，假设系统的传递函数模型为243237()4641s s G s s s s s ++=++++绘制出系统的阶跃响应曲线？2．对典型二阶系统222()2n n n G s s s ωζωω=++1）分别绘出2(/)n rad s ω=，ζ分别取0,0.25,0.5,1.0和2.0时的单位阶跃响应曲线，分析参数ζ对系统的影响，并计算ζ=0.25时的时域性能指标,,,,p r p s ss t t t e σ。
自动控制原理实验报告MATHLAB建模.doc
自动控制原理实验报告学院：机电工程学院班级：姓名：学号：指导老师：实验一：在MATLAB中创建系统模型一、实验目的：1．熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法。
2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。
二、实验过程：1.传递函数模型的建立1）多项式形式的传递函数①课本例2.23上机操作过程如下：②课本P62，习题2-3上机操作过程如下：2）零、极点形式的传递函数课本例2.24上机操作过程如下：3）分子、分母为因式乘积形式的传递函数课本例2.25上机操作过程如下：2.Simulink 建模①课本例题上机操作如下：设单位反馈系统的开环传递函数为：)1(1)(+=s s s G将其转换成Simulink 框图，输入为阶跃信号，它的Simulink 框图如下所示：② 比例环节1)(1=s G 和2)(1=s G 的SIMULINK 图形建模操作如下；比例环节1)(1=s G 的SIMULINK 图形如下图所示：比例环节2)(1=s G 的SIMULINK 图形3.课后练习用matlab求下列函数的拉氏变换（习题2-1），上机操作过程如图所示：实验二：在MATLAB中算特征根及绘制根轨迹图一、实验目的：1．掌握MATLAB下的根轨迹绘制方法；2．学会利用根轨迹进行系统分析。
二、实验过程：1）例3-21 试利用MATLAB函数求例3.1中k=2.k=20时系统的特征根，并分别判定稳定性。
上机操作过程如下：>> num=[2];den=conv([1 0],conv([0.1 1],[0.25 1]));g=tf(num,den);sys=feedback(g,1);>> pzmap(sys)p=pole(sys)p =-11.0314-1.4843 + 2.2470i-1.4843 - 2.2470i2）例3-22 二阶系统如图3.13所示，设Wn=1,试研究系统的单位阶跃响应与参数ξ的关系。
自动控制原理matlab实验报告
自动控制原理matlab实验报告1.由题意得：C(s)=R(s)*(11s+K)/(s2+12s+K)-N(S)/(S2+12S+k)该系统显然是稳定的。
为了减少扰动的影响，希望增益K>0。
扰动引起的稳态误差e ssn=1/K,现使扰动引起的稳态误差小于0.05，最大超调量小于0.1，则K的取值范围是：20<k<100。
实验中，选取K=20,25,30,40,100进行五次实验，实验结果记录如下：由表中数据可得，使扰动引起的稳态误差较小，且使单位阶跃输入下超调量也相对小的情况下，本系统应选取K=25。
实验中K取不同值时的响应如下：K=20 K=25K=30 K=40K=1002.C(s)=R(s)*Ka/(s2+k1s+Ka)-N(S)/(S2+k1S+ka)（1）在阶跃指令r(t)作用下，系统输出的超调量小于或等于10%； 由解得：代入σ=0.1，求出 ζ=0.59，取ζ=0.6。
因而，在满足σ%≤10%指标要求下，应选（2）在斜坡输入作用下，稳态误差达到最小； 令斜坡输入为r(t)=Bt,可得斜坡输入作用下的稳态误差：结合要求（1）可得此式表明K a 应取尽可能大（3）减小单位阶跃扰动的影响。
阶跃扰动作用下的稳态误差22)(ln 11σπ+=ζaa 1K 2.1K 2K =ζ=a 1ssrK BK K B e ==assrK B 2.1e =)s (sC )s (sEe n 0s ns ssn l i m l i m →→-==aa 12s 00s K 1s1K s K s 1s)s (N )s (G 1)s (G s l i m l i m -=++-=+-=→→%e100%21 / ζ- πζ - =σ可见，增大K a可以同时减小e ssn及e ssr。
在实际系统中，K a的选取必须受到限制，以使系统工作在线性区。
实验中选取以下几组数据进行仿真。
KA=100，K1=12 KA=576，K1=30KA=625，K1=30 KA=900，K1=40KA=1000，K1=45由上表及仿真图分析可知应取K a =1000，K 1=45.3. 此系统的特征方程为：s 4+8s 3+17s 2+(10+K 1)s+aK 1=0 由题目要求可得： 斜坡输入下的稳态误差：K K )a 64116(12600aK 126K 21111>--+><令斜坡输入为r(t)=At令稳态误差等于输入指令幅度的24%。
自动控制原理Matlab实验1(初步认识MATLAB和系统仿真)
《自动控制原理》课程实验报告实验名称初步认识MATLAB和系统仿真专业班级 ********* 学************号姓名**指导教师李离学院名称电气信息学院2012 年 11 月 5 日Lab1_1_1.m程序：y0=0.15;wn=sqrt(2);zeta=1/(2*sqrt(2));t=[0:0.1:10];c=(y0/sqrt(1-zeta^2));y=c*exp(-zeta*wn*t).*sin(wn*sqrt(1-zeta^2)*t+acos(zeta)); bu=c*exp(-zeta*wn*t);bl=-bu;plot(t,y,t,bu,'k--',t,bl,'k--'),gridxlabel('Time (sec)'),ylabel('y(t) (meters)')legend(['\omega_n=',num2str(wn),' \zeta=',num2str(zeta)])仿真结果：(1)零输入响应曲线理论分析:0<ζ<1,对于二阶响应，其瞬态响应应该是一个按照指数衰减的振荡过程，ζ越小，衰减越慢，该系统是欠阻尼系统。
从图中也可以看出，系统是零输入响应，是个震荡衰减的过程，符合理论判断。
Lab1_1_2.m程序：y0=0.15;wn=sqrt(2);zeta=1;t=[0:0.1:10];y=y0*(exp(-wn*t)+wn*t.*exp(-wn*t));plot(t,y),gridxlabel('Time (sec)'),ylabel('y(t) (meters)')legend(['\omega_n=',num2str(wn),' \zeta=',num2str(zeta)])仿真结果：(2)零输入响应曲线理论分析:在图中可以看到，随着时间的增加，响应在逐渐减小。
自动控制MATLAB实验报告
实验一 Matlab 使用方法和程序设计 热动102班 2010031294 唐旭一、实验目的1、理解 Matlab 软件使用的基本方法；2、理解 Matlab 的数据表示、基本运算和程序控制语句3、熟悉 Matlab 绘图命令及基本绘图控制4、学习使用 SIMULINK 进行系统仿真的方法二、实验内容：1、基本绘图命令（1）绘制余弦曲线 y=cos(t)，t ∈[0，2π] >> x=0:pi/15:2*pi;y=[cos(x)];plot(y)5101520253035-1-0.8-0.6-0.4-0.200.20.40.60.81（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t ∈[0，2π]>> x=0:pi/15:2*pi;y=[sin(x-0.5);cos(x-0.25)];plot(x,y)1234567-1-0.8-0.6-0.4-0.200.20.40.60.812、基本绘图控制绘制[0，4π]区间上的 y1=10sint 曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：x 轴和y 轴的显示范围（X 轴0：4Pi ；Y 轴：-10：10）、单位刻度大小相等、有网格线（3）标注控制：坐标轴名称（X 轴t ；Y 轴：y1）、标题（Sin 函数）指定位置创建说明性文本；>> x=linspace(0,4*pi,30);y1=10*sin(x); axis([0,4*pi,-10,10]); plot(x,y1,'r-.+') axis equal gridxlabel('t') ; ylabel('y1') ; title('Sin 函数');text(3,4,'热动102:唐旭') ; axis([0,4*pi,-10,10]);ty 1Sin 函数3、系统仿真PID 控制系统的结构如图所示，系统输入为斜坡输入，分别采用P 、PD 、PI 、PID 控制进行仿真，并分析P 、I 、D 对系统性能的影响。
《自动控制原理》MATLAB的使用分析与设计研究 仿真实验报告
兰州理工大学《自动控制原理》MATLAB分析与设计仿真实验报告院系：电信学院班级：姓名：学号：时间：2010年11月30日电气工程与信息工程学院《自动控制原理》MATLAB分析与设计仿真实验任务书（2009）一．仿真实验内容及要求：1．MATLAB软件要求学生通过课余时间自学掌握MATLAB软件的基本数值运算、基本符号运算、基本程序设计方法及常用的图形命令操作；熟悉MATLAB仿真集成环境Simulink的使用。
2．各章节实验内容及要求1）第三章线性系统的时域分析法• 对教材P136.3-5系统进行动态性能仿真，并与忽略闭环零点的系统动态性能进行比较，分析仿真结果；• 对教材P136.3-9系统的动态性能及稳态性能通过的仿真进行分析，说明不同控制器的作用；• 在MATLAB 环境下完成英文讲义P153.E3.3。
• 对英文讲义中的循序渐进实例“Disk Drive Read System ”，在100=a K 时，试采用微分反馈使系统的性能满足给定的设计指标。
2）第四章 线性系统的根轨迹法• 在MATLAB 环境下完成英文讲义P157.E4.5； • 利用MATLAB 绘制教材P181.4-5-(3）；• 在MATLAB 环境下选择完成教材第四章习题4-10或4-18，并对结果进行分析。
3）第五章 线性系统的频域分析法利用MATLAB 绘制本章作业中任意2个习题的频域特性曲线； 4）第六章 线性系统的校正利用MATLAB 选择设计本章作业中至少2个习题的控制器，并利用系统的单位阶跃响应说明所设计控制器的功能。
5）第七章 线性离散系统的分析与校正• 利用MATLAB 完成教材P383.7-20的最小拍系统设计及验证。
• 利用MATLAB 完成教材P385.7-25的控制器的设计及验证。
二．仿真实验时间安排及相关事宜1．依据课程教学大纲要求，仿真实验共6学时，教师可随课程进度安排上机时间，学生须在实验之前做好相应的准备，以确保在有限的机时内完成仿真实验要求的内容； 2．实验完成后按规定完成相关的仿真实验报告； 3．仿真实验报告请参照有关样本制作并打印装订；4．仿真实验报告必须在本学期第15学周结束之前上交授课教师。
自动控制原理MATLAB实验报告
实验一 典型环节得M ＡT ＬAB 仿真一、实验目得1.熟悉M ＡＴLAB 桌面与命令窗口,初步了解ＳIM ＵＬINK 功能模块得使用方法、 ２.通过观察典型环节在单位阶跃信号作用下得动态特性,加深对各典型环节响应曲线得理解。
3、定性了解各参数变化对典型环节动态特性得影响、二、实验原理1、比例环节得传递函数为其对应得模拟电路及SIMULI ＮK 图形如图１所示。
２.惯性环节得传递函数为uf C K R K R s C R R R Z Z s G 1,200,10012.021)(121121212===+-=+-=-=其对应得模拟电路及S ＩMULINK 图形如图2所示。
3.积分环节(I)得传递函数为其对应得模拟电路及SI ＭULINK 图形如图3所示。
4.微分环节(D)得传递函数为其对应得模拟电路及SIMULINK 图形如图4所示、 ﻩ５.比例+微分环节(PD)得传递函数为其对应得模拟电路及SIMUL ＩNK 图形如图５所示。
6.比例+积分环节(ＰI)得传递函数为其对应得模拟电路及SIMU ＬI ＮK 图形如图6所示。
图1 比例环节得模拟电路及SIMULINK 图形图3 积分环节得模拟电路及及SIMULINK 图形 图4 微分环节得模拟电路及及SIMULINK 图形图2惯性环节得模拟电路及SIMULINK 图形图5比例+微分环节得模拟电路及SIMULINK 图形曲线三、实验内容按下列各典型环节得传递函数,建立相应得SIMＵLINK仿真模型,观察并记录其单位阶跃响应波形。
①比例环节与;②惯性环节与③积分环节④微分环节⑤比例+微分环节(PD)与⑥比例+积分环节(PI)与四、实验报告记录各环节得单位阶跃响应波形,并分析参数对响应曲线得影响。
①比例环节:(如图7所示)(如图8所示)②惯性环节:(如图９所示) (如图１０所示)③积分环节: (如图11所示)④微分环节:(如图１2所示)⑤比例+微分环节(PＤ):(如图１3所示) (如图1４所示)⑥比例＋积分环节(PI):(如图15所示) (如图1６所示)实验二基于MATＬAB控制系统单位阶跃响应分析一、实验目得１、学会使用ＭATLＡB编程绘制控制系统得单位阶跃响应曲线。
自控实验-自动控制系统的MATLAB仿真分析
实验名称：自动控制系统的MATLAB仿真分析一、实验目的1.熟悉MATLAB在自动控制系统仿真中的应用；2.对自动控制系统进行仿真研究；3.掌握用MATLAB绘制自动控制系统根轨迹及对数频率特性的方法，掌握根据系统根轨迹及对数频率特性分析自动控制系统性能的方法。
二、实验设备1.计算机2.MATLAB软件三、实验内容1.用MATLAB提供的Simulink仿真软件工具对实验一中的各个典型环节及二阶系统进行阶跃响应仿真研究，将仿真获得的阶跃响应结果与模拟电路获得的阶跃响应结果进行比较。
（1）比例环节传递函数为200 ()51 G s=建立仿真模型，得到的输出结果如图所示：（2）积分环节传递函数为9.8 ()G ss=建立仿真模型，得到的输出结果如图所示：（3）一阶惯性环节传递函数为3.9 ()0.21G ss=+建立仿真模型，得到的输出结果如图所示：（4）比例积分环节传递函数为0.39781 ()0.102sG ss+=建立仿真模型，得到的输出结果如图所示：（5）比例微分环节传递函数为10 ()220s G ss=++建立仿真模型，得到的输出结果如图所示：（6）比例微分积分环节传递函数为51050 ()220sG ss s+=+++建立仿真模型，得到的输出结果如图所示：（7） 二阶系统的阶跃响应 ①0.325K ξ==传递函数为2()250()10250C s R s s s =++ 建立的仿真模型与阶跃响应仿真波形如下图所示：②0.510K ξ==传递函数为2()100()10100C s R s s s =++ 建立的仿真模型与阶跃响应仿真波形如下图所示：③0.75K ξ==传递函数为2()50()1050C s R s s s =++ 建立的仿真模型与阶跃响应仿真波形如下图所示：2. 单位负反馈系统的开环传递函数为：(1)()()(21)k s G s H s s s +=+仿真绘制K 从0~∞变化时的根轨迹，分析系统的稳定性。
自动控制原理Matlab仿真实验实验
源程序： numg=[1];deng=[500,0,0]; numc=[1,1];denc=[1,2]; [num1,den1]=series(numg,deng,numc,denc); [num,den]=cloop(num,den,-1); printsys(num,den) 运行结果： num/den = s + 1 ---------------------------500 s^3 + 1000 s^2 + 4 s + 4 3、 传递函数零、极点的求取 在命令窗口依次运行下面命令，并记录各命令运行后结果
示波器显示仿真结果：
(2)G2(s)= 1+1/2s Simulink 图形:
示波器显示仿真结果：
实验三 控制系统的时域分析
一、实验目的 学习利用 MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分 析系统的动态特性。 二、实验内容 (一) 稳定性 1、系统传函为 G(s) =
4
源程序： roots([1 7 3 5 2]) 运行结果： ans = - 6.6553 0.0327 + 0.8555i 0.0327 - 0.8555i - 0.4100 结论： 特征方程根具有正实部，因而系统是不稳定的 （二）阶跃响应 1. 二阶系统 G (s ) =
10 s + 2 s + 10
示波器显示仿真结果如下：
②惯性环节 G1(s)=1/s+1 和 G2(s)=1/0.5s+1。 Simulink 图形：
示波器显示仿真结果如下：
③积分环节 G1(s)=1/s。
Simulink 图形：
示波器显示仿真结果如下：
④微分环节 G1(s)= s。 Simulink 图形:
自动控制原理MATLAB仿真实验指导书(4个实验)
自动控制原理MATLAB仿真实验实验指导书电子信息工程教研室实验一典型环节的MA TLAB仿真一、实验目的1．熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法。
2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。
3．定性了解各参数变化对典型环节动态特性的影响。
二、SIMULINK的使用MATLAB中SIMULINK是一个用来对动态系统进行建模、仿真和分析的软件包。
利用SIMULINK功能模块可以快速的建立控制系统的模型，进行仿真和调试。
1．运行MA TLAB软件，在命令窗口栏“>>”提示符下键入simulink命令，按Enter键或在工具栏单击按钮，即可进入如图1-1所示的SIMULINK仿真环境下。
2．选择File菜单下New下的Model命令，新建一个simulink仿真环境常规模板。
图1-1 SIMULINK仿真界面图1-2 系统方框图3．在simulink仿真环境下，创建所需要的系统。
以图1-2所示的系统为例，说明基本设计步骤如下：1）进入线性系统模块库，构建传递函数。
点击simulink下的“Continuous”，再将右边窗口中“Transfer Fen”的图标用左键拖至新建的“untitled”窗口。
2）改变模块参数。
在simulink仿真环境“untitled”窗口中双击该图标，即可改变传递函数。
其中方括号内的数字分别为传递函数的分子、分母各次幂由高到低的系数，数字之间用空格隔开；设置完成后，选择OK，即完成该模块的设置。
3）建立其它传递函数模块。
按照上述方法，在不同的simulink的模块库中，建立系统所需的传递函数模块。
例：比例环节用“Math”右边窗口“Gain”的图标。
4）选取阶跃信号输入函数。
用鼠标点击simulink下的“Source”，将右边窗口中“Step”图标用左键拖至新建的“untitled”窗口，形成一个阶跃函数输入模块。
自动控制原理 - Matlab实验分析完整报告【优秀版】
利用MATLAB 进行自动控制原理的一些分析来自：我是痕痕的弟弟1、已知三阶系统开环传递函数为G （S ）=)232(2723+++s s s ，利用MATLAB 程序，画出系统的奈圭斯特图，求出相应的幅值裕量和相位裕量。
解： 程序如下：G=tf(3.5,[1,2,3,2]); %得到系统的传递函数 subplot(1,2,1);nyquist(G); %绘制奈圭斯特曲线gridxlabel('Real Axis')ylabel('Image Axis')[Gm,Pm,Weg,Wep]=margin(G) %求幅值和相角余度及对应的频率G_ c=feedback(G,1); %构造单位反馈系统subplot(1,2,2); %绘制单位阶跃响应曲线step(G_ c)gridxlabel('Time(secs)')ylabel('Amplitude')显示结果：Gm=1.1433 Pm=7.1688 Wcg=1.7323 Wcp=1.6541系统的奈圭斯特图如下（从MATLAB截图显示）：2、绘制二阶环节的伯特图。
解：MATLAB程序如下：figure('pos',[30 100 260 400],'color','w');axes('pos',[0.15 0.2 0.7 0.7]);wn=1w=[0,logspace(-2,2,200)]; %得到对数频率数组for zeta=[0.1 0.5 1 2] %分别绘制阻尼系数为0.1、0.5、1、2的二阶环节bode 图G=tf(1,[wn^-2 2*zeta/wn 1]); bode(G ,w); hold on end;grid程序运行后得到如下图（MATLAB 截图显示）：从图中可以看出，频率w 接近Wn 时产生谐振，阻尼比的大小确定谐振峰值的大小，阻尼比越小，谐振峰值越大。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
自动控制原理实验报告学院电子信息与电气工程学院实验一 MATLAB 及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点1、 系统的典型响应有哪些？2、 如何判断系统稳定性？3、 系统的动态性能指标有哪些？ 三、实验方法（一） 四种典型响应1、 阶跃响应：阶跃响应常用格式：1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。
2、),(Tn sys step ；表示时间范围0---Tn 。
3、),(T sys step ；表示时间范围向量T 指定。
4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。
2、 脉冲响应：脉冲函数在数学上的精确定义：0,0)(1)(0〉==⎰∞t x f dx x f其拉氏变换为：)()()()(1)(s G s f s G s Y s f ===所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式： ① )(sys impulse ； ②);,();,(T sys impulse Tn sys impulse③ ),(T sys impulse Y =（二） 分析系统稳定性 有以下三种方法：1、 利用pzmap 绘制连续系统的零极点图；2、 利用tf2zp 求出系统零极点；3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.四、实验内容 (一) 稳定性1． 系统传函为()27243645232345234+++++++++=s s s s s s s s s s G ，试判断其稳定性2． 用Matlab 求出253722)(2342++++++=s s s s s s s G 的极点。
%Matlab 计算程序num=[3 2 5 4 6];den=[1 3 4 2 7 2];G=tf(num,den);pzmap(G);p=roots(den)运行结果： p =-1.7680 + 1.2673i -1.7680 - 1.2673i 0.4176 + 1.1130i 0.4176 - 1.1130i -0.2991P ole-Zero MapReal AxisI m a g i n a r y A x i s-2-1.5-1-0.500.5-1.5-1-0.50.511.5图1-1 零极点分布图由计算结果可知，该系统的2个极点具有正实部，故系统不稳定。
%求取极点num=[1 2 2];den=[1 7 3 5 2];p=roots(den)运行结果： p =-6.6553 0.0327 + 0.8555i 0.0327 - 0.8555i -0.4100故253722)(2342++++++=s s s s s s s G 的极点s1=-6.6553 , s2=0.0327 + 0.8555i ,s3= 0.0327 - 0.8555i , s4=-0.41（二）阶跃响应1. 二阶系统()102102++=s s s G1）键入程序，观察并记录单位阶跃响应曲线2）计算系统的闭环根、阻尼比、无阻尼振荡频率，并记录 3）记录实际测取的峰值大小、峰值时间及过渡过程时间，并填表：由图1-3及其相关理论知识可填下表：3//πωπ==d p t =1.04724）修改参数，分别实现1=ζ和2=ζ的响应曲线，并记录5）修改参数，分别写出程序实现0121w w n =和022w w n =的响应曲线，并记录%单位阶跃响应曲线num=[10];den=[1 2 10];step(num,den);title('Step Response of G(s)=10/(s^2+2s+10)');4.52%(00.9)3.55%n s nt ζωζζω⎧∆=⎪⎪=<<⎨⎪∆=⎪⎩01234560.20.40.60.811.21.4Step Response of G(s)=10/(s 2+2s+10)Time (sec)A m p l i t u d e图1-2 二阶系统()102102++=s s s G 单位阶跃响应曲线%计算系统的闭环根、阻尼比、无阻尼振荡频率 num=[10];den=[1 2 10];G=tf(num,den); [wn,z,p]=damp(G)运行结果： wn =3.1623 3.1623 z =0.3162 0.3162 p =-1.0000 + 3.0000i -1.0000 - 3.0000i由上面的计算结果得系统的闭环根s= -1±3i ，阻尼比=ς3162.0、无阻尼振荡频率1623.3=n ω实验二 MATLAB 及仿真实验（控制系统的根轨迹分析）一 实验目的1．利用计算机完成控制系统的根轨迹作图 2．了解控制系统根轨迹图的一般规律 3．利用根轨迹图进行系统分析 二 预习要点1. 预习什么是系统根轨迹？2. 闭环系统根轨迹绘制规则。
三 实验方法（一） 方法：当系统中的开环增益k 从0到变化时，闭环特征方程的根在复平面上的一组曲线为根轨迹。
设系统的开环传函为：)()()(0s Q s N ks G =，则系统的闭环特征方程为：0)()(1)(10=+=+s Q s N ks G 根轨迹即是描述上面方程的根，随k 变化在复平面的分布。
（二） MATLAB 画根轨迹的函数常用格式：利用Matlab 绘制控制系统的根轨迹主要用pzmap ，rlocus ，rlocfind ，sgrid 函数。
1、零极点图绘制❑ [p,z]=pzmap(a,b,c,d)：返回状态空间描述系统的极点矢量和零点矢量，而不在屏幕上绘制出零极点图。
❑ [p,z]=pzmap(num,den)：返回传递函数描述系统的极点矢量和零点矢量，而不在屏幕上绘制出零极点图。
❑ pzmap(a,b,c,d)或pzmap(num,den)：不带输出参数项，则直接在s 复平面上绘制出系统对应的零极点位置，极点用×表示，零点用o 表示。
❑ pzmap(p,z)：根据系统已知的零极点列向量或行向量直接在s 复平面上绘制出对应的零极点位置，极点用×表示，零点用o 表示。
2、根轨迹图绘制❑ rlocus(a,b,c,d)或者rlocus(num,den)：根据SISO 开环系统的状态空间描述模型和传递函数模型，直接在屏幕上绘制出系统的根轨迹图。
开环增益的值从零到无穷大变化。
❑ rlocus(a,b,c,d,k)或rlocus(num,den,k)： 通过指定开环增益k 的变化范围来绘制系统的根轨迹图。
❑ r=rlocus(num,den,k) 或者[r,k]=rlocus(num,den) ：不在屏幕上直接绘出系统的根轨迹图，而根据开环增益变化矢量k ，返回闭环系统特征方程1＋k*num(s)/den(s)=0的根r ，它有length(k)行，length(den)-1列，每行对应某个k 值时的所有闭环极点。
或者同时返回k 与r 。
❑ 若给出传递函数描述系统的分子项num 为负，则利用rlocus 函数绘制的是系统的零度根轨迹。
（正反馈系统或非最小相位系统） 3、rlocfind()函数❑ [k,p]=rlocfind(a,b,c,d)或者[k,p]=rlocfind(num,den) 它要求在屏幕上先已经绘制好有关的根轨迹图。
然后，此命令将产生一个光标以用来选择希望的闭环极点。
命令执行结果：k 为对应选择点处根轨迹开环增益；p 为此点处的系统闭环特征根。
❑ 不带输出参数项[k,p]时，同样可以执行，只是此时只将k 的值返回到缺省变量ans中。
4、sgrid()函数❑ sgrid ：在现存的屏幕根轨迹或零极点图上绘制出自然振荡频率wn 、阻尼比矢量z对应的格线。
❑ sgrid(‘new ’)：是先清屏，再画格线。
❑ sgrid(z,wn)：则绘制由用户指定的阻尼比矢量z 、自然振荡频率wn 的格线。
四 实验内容 1．()()()21++=s s s k s G g 要求：二、 记录根轨迹的起点、终点与根轨迹的条数； 三、 确定根轨迹的分离点与相应的根轨迹增益； 四、 确定临界稳定时的根轨迹增益gL k%Matlab 计算程序z=[];p=[0 -1 -2];k=1;G=zpk(z,p,k);figure(1);pzmap(G) figure(2);rlocus(G)title('实验2.1所作曲线');（a ）由图2-2知，起点分别为0，-1，-2，终点为无穷远处，共三条根轨迹. (b) 结合图2-3和图2-5得分离点d=-0.4226，相应的根轨迹增益k=-0.3849. (c) 结合图2-3和图2-4得临界稳定时的根轨迹增益gL k =6.01P ole-Zero MapReal AxisI m a g i n a r y A x i s图2-1 零、极点分布图-7-6-5-4-3-2-1012-5-4-3-2-1012345实验2.1所作曲线Real AxisI m a g i n a r y A x i s图2-2 根轨迹图-7-6-5-4-3-2-112-5-4-3-2-1012345实验2.1所作曲线Real AxisI m a g i n a r y A x i s图2-3 根轨迹图（2）%求临界稳定时的根轨迹增益Kglz=[];p=[0 -1 -2];k=1;G=zpk(z,p,k); rlocus(G)title('实验2.1 临界稳定时的根轨迹增益Kgl'); [k,p]=rlocfind(G)运行结果：Select a point in the graphics windowselected_point =0.0059 + 1.4130i k =6.0139 p =-3.0013 0.0006 + 1.4155i 0.0006 - 1.4155i-5-4-3-2-1012345实验2.1 临界稳定时的根轨迹增益KglReal AxisI m a g i n a r y A x i s图2-4 根轨迹图（3）实验三 MATLAB及仿真实验（控制系统的频域分析）一实验目的1. 利用计算机作出开环系统的波特图2. 观察记录控制系统的开环频率特性3. 控制系统的开环频率特性分析二预习要点1.预习Bode图和Nyquist图的画法；2.映射定理的内容；3.Nyquist稳定性判据内容。