第五章穩(wěn)定性定義

1、1841年年Liouville證明了證明了Riccati方程方程:2( )( )( )dyr x yp x yq xdx解的存在，但不能用公式求解，所以解的存在，但不能用公式求解，所以微分方程研究的主流發(fā)生了變化，不微分方程研究的主流發(fā)生了變化，不解方程去判斷解的形式，即定性理論解方程去判斷解的形式，即定性理論和穩(wěn)定性理論。盡管是很古老的學科，和穩(wěn)定性理論。盡管是很古老的學科，但 這 里 還 有 很 多 問 題 需 要 研 究但 這 里 還 有 很 多 問 題 需 要 研 究 . . 非線性微分方程非線性微分方程 實際問題中所研究的對象往往是非常復雜的，需要非線性微分方程(組)來描述，非線性方

2、程能求出解析解的很少，需要進行數(shù)值計算或理論分析。微分方程的研究內(nèi)容微分方程的研究內(nèi)容 求解：解析解、近似解、數(shù)值解求解：解析解、近似解、數(shù)值解 基本理論：解的存在惟一性、連續(xù)性基本理論：解的存在惟一性、連續(xù)性 定性穩(wěn)定性：時間趨于無窮時解的性態(tài)定性穩(wěn)定性：時間趨于無窮時解的性態(tài) 分支理論：解性態(tài)發(fā)生改變的一些參數(shù)值分支理論：解性態(tài)發(fā)生改變的一些參數(shù)值 本章介紹非線性微分方程的基本研究辦法，其出發(fā)本章介紹非線性微分方程的基本研究辦法，其出發(fā)點是在無法求出解析解的情況下通過方程本身的形點是在無法求出解析解的情況下通過方程本身的形式來分析時間趨于無窮時解的性態(tài)。式來分析時間趨于無窮時解的性態(tài)。 不

