求解微分方程

2021-02-08 算法與數學之美

來源:http://blog.csdn.net/yunduanmuxue/article/details/11924089


一般地,含有未知函數及未知函數的導數或微分的方程稱為微分方程。微分方程中出現的未知函數的最高階導數的階數稱為微分方程的


按照不同的分類標準,微分方程可以分為線性或非線性,齊次或非齊次。


一般地,微分方程的不含有任意常數的解稱為微分方程的特解,含有相互獨立的任意常數,且任意常數的個數與微分方程階數相等的解稱為微分方程的通解一般解)。


下面介紹微分方程的求解方法。


一、一階微分方程


一階微分方程具有如下一般形式:



1.可分離變量

這類方程可以化為如下形式:



設  ,可以通過下式求解


如果  ,則易知  也是方程的解。


2.齊次方程

形如



的一階微分方程稱為齊次微分方程,簡稱齊次方程

通過變量替換,可以將這類方程化為可分離變量的方程來求解。

令 

 或 

其中  是新的未知函數,則有


將式代入原方程,得


分離變量,兩邊積分,有


求出積分後,再將  回代,便得到方程的解。


3.可化為齊次方程的微分方程


有些問題本身雖然不是齊次的,但是通過適當變換,可以化為齊次方程。

例如,對於形如


的方程,先求出兩條曲線


的交點  ,於是作平移變換


這時,原方程可以化為齊次方程



4.一階線性微分方程


形如


的方程稱為一階線性微分方程

當  時,上式變為


這個方程稱為一階齊次線性微分方程

一階齊次線性微分方程是可分離變量的方程,由上面的方法可以得到方程的通解為


其中  是任意常數。

下面再討論一階齊次非線性方程的通解。

將方程變形為


兩邊積分,得


若記  ,則



與齊次方程的通解相比較,易見其表達形式一致,只需將  換為函數  。由此可以引出求解一階非齊次線性微分方程的常數變異法。即在求出齊次方程通解後,將通解中的常數  變異為待定函數  ,並設一階非齊次方程的通解為


求導,得


將  和  代入原方程,得


積分,得


從而得到一階非齊次線性微分方程的通解為


上述公式可以寫成


從中可以看出,一階非齊次線性微分方程的通解是對應的齊次線性方程的通解與其本身的一個特解之和。這個結論對高階非齊次線性方程亦成立。


5.伯努利方程


形如


的方程稱為伯努利方程。其中  為常數,且  。

伯努利方程是一類非線性方程,但是通過適當的變換,就可以把它轉化為線性方程。

方程兩端同時除以  ,得


於是,令  ,就得到關於變量  的一階線性方程


利用線性方程的求解方法求出通解後,再回代原變量,便可得到伯努利方程的通解




二、特殊二階微分方程


本節介紹三種可以通過化簡來求解的二階微分方程形式。


1.  型


這是最簡單的二階微分方程,求解方法是逐次積分,得到


這類方程的解法,可以推廣到  階微分方程


只要連續積分  次,就可以得到這個方程含有  個任意常數的通解。


2.  型


這類方程的特點是不顯含未知函數  ,求解方法是:

令  ,則  ,原方程化為以  為未知函數的一階微分方程


設其通解為


然後回代變量,又可以得到一個一階微分方程


對它進行積分,可以得到方程的通解



3.   型


這類方程的特點是不顯含未知函數  ,求解方法是:

把  暫時看成自變量,並作變換  ,於是,由符合函數求導法則,有


這樣將原方程化為


這是一個關於變量  的一階微分方程,設它的通解為


這是可分離變量的方程,對其積分得到原方程的通解



三、二階線性微分方程


二階線性微分方程的一般形式是


函數  成為方程的自由項

當  時,方程變為


後者稱為二階齊次線性微分方程,前者稱為二階非齊次線性微分方程


下面介紹幾個重要的定理。


定理一:如果函數 與 是方程 (2) 的兩個解,則


