相平面分析法是研究二阶系统的一种图像方法。
他的主要思想是在相平面上绘制出与初始条件有关的运动轨迹,然后来分析这些轨迹的特征。
优点:
缺点:
系统的一族相平面轨迹叫做系统的相图。
有如下二阶系统,其中,以x1和x2作为坐标轴的二维平面叫做状态空间。
x ˙ 1 = f 1 ( x 1 , x 2 ) \dot{x}_{1}=f_{1}(x_{1},x_{2}) x˙1=f1(x1,x2)
x ˙ 2 = f 2 ( x 1 , x 2 ) \dot{x}_{2}=f_{2}(x_{1},x_{2}) x˙2=f2(x1,x2)
给定初始条件x(0),当时间t从0变化到无穷时,x(t)将在平面上绘制出一条曲线,叫做相平面轨迹。
一类重要的二阶系统微分方程形式是:
x ¨ + f ( x , x ˙ ) = 0 \ddot{x}+f(x,\dot{x})=0 x¨+f(x,x˙)=0
x 1 = x , x 2 = x ˙ x_{1}=x, x_{2}=\dot{x} x1=x,x2=x˙
利用上面这个变量代换,将这个二阶微分方程化成二元一次微分方程组。
x 1 ˙ = x 2 \dot{x_{1}}=x_{2} x1˙=x2
x 2 ˙ = − f ( x 1 , x 2 ) \dot{x_{2}}=-f(x_{1},x_{2}) x2˙=−f(x1,x2)
奇点就是相平面上的平衡点,而平衡点的定义是系统的状态能够永远保持稳定的点,也就是说,奇点是 x一阶导数为0的点。
f 1 ( x 1 , x 2 ) = 0 f_{1}(x_{1},x_{2})=0 f1(x1,x2)=0
f 2 ( x 1 , x 2 ) = 0 f_{2}(x_{1},x_{2})=0 f2(x1,x2)=0
求解以上两式即可得到平衡点。
非线性系统常常有多个孤立的奇点。
考察相轨迹的斜率:
d x 2 d x 1 = f 2 ( x 1 , x 2 ) f 1 ( x 1 , x 2 ) \frac{dx_{2}}{dx_{1}}=\frac{f_{2}(x_{1},x_{2})}{f_{1}(x_{1},x_{2})} dx1dx2=f1(x1,x2)f2(x1,x2)
由奇点处的定义可以知道,在奇点,这个式子变成了0/0形式的,也就是说,此时的斜率是不确定的。因此许多轨迹在这一点相交。
通过对称性可以简化系统的分析过程;
相图的对称性也可以通过分析相轨迹的斜率的对称性来判断。
一般微分方程的斜率形式如下:
d x 2 d x 1 = − f ( x 1 , x 2 ) x ˙ \frac{dx_{2}}{dx_{1}}=-\frac{f(x_{1},x_{2})}{\dot{x}} dx1dx2=−x˙f(x1,x2)
分析法中包含了系统微分方程的解析解。主要应用于一些特殊的非线性系统,比如分段线性系统。
在分析法中包含了两种求解方式。
求解系统的微分方程,分别解出x1和x2关于时间t的函数。然后消去时间t变量,从而构造出一个仅含x1和x2的方程:
g ( x 1 , x 2 , c ) = 0 g(x_{1},x_{2},c)=0 g(x1,x2,c)=0
其中常数c表示初始条件;
根据不同的初始条件绘制出一族相轨迹曲线,即可得到相图。
直接消去时间变量t。
比如一个质量弹簧系统:
x ¨ + x = 0 \ddot{x}+x=0 x¨+x=0
在方程中利用一个代换:
x ¨ = d x ˙ d x d x d t = x ˙ d x ˙ d x \ddot{x}=\frac{d\dot{x}}{dx} \frac{dx}{dt}=\dot{x}\frac{d\dot{x}}{dx} x¨=dxdx˙dtdx=x˙dxdx˙
带入可得:
x ˙ d x ˙ d x = − x \dot{x}\frac{d\dot{x}}{dx}=-x x˙dxdx˙=−x
这样,原方程化成了一个可分离变量的微分方程,因此可直接求解出x1与x2之间的函数关系。这里说成x1与x2,实际上也是x和x的导数:
x ˙ 2 + x 2 = x 0 2 \dot{x}^2+x^2=x_{0}^2 x˙2+x2=x02
在分析分段线性控制系统时,可分别绘制每段控制方程的相轨迹,然后简单的将每段轨迹同时绘制在一个相平面中即可。
对于无法进行分析求解的系统,可以利用等倾线法来绘制其相轨迹。
等倾线是由一组具有给定切线斜率的点的轨迹组成。
若相轨迹曲线上任意一点的斜率为 α \alpha α,则有下式成立:
d x 2 d x 1 = f 2 ( x 1 , x 2 ) f 1 ( x 1 , x 2 ) = α \frac{dx_{2}}{dx_{1}}=\frac{f_{2}(x_{1},x_{2})}{f_{1}(x_{1},x_{2})} = \alpha dx1dx2=f1(x1,x2)f2(x1,x2)=α
因此也可以得出具有斜率a的等倾线方程。
f 2 ( x 1 , x 2 ) = α f 1 ( x 1 , x 2 ) f_{2}(x_{1},x_{2})=\alpha f_{1}(x_{1},x_{2}) f2(x1,x2)=αf1(x1,x2)
下面给出一个简单的例子:
考虑质量弹簧系统,系统方程前文已经给出,其相轨迹上任一点处的斜率为:
d x 2 d x 1 = − x 1 x 2 = α \frac{dx_{2}}{dx_{1}}=-\frac{x_{1}}{x_{2}} = \alpha dx1dx2=−x2x1=α
因此等倾线方程为:
x 1 + α x 2 = 0 x_{1}+\alpha x_{2}=0 x1+αx2=0
也就是说,这个等倾线是一条直线。在这个直线上面绘制一系列的斜率为a的小短线,并且当a取不同值的时候,在同样绘制,最后得出相轨迹的切线方向场。绘制出系统的相轨迹如下图所示。
利用近似关系可以得到下式:
Δ x ≈ x ˙ Δ t \Delta x \approx \dot{x}\Delta t Δx≈x˙Δt
因此利用积分法可以得到计算时间的公式:
t − t 0 = ∫ x 0 x ( 1 / x ˙ ) d x t-t_{0} = \int_{x_{0}}^{x}(1/ \dot{x})dx t−t0=∫x0x(1/x˙)dx
一般形式的二阶线性系统总可以转化成下面的微分方程的形式:
x ¨ + a x ˙ + b x = 0 \ddot{x}+a \dot{x}+bx = 0 x¨+ax˙+bx=0
为了获取线性系统的相轨迹,需要求解上面的微分方程,得到x(t)的时间方程:
x ( t ) = k 1 e λ 1 t + k 2 e λ 2 t x(t)=k_{1}e^{\lambda_{1}t}+k_{2}e^{\lambda_{2}t} x(t)=k1eλ1t+k2eλ2t
其中 λ \lambda λ由微分方程的特征方程来确定。并且根据 λ 1 \lambda_{1} λ1和 λ 2 \lambda_{2} λ2的不同,奇点又分为以下几种情况:
因此可以看出,线性系统的稳定性特征是由其奇点的特性唯一确定的,但是非线性系统并不是这样。
非线性系统的局部特性可以通过一个线性系统来近似逼近。
对于非线性系统,如果奇点不在原点,那么可以通过把原点状态和奇点状态之间的差异定义为一个新的状态变量,就可以把奇点转化到原点。因此,可以简单的认为非线性系统有在原点的奇点。
由上文,非线性系统:
x ˙ 1 = f 1 ( x 1 , x 2 ) \dot{x}_{1}=f_{1}(x_{1},x_{2}) x˙1=f1(x1,x2)
x ˙ 2 = f 2 ( x 1 , x 2 ) \dot{x}_{2}=f_{2}(x_{1},x_{2}) x˙2=f2(x1,x2)
利用Taylor展开,可以得到:
x 1 ˙ = a x 1 + b x 2 + g 1 ( x 1 , x 2 ) \dot{x_{1}}=ax_{1}+bx_{2}+g_{1}(x_{1},x_{2}) x1˙=ax1+bx2+g1(x1,x2)
x 2 ˙ = c x 1 + d x 2 + g 2 ( x 1 , x 2 ) \dot{x_{2}}=cx_{1}+dx_{2}+g_{2}(x_{1},x_{2}) x2˙=cx1+dx2+g2(x1,x2)
其中个g1和g2是高阶项,如果忽略他们,那么该系统就变成了一个线性系统。
极限环是非线性系统特有的性质,它具有两个条件:封闭的,独立的。
根据相轨迹的运动方式,可以分为三种类型的极限环:
有以下三个定理可以用来判断系统极限环的存在性。
在介绍下面的定理之前,做如下规定:
N:表示极限环包围的所有节点、中心点和焦点数量。
S:表示极限环包围的所有鞍点数量。
1. Poincare定理
如果一个二阶自制系统的极限环存在,则有 N=S+1.
推论:极限环一定包含至少一个平衡点。
2. Poincare-Bendixson定理
如果一个二阶自制系统的相轨迹总是维持在一个有限的区域内部,那么下面说法中必定有一个是正确的:
(a)轨迹趋向于一个平衡点
(b)轨迹趋向于一个对称的稳定极限环
(c)轨迹自身是一个极限环
3. Bendixson定理
∂ f 1 ∂ x 1 + ∂ f 2 ∂ x 2 \frac{\partial f_{1}}{\partial x_{1}} + \frac{\partial f_{2}}{\partial x_{2}} ∂x1∂f1+∂x2∂f2
对于非线性系统,在相平面的一个区域内,如果上式不为0也不变号,那么在这个区域内不存在极限环。