前兩天看到 ei 的部落格裡面有這個,於是去翻了一下**
一些定義:
齊次線性微分方程
如果乙個函式 \(y(x)\) 滿足:
\[a_0y+a_1y'+a_2y''+\cdots+a_ny^=c(x)
\]則稱其滿足乙個線性微分方程。特別地,如果 \(c(x)=0\),則稱其滿足乙個齊次線性微分方程(linear homogeneous differential equation)。
根據所有的係數 \(\\) 是否為常數還可以進一步劃分為常係數齊次線性微分方程和變係數齊次線性微分方程。
p-遞迴數列
如果乙個數列 \(a_n\) 滿足
\(p_0(n)a_n+p_1(n)a_+p_2(n)a_+\cdots+p_d(n)a_=0\)
則稱其為 p-遞迴(p-recursive)的。其中所有的 \(p_i(x)\) 均為多項式。
微分有限
如果乙個函式 \(f(x)\)(不一定是多項式)滿足乙個齊次線性微分方程,則稱其為微分有限(d-finite)的。
關於微分有限函式我們有下述定理:
定理一記函式 \(f(x)\) 的級數展開為 \(f(x)=\sum_a_nx^n\),那麼 \(f(x)\) 是微分有限的當且僅當其係數 \(\\) 是 p-遞迴的。
定理二如果 \(f(x)\) 和 \(g(x)\) 都是微分有限的,那麼 \(\alpha f(x)+\beta g(x),f(x)g(x)\) 也是微分有限的。
定理三如果 \(f(x)\) 微分有限,\(g(x)\) 是代數函式(即只包含加減乘除乘方開方的函式),那麼 \(f(g(x))\) 是微分有限的。
以及一些簡單的例子:
冪函式都是微分有限的:
\[y=x^r\rightarrow ry-xy'=0
\]指數函式是微分有限的:
\[y=a^x\rightarrow y\ln a-y'=0
\]廣義超幾何函式是微分有限的:
不會證,告辭。
以及在資訊學競賽中的應用:所有化簡為多項式的函式都可以使用整式遞推的方式快速求出某一項的係數。同時可以彰顯人類智慧型
非線性微分方程有限差分解法
一般的,對於線性偏微分方程其有限差分解法比較容易編寫.對於非線性微分方程的數值求解而言,一般的,如果採用有限差分的方法進行離散,那麼離散後方程中就會出現非線性項.此時可以將非線性項看作是新的變數,利用已有的對 f x f x f x 迭代解法給出非線性項的近似解.然後再進入下一輪迴圈.總結下,對於非...
有限差分法Eluer演算法(求解常微分方程)
有限差分法c 實現 幾種eluer演算法 在matlab 裡面寫乙個程式 要求用 隱式尤拉法 backward euler 去解決常微分方程。下面是兩個例題,1.x 2x 還給出準確值是 x e 2t 要求 求出t 0,這是乙個初始值,然後算出在區間 0,5 的值。還給出 步長 step size ...
多元函式微分學
無窮小 當x 0時,a x 趨向於0,則說明a x 是無窮小。無窮下一般用a表示,高階無窮小用o表示。微分 0y ax o 0x 其中o 0x 這個符號,說明o是0x的高階無窮小。在極限的運算過程中,f x 在x0處的極限為a,則f x a 多元函式的概念 z f x,y 極限 從各個方向向 x0,...