3、求解微分方程而通過方程右端函數(shù)的不求解微分方程而通過方程右端函數(shù)的信息探討時間趨于無窮時解的性態(tài)信息探討時間趨于無窮時解的性態(tài)22280(1)(1)(1sin ()(0)( ),lim ( )?tdxxdtdxxxdtdxxxxxtdtxxxx tx t 例滿足的解為00002222(0),1,(0),1 ( /1)0,0,0(),lim( )?(),lim( )?rtrtttdxrxxxxx edtdxxkrxxxxdtkk xexkxxkxdxyy xyx tdtdyxx xyy tdt 幾個例子：， 利用極坐利用極坐標將方程標將方程組組2222(1)(1)dxyx xydtdyxy xy

4、dt 化為化為2(1)1drrrdtddt ),;(),;(),;();(,) 1 ();(2121221121nnnnnyyytgyyytgyyytgytgyyyyytgdtdy其其中中 ),(),(),();(,)2()(2121221121nnnnnyyygyyygyyygytgyyyyygdtdy其其中中非自治系統(tǒng)或非定常系統(tǒng)自治系統(tǒng)或定常系統(tǒng)2121);();(yyLytgytg 0 h 00)();(ytyytgdtdyhtt |0);(max),min(),(ytgMMbahRyt .);(ytgdtdy );(ytgbyyattR 00,|:|（1）在）在n+1維空間的區(qū)域維空

5、間的區(qū)域 上連續(xù)；上連續(xù)；（2）在）在R上關于上關于 滿足李普希茲條件，即存在滿足李普希茲條件，即存在L 0,使對使對y),(2yt，有，有R上任意兩點上任意兩點, ),(1yt);(ytgdtdy )(ty 思思路路)(tyx );(dtdxtfx 0)0;(,)(;()(;()();();( tfttgtxtgdttdytgxtf顯顯然然有有 );(ytgdtdy )(ty );( xtfdtdx 0 x 1892()年，李雅普諾夫就如何判斷系統(tǒng)穩(wěn)定性的問題，歸納成兩種方法簡稱第一法和第二法。第一法是通過求解系統(tǒng)的微分方程，然后根據(jù)解的性質(zhì)來判斷系統(tǒng)的穩(wěn)定性，同時，他還指出非線性系統(tǒng)在工作

6、點附近的一定范圍內(nèi)可以用線性化了的微分方程來近似地加以描述。如果線性化的特征方程式的根全部是負實數(shù)根，或者是具有負實數(shù)部分的復根，則該系統(tǒng)在工作點附近周圍是穩(wěn)定的，否則便是不穩(wěn)定的。 t 李氏第二法（亦稱直接法）的特點是不必求解系統(tǒng)的微分方程就可以對系統(tǒng)的穩(wěn)定性進行分析和判斷。它是從能量的觀點出發(fā)得來的他指出：若系統(tǒng)有一個平衡點，則當時，系統(tǒng)運動到平衡點時，則系統(tǒng)積蓄的能量必達到一個極小值。由此，李雅普諾夫創(chuàng)造了一個輔助函數(shù)，可以用它來衡量系統(tǒng)積蓄的能量，但它并非是一個真正的能量函數(shù)。只要這一函數(shù)符合李雅普諾夫提出的穩(wěn)定性理論準則就能用來判斷系統(tǒng)的穩(wěn)定性。因此應用李氏 o 的的。李李雅雅普普諾

7、諾夫夫意意義義下下穩(wěn)穩(wěn)定定是是的的零零解解則則稱稱方方程程組組均均有有確確定定的的解解的的由由初初始始條條件件方方程程組組時時滿滿足足使使當當任任一一給給定定的的若若對對 0(*),(t )( )()(*), 0),(, 0000000 xttxtxxtxxxt (*)0)0 ;();(tfxtfdtdx1/221niixx穩(wěn)穩(wěn)定定的的。李李雅雅普普諾諾夫夫意意義義下下漸漸近近是是則則稱稱零零解解均均有有確確定定的的解解滿滿足足初初始始條條件件時時使使當當且且穩(wěn)穩(wěn)定定的的零零解解若若方方程程組組 00)(lim )()(, 0,0(*)t00000 xtxtxxtxxx o0漸進穩(wěn)定漸進穩(wěn)定=

8、 =穩(wěn)定穩(wěn)定+ +吸吸引引 o 稱稱為為不不穩(wěn)穩(wěn)定定的的。的的零零解解則則方方程程組組使使得得某某個個至至少少確確定定的的解解使使由由初初始始條條件件滿滿足足總總有有一一個個怎怎樣樣小小無無論論若若對對某某個個給給定定的的0(*) , )( ,),()(,01010000 xtxtttxxtxxx .0,., 0)(lim )()(,0(*)00t00000穩(wěn)穩(wěn)定定的的稱稱全全局局為為全全局局漸漸近近穩(wěn)穩(wěn)定定的的或或簡簡則則稱稱零零解解即即若若穩(wěn)穩(wěn)定定域域為為全全空空間間定定域域或或吸吸引引域域稱稱為為漸漸近近穩(wěn)穩(wěn)則則域域均均有有的的解解滿滿足足初初始始條條件件時時當當且且僅僅當當域域且且漸漸

9、近近穩(wěn)穩(wěn)定定的的零零解解若若方方程程組組 xDtxtxxtxDxDx o .0(*), 0,(*)圍圍內(nèi)內(nèi)漸漸近近穩(wěn)穩(wěn)定定的的稱稱為為在在大大范范的的零零解解那那么么方方程程組組于于都都收收斂斂時時當當?shù)牡拿棵恳灰粋€個解解如如果果方方程程組組 xxt0001( ,0,)(0)010100rtdxxrxkxyxx txxxyxxdtkkxyxexyyy 一般解的穩(wěn)定例的零解 是漸近，性的定義類似，只需要做一個平例的零解 是穩(wěn)定而不是漸移變穩(wěn)定的換近穩(wěn)定的例例 用用Maple命令畫出下邊捕食被捕食系統(tǒng)的命令畫出下邊捕食被捕食系統(tǒng)的方向場及一些軌線圖方向場及一些軌線圖20.08,0.001.dxxx

10、ydtdyyxydt (圖圖5.1)用用Maple命令畫出的圖形命令畫出的圖形從計算機的模擬看出系統(tǒng)有多個周期解。從計算機的模擬看出系統(tǒng)有多個周期解。取取t t的變化范圍的變化范圍 -100t100, 100t100, 選取下面選取下面6 6組初始值組初始值 x(0)=1,y(0)=0, x(0)=0,y(0)=4, x(0)=20,y(0)=25,x(0)=1,y(0)=0, x(0)=0,y(0)=4, x(0)=20,y(0)=25,x(0)=40,y(0)=25, x(0)=60,y(0)=25, x(0)=80,y(0)=25.x(0)=40,y(0)=25, x(0)=60,y(0)=25, x(0)=80,y(0)=25.輸入Malpe命令如下DEtoolsphaseportrait(diff(x(t),t)=2*x(t)-0.08*x(t)*y(t), diff(y(t),t)=-y(t)+0.01*x(t)*y(t), x(t),y(t), t=-100.100, x(0)=1,y(0)=0, x(0)=0,y(0)=4, x(0)=20,y(0)=25, x(0)=40,y(0)=25, x(0)=60,y(0)=25, x(0)=80,y