写在前面

前面已经完成了控制系统的性能指标学习,从这节开始继续学习控制系统的分析方法,本文重点介绍分析方法概述和时域分析法。

一、控制系统的基本分析方法

  • 控制系统的基本分析方法包括:

    • 古典方法(经典控制理论):时域分析法、根轨迹法、频域分析法
    • 现代方法(现代控制理论):状态空间分析法
  • 利用上述方法分析系统的三大基本特性:

    • 能控能观性
    • 稳态性能
    • 动态性能

1.1 系统能控能观性

系统的动态性能与稳态性能前面已做介绍,这里介绍一下能控能观性。经典控制理论中并没有涉及这两个问题,因为经典控制理论讨论的是单入单出(SISO)系统输入输出的分析和综合问题,它的输入输出间的动态关系可以唯一的由传递函数来表示。

  • 能控性定义:
    对于一个 n n n阶系统 S S S,如果在有限的时间区间 t 0 ≤ t ≤ t a t_0\leq t ≤t_a t0tta 内,存在容许控制向量 u ( t ) u(t) u(t),能使系统从状态 x ( t 0 ) ≠ 0 x(t_0)≠0 x(t0)=0 转移到 x ( t a ) = 0 x(t_a)=0 x(ta)=0 ,则称状态 x ( t ) x(t) x(t) t 0 t_0 t0上能控。

  • 能观性定义:
    对于一个 n n n阶系统 S S S,如果对 t 0 t_0 t0 时刻,存在 t a t_a ta,即 t 0 < t a < ∞ t_0<t_a<\infty t0<ta<,根据 [ t 0 , t a ] [t_0,t_a] [t0,ta]上的 y ( t ) y(t) y(t)测量值能够唯一的确定系统在 t 0 t_0 t0 时刻的某初始状态 x 0 x_0 x0,则称 x 0 x_0 x0为系统在 [ t 0 , t a ] [t_0,t_a] [t0,ta]区间上的能观状态。

  • 系统能控能观性主要去解决两个问题:

    • ① 在有限时间内,控制作用能否使系统从初始状态转移到要求的状态?
    • ② 在有限时间内,能否通过系统输出的测量估计系统的初始状态?

    1.2 时域分析法

  • 时域分析法是根据系统的微分方程,以拉普拉斯变换作为数学工具,直接解出控制系统的时间响应。然后根据响应的表达式及其描述曲线来分析系统的控制性能,如稳定性、快速性、稳态精度等。

  • 时域法的特点:

    • ①直接在时间域中对系统进行分析校正,直观,准确。
    • ②可以提供系统时间响应的全部信息。
    • ③基于求解系统输出的解析解,比较烦琐。

  • 时域分析方法的基本假设

    • 系统的时间响应,不仅取决于系统本身的结构和参数,而且还与系统的初始状态以及加在系统上的外部作用信号有关。为了比较系统性能的优劣,对于外部作用信号和初始状态作典型化处理。

  • 时域法中部分动态指标的计算公式:

