惯性环节振荡环节二三阶系统的MATLAB仿真

第5章频率特性法频域分析法是一种图解分析法，可以根据系统的开环频率特性去判断闭环系统的性能，并能较方便地分析系统参量对系统性能的影响，从而指出改善系统性能的途径，已经发展成为一种实用的工程方法，其主要内容是：1）频率特性是线性定常系统在正弦函数作用下，稳态输出与输入的复数之比对频率的函数关系。

频率特性是传递函数的一种特殊形式，也是频域中的数学模型。

频率特性既可以根据系统的工作原理，应用机理分析法建立起来，也可以由系统的其它数学模型（传递函数、微分方程等）转换得到，或用实验法来确定。

2）在工程分析和设计中，通常把频率特性画成一些曲线。

频率特性图形因其采用的坐标不同而分为幅相特性(Nyquist图)、对数频率特性(Bode图)和对数幅相特性(Nichols图)等形式。

各种形式之间是互通的，每种形式有其特定的适用场合。

开环幅相特性在分析闭环系统的稳定性时比较直观，理论分析时经常采用；波德图可用渐近线近似地绘制，计算简单，绘图容易，在分析典型环节参数变化对系统性能的影响时最方便；由开环频率特性获取闭环频率指标时，则用对数幅相特性最直接。

3）开环对数频率特性曲线(波德图)是控制系统分析和设计的主要工具。

开环对数幅频特性L(ω)低频段的斜率表征了系统的型别(v)，其高度则表征了开环传递系数的大小，因而低频段表征系统稳态性能；L(ω)中频段的斜率、宽度以及幅值穿越频率，表征着系统的动态性能；高频段则表征了系统抗高频干扰的能力。

对于最小相位系统，幅频特性和相频特性之间存在着唯一的对应关系，根据对数幅频特性，可以唯一地确定相应的相频特性和传递函数。

4）奈奎斯特稳定性判据是利用系统的开环幅相频率特性G(jω)H(jω)曲线，又称奈氏曲线，是否包围GH平面中的(－l，j0)点来判断闭环系统的稳定性。

利用奈奎斯特稳定判据，可根据系统的开环频率特性来判断闭环系统的稳定性，并可定量地反映系统的相对稳定性，即稳定裕度。

稳定裕度通常用相角裕量和幅值裕量来表示。

一、 实验名称：典型环节的时域分析和频域分析二、实验目的：(1) 理解、掌握matlab 模拟典型环节的根本方法，包括：比例环节、积分环节、一阶微分环节、惯性环节和振荡环节等。

(2) 熟悉各种典型环节的阶跃响应曲线和频域响应曲线 (3) 理解参数变化对动态特性的影响三、 实验要求：(1) 一人一机，独立完成实验内容 。

(2) 根据实验结果完成实验报告，并用A4纸打印后上交。

四、 时间：2022年11月21日 五、 地点：信自楼234实验报告：一、比例环节的时域分析和频域分析 比例环节的传递函数：()G s k(1) 当k=1:3:10时，绘制系统的阶跃响应曲线，分析ｋ值的影响情况。

程序：for k=1:3:10;num=k;den=1;G=tf(num,den);figure(1);step(G); hold on; %翻开第1个图形窗口，绘制系统的阶跃响应曲线 endfigure(1); legend('k=1','k=4','k=7','k=10'); 曲线：结果分析：时域响应的结果就是把输入信号放大k 倍。

如图，输入信号为幅值为1的阶跃信号，因此，输出是幅值为k 的阶跃信号。

程序：for k=1:3:10;num=k;den=1;G=tf(num,den);figure(1);bode(G);hold on; %翻开第1个图形窗口，绘制系统的阶跃响应曲线 endfigure(1); legend('k=1','k=4','k=7','k=10');曲线：结果分析：比例环节对幅频有影响，输出信号的幅值为输入信号的20*lgk倍。

比例环节对相位没有影响，如图显示，相位特性为一条0度的程度线。

二、积分环节的时域分析和频域分析积分环节的传递函数：1 ()G ss=(1) 当k=1:3:10时，绘制系统()kG ss=的阶跃响应曲线，分析曲线特点。

实验一典型环节的模拟研究及阶跃响应分析1、比例环节可知比例环节的传递函数为一个常数：当Kp 分别为，1，2时，输入幅值为的正向阶跃信号，理论上依次输出幅值为，，的反向阶跃信号。

实验中，输出信号依次为幅值为，，的反向阶跃信号，相对误差分别为1.8%,2.2%,0.2%.在误差允许范围内可认为实际输出满足理论值。

2、 积分环节积分环节传递函数为：〔1〕T=0.1(0.033)时，C=1μf(0.33μf)，利用MATLAB ，模拟阶跃信号输入下的输出信号如图：与实验测得波形比较可知，实际与理论值较为吻合，理论上时的波形斜率近似为时的三倍，实际上为，在误差允许范围内可认为满足理论条件。

