控制工程(MATLAB)实验

《机电控制工程基础》

实验指导书

适用专业：机械设计制造及其自动化

机械电子工程

太原工业学院机械工程系

实验一 系统时间响应分析

实验课时数：2学时 实验性质：设计性实验 实验室名称：数字化实验室（机械工程系）

一、实验项目设计内容及要求

1. 实验目的

本实验的内容牵涉到教材的第3、4、5章的内容。本实验的主要目的是通过试验，能够使学生进一步理解和掌握系统时间响应分析的相关知识，同时也了解频率响应的特点及系统稳定性的充要条件。

2. 实验内容

完成一阶、二阶和三阶系统在单位脉冲和单位阶跃输入信号以及正弦信号作用下的响应，求取二阶系统的性能指标，记录试验结果并对此进行分析。

3. 实验要求

系统时间响应分析试验要求学生用MATLAB 软件的相应功能，编程实现一阶、二阶和三阶系统在几种典型输入信号（包括单位脉冲信号、单位阶跃信号、单位斜坡信号和正弦信号）作用下的响应，记录结果并进行分析处理：对一阶和二阶系统，要求用试验结果来分析系统特征参数对系统时间响应的影响；对二阶系统和三阶系统的相同输入信号对应的响应进行比较，得出结论。

4. 实验条件

利用数字化实验室的计算机，根据MATLAB 软件的功能进行简单的编程来进行试验。

二、具体要求及实验过程

1.系统的传递函数及其MATLAB 表达 (1)一阶系统 传递函数为：1

)(+=

Ts K

s G 传递函数的MATLAB 表达： num=[k]；den=[T,1]；G(s)=tf(num,den) (2)二阶系统

传递函数为：2

22

2)(n

n n

w s w s w s G ++=ξ

传递函数的MATLAB 表达： num=[2n w ]；den=[1，n w ξ2，2

n w ]；G(s)=tf(num,den)

（3）任意的高阶系统

传递函数为：n

n n n

m

m m m a s a s a s a b s b s b s b s G ++++++++=----11

101110)( 传递函数的MATLAB 表达：

num=[m m b b b b ,,,110- ]；den=[n

n a a a a ,,,110- ]；G(s)=tf(num,den) 若传递函数表示为：)())(()())(()(1010n

m p s p s p s z s z s z s K s G ------=

则传递函数的MATLAB 表达：

z=[m z z z ,,,10 ]；p=[n p p p ,,,10 ]；K=[K]；G(s)=zpk(z,p,k) 2.各种时间输入信号响应的表达

（1）单位脉冲信号响应：[y,x]=impulse[sys,t] （2）单位阶跃信号响应：[y,x]=step[sys,t] （3）任意输入信号响应：[y,x]=lsim[sys,u,t]

其中，y 为输出响应，x 为状态响应（可选）；sys 为建立的模型；t 为仿真时间区段（可选）

实验方案设计可参考教材相关内容，相应的M 程序可参考教材（杨叔子主编的《机械工程控制基础》第五版）提供的程序，在试验指导教师的辅导下掌握M 程序的内容和格式要求，并了解M 程序在MATLAB 软件中的加载和执行过程。

3.实验的具体内容

（1）完成一阶（选用不同的时间常数T ）、二阶系统（选择不同的阻尼比ξ和无阻尼固有频率n w ，而且阻尼比ξ要有欠阻、临界阻尼和过阻尼三种情况）在典型输入信号（单位脉冲、单位阶跃、正弦信号）作用下所对应的时间响应实验；

（2）完成二节系统性能指标的求取（设计的二阶系统必须是欠阻尼的二阶系统）。

（2）完成一稳定三阶系统的单位阶跃响应的实验；

（3）完成一个稳定的三阶系统和一个不稳定的三阶系统的单位脉冲响应，比较两响应曲线的差别并说明原因。

4.实验分析内容

（1）分析时间常数对一阶系统时间响应的影响；

（2）分析参数对二阶系统的时间响应的性能指标的影响；

（3）分析系统稳定性与系统特征值的关系；

（4）了解系统频率响应的特点。

三、实验报告

1、编写程序：

3、结果分析

科目：控制工程

姓名：颜人帅

班级学号:102012237

1:完成一阶系统，二阶系统在典型输入信号作用下所对应的时间响应实验 （1）一阶系统

传递函数为：1

)(+=

Ts K

s G ；其中K=10;T1=10,T2=20,T3=30; 1-1：一阶单位阶跃响应： 1、编写程序：

2、>> num1=[10];den1=[10 1];step(num1,den1) >> hold on

>> num2=num1;den2=[20 1];step(num2,den2) >> num3=num1;den3=[30 1];step(num3,den3) >> gtext('一阶系统单位阶跃响应') >> gtext(['图 1-1'])

gtext(['T1=10';'T2=20';'T3=30']) 运行结果：

结论：由图1-1 可知，随着时间常数T 的增大，一阶系统响应时间增大； 1-2：一阶单位脉冲响应：

num1=[10];den1=[10 1];impulse(num1,den1) >> hold on

>> num2=num1;den2=[20 1];impulse(num2,den2) >> num3=num1;den3=[30 1];impulse(num3,den3) >> gtext('一阶系统单位脉冲响应') >> gtext(['图 1-2'])