{ t r = π − β ω d t p = π 1 − ζ 2 ω n σ % = 3.5 ζ ω n \left\{ \begin{aligned} t_r&=\frac{\pi-\beta}{\omega_d}\\ t_p&=\frac{\pi}{\sqrt{1-\zeta^2}\omega_n}\\ \sigma\%&=\frac{3.5}{\zeta\omega_n} \end{aligned} \right. trtpσ%=ωdπβ=1ζ2 ωnπ=ζωn3.5

  • 时域法中系统稳定性的分析

    • 系统稳定的充要条件:
      • 系统所有闭环特征根均具有负的实部,或所有闭环特征根均位于左半s平面。

    • 劳斯(Routh)判据
      • 劳斯表第一列元素均大于零时系统稳定,否则系统不稳定且第一列元素符号改变的次数就是特征方程中正实部根的个数。

  • 时域法中系统稳态误差的计算:静态误差系数法

e s s = lim ⁡ s → 0 s Φ e ( s ) R ( s ) = lim ⁡ s → 0 s R ( s ) 1 1 + G 1 ( s ) H ( s ) = lim ⁡ s → 0 s R ( s ) 1 1 + K s v G 0 ( s ) e_{ss}=\underset{s\rightarrow0}{\lim}s\Phi_e(s)R(s)=\underset{s\rightarrow0}{\lim}sR(s)\frac{1}{1+G_1(s)H(s)}=\underset{s\rightarrow0}{\lim}sR(s)\frac{1}{1+\frac{K}{s^v}G_0(s)} ess=s0limsΦe(s)R(s)=s0limsR(s)1+G1(s)H(s)1=s0limsR(s)1+svKG0(s)1

r ( t ) = A ⋅ 1 ( t ) e s s p = lim ⁡ s → 0 s Φ e ( s ) R ( s ) = lim ⁡ s → 0 s ⋅ A s ⋅ 1 1 + G 1 ( s ) H ( s ) = A 1 + lim ⁡ s → 0 G 1 ( s ) H ( s ) = A 1 + K p r(t)=A\cdot1(t)e_{ssp}=\underset{s\rightarrow0}{\lim}s\Phi_e(s)R(s)=\underset{s\rightarrow0}{\lim}s\cdot\frac{A}{s}\cdot\frac{1}{1+G_1(s)H(s)}=\frac{A}{1+\underset{s\rightarrow0}{\lim}G_1(s)H(s)}=\frac{A}{1+K_p} r(t)=A1(t)essp=s0limsΦe(s)R(s)=s0limssA1+G1(s)H(s)1=1+s0limG1(s)H(s)A=1+KpA

r ( t ) = A ⋅ t e s s v = lim ⁡ s → 0 s Φ e ( s ) R ( s ) = lim ⁡ s → 0 s ⋅ A s 2 ⋅ 1 1 + G 1 ( s ) H ( s ) = A lim ⁡ s → 0 s G 1 ( s ) H ( s ) = A K v r(t)=A\cdot t e_{ssv}=\underset{s\rightarrow0}{\lim}s\Phi_e(s)R(s)=\underset{s\rightarrow0}{\lim}s\cdot\frac{A}{s^2}\cdot\frac{1}{1+G_1(s)H(s)}=\frac{A}{\underset{s\rightarrow0}{\lim}sG_1(s)H(s)}=\frac{A}{K_v} r(t)=Atessv=s0limsΦe(s)R(s)=s0limss2A1+G1(s)H(s)1=s0limsG1(s)H(s)A=KvA

r ( t ) = A 2 ⋅ t 2 e s s a = lim ⁡ s → 0 s Φ e ( s ) R ( s ) = lim ⁡ s → 0 s ⋅ A s 3 ⋅ 1 1 + G 1 ( s ) H ( s ) = A lim ⁡ s → 0 s 2 G 1 ( s ) H ( s ) = A K a r(t)=\frac{A}{2}\cdot t^2 e_{ssa}=\underset{s\rightarrow0}{\lim}s\Phi_e(s)R(s)=\underset{s\rightarrow0}{\lim}s\cdot\frac{A}{s^3}\cdot\frac{1}{1+G_1(s)H(s)}=\frac{A}{\underset{s\rightarrow0}{\lim}s^2G_1(s)H(s)}=\frac{A}{K_a} r(t)=2At2essa=s0limsΦe(s)R(s)=s0limss3A1+G1(s)H(s)1=s0lims2G1(s)H(s)A=KaA

其中,时域分析 K p K_p Kp是静态位置误差系数; K v K_v Kv是静态速度误差系数; K a K_a Ka是静态加速度误差系数。



人一旦受到责任感的驱使,就能创造出奇迹来。


喜欢我的分享,记得留下足迹!😐

12-09 11:54