也是方程 (2) 的解,其中  是任意常數。


這個性質表明其次線性方程的解符合疊加定理


補充一個定義:如果  是定義在區間  內的兩個函數。如果存在兩個不全為零的常數  ,使得在區間  內恆有


則稱這兩個函數在區間  內線性相關,否則稱為線性無關


定理二:如果函數 與 是方程 (2) 的兩個線性無關的特解,則


也是方程 (2) 的通解,其中  是任意常數。


定理三:設  是方程 (1) 的一個特解,而  是其對應的齊次方程 (2) 的通解,則


是二階非齊次線性微分方程 (1) 的通解。


定理四:設  分別是方程


的特解,則  是方程  的特解。


這個定理通常稱為非齊次線性微分方程的解的疊加定理


定理五:設  是方程 


的解,其中  為實值函數, 為純虛數,則  分別是方程


的解。


在方程 (1) 中,係數是隨  變化的,這類方程求解比較困難,下面介紹處理這類方程的兩種方法。


1.降階法


考慮二階齊次線性方程


設  是方程的一個已知的非零特解,作變量替換


其中  為待定係數,求  的一階導數和二階導數,得


將他們代入方程,得


各項係數是  的已知函數,因為  是原方程的解,所以,其中  的係數為 0 。

故上式可以化為


再做變量替換  ,得


分離變量


兩邊積分,得其通解


其中  為任意常數。

對  積分,得


其中  為任意常數。

代回原變量,就得到原方程的通解


這個公式稱為二階線性微分方程的劉維爾公式

綜上所述,對於二階齊次線性方程,如果已知其一個非零特解,作變量替換  就可將其姜維一階齊次線性方程,從而求得通解。

對於二階非齊次線性方程,若已知其對應的齊次方程的一個特解,作同樣的變量替換(因為這種變換並不影響方程的右端),也能使非齊次方程降一階。


2.常數變異法


設有二階非齊次線性方程


其中  在某區間上連續,如果其對應的齊次方程


的通解  已經求得,可通過如下方法求得其通解。

設非齊次方程具有形如


的特解,對  求導數,得



把特解代入原方程,可以得到確定  的一個方程,因為這裡有兩個未知函數,所以還需要添加另一個條件,為計算髮辮,我們補充如下條件:  ,這樣



代入原方程,並注意到  是齊次方程的解,經整理得


與補充條件聯立,得方程組


因為  線性無關,即  常數,所以


設  則有  ,所以上述方程有唯一解,解得


積分並取一個原函數,得


於是,所求特解為


所以,所求方程的通解為

  

為了促進交流,方便大家交換知識與經驗,本號開通了「算數交流群」,願意進群的讀者,請先加個人微信號MissAmath !添加個人微信號時,請備註「算數」二字,謝謝!


熱經典文章推薦:

矩陣的秩與行列式的幾何意義

微信紅包實現原理猜想


回復以下關鍵字獲取相關文章:


數據挖掘 | 機器學習 | 數學之美 | 遊戲算法 | 生活數學 | 排名算法|大型網站技術演進 | 數學名人 | 學科概述 | 計算機科學 | 搜尋引擎



據說好多人都不知道長按圖片也能關注,你知道嗎?


