版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、<p> 畢 業(yè) 論 文(設(shè)計(jì))</p><p> 論文題目 常微分方程及其matlab求解 </p><p><b> 開題報(bào)告</b></p><p><b> 目 錄</b></p><p><b> 摘要1<
2、/b></p><p><b> 關(guān)鍵字1</b></p><p><b> 引言1</b></p><p> 第一章 一階微分方程的初等解法1</p><p> 1.1 變量分離微分方程與變量代換1</p><p> 1.1.1變量分離微分方程
3、2</p><p> 1.1.2 可化為變量分離微分方程的類型2</p><p> 1.2線性分式方程3</p><p> 1.3線性方程與伯努利方程4</p><p> 1.4全微分方程與積分因子6</p><p> 1.4.1 全微分方程6</p><p> 1.4.
4、2 積分因子7</p><p> 1.5一階隱方程8</p><p> 1.5.1可解出y或x的方程的解法8</p><p> 1.5.2 不顯含y(或x)的方程10</p><p> 1.6 matlab 在一階方程中的應(yīng)用10</p><p> 第二章 高階微分方程的解法13</p&
5、gt;<p> 2.1 線性微分方程的一般理論13</p><p> 1 齊次線性微分方程的解的性質(zhì)與結(jié)構(gòu)13</p><p> 2.非奇次線性微分方程與常數(shù)變易法14</p><p> 2.2常系數(shù)線性微分方程的解法14</p><p> 1.復(fù)值函數(shù)與復(fù)值解15</p><p>
6、 2.3.非齊次線性微分方程的比較系數(shù)法18</p><p> 2.4用matlab解高階線性方程20</p><p> 第三章 線性微分方程組24</p><p> 3.1矩陣指數(shù)24</p><p> 3.2基解矩陣的計(jì)算25</p><p> 1.矩陣A存在n個(gè)線性無(wú)關(guān)的特征向量的情形2
7、5</p><p> 3.3用MATLAB解線性方程組29</p><p><b> 致 謝30</b></p><p><b> 參考文獻(xiàn)30</b></p><p> Ordinary Differential Equation And Carry On By MATLAB30&
8、lt;/p><p> Abstract30</p><p> Keywords30</p><p> 常微分方程及其matlab求解</p><p><b> 摘要:</b></p><p> 本文主要討論了一階常微分方程和高階常微分方程的相關(guān)解法問(wèn)題,并用matlab求解相關(guān)方程.文章
9、首先探討了一階常微分方程的解法,討論的主要類型有:變量可分離方程、可化為變量可分離方程的類型、齊次方程、一階線性微分方程;在解決這些類型的一階常微分方程時(shí),用到的方法有:變量分離法和一階線性方程的常數(shù)變易法.然后討論了高階常微分方程的解法的問(wèn)題,所討論的解法有:非齊線性方程的常數(shù)變易法、常系數(shù)齊線性方程的歐拉待定指數(shù)法和非齊線性方程的比較系數(shù)法.最后簡(jiǎn)單討論了線性微分方程組的解法和性質(zhì)。</p><p> 關(guān)鍵
10、字: 一階常微分方程 高階常微分方程 解法 matlab</p><p><b> 引言:</b></p><p> 常微分方程作為現(xiàn)代數(shù)學(xué)的一個(gè)重要分支,它的產(chǎn)生幾乎與微積分是同時(shí)代的,經(jīng)過(guò)歷史的演變,它已經(jīng)是各種應(yīng)用學(xué)科和數(shù)學(xué)理論研究都不可缺少的工具。隨著計(jì)算機(jī)技術(shù)的飛速發(fā)展,更使它的應(yīng)用滲透到力學(xué)、天文、物理等各個(gè)領(lǐng)域。</p><p&g
11、t; 遺憾的是,絕大多數(shù)微分方程定解問(wèn)題的解不能以實(shí)用的解析形式來(lái)表示,這就產(chǎn)生了理論與應(yīng)用的矛盾:一方面,人們建立了大量實(shí)用的數(shù)學(xué)模型,列出了反映客觀現(xiàn)象的微分方程;另一方面,人們又無(wú)法得到這些方程的準(zhǔn)確解以定量地描述客觀過(guò)程。從20世紀(jì)80年代以來(lái),世界各國(guó)所開發(fā)的數(shù)學(xué)類科技軟件多達(dá)幾十種,在我國(guó)流行的數(shù)學(xué)軟件主要有四種:MATLAB、Mathematica、Maple和MathCAD。其中,MATLAB有著其它幾種數(shù)學(xué)軟件無(wú)法比
12、擬的優(yōu)勢(shì)和適用面,近幾年,MATLAB已成為科學(xué)工作者首選的數(shù)學(xué)軟件。本文把微分方程求解和MATLAB有機(jī)的結(jié)合起來(lái),全面介紹了微分方程的求解在MATLAB中的實(shí)現(xiàn),使得讓數(shù)學(xué)基礎(chǔ)不深厚的讀者同樣能輕易利用MATLAB解決較高深的微分方程問(wèn)題。</p><p> 第一章 一階微分方程的初等解法</p><p> 1.1 變量分離微分方程與變量代換</p><p&
13、gt; 變量分離微分方程是一種最簡(jiǎn)單也是最基本的可用初等方法求解的微分方程類型,對(duì)一般的微分方程總是設(shè)法尋求適當(dāng)?shù)淖兞看鷵Q,將其化為變量分離微分方程來(lái)求解。</p><p> 1.1.1變量分離微分方程</p><p><b> 形如</b></p><p><b> (1.1)</b></p>&l
14、t;p> 的方程,稱為變量分離微分方程,其中分別是的連續(xù)函數(shù)。</p><p><b> 如果,方程改寫為:</b></p><p><b> (1.2)</b></p><p> 如果是方程的解,且,則在解的定義域內(nèi)滿足</p><p><b> (1.3)</b&
15、gt;</p><p><b> 兩邊積分得</b></p><p><b> (1.4)</b></p><p><b> 則是由隱函數(shù)方程</b></p><p><b> ?。?.5)</b></p><p> 所確定
16、的函數(shù),對(duì)一切允許值C,都將是方程的解,即(1.5)實(shí)際給出了方程的通解公式。</p><p> 如果有實(shí)根(k=1,2,3,…,m),則可以驗(yàn)證(k=1,2,3,…,m)也能是方程的解。有時(shí)可以通過(guò)擴(kuò)大常數(shù)C的取值范圍,使其包含于通解表達(dá)式中。</p><p> 1.1.2 可化為變量分離微分方程的類型</p><p><b> 對(duì)于方程</
17、b></p><p><b> 在變量代換下變成</b></p><p><b> ?。?.6)</b></p><p><b> 整理得:</b></p><p><b> (1.7)</b></p><p> 不難
18、看出,得到的方程為變量分離微分方程當(dāng)且僅當(dāng):</p><p><b> ?。?.8)</b></p><p> 由此我們可得出以下結(jié)論:</p><p> 變量代換能化為變量分離微分方程的一般類型為:</p><p><b> ?。?.9)</b></p><p> 其
19、中均為任意的連續(xù)函數(shù),在變量代換u下,方程化為變量分離方程</p><p><b> ?。?.10)</b></p><p> 求出(1.10)的解,用代回原變量,即求的原方程解.</p><p><b> 1.2線性分式方程</b></p><p><b> 我們把形如</b
20、></p><p><b> ?。?.11)</b></p><p> 的方程稱為線性分式方程,這里均為常數(shù)。</p><p> 當(dāng)時(shí),方程(1-11)是齊次方程,當(dāng)不全為零時(shí),如何化為某種已知的可解類型?</p><p><b> 當(dāng)?shù)那樾巍?lt;/b></p><p&
21、gt;<b> 設(shè),則方程可以寫成</b></p><p><b> (1.12)</b></p><p> 令,方程就化為變量分離微分方程</p><p><b> ?。?.13)</b></p><p><b> 當(dāng)時(shí)</b></p>
22、;<p><b> 解方程組</b></p><p><b> (1.14)</b></p><p><b> 設(shè)所求得的解為 。</b></p><p><b> 作坐標(biāo)平移變換</b></p><p><b> (1.
23、15)</b></p><p> 將方程(1.11)變成齊次方程</p><p><b> ?。?.16)</b></p><p> 經(jīng)過(guò)變換將方程(1-16)化為變量微分方程。求解所得變量微分方程后,逐步代回原來(lái)的變量,求得原方程的解。</p><p> 1.3線性方程與伯努利方程</p>
24、<p> 我們把一階線性方程通常寫成其標(biāo)準(zhǔn)形式:</p><p><b> (1.17)</b></p><p> 其中,為連續(xù)函數(shù),當(dāng)時(shí),方程成為:</p><p><b> (1.18)</b></p><p> 稱方程(1.18)為方程(1.17)對(duì)應(yīng)的齊次線性方程,而
25、稱(1.17)為非齊次線性方程。</p><p> 齊次線性方程(1.18)是變量分離微分方程,可求其通解為:</p><p><b> (1.19)</b></p><p> 為了求(1.17)的解,設(shè)想用兩個(gè)新的未知函數(shù) 的乘積表示原來(lái)的未知函數(shù),即:</p><p><b> ?。?.20)<
26、/b></p><p><b> 代入方程得</b></p><p><b> ?。?.21)</b></p><p><b> 將其整理得:</b></p><p><b> (1.22)</b></p><p>
27、設(shè)為齊次方程的解,則方程變成</p><p><b> (1.23)</b></p><p> 這是一個(gè)變量分離微分方程,很容易求得其通解為;</p><p><b> (1.24)</b></p><p> 最后得到非齊次線性方程的通解</p><p><b&
28、gt; ?。?.25)</b></p><p> 以上這種方法被稱為常數(shù)變易法。</p><p><b> 形如 </b></p><p><b> (1.26)</b></p><p> 的方程稱為伯努利方程,將方程(1.26)兩邊同時(shí)除以,方程變成</p>&l
29、t;p><b> 由于</b></p><p> 即關(guān)于y的導(dǎo)數(shù)恰好為的1-n(常數(shù))倍,于是</p><p><b> 方程化為</b></p><p><b> 令</b></p><p> 伯努利方程化為了線性方程</p><p>
30、 求得此線性方程的通解,代回原變量,就可得到伯努利方程的通解。此外,如果n>0,則y=0顯然也是伯努利方程的解。</p><p> 1.4全微分方程與積分因子</p><p> 全微分方程是另一種既簡(jiǎn)單又基本的可積方程類型,一般的方程可以通過(guò)求得其積分因子,乘以積分因子而化成全微分方程求解,這是一階微分方程初等積分方法的第二種有效途徑。</p><p>
31、 1.4.1 全微分方程</p><p> 我們討論一階對(duì)稱形方程</p><p> M(x,y)dx+N(x,y)dy=0. (1.27) </p><p> 而且約定M(x,y),N(x,y)在所考慮的單連
32、通區(qū)域G都具有一階連續(xù)偏導(dǎo)數(shù)。如果其左端恰好是一個(gè)全微分,即存在一個(gè)可微的二元函數(shù)u(x,y),使得</p><p> M(x,y)dx+N(x,y)dy=du(x,y), (1.28)</p><p> 則稱這個(gè)方程為全微分方程,u(x,y)稱為左端全微分的一個(gè)原函數(shù).</p><p> 根據(jù)
33、數(shù)學(xué)分析中關(guān)于線積分與路徑無(wú)關(guān)的幾個(gè)等價(jià)命題及全微分方程的定義,方程(1.27)為全微分方程的充分必要條件是:在所考慮的單連通區(qū)域G內(nèi)有</p><p><b> ?。?.29)</b></p><p> 如果方程(1.27)是全微分方程,設(shè)u(x,y)是使得條件(1.28)成立的可微函數(shù),則方程(1.27)就成為du(x,y)=0.設(shè)函數(shù)y=(x)是方程(1.27
34、)的解,</p><p><b> 則</b></p><p><b> (1.30)</b></p><p><b> 從而</b></p><p><b> (1.31)</b></p><p><b> 這
35、里C是某個(gè)常數(shù)。</b></p><p> 這表明,全微分方程(1.27)的解是由隱方程</p><p> u=(x,y)=C (1.32)</p><p> 所確定的隱函數(shù)。反之,對(duì)于保證隱方程(1.32)有解的任意常數(shù)C,方程(1.32)確定的隱函數(shù)是則,于是
36、</p><p><b> (1.33)</b></p><p> 這表明是全微分方程(1.27)的解。</p><p> 總而言之,全微分方程(1.27)的解是由隱函數(shù)方程(1.32)所確定的隱函數(shù),反之由方程(1.32)所確定的隱函數(shù)也是必然是全微分方程(1.27)的解。</p><p> 求出方程(1.27
37、)左端微分式的原函數(shù)u(x,y)成為解全微分方程的關(guān)鍵。利用與路徑無(wú)關(guān)的線積分化為定積分求這個(gè)原函數(shù)的方法,不失一般性,取,有</p><p><b> (1.34)</b></p><p> 最后得到全微分方程(1.27)的通解為</p><p> u(x,y)=C,這里C為任意常數(shù)。</p><p> 1.4
38、.2 積分因子</p><p> 下面我們看這樣一個(gè)簡(jiǎn)單的方程</p><p> ydx-xdy=0.</p><p> 它顯然不是全微分方程,但只要分別乘以下列因子</p><p><b> 方程就分別化為</b></p><p> 這些都是全微分方程,相應(yīng)得到方程的通解分別為<
39、/p><p> 由此可見,非全微分方程可以通過(guò)乘以某個(gè)非零因子而化為全微分方程,并且這樣的因子還不是唯一的,相應(yīng)于方程乘以不同的這樣的因子化成的全微分方程所得到的通解在形式上還可以是多樣的。</p><p> 如果存在連續(xù)可微的函數(shù)使得</p><p><b> ?。?.35)</b></p><p> 成為一個(gè)全微分
40、方程,則稱為方程(1.27)的一個(gè)積分因子.顯然,方程(1.35)與方程(1.27)是同解的.</p><p> 為了求方程(1.27)的積分因子,首先需了解作為方程(1.27)的因子所具備的條件.根據(jù)積分因子的定義及全微分方程的充要條件(1.33),為方程(1.27)的積分因子當(dāng)且僅當(dāng)</p><p><b> 即</b></p><p>
41、;<b> ?。?.36)</b></p><p> 這意味著要求一般方程(1.27)的積分因子,需要一個(gè)一階線性偏微分方程(1.36),這比解方程(1.27)本身反而困難得多.盡管如此,我們可以考慮求某些特殊形式的積分因子,以簡(jiǎn)化條件(1.36)使其求解成為可能.</p><p><b> 1.5一階隱方程</b></p>&
42、lt;p> 接著我們討論一階隱方程</p><p> . (1.37)</p><p> 如果可以將從方程中解出,求解方程就歸結(jié)到一個(gè)或幾個(gè)顯示微分方程的求解問(wèn)題.,但是如果難以從方程中解出,則可以采用引進(jìn)參數(shù)的辦法使之變成導(dǎo)數(shù)已解出的方程類型,這就是現(xiàn)在主要介紹的方法。</p><p> 1.5.
43、1可解出y或x的方程的解法</p><p><b> 1.首先討論形如</b></p><p><b> (1.38)</b></p><p> 的方程的解法,這里假設(shè)函數(shù)有連續(xù)的偏導(dǎo)數(shù)。</p><p> 引進(jìn)參數(shù),則(1.38)變?yōu)?lt;/p><p><b&g
44、t; ?。?.39)</b></p><p> 將(1.39)兩邊對(duì)x 求導(dǎo)數(shù),并以代入,得到</p><p><b> (1.40)</b></p><p> 方程(1.40)是關(guān)于x,p的一階微分方程,但它的導(dǎo)數(shù)已解出,于是我們可以按照前面的方法給出解。</p><p> 若已求得(1.40)的通
45、解的形式為:</p><p> , (1.41)</p><p><b> 代入(1.39)得</b></p><p><b> ?。?.42)</b></p><p> 若已求得(1.40)的通解的形式為:</p>
46、;<p><b> ?。?.43)</b></p><p> 則得到(1.38)的參數(shù)形式的通解為</p><p><b> (1.44)</b></p><p> 其中p為參數(shù),c為任意常數(shù)。</p><p> 若已求得(1.40)的通解的形式為:</p>&l
47、t;p><b> ?。?.45)</b></p><p> 則得到(1.38)的參數(shù)形式的通解為:</p><p><b> ?。?.46)</b></p><p> 其中p為參數(shù),c為任意常數(shù)。</p><p><b> 2.形如:</b></p>
48、<p><b> ?。?.47)</b></p><p> 的方程的求解方法與方程(1.38)的求解方法完全類似,這里假設(shè)函數(shù)有連續(xù)偏導(dǎo)數(shù)。</p><p> 引進(jìn)參數(shù)則(1.47)變?yōu)椋?lt;/p><p><b> ?。?.48)</b></p><p> 將式(1.48)兩邊對(duì)y
49、求導(dǎo)數(shù),然后以代入,得到</p><p><b> ?。?.49)</b></p><p> 方程(1.49)是關(guān)于y,p的一階微分方程,但是它的導(dǎo)數(shù)已解出,可按照之前的方法求解。設(shè)求得通解為:</p><p><b> ?。?.50)</b></p><p> 則(1.47)的通解為:<
50、/p><p><b> ?。?.51)</b></p><p> 1.5.2 不顯含y(或x)的方程</p><p><b> 3.現(xiàn)在討論形如</b></p><p> =0 (1.52)</p>&l
51、t;p><b> 的方程的解法。</b></p><p> 記,從幾何的觀點(diǎn)看,代表Oxp平面上的一條曲線,設(shè)把這曲線表示為適當(dāng)?shù)膮?shù)形式</p><p><b> ?。?.53)</b></p><p> 這里t為參數(shù),在注意到,沿方程(1.52)的任何一條積分曲線上,恒滿足基本條件,以(1.53)代入上式得
52、</p><p><b> ?。?.54)</b></p><p><b> 兩邊積分,得到</b></p><p><b> ?。?.55)</b></p><p> 于是得到方程(1.52)的參數(shù)形式的通解為</p><p><b>
53、 (1.56)</b></p><p><b> C為任意常數(shù)。</b></p><p> 1.6 matlab 在一階方程中的應(yīng)用</p><p> [設(shè)計(jì)舉例1] 求方程的通解和滿足初始條件的解.>></p><p> clear all;</p><p>
54、eq='(1+x*exp(x*y))*Dy=-(1+y*exp(x*y))'</p><p> sols=dsolve(eq,'x')</p><p> y=dsolve(eq,'y(0)=-1','x')</p><p><b> 運(yùn)行結(jié)果如下</b></p>
55、<p> sols=--(x*C1+x^2+lambertw(x*exp(-(C1+x))))/x</p><p> y=-(x^2+lambertw(x*wxp(-x^2)))/x</p><p> [設(shè)計(jì)舉例2]判斷方程為線性方程,繪出積分曲線簇.</p><p> >>maple(‘with(DEtools)’);</p&
56、gt;<p> maple(‘odeadvisor(x*diff(y(x),x))+y(x)=sin(x))’)</p><p> 上機(jī)運(yùn)行顯示[_ linear]表示方程確實(shí)為線性。再運(yùn)行如下程序:</p><p> eq='x*Dy+y=sin(x)';</p><p> y=dsolve(eq,'x')&l
57、t;/p><p> for i=-4:0.5:4</p><p> ezplot(subs(y,'C1',i),[-4,4,-5,5],5)</p><p><b> hold on</b></p><p><b> end</b></p><p><
58、;b> 運(yùn)行結(jié)果如下</b></p><p> y =-(C2 + cos(x))/x</p><p> [設(shè)計(jì)舉例3] 求單變量積分因子的函數(shù)與求形如的積分因子.</p><p> function mu=intfactor(M,N)</p><p> G=diff(M,’y’)-diff(N,;x’);<
59、/p><p> if diff(G/N,’y’)= =0 mu=exp(int(G/N,’x’));</p><p> elseif diff(G/M,’x’)= =0 mu=exp(int(-G/M,’y’));</p><p> else mu=’Empty’;</p><p><b> end</b>&l
60、t;/p><p><b> end</b></p><p> 以下面兩個(gè)方程為例調(diào)用這個(gè)函數(shù)求出積分因子:</p><p> >>syms x y</p><p> eq1=’exp(x)+3*y^2+2*x*y*Dy=0’;</p><p> M=’exp(x)+3*y^2’;
61、N=’2*x*y’;</p><p> mu1=intfactor(M,N)</p><p> eq2=’y+(exp(1/y)-x)*Dy=0’;</p><p> M=’y’;N=’exp(1/y)-x’;</p><p> mu2=intfactor(M,N)</p><p> 分別返回所求積分因子為m
62、n1=和mu2=1/.下面是求形如的積分因子的例子:</p><p> >>eq3=’y+(x-3*x^3*y^2)*Dy=0’;</p><p> M=’y’;N=’x-3*x^3*y^2’;</p><p> G=diff(M,’y’)-diff(N,’x’);</p><p> f=simple(G/y*N-x*M)
63、)根據(jù)運(yùn)行所得的f=-3/x/y,再運(yùn)行</p><p> >>mu=subs(exp(int(‘-3/u’,’u’)),’u’,’x*y’)</p><p> 求得積分因子為mu=1/(x*y) .</p><p> [設(shè)計(jì)舉例4]驗(yàn)證齊次方程有形如的積分因子.</p><p> 我們只要驗(yàn)證乘以原方程后會(huì)是全微分方程.
64、</p><p><b> >>clear</b></p><p> eq=sym(‘x^2*y-2*x*y+(3*x^2*y-3^3)*Dy=0’)</p><p> M=sym(‘x^2*y-2*x*y^2’);N=sym(‘3*x^2*y-x^3’);</p><p> mu=simple(1/
65、(‘x’*M+’y’*N))</p><p> eq1=simple(mu*eq)</p><p> 首先請(qǐng)注意到這里都是用sym()定義方程與函數(shù)表達(dá)式,因?yàn)橐獙?duì)它們進(jìn)行乘積的運(yùn)算.這里求出了用乘以原方程以后的新方程:</p><p> eq1=1/y-2/x+3/y*Dy-x/y^2*Dy=0.</p><p> 為了用Maple
66、的DEtools中的方程類型判斷指令odeadvisor,得先將它的形式加以改變:</p><p> >>maple(‘with(DEtools)’);</p><p> Maple(‘odeadvisor(1/y(x)-2/x+(3/y(x)-x/y(x)^2*diff(y(x),x)=0,[exact])’)</p><p> 答案是ans=[
67、_ exact].</p><p> [設(shè)計(jì)舉例5]解隱方程</p><p> 先試圖在MATLAB中直接用dsolve求解,</p><p> sols=dsolve(‘Dy+2*x*Dy-y=0’,’x’)</p><p> 結(jié)果返回一個(gè)形式相當(dāng)復(fù)雜的顯式解.對(duì)于某些隱方程,我們通過(guò)Maple接口,用dsolve指令加上參數(shù)imp
68、licit可能求出方程的參數(shù)形式的解.</p><p> >>clear all</p><p> sols=maple(‘dsolve(diff(y(x),x)^3+2*x*diff(y(x),)-y(x)=0,implicit)’)</p><p><b> 返回的結(jié)果為</b></p><p>
69、 sols=[x(_T)=1/_T^2*(-3/4*_T^4+_C1),</p><p> y=(_T)=_T^3+2/_T*(-3/4*_T^4+_C1)],y(x)=0</p><p> 有時(shí)也可以再加上參數(shù)parametric制定返回參數(shù)形式的解.</p><p> >>sols=maple(‘dsolve(diff(y(x),x)^3+2*
70、x*diff(y(x),x)-y(x)=0,</p><p> parametric,implicit)’)</p><p> 返回結(jié)果形式上略有不同:</p><p> sols=[x(_T)=1/4*(-3*_T^4+2*_C1)/_T^2,</p><p> y(_T)=1/2*(-_T^4+2*_C1)/_T].</p&
71、gt;<p> 第二章 高階微分方程的解法</p><p> 2.1 線性微分方程的一般理論</p><p> 我們討論如下的n階線性微分方程</p><p><b> (2.1)</b></p><p> 其中都是區(qū)間上的連續(xù)函數(shù)。</p><p><b>
72、 如果,則方程變?yōu)椋?lt;/b></p><p><b> (2.2)</b></p><p> 我們稱之為n階齊次線性微分方程,而稱(2.1)為n階非齊次線性微分方程。</p><p> 定理1 如果都是區(qū)間上的連續(xù)函數(shù),則對(duì)于任一及任意的,方程(2.1)存在唯一解,定義于區(qū)間上,且滿足初值條件</p><
73、p><b> ?。?.3)</b></p><p> 1 齊次線性微分方程的解的性質(zhì)與結(jié)構(gòu)</p><p> 首先討論齊次線性微分方程</p><p> 定理2 如果是方程2.2的k個(gè)解,則它們的線性組合也是2.2的解,這里是任意常數(shù)。</p><p> 特別地,當(dāng)時(shí),即方程2.2有解</p>
74、<p><b> ?。?.4)</b></p><p> 它含有n個(gè)任意常數(shù),我們需要知道2.4在什么條件下能夠成為2.2的通解,為了討論需要,我們引進(jìn)函數(shù)線性相關(guān)和線性無(wú)關(guān)及朗斯基行列式等概念。</p><p> 如果存在不全為零的常數(shù),使得恒等式</p><p> 對(duì)于所有都成立,稱這些函數(shù)是線性相關(guān)的,否則就稱線性無(wú)關(guān)。
75、</p><p> 由定義在區(qū)間上的k個(gè)可微k-1次的函數(shù)所組成的行列式</p><p><b> 稱為朗斯基行列式。</b></p><p> 定理 3 若函數(shù)在區(qū)間上線性相關(guān),則在[a,b]上它們的朗斯基行列式。否則</p><p> 定理 4 n階齊次線性微分方程2.2一定存在n個(gè)線性無(wú)關(guān)的解。</
76、p><p> 定理5 如果 是方程2.2的n個(gè)線性無(wú)關(guān)的解,則方程2.2的通解可表示為:</p><p><b> (2.5)</b></p><p><b> 其中是任意常數(shù)。</b></p><p> 2.非奇次線性微分方程與常數(shù)變易法</p><p> 考慮n階
77、非齊次線性微分方程</p><p><b> ?。?.1)</b></p><p> 易見(2.2)是它的特殊形式,可指出兩者之間解的性質(zhì)和結(jié)構(gòu)有十分密切的聯(lián)系,可直接驗(yàn)證如下性質(zhì)</p><p> 性質(zhì)1 如果方程(2.2)的解,則也是(2.1)的解</p><p> 性質(zhì)2 方程(2.1)的任意兩個(gè)解之差必是(
78、2.2)的解</p><p> 其次我們有如下的定理</p><p> 定理6 設(shè)為方程(2.2)的基本解組,而是方程(2.1)的某一解,則方程(2.1)的通解可表為:</p><p><b> ?。?.6)</b></p><p><b> 其中 為任意常數(shù)</b></p>&
79、lt;p> 2.2常系數(shù)線性微分方程的解法</p><p> 1.復(fù)值函數(shù)與復(fù)值解</p><p> 如果對(duì)于區(qū)間中的每一實(shí)數(shù)t,有復(fù)數(shù)與它對(duì)應(yīng),其中是在區(qū)間上定義的實(shí)數(shù),是虛數(shù)單位,我們稱為在區(qū)間上給定了一個(gè)復(fù)值函數(shù)。</p><p> 在討論常系數(shù)線性微分方程時(shí),函數(shù)將起著重要的作用,這里K是復(fù)數(shù)常數(shù),設(shè)是任一復(fù)數(shù),這里是實(shí)數(shù),而t為實(shí)變量。<
80、;/p><p><b> 我們定義 </b></p><p><b> (2.7)</b></p><p><b> 可推得</b></p><p><b> ?。?.8a)</b></p><p><b> ?。?.8
81、b)</b></p><p> 如果以表示復(fù)數(shù)的共軛復(fù)數(shù),那么容易得:</p><p><b> ?。?.9)</b></p><p> 定理7 如果方程(2.2)中所有系數(shù)都是實(shí)值函數(shù),而 是方程的復(fù)值解,則的實(shí)部、虛部和共軛復(fù)值解函數(shù)也都是方程(2.2)的解</p><p> 2.2設(shè)齊次線性微分方
82、程中所有系數(shù)都是常數(shù),即方程有如下形狀</p><p><b> ?。?.10)</b></p><p> 其中為n階常系數(shù)齊次線性微分方程.正如我們前面所指出的,它的求解問(wèn)題可以歸結(jié)為代數(shù)方程求根問(wèn)題,現(xiàn)在就來(lái)具體討論方程(2.10)的解法.按照一般的理論,為了求方程(4.10)的通解,只需要求出它的基本解組.下面介紹求(2.10)的基本解組的特征根法</p
83、><p> 回顧一階常系數(shù)齊次線性微分方程</p><p><b> ,</b></p><p> 我們知道它有形如的解,且它的通解就是.這啟示我們對(duì)于方程(2.10)也去試求指數(shù)函數(shù)形式的解</p><p> , (2.11)</p><p&g
84、t; 其中是待定常數(shù),可以是實(shí)的,也可以是復(fù)的.</p><p><b> 注意到</b></p><p> 其中是的n次多項(xiàng)式.易知,(2.11)為方程(2.10)的解的充要條件是是代數(shù)方程</p><p><b> (2.12)</b></p><p> 的根.因此,方程(2.12)將
85、起著預(yù)示方程(2.10)的解的特性的作用,我們稱它為方程(2.10)的特征方程,它的根就稱為特征根.下面根據(jù)特征根的不同情況分別進(jìn)行討論.</p><p><b> 特征根是單根的情形</b></p><p> 設(shè)是特征方程(2.12)的n個(gè)彼此不相等的根,則相應(yīng)地方程(2.10)有如下n個(gè)解:</p><p><b> ?。?.
86、13)</b></p><p> 我們指出這n個(gè)解在區(qū)間上線性無(wú)關(guān),從而組成方程的基本解組.事實(shí)上,這時(shí)</p><p> 如果均為實(shí)數(shù),則2.13是方程的線性無(wú)關(guān)的實(shí)數(shù)解。</p><p> 如果特征方程有復(fù)根,復(fù)根將成對(duì)出現(xiàn)。設(shè),則,也是特征根,因而與這對(duì)共軛復(fù)根對(duì)應(yīng)的,方程2.10兩個(gè)復(fù)值解</p><p><b
87、> (2.14a)</b></p><p><b> ?。?.14b)</b></p><p> 根據(jù)定理7,它們的實(shí)部和虛部也是方程的根,這樣一來(lái),對(duì)應(yīng)于特征方程的一對(duì)共軛復(fù)根,我們可求得方程2.10的兩個(gè)實(shí)根</p><p> ?。?)特征根有重根的情形</p><p> 設(shè)特征方程有k重根,
88、則如所周知</p><p> , (2.15)</p><p> 先設(shè),即特征方程有因子,于是</p><p> 也就是特征方程的形狀為</p><p><b> ?。?.16)</b></p><p> 而對(duì)應(yīng)的2.10變?yōu)?lt;/p>
89、<p><b> ?。?.17)</b></p><p> 易見它有k個(gè)解,而且它們是線性無(wú)關(guān)的,這樣一來(lái),特征方程的k重零根就對(duì)應(yīng)于方程(2.10)的k個(gè)線性無(wú)關(guān)解。</p><p> 如果,我們作變量變換,注意到</p><p><b> ?。?.18)</b></p><p>
90、;<b> 可得</b></p><p><b> (2.19)</b></p><p> 于是(2.10)化為:</p><p> 其中仍為常數(shù),而相應(yīng)的特征方程為</p><p><b> ?。?.20)</b></p><p><b&
91、gt; 直接計(jì)算得</b></p><p><b> ,</b></p><p><b> ?。?.21)</b></p><p> 可見(2.12)的根,對(duì)應(yīng)于(2.20)的根,而且重?cái)?shù)相同,這樣問(wèn)題就化為前面已經(jīng)討論過(guò)的情形了。</p><p> 我們知道,方程(2.20)的
92、對(duì)應(yīng)于方程(2.19)的個(gè)解,因而對(duì)應(yīng)于特征方程(2.12)的,方程(2.10)有個(gè)解 </p><p><b> ?。?.22)</b></p><p> 同樣,假設(shè)特征方程(2.10)的其他根的重?cái)?shù)依次為,而且,則方程(2.10)對(duì)應(yīng)地有解</p><p><b> ?。?.23)</b></p>&
93、lt;p> 這樣(2.22)和(2.23)就組成全部2.10的全部解。</p><p> 對(duì)于特征方程有重復(fù)根的情況,譬如假設(shè)是k重特征根,則也是k重特征根,仿(1)一樣處理,我們將得到方程(2.10)的2k個(gè)實(shí)值解</p><p> 2.3.非齊次線性微分方程的比較系數(shù)法</p><p> 現(xiàn)在討論常系數(shù)非齊次線性微分方程</p>&l
94、t;p><b> ?。?.24)</b></p><p> 的求解問(wèn)題,這里是常數(shù)。</p><p> 本來(lái)有了前面的討論結(jié)果這一問(wèn)題已經(jīng)可以解決了,但是求解往往是比較繁瑣的,而且必須經(jīng)過(guò)積分運(yùn)算,下面介紹當(dāng)具有某些特殊形狀時(shí)所適用的方法—比較系數(shù)法</p><p><b> 類型Ⅰ</b></p>
95、<p> 設(shè),其中及為實(shí)常數(shù),那么方程(2.24)有形如</p><p><b> ?。?.25)</b></p><p> 的特解,其中k為特征方程F()=0的根的重?cái)?shù)(單根相當(dāng)于k=1;當(dāng)不是特征根時(shí),取k=0),而是待定常數(shù),可以通過(guò)比較系數(shù)來(lái)確定.</p><p><b> 如果,則此時(shí)</b>
96、</p><p> 現(xiàn)在再分兩種情形討論.</p><p> 在不是特征根的情形,即,因而,這時(shí)取k=0,以代入方程(2.24),并比較t的同次冪的系數(shù),得到常數(shù)必須滿足的方程</p><p><b> ?。?.26)</b></p><p> 注意到,這些待定常數(shù)可以從方程組(2.26)唯一地逐個(gè)確定出來(lái).<
97、;/p><p> (b)在是k重特征根的情形,,而,也就是.這時(shí)相應(yīng)地,方程(2.24)將為</p><p><b> (2.27)</b></p><p> 令,則方程(2.27)化為</p><p><b> ?。?.28)</b></p><p> 對(duì)方程()來(lái)說(shuō),
98、由于已不是它的特征根.因此,由(1)知它有形如的特解,因而方程(2.28)有特解滿足</p><p> 這表明是t的m+k次多項(xiàng)式,其中t的冪次的項(xiàng)帶有任意常數(shù).但因我們只需要知道一個(gè)特解就夠了.我們特別地取這些任意常數(shù)均為零,于是我們得到方程(2.28)的一個(gè)特解</p><p> 這里是已確定了的常數(shù).</p><p> ②如果,則此時(shí)可像 做法一樣
99、,作變量變換,將方程(2.24)化為</p><p><b> ?。?.29)</b></p><p> 其中都是常數(shù).并且特征方程(2.12)的根對(duì)應(yīng)于方程(2.29)的特征方程的零根,并且重?cái)?shù)也相同.因此,利用上面的結(jié)果就有如下結(jié)論:</p><p> 在不是特征方程(2.12)的根的情形,方程(2.29)有特解</p>
100、<p> 從而方程(2.24)有特解</p><p><b> ;</b></p><p> 在是特征方程(2.12)的k重根的情形,方程(2.29)有特解,從而方程(2.24)有特解</p><p><b> 類型Ⅱ</b></p><p> 設(shè),其中為常數(shù),而A(t),B(t
101、)是帶實(shí)數(shù)系數(shù)的t的多項(xiàng)式,其中一個(gè)的次數(shù)為m,而另一個(gè)的次數(shù)不超過(guò)m,那么我們有如下結(jié)論:方程(2.24)有形如</p><p><b> ?。?.30)</b></p><p> 的特解,這里k為特征方程的根的重?cái)?shù),而P(t),Q(t)均為待定的帶實(shí)系數(shù)的次數(shù)不高于m的t的多項(xiàng)式,可以通過(guò)比較系數(shù)的方法來(lái)確定.</p><p> 事實(shí)上
102、,回顧一下類型Ⅰ的討論過(guò)程,易見當(dāng)不是實(shí)數(shù),而是復(fù)數(shù)時(shí),有關(guān)結(jié)論仍然正確.現(xiàn)將f(t)表為指數(shù)形式</p><p> 根據(jù)非齊次線性微分方程的疊加原理,方程</p><p><b> 與</b></p><p> 的解之和必為方程(2.24)的解.</p><p> 注意到,易知,若為的解,則必為的解.因此,直接
103、利用類型Ⅰ的結(jié)果,可知方程(2.24)有解形如</p><p><b> =,</b></p><p> 其中D(t)為t的m次多項(xiàng)式,而P(t)=2Re{D(t)},Q(t)=2lm{D(t)}.</p><p> 2.4用matlab解高階線性方程</p><p> [設(shè)計(jì)舉例1]求方程的通解及其滿足初始條件
104、的特解這里為常數(shù)。</p><p><b> 程序如下:</b></p><p><b> clear</b></p><p> syms m k t0 x0 v0</p><p> eqn='m*D2x+k*x=0'</p><p> genso
105、l=dsolve(eqn,'t')</p><p> sol=simple(dsolve(eqn,'x(t0)=x0','Dx(t0)=v0'))</p><p> 運(yùn)行結(jié)果保留了相應(yīng)的文字系數(shù)和參數(shù)m, k, t0, x0, v0.如下:</p><p> gensol = C1*sin((k/m)^(1/2)
106、*t)+C2*cos((k/m)^(1/2)*t)</p><p> sol=(x0*(k/m)^(1/2)*cos((k/m)^(1/2)*t-t0*(k/m)^(1/2))+v0*sin((k/m)^(1/2)*t-t0*(k/m)^(1/2)))/(k/m)^(1/2)</p><p> [設(shè)計(jì)舉例2]求方程的通解</p><p><b> 程
107、序?yàn)椋?lt;/b></p><p> % 編輯一個(gè)函數(shù)文件取名為hlgsol.m保存在當(dāng)前工作路徑</p><p> function sol=hlgsol(lamda,lamdnum,m,C) </p><p><b> syms t</b></p><p><b> phi=[];</
108、b></p><p> for j=1:lamdnum</p><p> for s=1:m(j)</p><p> if imag(lamda(j))==0 </p><p> phi=[phi;t^(s-1)*exp
109、(lamda(j)*t)]; </p><p> elseif imag(lamda(j))>0 </p><p> phi=[phi;t^(s-1)*cos(imag(sym(lamda(j)))*t)*exp(real(sym
110、(lamda(j)))*t)];</p><p> else phi=[phi;t^(s-1)*sin(-imag(sym(lamda(j)))*t)*exp(real(sym(lamda(j)))*t)];</p><p><b> end </b></p><p><b> end</b></p>
111、<p><b> end </b></p><p> sol=C*phi;</p><p><b> return</b></p><p> % 通過(guò)下面的命令腳本程序調(diào)用函數(shù)M文件hlgsol.m求高階齊線性方程通解</p><p><b> clear</b&
112、gt;</p><p> coef=sym('[1,5,12,18,18,12,5,1]');</p><p> lambda=sym('lambda'); n=length(coef)-1; chareq='0';</p><p> for j=1:n+1 chareq=chareq+coef(j)*lambd
113、a^(n+1-j); end</p><p> solvlamd=[solve(chareq,lambda)]</p><p> % 求出了所有特征根,然后運(yùn)行 </p><p> lamda=[-1/2+1/2*i*3^(1/2),-1/2-1/2*i*3^(1/2),-1];</p><p> lamdnum=3; m=[2,2,
114、3];</p><p> syms C1 C2 C3 C4 C5 C6 C7</p><p> C=[C1,C2,C3,C4,C5,C6,C7];</p><p> sols=hlgsol(lamda,lamdnum,m,C)</p><p><b> 運(yùn)行結(jié)果如下:</b></p><p&g
115、t; solvlamd =</p><p><b> -1</b></p><p><b> -1</b></p><p><b> -1</b></p><p> - (3^(1/2)*i)/2 - 1/2</p><p> - (3^(1
116、/2)*i)/2 - 1/2</p><p> (3^(1/2)*i)/2 - 1/2</p><p> (3^(1/2)*i)/2 - 1/2</p><p><b> sols =</b></p><p> C5/exp(t) + (C6*t)/exp(t) + (C1*cos((3^(1/2)*t)/2))
117、/exp(t/2) + (C3*sin((3^(1/2)*t)/2))/exp(t/2) + (C7*t^2)/exp(t) + (C2*t*cos((3^(1/2)*t)/2))/exp(t/2) + (C4*t*sin((3^(1/2)*t)/2))/exp(t/2)</p><p> [設(shè)計(jì)舉例3] 解歐拉方程</p><p> 程序如下% 在編輯器編輯一個(gè)取名為eulachar
118、sol.m的函數(shù)文件用于構(gòu)造特征方程并求出其解</p><p> function ksol=eulacharsol(coef) </p><p> n=length(coef)-1;</p><p><b> sumer=0;</b><
119、;/p><p><b> for j=0:n</b></p><p> product=1;</p><p> for i=0:j-1</p><p> product=product*(sym('k')-i);</p><p><b> end</b>
120、</p><p> sumer=sumer+coef(n+1-j)*product;</p><p><b> end</b></p><p> ksol=solve(sumer,'k');</p><p><b> return</b></p><p&g
121、t; % 在編輯器編輯一個(gè)取名為eulagensol.m的函數(shù)文件最后調(diào)用其求出歐拉方程的通解,兩個(gè)文件保存在當(dāng)前工作路徑.</p><p> function sols=eulagensol(kroot,knum,kmul,C)</p><p><b> syms x</b></p><p><b> phi=[];<
122、/b></p><p> for j=1:knum</p><p> for s=1:kmul(j)</p><p> if imag(kroot(j))==0 phi=[phi;(log(x))^(s-1)*x^sym(kroot(j))];</p><p> elseif imag(kroot(j))>0 </p
123、><p> phi=[phi;x^real(sym(kroot(j)))*log(x)^(s-1)*cos(imag(sym(kroot(j)))*log(x))];</p><p> else phi=[phi;x^real(sym(kroot(j)))*log(x)^(s-1)*sin(-imag(sym(kroot(j)))*log(x))]; </p><p&
124、gt; end </p><p><b> end;</b></p><p><b> end </b></p><p> sols=C*phi; </p><p><b> return&
125、lt;/b></p><p> % 將以下程序段在編輯器中編輯好,并將其命名為eulafch.m</p><p><b> clear</b></p><p> coef=input('請(qǐng)輸入歐拉方程由高階到低階順序的系數(shù)向量: '); % [1 13 47 48 9 -3 7 -8]</p>
126、<p> ksol=eulacharsol(coef)</p><p> % 求解特征方程,以下根據(jù)運(yùn)行過(guò)程的提示輸入從這里的結(jié)果得到的所有不同特征根,個(gè)數(shù)</p><p> % 并計(jì)算出它們的重?cái)?shù).</p><p> kroot=input('請(qǐng)輸入不同特征根構(gòu)成的向量kroot: '); % [
127、 1/2-1/2*i*3^(1/2), 1/2+1/2*i*3^(1/2),2]</p><p> knum=length(kroot); </p><p> kmul=input('請(qǐng)輸入不同特征根重?cái)?shù)對(duì)應(yīng)向量kroot的順序構(gòu)成的向量kmul: '); % [2,2,3]</p><p> C=input('請(qǐng)輸入任意
128、常數(shù)C1,C2,...構(gòu)成的向量[sym(''C1''),sym(''C2''),...,sym(''Cn'')]: ');</p><p> sols=eulagensol(kroot,knum,kmul,C)</p><p><b> 運(yùn)行結(jié)果如下</b>&
129、lt;/p><p> 請(qǐng)輸入歐拉方程由高階到低階順序的系數(shù)向量: [1 13 47 48 9 -3 7 -8]</p><p><b> ksol =</b></p><p><b> 2</b></p><p><b> 2</b></p><p>
130、;<b> 2</b></p><p> 1/2 - (3^(1/2)*i)/2</p><p> 1/2 - (3^(1/2)*i)/2</p><p> (3^(1/2)*i)/2 + 1/2</p><p> (3^(1/2)*i)/2 + 1/2</p><p> 請(qǐng)輸入不同特
131、征根構(gòu)成的向量kroot: [ 1/2-1/2*i*3^(1/2), 1/2+1/2*i*3^(1/2),2]</p><p> 請(qǐng)輸入不同特征根重?cái)?shù)對(duì)應(yīng)向量kroot的順序構(gòu)成的向量kmul: [2,2,3]</p><p> 請(qǐng)輸入任意常數(shù)C1,C2,...構(gòu)成的向量[sym('C1'),sym('C2'),...,sym('Cn')
132、]: [1]</p><p><b> sols =</b></p><p><b> 2</b></p><p> x^(1/2)*sin((3^(1/2)*log(x))/2)</p><p> x^(1/2)*sin((3^(1/2)*log(x))/2)*log(x)</p&g
133、t;<p> x^(1/2)*cos((3^(1/2)*log(x))/2)</p><p> x^(1/2)*cos((3^(1/2)*log(x))/2)*log(x)</p><p><b> x^2</b></p><p> x^2*log(x)</p><p> x^2*log(x)^2
134、</p><p> 第三章 線性微分方程組</p><p> 先考慮常系數(shù)齊次線性方程組</p><p><b> ?。?.1)</b></p><p> 的基解矩陣的求法,這里矩陣A是的常數(shù)矩陣.對(duì)于非齊次線性方程組</p><p> 的解,可以應(yīng)用常數(shù)變異法求出.</p>
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 眾賞文庫(kù)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 用matlab解常微分方程
- 課程設(shè)計(jì)-常微分方程的求解
- 常微分方程的解法,建模,matlab算法
- matlab求解微分方程(組)及偏微分方程(組)
- 常微分方程
- 第五講 matlab求解微分方程
- 畢業(yè)論文----常微分方程積分因子法的求解
- 16 常微分方程
- 常微分方程教案
- 常微分方程答案
- 常微分方程試題
- 基于LS-SVM的常微分方程求解.pdf
- 階常微分方程初等解法畢業(yè)
- 常微分方程畢業(yè)論文doc
- 《matlab求解偏微分方程常見問(wèn)題》
- matlab求解延遲微分方程的注意事項(xiàng)
- 常微分方程模擬試題
- 實(shí)驗(yàn)二微分方程與差分方程模型matlab求解
- 常微分方程答案習(xí)題
- 常微分方程課程總結(jié)
評(píng)論
0/150
提交評(píng)論