| 微分方程 |
| --- |
| Navier–Stokes differential equations used to simulate airflow around an obstruction模擬有阻礙下空氣流動的纳维-斯托克斯方程 |
| 領域

| * 自然科学
* 工程学
 |
| --- |
| * 天文學
* 物理学
* 化學
* 生物学
* 地质学
 |
| 应用数学 |
| * 连续介质力学
* 混沌理论
* 动力系统
 |
| 社会科学 |
| * 经济学
* 人口動力學（英语：）
 |

 |
| 
分類 |
| 種類

|  |
| --- |
| 
* 常微分
* 偏微分
* 微分代數（英语：）
* 積分-微分（英语：）
* 分数
* 线性
* 非線性

 |
| 
依變數種類 |
| * 自变量和因变量

* 自治
* 耦合 / 解耦合
* 全微分
* 齊次（英语：） / 非齊次

 |
| 
特徵 |
| 
* 階數
* 微分算子

* 表示法（英语：）
 |

 |
| 和過程的關係
* 差分 (離散下的類比)

* 隨機
	+ 随机偏微分
* 时滞

 |
| 
解 |
| 存在和唯一
* 柯西-利普希茨定理
* 皮亚诺存在性定理
* Carathéodory存在性定理（英语：）
* Cauchy–Kowalevski定理（英语：）
 |
| 通用主題
* 朗斯基行列式
* 相圖
* 相空間
* 李雅普诺夫 / 漸進 / 指數穩定
* 收斂速度
* 級數 / 積分解
* 數值積分
* 狄拉克δ函数
 |
| 解法
* 特征线
* 欧拉
* 指數響應公式（英语：）
* 有限差分 (克兰克－尼科尔森)
* 有限单元
	+ 無限單元
* 有限體積
* 伽辽金
	+ Petrov–Galerkin（英语：）
* 积分因子
* 积分变换
* 摄动
* 龙格－库塔

* 分離變數
* 待定係數
* 參數變換
 |
| 人物
* 艾萨克·牛顿
* 戈特弗里德·莱布尼茨
* 萊昂哈德·歐拉
* 埃米尔·皮卡
* Józef Maria Hoene-Wroński（英语：）
* Ernst Lindelöf（英语：）
* 鲁道夫·利普希茨
* 奧古斯丁-路易·柯西
* 約翰·克蘭克
* 菲利斯·尼科爾森
* 卡爾·龍格
* 馬丁·威廉·庫塔

 |
| * 
* 论
* 编
 |

在数学和计算机科学中，**欧拉方法**（英語：），是一种一阶数值方法，用以对给定初值的常微分方程求解。

欧拉方法是常微分方程數值方法（英语：）中最基本的显式方法；也是一个一阶方法，意味着其局部截断误差正比于步长的平方，并且其全局截断误差正比于步长。

非正式的几何描述
--------

欧拉方法的图示。待求的曲线为蓝色，它的折線近似为红色。
考虑计算這樣的一个未知曲線的形状：它具有给定的起点并且满足一个给定的微分方程。 这里，所谓"微分方程"可以看作能够通过曲线上任意点的位置而计算出这一点的切线斜率的公式。

思路是，一开始只知道曲線的起点（假设为A\_{0}），曲線其他部份是未知的，不過通过微分方程，A\_{0}的斜率可以被计算出来，也就得到了切线。

顺着切线向前走一小步到点A\_{1}。如果假设A\_{1}是曲线上的一点（实际上通常不是），那么同样的道理就可以确定下一条切线，依此类推。在经过几步之后，一条折线（英语：）A\_{0}A\_{1}A\_{2}A\_{3}\dots 就被计算出来了。一般情况下，这条折线与原先的未知曲线偏离不远，并且任意小的误差都可以通过减少步长来得到。

欧拉方法的推导
-------

