稳定性理论

微分方程的稳定性理论简介一阶方程的平衡点及稳定性设有微分方程()()tfxx（1）右端方程不显含自变量t，称为自治方程。代数方程()0fx的实根0xx称为方程（1）的平衡点（或齐点）它也是方程（1）的解（齐解）。如果存在某个邻域，使方程（1）的解()xt从这个邻域内的某个(0)x出发，满足0lim()txtx（3）则称平衡点0x是稳定的（稳定性理论中称渐近稳定）；否则，称0x是不稳定的(不渐近稳定)判断平衡点0x是否稳定点通常有两种方法。利用定义即（3）式称间接法。不求方程（1）的解()xt，因而不利用（3）式的方法称直接法。下面介绍直接法。将()fx在0x点做Taylor展开，只取一次项，方程（1）近似为'00()xtfxxx（）（）（4）（4）称为（1）的近似方程，0x也是方程（4）的平衡点。关于0x点稳定性有如下结论：若'0fx（）0,则0x对于方程（4）和（1）都是稳定的；若'0fx（）0，则0x对于方程（4）和（1）都是不稳定的。0x对于方程（4）的稳定性很容易由定义（3）式证明，因为若记'0()fxa，则（4）的一般解是0()atxtcex其中c是由初始条件决定的常数，显然，当0a时（3）式成立。二阶方程的平衡点和稳定性二阶方程可用两个一阶方程表示为112212()(,)()(,)xtfxxxtgxx（6）右端不显含t，是自治方程。代数方程组1212(,)0(,)0fxxgxx（7）的实根011xx，022xx称为方程（6）的平衡点，记做00012(,)Pxx。如果存在某个邻域，使方程（6）的解1()xt，2()xt从这个邻域内的某个12((0),(0))xx出发，满足011lim()txtx，022lim()txtx（8）则称平衡点0P是稳定的（渐近稳定）；否则，称0P是不稳定的（不渐近稳定）。为了用直接法讨论方程的（6）的平衡点的稳定性，先看线性常系数方程1112221122()()xtaxaxxtbxbx（9）系数矩阵记做11aAb22ab（10）为研究方程（9）的唯一平衡点0P（0，0）的稳定性，假定A的行列式det0A(11)0P（0，0）的稳定性由（9）的特征方程det()0AI（12）的根（特征根）决定。方程（12）可以写成更加清晰的形式2120()detpqpabqA（13）特征根记作1，2，则1，221(4)2ppq（14）方程（9）的一般解具有形式121212()ttcece或121212()ttcecte，1,c，2,c为任意常数。按照稳定性的定义（8）式可知，当1，2为负数或有负实部时0P（0，0）是稳定平衡点；而当1，2有一个为正数或有正实部时0P（0，0）不是稳定平衡点。在条件（11）下1，2不可能为零。微分方程稳定性理论将平衡点分为结点、焦点、鞍点、中心等类型，完全由特征根1，2或相应的,pq取值决定。下表简明地给出了这些结果，表中最后一列指按照定义（8）式得到的关于稳定性的结论。1，2p，q平衡点类型稳定性12020,0,4pqpq稳定结点稳定12020,0,4pqpq不稳定结点不稳定1200q鞍点不稳定12020,0,4pqpq稳定退化结点稳定12020,0,4pqpq不稳定退化结点不稳定1,2i，020,0,4pqpq稳定焦点稳定1,2i，020,0,4pqpq不稳定焦点不稳定1,2i，00,0pq中心不稳定表1由特征方程决定的平衡点的类型和稳定性由表1可以看出，根据特征方程的的系数,pq的正负很容易判断平衡点的稳定性，准则如下：若0,0pq（15）则平衡点稳定；若0p或0q（16）则平衡点不稳定。以上是对线性方程（9）的平衡点0P（0，0）稳定性的结论，对于一般的非线性方程（6），可以用近似线性方法判断其平衡点00012(,)Pxx的稳定性。在0P点将12(,)fxx和12(,)gxx作Taylor展开，只取一次项，得（6）的近似线性方程1212000000112111222000000212111222()(,)()(,)()()(,)()(,)()xxxxxtfxxxxfxxxxxtgxxxxgxxxx（17）系数矩阵记作11xxfAg2002012(,)xxPxxfg（18）特征系数为110()xxPpfg，detqA（19）显然，0P点对于方程（17）的稳定性由表1或在准则（15），（16）决定，而且已经证明了如下结论：若方程（17）的特征根不为零或实部不为零，则0P点对于方程（6）的稳定性与对于近似方程（17）的稳定性相同，即由准则（15），（16）决定。最后，提出以下几点值得注意：1.平衡点及其稳定性的概念只是对自治方程（1），（6）而言才有意义。2.非线性方程（1），（6）的平衡点的稳定性，与相应的近似线性方程（4），（17）的平衡点的稳定性一致，是在非临界情况下（即0a或,0pq）得到的，在临界情况下（即0a或,0pq）二者可以不一致。3.在讨论平衡点稳定性时，对初始点的要求是存在一个邻域，这是局部稳定的定义。如果要求对任意的初始点，（3），（8）式成立，称为全局稳定。对于线性方程，局部稳定与全局稳定是等价的，对于非线性方程，二者不同。4.对于临界情况，和非线性方程的全局稳定，可以利用相轨线分析方法讨论。微分方程的定性理论对某些实际问题，建模的主要目的并不是要寻求动态过程每个瞬间的性态，而是研究某种意义下稳定状态的特征，特别是当时间充分长后动态过程的变化趋势。为分析这种稳定与不稳定的规律，常常不需要求解微分方程，而是利用微分方程的稳定性理论，直接研究平衡状态的稳定性即可。因此，常微分方程的定性和稳定性理论已成为数学建模必备的基础理论知识。动力学体系、自治系统与非自治系统1.基本概念考虑微分方程的初值问题(,)dXFtXdt(1)00()XtX（2）其中1,2,(...,)TnXxxx(,)FtX=11,21,1,(,...),(,...),...,(,...)Tnnnnftxxftxxftxx以下我们都假设(,)FtX是t、X的函数，且保证解的唯一性，即(,)FtX对X满足利普希茨（Lipschitz）条件：存在L，使|(,)(,)|||FtXFtXLXX，于是初值问题（1）、（2）存在唯一的解00(;,)XXttX,(3)设方程组（1）表示某一运动系统，其中自变量t视为时间，而X是在n维空间nR中质点运动时点的坐标1,2(,...,)nxxx。这时解（3）称为运动系统（1）在时刻t质点通过点0,0()tX的一个运动。在把时间t当做参数的这种解释下，称（1）是一个动力系统，称空间nR为相空间。参数方程（3）在相空间中确定的曲线称为相轨线，简称轨线。以下只考虑2n的情形，这时（1）、（2）变为(,,)(,,)dxPtxydtdyQtxydt（4）它满足初值条件00,00()()xtxyty的解为()()xxtyyt（5）方程组（4）是二维动力系统，Oxy平面就称为动力系统（4）的相平面。以t为参数，解（5）在相平面上所描绘的曲线就是相轨线或轨线。如果方程组(4)的右端函数显含自变量t，则称它为非自治系统，相应地，把右端函数不显含t的方程组(,)(,)dxPxydtdyQxydt（6）称为自治系统（或定常系统）。例1考虑自治系统dxydtdyxdt显然，cos,sinxtyt是方程组满足初值条件(0)1,(0)0xy的解。它在三维空间{(,,)}txy中表示的曲线是一条螺旋线。如果上述方程是描述质点在平面上运动的的动力系统，则以上的解在平面的轨线是一个圆221xy，它是上述曲线在Oxy平面上的投影，当t增加时，轨线的方向如图所示，它表明当时刻0t时经过点(1,0)的质点做逆时针方向的周期性运动。2.自治系统相轨线的基本性质假设自治系统（6）的右端函数在相平面2R满足存在唯一性定理条件，则它的轨线有以下基本性质：性质1设()xxt，()yyt是自治系统（6）的一个解，则(),()xxtcyytc也是(6)的解，其中c是任意常数。性质2自治系统（6）经过相平面上任意一点0,0()xy存在唯一的一条轨线。