1. 程式人生 > 其它 >求解微分方程的一些方法

求解微分方程的一些方法

當年已經學過了,可是忘光了。從知乎上找到了一個課程,可是和之前老師講的不一樣,在這裡說明一下。

求解微分方程,是解一個含有微分的方程。因為含有微分,它和一般的方程可不一樣,求解的結果裡會具有一個常數\(C\)。若想要去掉這個常數\(C\),需要附加條件。這個附加條件表現為:

\[y'(x_1)=e_1,\\ y(x_2)=e_2 \]

假若\(x_1=x_2\),稱這個附加條件下的問題為初值問題。反之,則稱為條件值問題。一般遇見的都是初值問題。

在微分形式以及它的變種中,初值條件僅僅為:

\[y(x_0)=y_0 \]

要解決初值問題,本質上不需要尋找額外的方法。只要完成了求解,再代入初值即可解決初值問題。當然,或許存在額外的解法

大抵來說,這個教程的內容是:將微分方程分為幾類,在這之後,每一類都有自己的獨特解法。

微分方程的分類與計算

  1. 標準形式
\[y'=f(x,y) \]

當然,這只是一個範例。如果標準形式也存在著解法,我們就沒有必要去討論不同形式下的解法了。

  1. 微分形式
\[\frac{dy}{dx}=-\frac{M(x,y)}{N(x,y)}\\ \Rightarrow M(x,y)dx+N(x,y)dy=0 \]

同上標準形式,這只是一個範例。

  1. 可分離變數的形式
\[M(x)dx+N(y)dy=0 \]

直接進行積分,即可求解。這是求解最簡單的一個形式。

\[\int M(x)dx+\int N(y)dy=\int 0=C \]

當然,它存在著求解初值問題的額外方法:

\[假設初值條件:\\ y'(x_0)=y_0,\\ 那麼,可以求解以: \int_{x_0}^xM(x)+\int_{y_0}^yN(y)=0 \]
  1. 齊次方程
\[對於y'=f(x,y),有:f(tx,ty)=f(x,y) \]

假設一個齊次方程:

\[\frac{dy}{dx}=f(x,y) \]

由於不是一個可分離變數的方程,顯然不能夠直接求解。
由於y是x的函式\((y=y(x))\),顯然這個形式可以變化。比如,\(y=x\cdot y(x)\),不過這樣會在符號的使用上引發問題,所以改寫為\(y=x\cdot v(x)\)
為什麼要改寫?因為是齊次方程,\(f(tx,ty)=f(x,y)\)

,如果\(y=x\cdot v(x)\),那麼有:

\[\frac{dy}{dx}=f(x,y)\\ \Rightarrow \frac{dy}{dx}=f(x,xv)=f(1,v)=f(v)\\ 又\because \frac{dy}{dx}=v+x\frac{dv}{dx} \Rightarrow v+x\frac{dv}{dx}=f(v), x\frac{dv}{dx}=f(v)-v=g(v), \frac{dv}{dx}=\frac{g(v)}{x} \]

最後會變化為可分離變數的形式。

  1. 恰當方程
\[對於M(x,y)dx+N(x,y)dy=0,有:\exists F(x,y),使得: \begin{cases} \frac{\partial F(x,y)}{\partial x}=M(x,y)\\ \frac{\partial F(x,y)}{\partial y}=N(x,y) \end{cases} \]\[\therefore M(x,y)dx+N(x,y)dy=0\Rightarrow dF(x,y)=0 \]

首先,我們有著一個確認恰當方程的方法。如下:

\[\frac{\partial M(x,y)}{\partial y}= \frac{\partial N(x,y)}{\partial x} \]

如果符合上式,那麼這就是一個恰當方程。

接下來,我們就可以根據下式確定\(F(x,y)\)

\[\begin{cases} \frac{\partial F(x,y)}{\partial x}=M(x,y)\\ \frac{\partial F(x,y)}{\partial y}=N(x,y) \end{cases} \]

確定後,原式可以變化為:

\[dF(x)=0, \int dF(x)=\int 0=0, \therefore F(x)=C(C為任意常數) \]

這樣就直接得到了對應的隱式解。從這個隱式解,或許可以得到顯式解。

額外的情況,即使原方程不是恰當方程,可以將其變化為恰當方程。具體方式為:

\[M(x,y)dx+N(x,y)dy=0 \Rightarrow I(x,y)[M(x,y)dx+N(x,y)dy]=0 \]

在一些情況下,我們可以通過一些固定的方式來尋找這個\(I(x,y)\),如下:

假若有:

\[(\frac{\partial M}{dy}-\frac{\partial N}{dx})=N\cdot g(x) \quad or\quad M\cdot h(y)\\ \]

那麼有:

\[I(x,y)=e^{\int g(x)dx \quad or \quad -\int h(y)dy} \]

假若有:

\[\begin{cases} M(x,y)=yf(xy)\\ N(x,y)=xg(xy) \end{cases} \]

那麼有:

\[I(x,y)=\frac{1}{xM-yN} \]
  1. 線性方程
\[y'+p(x)y=q(x) \]

所有的線性方程,都可以變化為恰當方程,且為:

\[(\frac{\partial M(x,y)}{\partial y}- \frac{\partial N(x,y)}{\partial x})= N \cdot g(x) \]
  1. 伯努利方程
\[y'+p(x)y=q(x)y^n \]

\(令z=y^{1-n}\),即可轉化為線性方程。