图示为方程{\displaystyle y'=y,y(0)=1}的数值积分。蓝色为欧拉法，绿色为中点法（英语：），红色为精确解{\displaystyle y=e^{t}}。所用步长為h=1.0。
图示为同一个方程在步长h=0.25时的结果。可以看出中点法比欧拉法收敛更快。
以以下微分方程為例

{\displaystyle y'(t)=f(t,y(t)),\qquad y(t\_{0})=y\_{0},}
希望用 *y* 在點 (t0,y(t0)) 附近的線性近似來得到其近似解（也就是 *y* 的泰勒展開式的前二項）。利用時間 *t*n 時的數值，若用單步的欧拉方法，可得到時間 *t*n+1 = *t*n + *h* 時的近似值如下：

y\_{{n+1}}=y\_{n}+hf(t\_{n},y\_{n}).\qquad \qquad 
欧拉方法是一種顯型方法，也就是說 y\_{n+1} 的解是 y\_{i}, i\leq n 的顯函數。

欧拉方法可以求解一階的微分方程，而任何N階的微分方程都可以表示成一階的微分方程。

對於微分方程

{\displaystyle y^{(N)}(t)=f(t,y(t),y'(t),\ldots ,y^{(N-1)}(t))}
可以通過新設輔助變量 {\displaystyle z\_{1}(t)=y(t),z\_{2}(t)=y'(t),\ldots ,z\_{N}(t)=y^{(N-1)}(t)}，得到以下的等價方程

{\displaystyle \mathbf {z} '(t)={\begin{pmatrix}z\_{1}'(t)\\\vdots \\z\_{N-1}'(t)\\z\_{N}'(t)\end{pmatrix}}={\begin{pmatrix}y'(t)\\\vdots \\y^{(N-1)}(t)\\y^{(N)}(t)\end{pmatrix}}={\begin{pmatrix}z\_{2}(t)\\\vdots \\z\_{N}(t)\\f(t,z\_{1}(t),\ldots ,z\_{N}(t))\end{pmatrix}}}
這是一個以{\displaystyle \mathbf {z} (t)}為變量的一階系統，因此可以用歐拉法求解，也可以使用其他的一階數值方法。

应用例题
----

设微分方程为 y'=y ，初始值为 y(0)=1，试用欧拉方法求 y\_{3}的近似值，步长为 h=1。

欧拉法為：

{\displaystyle y\_{n+1}=y\_{n}+hf(t\_{n},y\_{n}).}
首先求f(t\_{0},y\_{0})（当n = 0），f的定義為{\displaystyle f(t,y)=y}，因此有

{\displaystyle f(t\_{0},y\_{0})=f(0,1)=1.}
透過以上步驟，求得解曲線在点(0,1)的切线斜率。回顾直線斜率的定义：y变化量和t变化量的比值，亦記作{\displaystyle \Delta y/\Delta t}。 

接著是

{\displaystyle h\cdot f(y\_{0})=1\cdot 1=1.}
{\displaystyle y\_{0}+hf(y\_{0})=y\_{1}=1+1\cdot 1=2.}
重复以上步骤求出y\_{2} 和 y\_{3}的值。

{\displaystyle y\_{2}=y\_{1}+hf(y\_{1})=2+1\cdot 2=4}
{\displaystyle y\_{3}=y\_{2}+hf(y\_{2})=4+1\cdot 4=8}
由于欧拉法属于递归算法，把運算整理成表格也許有助於避免計算錯誤。

| y\_n | t\_{n} | y'(t) | h | dy | y\_{n+1} |
| --- | --- | --- | --- | --- | --- |
| 1 | 0 | 1 | 1 | 1 | 2 |
| 2 | 1 | 2 | 1 | 2 | 4 |
| 4 | 2 | 4 | 1 | 4 | 8 |

局部截尾误差
------

欧拉法的局部截尾误差（Local truncation error, LTE）是指在实施一次欧拉法所产生的误差，是指经过一步的数值解y\_{1}与在t\_{1}=t\_{0}+h时精确解的误差。数值解y\_{1}由以下给出：

y\_{1}=y\_{0}+hf(t\_{0},y\_{0}).\quad 
对于精确解，使用泰勒级数展开给出：

y(t\_{0}+h)=y(t\_{0})+hy'(t\_{0})+{\frac {1}{2}}h^{2}y''(t\_{0})+O(h^{3}).
欧拉法的局部截尾误差为：

{\mathrm {LTE}}=y(t\_{0}+h)-y\_{1}={\frac {1}{2}}h^{2}y''(t\_{0})+O(h^{3}).
当y拥有三阶有界导数时，这个结果是成立的。

结果显示：当步长h很小时，局部截尾误差近似与 h^{2} 成比例。也就是说，欧拉法没有其他的高阶方法如龙格－库塔法和线性多步法（英语：）精确，这些方法的局部截尾误差与h^{p}（*p*>2）成比例。

全局截尾误差
------

全局截尾误差（Global truncation error, GTE）是指在一个固定时间t时的误差，但是很多步之后该方法需要以从初始时间到达该时间来计算。全局截尾误差可以看做是一个每一步的局部截尾误差的累积效应。 经过的步骤數為(t-t\_{0})/h，而每步的误差则正比于h^{2}。因此，可以预期全局截尾误差是正比于h的。

这个直观的推测可以被嚴謹地證明。如果解y存在二阶有界导数，并且f關於y是利普希茨连续的，那么全局截尾误差是有界的：

|{\text{GTE}}|\leq {\frac {hM}{2L}}(e^{{L(t-t\_{0})}}-1)\qquad \qquad 
其中 M 是在给定区间内y的二阶导数的上界，L 是f的利普希茨常数。

这种精确的形式其实是没有什么意义的，通常情况下这个上界都會嚴重高估了欧拉法所造成的实际误差。重要的是，这顯示了全局截尾误差是近似正比于h的，所以欧拉法被稱为是一阶的。

相關條目
----

* 克兰克－尼科尔森方法
* 梯度下降法，也是用有限步進行，計算函數的最小值。
* 龍格－庫塔法列表（英语：）
* 線性多步法（英语：）
* 數值積分（計算定積分）
* 常微分方程數值方法（英语：）

註腳
--

1. ↑  命名自它的发明者萊昂哈德·歐拉
2. ↑  即初值問題
3. ↑  即每步误差
4. ↑  欧拉方法经常应用于作为构建一些更复杂方法的基础，例如预估-校正方法。

参考资料
----

1. ↑  Butcher 2003，第3頁；Hairer, Nørsett & Wanner 1993，第2頁
2. ↑  Butcher 2003，第60頁
3. ↑  Atkinson 1989，第344頁
4. ↑  Butcher 2003，第49頁
5. ↑  Atkinson 1989，第346頁；Lakoba 2012，公式 (1.16)
6. ↑  Iserles 1996，第7頁
7. ↑  Butcher 2003，第63頁

参考文献
----

* Atkinson, Kendall A.  2nd, New York: John Wiley & Sons, 1989, ISBN 978-0-471-50023-0 .
* Ascher, Uri M.; Petzold, Linda R. Philadelphia: Society for Industrial and Applied Mathematics, 1998, ISBN 978-0-89871-412-8 .
* Butcher, John C. New York: John Wiley & Sons, 2003, ISBN 978-0-471-96758-3 .
* Hairer, Ernst; Nørsett, Syvert Paul; Wanner, Gerhard, Berlin, New York: Springer-Verlag, 1993, ISBN 978-3-540-56670-0 .
* Iserles, Arieh, Cambridge University Press, 1996, ISBN 978-0-521-55655-2 .
* Lakoba, Taras I.  (PDF) (Lecture notes for MATH334, University of Vermont), 2012  [2016-01-02], （原始内容存档 (PDF)于2012-07-12） .

|  |  |  |  |  |  |
| --- | --- | --- | --- | --- | --- |
| 

| * 
* 论
* 编
莱昂哈德·欧拉 |
| --- |
|  |
| 
* 歐拉-拉格朗日方程
* 歐拉-洛特卡方程（英语：）
* 欧拉-麦克劳林求和公式
* 欧拉-丸山法
* 歐拉-馬斯刻若尼常數
* 歐拉-泊松-達布方程（英语：）
* 歐拉-羅德里格斯公式（英语：）
* 欧拉-特里科米方程
* 歐拉連分數公式（英语：）
* 欧拉临界负载
* 欧拉公式
* 欧拉四平方和恒等式
* 欧拉恒等式
* 歐拉泵和渦輪方程（英语：）
* 欧拉旋转定理
* 欧拉猜想
* 欧拉定理 (数论)
* 欧拉方程 (流体动力学)
* 歐拉函數 (複變函數)
* 欧拉方法
* 欧拉数
* 欧拉数 (物理学)
* 欧拉-伯努力栋梁方程
* 以他命名的事物（英语：）

 |
|  |
| 
* 分类 **分類**

 |

 |