簡化您的工作流程：搜尋 miniwebtool。
添加插件
相關工具
絕對值不等式求解器伯努利微分方程式求解器特徵值和特徵向量計算機一階常微分方程求解器矩陣對角化計算機相關變化率求解器龍格-庫塔RK4方法計算機二階常微分方程求解器線性方程組求解器非線性方程組求解器
主頁 > 數學 > 微積分 > 常微分方程組求解器
 

常微分方程組求解器

以符號和數值方式求解常微分方程組 x' = Ax。自動對平衡點進行分類（鞍點、節點、螺旋點、中心點），逐步推導特徵值和特徵向量，寫出閉式通解和特解，並繪製帶有動態軌跡的互動式相圖 — 適用於 2×2、3×3 線性以及非線性 2D 系統。

常微分方程組求解器
快速範例：
穩定節點 中心點（振盪器） 穩定螺旋點 鞍點 重根（退化） 單擺 Van der Pol 振盪器 3×3 三階系統
求解 \(\dfrac{d}{dt}\!\begin{pmatrix}x\\y\end{pmatrix} = A \begin{pmatrix}x\\y\end{pmatrix}\)，其中 \(A\) 是具有實數項的 2×2 矩陣。
[
[
支援小數（例如 2.5）、負數（-1）和簡單分數（3/4）。
求解 \(\dfrac{d}{dt}\!\begin{pmatrix}x_1\\x_2\\x_3\end{pmatrix} = A \begin{pmatrix}x_1\\x_2\\x_3\end{pmatrix}\) — 獲取特徵值、特徵向量和閉合形式通解。
[
[
[
預設值為 \( s^3 + 6s^2 + 11s + 6 \) 的伴隨矩陣，特徵值為 \(-1, -2, -3\)。
數值求解一般 2D 非線性系統：\(x' = f(x, y)\)，\(y' = g(x, y)\)。可以使用 sin、cos、exp、sqrt 等函數及常數 pi、e。
不計算閉合形式解；而是通過高精度 RK4 積分產生軌跡和相圖。
x(0)
y(0)
時間跨度 T

Embed 常微分方程組求解器 Widget

常微分方程組求解器

這款常微分方程組求解器是一個多功能的微分方程工具箱，適用於耦合的線性與非線性系統。只需輸入 2×2 或 3×3 係數矩陣，該工具即可執行完整的特徵值/特徵向量分析，以 LaTeX 格式寫出閉合形式的通解和特解，將原點平衡點分類為鞍點、節點、螺旋點或中心點，並繪製帶有動畫軌跡的互動式相圖。對於非線性平面系統，您可以輸入任意右側項 \(f(x,y)\) 和 \(g(x,y)\)，該計算機將產生高精度的 RK4 相圖。

什麼是常微分方程組？

常微分方程組通過其導數將多個單變量（通常是時間 \(t\)）的未知函數耦合在一起。以最簡潔的形式表示為：

\[ \mathbf{x}'(t) = \mathbf{F}(t, \mathbf{x}(t)) , \qquad \mathbf{x} \in \mathbb{R}^n \]

當 \(\mathbf{F}(t, \mathbf{x}) = A\mathbf{x}\) 且 \(A\) 為常數矩陣時，該系統稱為線性且自治的。這正是理論最優美之處：整個長期行為完全由矩陣 \(A\) 的特徵值決定。

線性系統的特徵值求解步驟

對於 \(\mathbf{x}' = A\mathbf{x}\)，標準方法如下：

  1. 計算特徵多項式 \(\det(\lambda I - A) = 0\)。
  2. 求解特徵值 \(\lambda_1, \lambda_2, \dots\)。
  3. 對於每個特徵值，通過求解 \((A - \lambda I) v = 0\) 找到特徵向量 \(v\)。
  4. 將通解組合成線性組合：\(\mathbf{x}(t) = c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t} + \cdots\)。
  5. 將初始條件 \(\mathbf{x}(0)\) 代入通解以確定常數 \(c_i\)。

2×2 系統的三種情況

特徵值 通解 相圖特徵
實數相異 \(\lambda_1 \ne \lambda_2\) \(c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t}\) 正負號不同為鞍點；否則為節點
共軛複數 \(\alpha \pm i\beta\) \(e^{\alpha t}[c_1(p\cos\beta t - q\sin\beta t) + c_2(p\sin\beta t + q\cos\beta t)]\) 螺旋點 (\(\alpha \ne 0\)) 或中心點 (\(\alpha = 0\))
重根 \(\lambda_1 = \lambda_2 = \lambda\) \(c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t}\) 退化節點

跡-行列式平面 (Trace-Determinant Plane)

對於一個跡為 \(T = a_{11} + a_{22}\) 且行列式為 \(D = a_{11} a_{22} - a_{12} a_{21}\) 的 2×2 矩陣，整個分類可以歸納在一個圖表中：

\[ \begin{array}{l} D < 0 \Rightarrow \text{鞍點} \\ D > 0, \; T^2 - 4D > 0 \Rightarrow \text{節點 (穩定，若 } T < 0\text{)} \\ D > 0, \; T^2 - 4D < 0 \Rightarrow \text{螺旋點 (穩定，若 } T < 0\text{), 中心點，若 } T = 0 \\ T^2 - 4D = 0 \Rightarrow \text{退化節點 (重根特徵值)} \end{array} \]

這就是為什麼結果面板會突出顯示 \(T\)、\(D\) 和 \(\Delta = T^2 - 4D\) — 這三個數字足以確定平衡點的名稱。

非線性系統與相圖

大多數現實世界的 ODE 都是非線性的，沒有閉合形式解。該工具通過四階 Runge–Kutta (RK4) 方法對方程進行數值積分，該方法具有 \(O(h^5)\) 的局部截斷誤差，是處理平滑向量場的標準方法。

\[ \begin{aligned} k_1 &= f(t_n, y_n) \\ k_2 &= f(t_n + h/2, y_n + h k_1 / 2) \\ k_3 &= f(t_n + h/2, y_n + h k_2 / 2) \\ k_4 &= f(t_n + h, y_n + h k_3) \\ y_{n+1} &= y_n + \tfrac{h}{6}(k_1 + 2k_2 + 2k_3 + k_4) \end{aligned} \]

相圖疊加了：

如何使用此計算機

  1. 透過表單頂部的標籤選擇模式 — 2×2 線性、3×3 線性或 2D 非線性。
  2. 填寫係數或方程。點擊任何「快速範例」以預填典型系統（穩定節點、中心點、鞍點、單擺、Van der Pol 振盪器等）。
  3. 輸入初始條件 \((x_0, y_0)\) 和時間跨度 \(T\)。對於振盪器，典型的 \(T\) 值為 6–20，對於快速衰減的穩定系統，典型的 \(T\) 值為 3–6。
  4. 點擊求解。完整的結果頁面將會出現，包含分類、特徵值、特徵向量、閉合形式解（線性模式）、動畫相圖和時間序列圖。
  5. 如果您想再次觀看運行器遍歷 IC 曲線，可以使用相圖下方的按鈕重新播放軌跡。

實例分析 — 阻尼諧振子

阻尼諧振子 \(\ddot{x} + 2\zeta \omega \dot{x} + \omega^2 x = 0\) 可以通過令 \(y = \dot{x}\) 重寫為 2D 系統：

\[ \begin{pmatrix} \dot{x} \\ \dot{y} \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -\omega^2 & -2\zeta\omega \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} \]

對於 \(\omega = 1\) 且 \(\zeta = 0.2\)（欠阻尼）的情況，矩陣為 \(A = \begin{pmatrix} 0 & 1 \\ -1 & -0.4 \end{pmatrix}\)。跡 \(T = -0.4\)，行列式 \(D = 1\)，判別式 \(\Delta = 0.16 - 4 = -3.84 < 0\)，因此我們得到一個穩定螺旋點，特徵值為 \(-0.2 \pm 0.9798\,i\)。軌跡向原點螺旋收斂，時間序列顯示指數衰減的正弦曲線。

應用領域

提示與注意事項

常見問題解答

什麼是常微分方程組？

常微分方程組 (ODEs) 是一組耦合方程，它們將多個未知函數的導數與單個自變量（通常是時間）聯繫起來。經典形式是 \( \mathbf{x}'(t) = F(t, \mathbf{x}(t)) \)，其中 \( \mathbf{x} \) 是狀態向量，\(F\) 是向量場。線性系統可以簡潔地寫為 \( \mathbf{x}' = A\mathbf{x} + \mathbf{b} \)，其行為幾乎完全由係數矩陣 \(A\) 的特徵值決定。

特徵值如何對 2×2 線性系統的平衡點進行分類？

對於 2×2 系統 \( \mathbf{x}' = A\mathbf{x} \)，原點由 \(A\) 的跡 \(T\) 和行列式 \(D\) 分類：\(D < 0\) 產生鞍點（不穩定）；\(D > 0\) 且 \(T^2 > 4D\) 產生節點（\(T < 0\) 時穩定，\(T > 0\) 時不穩定）；\(D > 0\) 且 \(T^2 < 4D\) 產生螺旋點（\(T < 0\) 時穩定，\(T > 0\) 時不穩定，\(T = 0\) 時為純中心點）。臨界情況 \(T^2 = 4D\) 產生退化節點。

當特徵值為複數時，閉合形式解長什麼樣？

如果 \(A\) 具有共軛複數特徵值 \( \alpha \pm i\beta \)，其複數特徵向量為 \( v = p + iq \)，則實數通解為 \( \mathbf{x}(t) = e^{\alpha t} \left[ c_1 (p \cos\beta t - q \sin\beta t) + c_2 (p \sin\beta t + q \cos\beta t) \right] \)。指數 \(e^{\alpha t}\) 控制振幅（增長、衰減或常數），而正弦和餘弦處理旋轉部分。

當矩陣有重根特徵值時會發生什麼？

如果矩陣具有重根特徵值 \(\lambda\) 但只有一個線性無關的特徵向量 \(v\)，則還需要一個滿足 \( (A - \lambda I) w = v \) 的廣義特徵向量 \(w\)。通解形式為 \( \mathbf{x}(t) = c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t} \)。如果特徵空間恰好是二維的，則矩陣在該不變子空間上是單位矩陣的純量倍數，解簡化為 \( \mathbf{x}(t) = (c_1 v_1 + c_2 v_2) e^{\lambda t} \)。

此工具可以符號化求解非線性系統嗎？

非線性模式使用四階 Runge–Kutta (RK4) 積分器進行數值求解並繪製相圖。大多數非線性系統沒有閉合形式解，因此這是標準做法。您仍可以通過線性化讀取平衡點附近的局部行為（2×2 線性模式可處理此部分）— 只需計算固定點處的 Jacobian 矩陣並將其作為 \(A\) 輸入即可。

什麼是相圖？

相圖是 \(x\)–\(y\) 平面中 2D 系統解的幾何圖形。每個解描繪出一條稱為軌跡的曲線，軌跡集合與向量場箭頭共同揭示了定性行為：解是向內螺旋、向外鞍開、持續振盪還是趨於平衡點。相圖讓系統的全局結構一目了然。

延伸閱讀

引用此內容、頁面或工具為：

"常微分方程組求解器" 於 https://MiniWebtool.com/zh-tw/常微分方程組求解器/，來自 MiniWebtool，https://MiniWebtool.com/

由 miniwebtool 團隊開發。更新日期：2026年4月23日

您還可以嘗試我們的 AI數學解題器 GPT，通過自然語言問答解決您的數學問題。

微積分:

常用工具:

隨機撲克牌產生器分數計算機隨機信用卡生成器真心話大冒險產生器比例計算機相對標準偏差計算機百分比增加計算機斜邊計算機最簡分數計算機🎮 遊戲靈敏度轉換器百分比增長率計算機kpa到psi轉換器百分比減少計算機毛利率計算機質數分解計算機羅馬數字轉換器百分比折扣計算機HEX計算機MAC地址查找圓計算機kg到lbs轉換器Instagram用戶ID查詢磅轉公斤轉換器隨機字母生成器校正鈣計算機圖片分割器年齡計算機查找並替換文字隨機名稱生成器隨機餐點產生器坡度與傾斜度計算機ANC計算機加價計算機調整影片速度最小公倍數計算機太陽、月亮與上升星座計算機 🌞🌙✨年份天數計算機 - 今天是今年的第幾天影片轉圖片擷取器凱薩密碼工具商和餘數計算機反向文字百分比變化計算機分數百分比轉換器CAGR計算機Bar to PSI 轉換器簡單利息計算機AI內容檢測器定期存款計算機SRT時間偏移樂透號碼生成器因子計算機隨機日期生成器📅 日期計算機百分比誤差計算機分數到小數計算機🌡️ 體感溫度計算機對數計算機可整除測試計算機隨機選擇器合併影片音訊分割器剪刀石頭布產生器棒球打擊率計算機標準偏差計算機 - 高精度影片壓縮器賓果卡生成器隨機數學題產生器psi到kpa轉換器平均值計算機Facebook用戶ID查詢FPS 轉換器畢達哥拉斯定理計算機YouTube頻道統計小字體生成器 ⁽ᶜᵒᵖʸ ⁿ ᵖᵃˢᵗᵉ⁾直角三角形計算機二次公式計算機複利計算機橢圓 周長計算機SRT轉換為TXT工具文字重複工具ERA計算機壓力轉換器小數到分數計算機⏱️ 小時計算機比率與百分比計算機總和計算機考拉茲猜想計算機OPS計算機百分比到ppm轉換器隨機英文單字產生器每局被上壘率 WHIP 計算機輪胎尺寸計算機純利潤計算機MAC地址產生器心算訓練器填字遊戲製作器最大公因子計算機跑步配速計算機HEX轉換器組合計算機公因子計算機迷宮產生器隨機錦標賽對陣生成器兩點間距離計算機我的幸運數字是什麼摩爾斯電碼生成器比率計算機速度計算機斜率計算機克到磅轉換器積分計算機保齡球計分計算機質數檢查器數獨產生器與求解器GUID / UUID 生成器刪除線文字產生器愛情兼容性計算機文件大小轉換器條碼產生器cpm計算機常規時間到小數時間轉換器圖片打碼工具土星回歸計算機汽車折舊計算機演講時間計算機騎行速度計算機汽車貸款計算機影片分割器血壓解讀器磅到克轉換器星期幾計算機自行車齒輪比計算機上壘率計算機步數距離計算機投球命中率計算機分數簡化鋼筋計算機隨機字符串生成器線規計算機中位數計算機科學記數法計算機GIF 轉 MP4 轉換器百分比計算機速度轉換器樓梯計算機行數統計工具發音音標轉換器PSI 轉 Bar 轉換器獲取字符串長度隨機顏色生成器高斯分布產生器日曆排序數字翻轉圖片階乘計算機💧 露點計算機體型計算機AI標點符號添加器弧長計算機厘米到英尺和英寸轉換器碳水化合物計算機YouTube收益估算器多分數計算機YouTube Shorts 收益計算機淨利潤率計算機為影片新增浮水印燃油費用計算機旋轉影片橢圓面積計算機進製轉換器半衰期計算機慣性矩計算機數字轉分數轉換器變異係數計算機預期壽命計算機工時計算機校正鈉計算機隨機化數字SRT合併工具圖片壓縮器角度轉換器體積轉換器ppm到百分比轉換器棒球長打率計算機線性迴歸計算機鋼材重量計算機債務資產比率計算機🌐 時區轉換器比較分數計算機移除標點符號線上工具ASCII碼表MP3循環播放器圖片旋轉器WHtR計算機刪除空行現值計算機隨機時間生成器骰子機率計算機分壓器計算機複數計算機擊球率計算機淨得分率計算機K/D 比率計算機Elo 等級分計算機心率恢復計算機健行時間計算機划船配速計算機跑步年齡分級計算機FTP功率區間計算機Beep 測試計算機ACFT 分數計算機Wilks 與 DOTS 計算機輪圈偏移量計算機輪胎負載指數與速度等級查詢每英里成本計算機租賃買斷計算機辛烷值混合計算機二衝程機油混合計算機引擎排氣量計算機沙發過門計算機木柴考得計算機空氣清淨機cadr計算機除濕機尺寸計算機吊扇尺寸計算機窗簾尺寸計算機地毯尺寸計算機掛畫高度計算機電視安裝高度計算機電視尺寸計算機池塘容積與襯墊計算機泳池鹽量計算機泳池容積計算機熱水器尺寸計算機環氧樹脂計算機欄杆立柱間距計算機踢腳板與飾條計算機壁板計算機露臺木器染色劑計算器草籽計算機草皮計算機瀝青計算機立方碼計算機天線長度計算機電線管填充率計算機串聯並聯電容計算機感抗計算機房間照明計算機勒克斯轉流明計算機流明轉瓦特轉換器發電機功率計算機mah到wh轉換器三相功率計算機kVA計算機安培轉瓦特計算機瓦特轉安培計算機串聯電阻計算機摩擦力計算機斜面計算機機械利益計算機聲速計算機波速計算機浮力計算機終端速度計算機德布羅意波長計算機光子能量計算機emc2計算機時間膨脹計算機克卜勒第三定律計算機逃逸速度計算機萬有引力計算機比爾-朗伯定律計算機能斯特方程式計算機滲透壓計算機沸點上升計算機凝固點下降計算機百分比組成計算機當量濃度計算機重量莫耳濃度計算機pKa到Ka轉換器亨德森-哈塞爾巴爾赫計算機理論產量計算機限量試劑計算機電子組態計算機互動式元素週期表AI教案產生器AI測驗產生器引用產生器 (APA/MLA/Chicago)出席率計算機AP分數計算機ACT分數計算機SAT分數計算機百分比轉CGPA轉換器CGPA轉百分比轉換器簡易評分器 EZ Grader養育孩子成本計算機寶寶奶量計算機尿布尺寸計算機嬰兒名字產生器寶寶眼睛顏色預測器兒童BMI百分位計算機兒童身高預測器hcg翻倍時間計算機試管嬰兒預產期計算機著床計算機中國生男生女預測器ISO 8601 日期格式化工具儒略日轉換器午睡計算機月相計算機日出日落計算機世界時鐘日期至羅馬數字轉換器退休倒數計時戒癮天數計算機半歲生日計算機紀念日計算機週數計算機小費分配計算機電子郵件行銷ROI計算機每筆潛在客戶成本計算機營運資本計算機隨機RPG角色生成器