下载文档原格式
自动控制原理实验报告 ——控制系统的阶跃响应09021209 侯竟骁一、实验目的 1、观察学习控制系统的单位阶跃响应； 2、记录单位阶跃响应曲线； 3、掌握时间响应分析的一般方法。
二、实验步骤 1、开机执行程序c:\ml\bin\matlab-s.exe（用鼠标双击图标）进入MATLAB 命令窗口：“Command Windows ”。
2、建立系统模型 在MATLAB 命令窗口上，以立即命令方式建立系统的传递函数。
在MATLAB 下，系统传递函数有三种描述方式，在实验中只用到多项式模型和零点极点模型 多项式模型)()()(s s s den num G =式中“num(s)”表示分子多项式的系数，“den(s)”表示分母多项式的系数，全部按照复自变量s 的降幂排列，以行向量的方式输入。
例如，程序为 num=[0 1 3]; 分子多项式系数 den=[1 2 2 1]; 分母多项式系数 printsys(num,den); 构造传递函数G(s)并显示 零点极点模型∏∏--=ni mj s s s )()()(p z k G式中，k 为增益值，z j 为第j 个零点值，p i 为第i 个零点值。
例如，程序为 k=2; 赋增益值，标量 z=[1]; 赋零点值，向量 p=[-1 2 -3]; 赋极点值，向量 [num,den]=zp2tf(z,p,k); 零点极点模型转换成多项式模型 printsys(num,den); 构造传递函数G(s)并显示给定系统传递函数)(s G 的多项式模型，求系统的单位脉冲响应。
传递函数为)()()(s s s den num G =式中，num (s)为系统传递函数)(s G 的分子多项式系数向量，den (s)为系统传递函数)(s G 的分母多项式系数向量。
函数格式1：给定num 、den 求系统的阶跃响应。
时间向量t 的范围自动设定。
函数格式2：时间向量t 的范围可以由人工给定。
实验一 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的控制系统频域及根轨迹分析二、实验目的：（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 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点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 分析一、 实验目的1)、掌握如何使用Matlab 进行系统的时域分析 2)、掌握如何使用Matlab 进行系统的频域分析 3)、掌握如何使用Matlab 进行系统的根轨迹分析 4)、掌握如何使用Matlab 进行系统的稳定性分析 5)、掌握使用Bode 图法进行控制系统设计的方法 二、 实验内容 1、时域分析法根据下面传递函数模型：绘制其单位阶跃响应曲线并从图上读取最大超调量，绘制系统的单位脉冲响应、零输入响应曲线。
1)、某单位负反馈系统传递函数为：8106)65(5)(232+++++=s s s s s s Gt (seconds)c (t )t (seconds)c (t )结论：2)、典型二阶系统传递函数为：2222)(nn nc s s s G ωξωω++= 当ζ=0.7,ωn 取2、4、6、8、10、12的单位阶跃响应。
Step ResponseTime (seconds)00.51 1.52 2.53 3.54结论：3)、典型二阶系统传递函数为：2222)(nn nc s s s G ωξωω++= 当ωn =6,ζ取0.2、0.4、0.6、0.8、1.0、1.5、2.0的单位阶跃响应。
Time (seconds)结论：2、频率分析法根据下面传递函数模型，绘制出系统的频率响应曲线，包括Bode 图和Nyquist 图，并从图上读取相角交接频率、截止频率，并求出幅值裕度和相角裕度。
1)、典型二阶系统传递函数为：2222)(nn nc s s s G ωξωω++= 当ζ=0.7,ωn 取2)、4)、6)、8、1)0、1)2)的伯德图和奈奎斯特图。
Wn=2M a g n i t u d e (d B )10101010P h a s e (d e g )Bode DiagramGm = Inf dB (at Inf rad/s) , Pm = 164 deg (at 0.4 rad/s)Frequency (rad/s)Real AxisI m a g i n a r y A x i sWn=4M a g n i t u d e (d B )10101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i sWn=6M a g n i t u d e (d B )1010101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i sWn=8M a g n i t u d e (d B )1010101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i sWn=10M a g n i t u d e (d B )1010101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i sWn=12M a g n i t u d e (d B )1010101010P h a s e (d e g )Bode DiagramFrequency (rad/s)Real AxisI m a g i n a r y A x i s2)、典型二阶系统传递函数为：2222)(nn nc s s s G ωξωω++= 当ωn =6,ζ取0.2、0.4、0.6、0.8、1.0、1.5、2.0的伯德图和奈奎斯特图。
《自动控制原理》Matlab求解控制系统时域响应实验一、实验目的(1)学习控制系统在阶跃信号、脉冲信号、速度函数、加速度函数等不同输入信号下的响应(2)通过计算机的仿真图形观测二阶控制系统的时域响应曲线。
(3)改变ξ和ωn，观测参数变化时对典型环节时域响应的影响。
二、实验仪器Matlab7.0，计算机三、实验原理1、求解连续系统时域响应的MATLAB函数Impulse函数计算/绘制连续系统的单位脉冲响应句法：step函数lsim函数计算/绘制连续系统在0～t时间内的对任意输入u的响应。
上述函数句法中带输出变量引用时，得到系统的时域响应数据而不直接显示响应曲线，此时需通过绘图命令plot查看响应曲线。
四、实验内容及步骤（分析：无阻尼固有频率一定，阻尼比越小，系统超调量越大，振荡越厉害）五、实验原始数据记录与数据处理当ζ=0.7,ωn取2、4、6、8、10、12时的单位阶跃响应曲线。
六、实验结果与分析讨论无阻尼固有频率一定，阻尼比越小，系统超调量越大，振荡越厉害。
阻尼比一定，无阻尼固有频率越大，系统能更快达到稳定值，响应的快速性越好七、结论控制系统在阶跃信号、脉冲信号、速度函数、加速度函数等不同输入信号下的响应不同，无阻尼固有频率一定，阻尼比越小，系统超调量越大，振荡越厉害。
阻尼比一定，无阻尼固有频率越大，系统能更快达到稳定值，响应的快速性越好。
八、实验心得体会（可略）通过这次实验我学到了很多，首先是对单位脉冲响应，单位加速度响应，单位速度响应，单位阶跃响应有了一定的认知，对matlab软件有了更深的认识，会使用的功能越来越多了，希望在今后的学习生涯中能学到更多的知识。
仿真实验○一：控制系统的时域分析一、实验目的：1．观察控制系统的时域响应；2．记录单位阶跃响应曲线；3．掌握时间响应分析的一般方法；4．初步了解控制系统的调节过程。
二、实验步骤：1．开机进入Matlab6.1运行界面。
2．Matlab指令窗："Command Window". 运行指令：con_sys; 进入本次实验主界面。
3．分别双击上图中的三个按键，依次完成实验内容。
4．本次实验的相关Matlab函数：tf([num],[den])可输入一传递函数。
step(G,t)在时间范围t秒内，画出阶跃响应图。
三、实验内容：1、观察一阶系统G=1/(T+s) 的时域响应：取不同的时间常数T，分别观察该系统的脉冲响应、阶跃响应、斜坡响应以及单位加速度响应。
结论：时间常数越小，响应越迅速。
2、二阶系统的时域性能分析：（1）调节时间滑块，使阶跃响应最终出现稳定值。
（2）结合系统的零极点图，观察自然频率与阻尼比对极点位置的影响。
（3）结合时域响应图，观察自然频率与阻尼比对阶跃响应的影响。
结论：阻尼比越小，极点越靠近虚轴，超调量减小，但响应速度变慢。
自然频率减小，极点靠近虚轴，响应速度减小，超调几乎不变。
（4）调节自然频率与阻尼比，要求：Tr<0.56s ，Tp<1.29s，Ts<5.46，超调不大于5％.记录下满足上述要求的自然频率与阻尼比。
调节完成之后的响应曲线如图。
此时自然频率为14.5872rad/sec,阻尼比为0.77456。
各项参数完全满足要求。
3、结合《自动控制原理》一书，Page 135，题3_10. 分别观察比例_微分与测速反馈对二阶系统性能的改善。
（1）.按原始的调节参数输入，调节时间滑块，使阶跃响应最终出现稳定值。
（2）采用不同的G输入，观察各项性能指数。
结论：增大分母中间的参数，相当于增大系统阻尼比，从而减小超调量(3).分别取不同的K3，观察比例_微分控制对系统性能的改善。
自动控制原理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 使用方法和程序设计 热动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 对系统性能的影响。
一、实验目的：一、实验目的：1、 目的：用MATLAB 的控制系统工具箱中线性系统时域响应分析模块的功能，更好的借助MATLAB 和控制工具箱非常方便的进行自动控制系统的时域分析、频域分析和系统校正。
熟悉MA TLAB 平台的使用平台的使用2、 手段：MA TLAB 应用软件应用软件3、 效果：出现系统的单位阶跃响应、出现系统的单位脉冲响应、系统的零输入响应等的效果图，帮助更好的分析模块功能的效果图，帮助更好的分析模块功能二、相关函数：二、相关函数：Step （num,den ）num----分子的系数；分子的系数；den-----分母的系数分母的系数三、实验内容：三、实验内容：a) 绘制时间常数T=0.5s ，1s ，2s 时惯性环节的单位阶跃响应曲线族。
时惯性环节的单位阶跃响应曲线族。
解：惯性环节的传递函数为G(s)=1/Ts+1在MA TLAB 窗口输入下列命令窗口输入下列命令num=1;den1=[0.5,1];den2=[1,1];den3=[2,1];[y1,x,t1]=step(num,den1);[y2,x,t2]=step(num,den2);[y3,x,t3]=step(num,den3);Plot(t1,y1,Plot(t1,y1,’’-b’b’,t2,y2,,t2,y2,,t2,y2,’’-r -r’’,t3,y3,,t3,y3,’’-g -g’’)Title(Title(‘‘Step Response Step Response’’);Xlabel(Xlabel(‘‘Time(secs)Time(secs)’’);Ylabel(Ylabel(‘‘Amplitude Amplitude’’)效果图如下：2、二阶系统的传递函数为G(s)= wn*wn/s*s+2 ζw ns+wn*wn 用时域分析模块中的函数可以方便的分析无阻尼自然振荡角频率Wn 、阻尼比ζ对系统暂态响应性能的影响。
实验名称：自动控制系统的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在自动控制理论中的应用二、实验目的熟悉并掌握MATLAB在自动控制理论中的数学计算和绘图功能、simulink仿真功能等。
三、实验内容（所有）1、传递函数的描述法2、自动控制系统结构框图的模型表示3、线性系统的时域分析4、线性系统的频域分析5、线性系统的根轨迹分析6、状态空间描述法四、实验步骤（部分）【实验一】实验目的：观察二阶振荡环节中，参数ζ和Wn分别变化时对输出波形的影响。
实验内容：二阶标准传递函数：（1）令Wn不变，ζ取不同的值。
（0<ζ<1）练习：令Wn=5不变，ζ等于0.2和0.707结论：Wn相同，ζ等于0.707时响应更快（2）令ζ不变，Wn取不同的值。
练习：令ζ=0.25，Wn等于1和10结论：ζ相同，Wn越大响应越快实验代码：>> num=[25];>> den=[1,2,25];>> G1=tf(num,den);>> den2=[1,7.07,25];>> G2=tf(num,den2);>> num2=[1];>> den3=[1,0.5,1];>> den4=[1,5,100];>> G3=tf(num2,den3);>> num3=[100];>> G4=tf(num3,den4);>> step(G1);hold on>> step(G2);hold on>> step(G3);hold on>> step(G4);hold on实验结果：ζ=0.2ζ=0.707Wn=10Wn=1实验结论：我们可以很直观的看到，当Wn相同，ζ等于0.707时比ζ等于0.2时响应更快；ζ相同，Wn越大响应越快。
但是因为ζ范围是0到1，而ζ的取值到底是怎么样影响系统输出的，是否是越大响应越快，就可以通过下面一个实验来进行验证。
《自动控制原理》Matlab求解控制系统数学模型实验一、实验目的(1)熟练运用matlab软件，求解控制系统数学模型(2)掌握传递函数在matlab中的表达方法(3)掌握matlab求解拉氏变换和反变换(4)掌握matlab求系统极值点和零点判断系统稳定性二、实验仪器装配Matlab7.0的计算机三、实验原理传递函数在matlab中的表达方法控制系统的传递函数模型为：在MATLAB中，分子/分母多项式通过其系数行向量表示，即：num = [b0 b1 … bm]den = [a0 a1 … an]此时，系统的传递函数模型用tf函数生成，句法为：sys=tf(num, den)其中，sys为系统传递函数。
如：num = [1 5 0 2]; den = [2 3 15 8];则：sys=tf(num, den)输出为：Transfer function:传递函数的转换[num,den]=zp2tf(z,p,k)[z,p,k]=tf2zp(num,den)实际系统往往由多个环节通过串联、并联及反馈方式互连构成。
MATLAB提供的三个用于计算串联、并联及反馈连接形成的新系统模型的函数。
四、实验内容及步骤2、用MATLAB展求拉氏变换和反变换在MATLAB中，多项式通过系数行向量表示，系数按降序排列如要输入多项式：x4-12x3+25x+126>> p=[1 -12 0 25 126]-p = 1 -12 0 25 1263、连续系统稳定性分析的MATLAB函数roots函数：求多项式的根句法： r=roots(p)其中，r为由多项式根组成的列向量。
➢pole函数：计算系统的极点句法： p=pole(sys)其中，p为由极点组成的列向量zero函数：计算系统的零点句法： r=zero(sys) 或 [z, k]=zero(sys)其中，r为由多项式根组成的列向量。
k为零极点增益模型之增益pzmap函数：绘制零极点分布图句法： pzmap(sys) 或 [p,z] = pzmap(sys)五、实验原始数据记录与数据处理在MATLAB中，多项式通过系数行向量表示，系数按降序排列如要输入多项式：x4-12x3+25x+126>> p=[1 -12 0 25 126]-p = 1 -12 0 25 126六、实验结果与分析讨论七、结论掌握 MATLAB命令窗口的基本操作;掌握MATLAB 建立控制系统数学模型的命令及模型相互转换的方法;掌握了使用各种函数命令建立控制系统数学模型.八、实验心得体会（可略）通过该试验我们熟悉 MATLAB 实验环境，掌握 MATLAB命令窗口的基本操作;掌握MATLAB 建立控制系统数学模型的命令及模型相互转换的方法;掌握了使用各种函数命令建立控制系统数学模型:完成实验的范例题和自我实践，并记录结果;编写M文件程序，完成简单连接的模型等效传递函数，并求出相应的零极点。
利用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 时产生谐振，阻尼比的大小确定谐振峰值的大小，阻尼比越小，谐振峰值越大。
