自动控制原理结课论文论文题目：时域分析的Matlab实现时域分析的Matlab实现摘要分析和设计系统的首要工作是确定系统的数学模型。

一旦建立了合理的、便于分析的数学模型，就可以对已组成的控制系统进行分析，从而得出系统性能的改进方法。

经典控制理论中，常用时域分析法、根轨迹法或频率分析法来分析控制系统的性能。

本文采用MATLAB 语言编程实现了高阶系统时域分析，分析了其稳定性、快速性、准确性，并应用实例验证了其有效性。

[关键词] 时域分析高阶系统MATLAB 实现目录一、引言 (1)二、时域分析基础理论 (1)（一）典型输入信号和时域性能指标 (2)1、典型输入信号 (2)2、时域性能指标 (4)（二）一阶系统的时域分析 (5)1、单位阶跃响应 (5)2、单位斜坡响应 (7)3、单位脉冲响应 (7)（三）高阶系统的时域分析 (8)三、基于MATLAB实现高阶系统的时域分析 (10)四、高阶系统时域分析的MATLAB 实现 (11)（一）应用经典法求解 (12)（二）MATLAB实现 (12)1、系统稳定性分析 (13)2、系统的快速性分析 (16)3、系统的准确性分析 (16)（三）应用MATLAB分析系统的动态特性 (17)五、结论 (19)参考文献 (20)时域分析的Matlab实现一、引言信号与系统的分析在自动控制领域有十分重要的作用。

进行分析时，一般先抽象为数学模型，然后讨论系统本身的初始状态以及不同激励时的响应。

对于高阶的微分方程，由于计算量庞大，人工计算难于实现。

经典控制理论对高阶系统进行时域分析通常采用拉氏反变换的方法求系统响应，系统阶次越高，进行拉氏反变换的困难就越大，因此，用经典法对高阶系统进行时域分析是一件较困难的事。

本文采用MATLAB 语言编程，设计了对高阶系统进行时域性能辅助分析程序，充分发挥了MATLAB 人机交互性好、函数调用方便、数学运算与绘图功能强大的优势，使分析效率和准确性大为提高。

二、时域分析基础理论那什么是时域分析呢？时域分析是指控制系统在一定的输入下，根据输出量的时域表达式，分析系统的稳定性、瞬态和稳态性能。

由于时域分析是直接在时间域中对系统进行分析的方法，所以时域分析具有直观和准确的优点。

系统输出量的时域表示可由微分方程得到，也可由传递函数得到。

在初值为零时，一般都利用传递函数进行研究，用传递函数间接的评价系统的性能指标。

具体是根据闭环系统传递函数的极点和零点来分析系统的性能。

此时也称为复频域分析。

（一）典型输入信号和时域性能指标1、典型输入信号控制系统的输出响应是系统数学模型的解。

系统的输出响应不仅取决于系统本身的结构参数、初始状态，而且和输入信号的形式有关。

初始状态可以作统一规定，如规定为零初始状态。

如再将输入信号规定为统一的形式，则系统响应由系统本身的结构、参数来确定，因而更便于对各种系统进行比较和研究。

自动控制系统常用的典型输入信号有下面几种形式：1.阶跃函数 定义为⎩⎨⎧<≥=000u(t) t t U式中U 是常数，称为阶跃函数的阶跃值。

U=1的阶跃函数称为单位阶跃函数，记为1(t)。

如图2-1所示。

单位阶跃函数的拉氏变换为1/s 。

在t=0处的阶跃信号，相当于一个不变的信号突然加到系统上，如指令的突然转换、电源的突然接通、负荷的突变等，都可视为阶跃作用。

2.斜坡函数 定义为⎩⎨⎧<≥=000u(t)t t Ut这种函数相当于随动系统中加入一个按恒速变化的位置信号，恒速度为U 。

当U=1时，称为单位斜坡函数，如图2-2所示。

单位斜坡函数的拉氏变换为 1/s 2。

3.抛物线函数 定义为⎪⎩⎪⎨⎧<≥=00021u(t)2t t Ut这种函数相当于系统中加入一个按加速度变化的位置信号，加速度为U 。

当U=1时，称为单位抛物线函数，如图2-3所示。

单位抛物线函数的拉氏变换为1/s 3。