3、 惯性环节惯性环节传递函数为：K = R f /R 1,T = R f C,（1） 保持K = R f /R 1= 1不变，观测秒，秒〔既R 1 = 100K,C = 1μf ，μf 〕时的输出波形。

利用matlab 仿真得到理论波形如下：时t s 〔5%〕理论值为300ms,实际测得t s =400ms 相对误差为：〔400-300〕/300=33.3%,读数误差较大。

K 理论值为1，实验值，相对误差为〔〕/2.28=7%与理论值较为接近。

时t s 〔5%〕理论值为30ms,实际测得t s =40ms 相对误差为：〔40-30〕/30=33.3% 由于ts 较小，所以读数时误差较大。

K 理论值为1，实验值，相对误差为〔〕/2.28=7%与理论值较为接近（2） 保持T = R f s 不变，分别观测K = 1，2时的输出波形。

K=1时波形即为〔1〕中时波形K=2时，利用matlab 仿真得到如下结果：t s 〔5%〕理论值为300ms,实际测得t s =400ms相对误差为：〔400-300〕/300=33.3% 读数误差较大K 理论值为2，实验值， 相对误差为〔〕/2=5.7%if i o R RU U -=1TS K)s (R )s (C +-=与理论值较为接近。

典型环节时域特性的仿真实验1、通过观察典型环节在单位阶跃信号作用下的动态特性，熟悉各种典型环节的响应曲线。

2 、定性了解各参数变化对典型环节动态特性的影响。

3、初步了解Matlab中Simulink的使用方法。

研究典型环节(比例、积分、微分、惯性、二阶)在阶跃输入信号及白噪声干扰信号输入的响应。

1.1 运行Matlab，在命令窗口“Command Window”下键入“Simulink”后回车，则打开相应的系统模型库；或者点击菜单上的“Simulink”图标，进入系统仿真模型库。

然后点击左上角“创建新文件图标”，打开模型编辑窗口。

1.2 调出模块在系统仿真模型库中，把要求的模块都放置在模型编辑窗口里面。

从信号源模块包(Sources)中拖出1个阶跃信号(step)和1个白噪声信号发生器 (band-limited white noise) ；从输出模块包(Sinks)中拖出1个示波器(Scope)；从连续系统典型环节模块包(Continuous) 中拖出1个微分环节(Derivative)和3个传函环节(Transfer Fcn)；从数学运算模块包(Math Operations)中拖出1个比例环节(Gain)和1个加法器 (Sum) ；从信号与系统模块包(Signals Routing) 拖出1个汇流排(Mux)；所有模块都放置在模型编辑窗口里面。

1.3 模块参数设置双击打开3个传函环节(Transfer Fcn)，通过设定参数 (参照图1的数据)，分别构成积分、惯性和二阶环节；打开比例环节，设定比例增益为2；打开白噪声信号发生器，设定功率(Noise power)为0.0001，采样时间(Sample time) 为0.05。

1.4 模块连接将各模块连接成如图1所示的仿真模型系统。

图1仿真模型系统22.1 双击Scope打开示波器，点击按钮“”启动仿真，画出输入信号波形图。

2.2 将比例环节的输出端接到汇流排(如图1所示)，打开示波器, 点击按钮“”启动仿真，观察比例环节的阶跃响应及对白噪声信号是否敏感，然后画出波形图。

实验一 一阶系统及二阶系统时域特性MatLab 仿真实验一、实验目的1、使学生通过实验中的系统设计及理论分析方法，帮助学生进一步理解自动控制系统的设计与分析方法。

2、熟悉仿真分析软件。

3、利用Matlab 对一、二阶系统进行时域分析。

4、掌握一阶系统的时域特性，理解常数T 对系统性能的影响。

5、掌握二阶系统的时域特性，理解二阶系统重要参数对系统性能的影响。

二、实验设备计算机和Matlab 仿真软件。

三、实验内容1、一阶系统时域特性 一阶系统11)(+=Ts s G ，影响系统特性的参数是其时间常数T ，T 越大，系统的惯性越大，系统响应越慢。

Matlab 编程仿真T=0.4，1.2，2.0，2.8，3.6，4.4系统单位阶跃响应。

2、二阶系统时域特性a 、二阶线性系统 16416)(2++=s s s G 单位脉冲响应、单位阶跃响应、单位正弦输入响应的 Matlab 仿真。

b 、下图为具有一微分负反馈的位置随动系统框图，求出系统的闭环传递函数，根据系统瞬态性能指标的定义利用Matlab 分别计算微分反馈时间常数τ为0，0.0125，0.025时系统的上升时间、峰值时间、最大超调量和调整时间。

C 、二阶线性系统3612362++s s ξ，当ξ为0.1，0.2，0.5，0.7，1.0，2.0时，完成单位阶跃响应的Matlab 仿真，分析ξ值对系统响应性能指标的影响。

四、实验要求1、进入机房，学生要严格遵守实验室规定。

2、学生独立完成上述实验，出现问题，教师引导学生独立分析和解决问题。

3、完成相关实验内容，记录程序，观察记录响应曲线，响应曲线及性能指标进行比较，进行实验分析4、分析系统的动态特性。

5、并撰写实验报告，按时提交实验报告。

五、Matlab 编程仿真并进行实验分析1、一阶系统时域特性实验代码：运行曲线：实验分析：由上图分析可知，一阶系统时间常数越大，图像图线越晚达到常值输出，即时间常数T影响系统参数，时间常数越大，系统的惯性越大，系统响应越慢。

基于MATLAB 的自动控制原理实验仿真系统的设计基于MATLAB 的自动控制原理实验仿真系统的设计基于MATLAB 的自动控制原理实验仿真系统的设计田晴，张茁（河北联合大学电气工程学院，河北唐山063000）摘要：MATLAB的图形用户界面（GUI）为基于对象的可视化编程，本文以此为基础，进行了自动控制原理实验仿真系统的软件设计，仿真实验系统囊括了控制理论的大部分要点，较实验室传统性实验更全面、具体。

应用GUIDE的设计，该系统操作简单，知识点表现清晰明了，学生能够在轻松的环境下，且不受地域时域的限制，掌握控制理论知识。

关键词：自动控制原理；实验仿真系统；MATLAB；GUIDE基金项目：河北联合大学教育教学改革项目，项目编号：Y1340-10一、引言《自动控制原理》是自动化专业的基础课程，是控制科学与工程学科的一门方法论课程，主要培养学生掌握控制系统的分析和设计方法，其内容之多，理论性之强，决定了课程学习的难度。

而实验课作为课堂教学的辅助内容，是培养学生自主性和创新性的重要环节。

目前实验室的实验教学采用模拟电路实验台，将集成电路模块进行连线，形成典型系统，通过示波器观察响应曲线。

传统性实验训练了学生对以传递函数为核心的控制系统与模拟电路系统之间的联系的认识和实践能力，但也存在其局限性：（1）价格昂贵，占地很大，因为实验台有限，只能几个同学共用一个实验台，难以满足几百学生实验教学的需要；（2 ）同一个实验，教师要对学生分拨讲解，重复进行，浪费人力、物力；（3 ）由于实验设备的长期工作，造成电容积分饱和，致使出现实验误差；（4 ）实验设备高度集成，操作性复杂，参数变化有限，局限了综合性设计性实验的开展；（5 ）实验室难以做到全开放性，学生实验受地域和时域的限制。

因此，研制实验仿真系统是解决上述问题的有效措施。

二、软件的总体结构设计MATLAB的图形用户界面（GUI）可实现可视化编程，不仅形象生动、互动友善、操作灵活，而且为人们提供了定性定量结合、局域全域结合、时域频域结合、模拟数字结合的数据探索、科学分析的仿真平台。

河南机电高等专科学校《自动控制原理》实验指导书专业：电气自动化技术、计算机控制技术生产过程自动化技术等吴君晓编2008年9月目录实验一 (2)实验二 (4)实验三 (6)实验四 (8)实验五 (10)实验六 (12)实验七 (14)实验八 (15)实验九 (17)实验一建立MATLAB环境下控制系统数学模型一. 实验目的1.熟悉MATLAB实验环境，掌握MATLAB命令窗口的基本操作。

2.掌握MATLAB建立控制系统数学模型的命令及模型相互转换的方法。

3.掌握使用MATLAB命令化简模型基本连接的方法。

二、实验设备和仪器1．计算机2. MATLAB软件三、实验原理控制系统常用的数学模型有四种：传递函数模型（tf对象）、零极点增益模型（zpk对象）、结构框图模型和状态空间模型（ss对象）。

经典控制理论中数学模型一般使用前三种模型，状态空间模型属于现代控制理论范畴。

1.传递函数模型（也称为多项式模型）连续系统的传递函数模型为：在MATLAB中用分子、分母多项式系数按s的降幂次序构成两个向量：num = [ b0 , b1 ,…, bm ] ，den = [ a0 , a1 ,…, an]。

用函数tf ( )来建立控制系统的传递函数模型，其命令调用格式为：G = tf ( num , den )注意：对于已知的多项式模型传递函数，其分子、分母多项式系数两个向量可分别用G.num{1}与G.den{1}命令求出。

2.零极点增益模型零极点模型是是分别对原传递函数的分子、分母进行因式分解，以获得系统的零点和极点的表示形式。

式中，K为系统增益，z1，z2，…，z m为系统零点，p1，p2，…，p n为系统极点。

在MATLAB 中，用向量z，p，k构成矢量组[ z, p, k ]表示系统。

即z = [ z1, z2 ,…,z m ]，p = [ p1, p2,…, p n ]，k = [ k ],用函数命令zpk ( )来建立系统的零极点增益模型，其函数调用格式为：G = zpk ( z, p, k )3.控制系统模型间的相互转换零极点模型转换为多项式模型: G＝zpk(G)多项式模型转化为零极点模型: G＝tf(G)4.系统反馈连接之后的等效传递函数两个环节反馈连接后，其等效传递函数可用feedback ( )函数求得。