相關焦點

  • 先從求解符號微分方程講起
    本次公眾號將介紹符號微分方程的求解方法,包括求解符號微分方程的通解和特解。在MATLAB中,求解符號微分方程通解的命令格式為:y=dsolve('equation','x');其中,equation為符號微分方程,x為符號自變量。利用該命令可以非常方便的求出符號微分方程的解析通解。舉例說明:求微分方程y'=x+1的通解y=f(x)。
  • 隱式常微分方程擬合求解
    但如果微分方程是隱函數形式,則不能直接求解,必須先進行轉換成顯示形式後,方能進行求解。下面以三個實例分別介紹一般常微分方程擬合及不同形式的隱式常微分方程擬合。1.案例之一:一般正常微分方程(顯示微分方程)求解已知常微分方程及相關數據如下:t0,0.5,1.5,3.0,6.0,12.0,18.0,24.0,36.0y0,66.8041,40.0633,22.8622,8.2312,3.1214,1.8126,1.0809,0.7138其中x系列數據缺失。
  • 用Python數值求解偏微分方程
    然而,只有很少的微分方程可以解析求解,尤其對於偏微分方程,能解析求解的種類更是寥寥可數。更多的微分方程可以採用數值法進行求解,只要精度足夠高,就可以滿足科學和工程上的需求。Python的Numpy、Scipy工具包可以很好地實現此功能,Matplotlib工具包則可以將求解結果畫為非常直觀的圖形。接下來,我們先以常微分方程的數值求解為例,引入差分的思想,再將其推廣到偏微分方程中。
  • 偏微分方程:作用、分析與數值求解
    報告題目:偏微分方程:作用、分析與數值求解 報 告 人:江松 研究員 北京應用物理與計算數學研究所 報告時間:2020年9月27日9:00 報告地點:數學學院二樓報告廳 校內聯繫人:張然zhangran@jlu.edu.cn 報告摘要: 科學與工程技術
  • 簡諧振動二階微分方程求解過程
    很多同學對於簡諧振動的表達形式,尤其是關於角頻率的定義,與振動的二次微分方程的關係不是很清楚。
  • 2019數學建模國賽|Matlab 求解微分方程(組)
    2.函數 dsolve 求解的是常微分方程的精確解法,也稱為常微分方程的符號解.但是,有大量的常微分方程雖然從理論上講,其解是存在的,但我們卻無法求出其解析解,此時,我們需要尋求方程的數值解,在求常微分方程數值解方面,MATLAB 具有豐富的函數,將其統稱為 solver,其一般格式為:           [T,Y]=solver(odefun,tspan,y0)
  • 《常見一階微分方程》類型及其一般求解思路與步驟
    一、《高等數學》一階微分方程分類:第一類:可分離變量的微分方程及其分離變量的求解方法,包括齊次微分方程(換元法)。 第二類:一階線性微分方程,其中齊次線性微分方程的求解歸結為可分離變量的微分方程;而非齊次線性微分方程基於常數變易法,或稱為待定函數法,直接得到非齊次線性微分方程的通解或者基於線性微分方程解的結構求得其一個特解來求通解:非齊次線性微分方程的特解
  • 高等數學《常微分方程》內容總結、題型求解方法與典型題
    一、求解一階微分方程的基本思路1.改寫結構,對比標準可求解類型適當變換微分方程描述形式,比對標準類型方程結構。●伯努利方程:通過兩端同時除以yn,將方程轉換為一階線性微分方程求解.●全微分方程:它的判定和求解方法,使用曲線積分相關的理論與方法求解.
  • 薛定宇教授大講堂(卷Ⅴ):MATLAB微分方程求解|文末贈書
    2.2.4 Legendre微分方程與 Legendre函數  192.2.5 Airy函數  202.3常係數線性微分方程的求解  212.3.1線性常係數微分方程解析解的數學描述  212.3.2基於 Laplace變換的求解方法 222.3.3非齊次微分方程的求解  242.3.4非零初值的微分方程求解
  • MATLAB常微分方程數值求解
    有一類常微分方程,其解的分量有的變化很快,有的變化很慢,且相差懸殊,這就是所謂的剛性問題 (Stiff) 。對於剛性問題,數值解算法必須取很小步長才能獲得滿意的結果,導致計算量會大大增加。解決剛性問題需要有專門方法。
  • AI已能求解微分方程,數學是這樣一步步「淪陷」的
    今天我們就簡單梳理一下機器學習解方程的近些年最新進展。首先,先得說明,所謂的AI解方程,解的可不是簡單的數值方程!人家正經求解的是面對實際工程問題的偏微分方程,又稱PDE。什麼是偏微分方程呢?,那他就被叫做偏微分方程。簡稱PDE。大自然中上到宇宙行星運動,下到微觀粒子的相互作用都是用PDE進行描述的。而我們介紹的就是用AI來求解這些微分方程。
  • 「每周一識」一階非齊次線性微分方程求解及應用舉例
    本文介紹一階非齊次線性微分方程的通解的應用、特解求解舉例,以及二階微分方程可用該通解求解的情形。一、方程通解公式一階非齊次線性微分方程的解析式為:y'+p(x)=q(x),則其通解表達式如下:y=e^[-∫p(x)]dx{∫q(x)*e^[∫p(x)dx]dx+c}.
  • 描述物質運動變化的數學學科:常微分方程、偏微分方程
    隨著分析學對函數引入微分運算,表示未知函數的導數以及自變量之間的關係的方程進入數學家的視野,這就是微分方程。微分方程的形成與發展與力學、天文學、物理學等科學技術的發展密切相關。因為在現實的世界中,物質的運動及其變化規律在數學上是用函數關係來描述的,這意味著問題的解決就是要去尋求滿足某些條件的函數,而這類問題就轉換為微分方程的求解問題。
  • MATLAB常微分方程
    4.8.1  常微分方程眾所周知,對一些典型的常微分方程,能求解出它們的一般表達式,並用初始條件確定表達式中的任意常數。但實際中存在有這種解析解的常微分方程的範圍十分狹窄,往往只局限在線性常係數微分方程(含方程組),以及少數的線性變係數方程。對於更加廣泛的、非線性的一般的常微分方程,通常不存在初等函數解析解。
  • MATLAB偏微分方程
    4.8.2  偏微分方程在自然科學的很多領域內,都會遇到微分方程初值問題,特別是偏微分方程,它的定解問題是描述自然界及科學現象的最重要的工具
  • 《歐拉方程及微分方程建模》思路與方法
    一、歐拉方程及其求解方法具有結構即於是,就可以通過常係數線性微分方程的求解方法求該方程的通解了。 【注】歐拉方程其實就是一種線性微分方程的結構,只不過不具有直接的顯性結果,需要換元變換得到。 三、微分方程模型求解實際問題的基本步驟(1) 確定模型類型:注意到實際問題中與數學中的導數相關的常用詞語。
  • [基礎理論]偏微分方程的類型
    偏微分方程(PDE)是真實世界常見現象的一種數學語言描述,二階偏微分方程始終是重要的研究對象。
  • 常微分方程中的重要方程:黎卡提方程(一階二次非線性微分方程)
    前面我們了解了什麼是一階線性微分方程,可分離變量微分方程,以及齊次微分方程,本篇講升上一個高度,一階微分方程中的二次微分方程義大利數學家在17世紀提出了著名的「黎卡提方程」,這個方程看上去挺簡單的,但分析起來相當複雜
  • Matlab/Simulink建模詳解:一階時變偏微分方程的求解
    Matlab/Simulink建模詳解:一階時變偏微分方程的求解 這一次日篤小編來教大家如何在simulnk裡面,求解偏微分方程(Partial Differential Equation-PDE)。跟ODE不同的是,PDE的變量不僅是時間的函數,而且也是空間的函數。舉個簡單的例子:炒菜的鍋的鍋壁溫度,溫度這個變量不僅僅是時間的參數(不同時候,鍋壁同一個點的溫度不一樣),而且也是空間的參數(同一時刻,鍋壁不同位置的溫度不一樣)。求解PDE的時候可以直接求解,也可以先轉化成ODE,再用求解ODE的方法求解。這裡,小編先將PDE轉化成ODE再求解。
  • 常微分方程真題講解
    1.理解常微分方程的概念,理解常微分方程的階、解、通解、初始條件和特解的概念。(二)二階常係數線性微分方程      2.會求解二階常係數齊次線性微分方程。注意:熟記通解的結構公式對於準確快速答題至關重要,更多一階微分方程內容請查看往期推送:常微分方程考點(1):一階常微分方程.