4.单位脉冲函数δ(t) 定义为⎪⎩⎪⎨⎧=⎩⎨⎧≠=∞==⎰∞∞-0)(000)(u(t)dt t t t t δδ图2-1 单位阶跃函数 图2-2 斜坡函数图2-3 抛物线函数 图2-4 单位脉冲函数 单位脉冲函数的积分面积是1。

单位脉冲函数如图2-4所示。

其拉氏变换为1。

单位脉冲函数在现实中是不存在的，它只有数学上的意义。

在系统分析中，它是一个重要的数学工具。

此外，在实际中有很多信号与脉冲信号相似，如脉冲电压信号、冲击力、阵风等。

5.正弦函数 定义为t A t u ωsin )(=其中A 为振幅，ω为角频率。

其拉氏变换为22ωω+s A 。

ttt t用正弦函数作输入信号，可以求得系统对不同频率的正弦输入函数的稳态响应，由此可以间接判断系统的性能。

2、时域性能指标时域中评价系统的暂态性能，通常以系统对单位阶跃输入信号的暂态响应为依据。

这时系统的暂态响应曲线称为单位阶跃响应或单位过渡特性，典型的响应曲线如图2-5所示。

为了评价系统的暂态性能，规定如下指标：图2-5 单位阶跃输入信号下的暂态响应1.延迟时间t d指输出响应第一次达到稳态值50%所需的时间。

2.上升时间t r指输出响应从稳态值的10%上升到90%所需的时间。

对有振荡的系统，则取响应从零到第一次达到稳态值所需的时间。

3.峰值时间t p指输出响应超过稳态值而达到第一个峰值（即y(t p)）所需的时间。

4.调节时间t s指当输出量y(t)和稳态值y(∞)之间的偏差达到允许范围（一般取2%或5%）以后不再超过此值所需的最短时间。

5.最大超调量（或称超调量）σp% 指暂态过程中输出响应的最大值超过稳态值的百分数。

即%100)()]y(-)[y(t %p p ⨯∞∞=y σ6.稳态误差e ss 指系统输出实际值与希望值之差。

在上述几项指标中，峰值时间t p 、、上升时间t r 和延迟时间t d 均表征系统响应初始阶段的快慢；调节时间t s 表征系统过渡过程（暂态过程）的持续时间，从总体上反映了系统的快速性；而超调量σp %标志暂态过程的稳定性；稳态误差反映系统复现输入信号的最终精度。

（二）一阶系统的时域分析凡是可用一阶微分方程描述的系统称一阶系统。

一阶系统的传递函数为1Ts 1G(s)+=式中T 称为时间常数，它是表征系统惯性的一个重要参数。

所以一阶系统是一个非周R(s)期的惯性环节。

图2- 6为一阶系统的结构图。

图2-6 一阶系统的结构图下面分析在三种不同的典型输入信号作用下一阶系统的时域分析。

1、单位阶跃响应当输入信号u(t)=1(t)时，U(s)=1/s ，系统输出量的拉氏变换为11)1s(Ts 1Y (s) +-=+=TS T s对上式取拉氏反变换，得单位阶跃响应为)0(1Y (t)≥-=-t e T t由此可见，一阶系统的阶跃响应是一条初始值为0，按指数规律上升到稳态值1的曲线，见图2-7。

由系统的输出响应可得到如下的性能。

图2-7 一阶系统的阶跃响应曲线1.由于Y(t)的终值为1，因此系统稳态误差为0。

2.当t=T 时，Y(T)=0.632。

这表明当系统的单位阶跃响应达到稳态值的63.2%时的时间，就是该系统的时间常数T 。

单位阶跃响应曲线的初始斜率为T e T t Tt t 11dt dy(t)00===-=这表明一阶系统的单位阶跃响应如果以初始速度上升到稳态值1，所需的时间恰好等于T 。

3.根据暂态性能指标的定义可以求得调节时间为 t s =3T(s) (±5%的误差带)t s =4T(s) (±2%的误差带)延迟时间为 t d =0.69T(s)上升时间为 t r =2.20T(s)峰值时间和超调量都为0。

2、单位斜坡响应当输入信号u(t)=t 时，U(s)=1/s 2，系统输出量的拉氏变换为)0(11)1(Ts s 1Y(s)222≥++-=+=t TS T s T s对上式取拉氏反变换，得单位斜坡响应为0)(Y (t)≥+-=-t Te T t T t其中（t-T ）为稳态分量，Te -t/T 为暂态分量。

单位斜坡响应曲线如2-8。

图2-8 单位斜坡响应曲线由一阶系统单位斜坡响应可分析出，系统存在稳态误差。

因为u(t)=t ，输出稳态为t –T,所以稳态误差为e ss =t – (t – T)=T 。

从提高斜坡响应的精度来看，要求一阶系统的时间常数T 要小。

3、单位脉冲响应当u(t)=δ(t)时，系统的输出响应为该系统的脉冲响应。

因为L[δ(t)]=1，一阶系统的脉冲响应的拉氏变换为T s T s G /1/1)(Y(s) +== 对应单位脉冲响应为)0(1Y (t) ≥=-t e T Tt单位脉冲响应曲线如图2-9。

时间常数T 越小，系统响应速度越快。

图2-9 脉冲响应曲线 （三）高阶系统的时域分析设高阶系统的传递函数可表示为：)(...... )(11101110m n a s a s a s a b s b s b s b s n n n n m m m m ≥++++++=----φ设闭环传递函数的零点为-z 1,-z 2,…,-z m ，极点为-p 1,-p 2,…, -p n ，则闭环传递函数可表示为：)())...()(())...()(( )(2121m n p s p s p s z s z s z s K s Y n m ≥++++++=当输入信号为单位阶跃信号时，输出信号为：∏∏∏===++++=r k nk nk k q j j m i i s s p s s z s K s Y 12211)2()()()(ωωζ式中n =q +2r ，而q 为闭环实极点的个数，r 为闭环共轭复数极点的对数。

用部分分式展开得∑∑==++-+++++=r 1k 222q 1j j 021)(p s s Y (s)nk nk k k nk k nk k K j s s C s B A A ωωζζωωζ对上式取反拉氏变换得：)0(1sin 1cos A y(t)12r 1k 2q 1j 0≥-+-++=∑∑∑=-=-=-t te C t e B eA r k k nk t K k nk t K t p j nk k nk k j ζωζωωζωζ由上式分析可知，高阶系统的暂态响应是一阶惯性环节和二阶振荡响应分量的合成。

系统的响应不仅和k ζ、ωnk 有关，还和闭环零点及系数Aj 、B k 、C k 的大小有关。

这些系数的大小和闭环系统的所有的极点和零点有关，所以单位阶跃响应取决于高阶系统闭环零极点的分布情况。

从分析高阶系统单位阶跃响应表达式可以得到如下结论：1. 高阶系统暂态响应各分量衰减的快慢由j p -和k ζ、ωnk 决定，即由闭环极点在s 平面左半边离虚轴的距离决定。

闭环极点离虚轴越远，相应的指数分量衰减的越快，对系统暂态分量的影响越小；反之，闭环极点离虚轴越近，相应的指数分量衰减的越慢，系统暂态分量的影响越大。