簡化您的工作流程:搜尋 miniwebtool。
添加插件
相關工具
伯努利微分方程式求解器歐拉方法計算機一階常微分方程求解器感抗計算機運動學計算機龍格-庫塔RK4方法計算機二階常微分方程求解器
主頁 > 數學 > 微積分 > 常微分方程組求解器
 

常微分方程組求解器

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

常微分方程組求解器
求解 \(\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)\)。可以使用 sincosexpsqrt 等函數及常數 pie
不計算閉合形式解;而是通過高精度 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轉換器🎮 遊戲靈敏度轉換器百分比減少計算機質數分解計算機毛利率計算機百分比折扣計算機MAC地址查找HEX計算機羅馬數字轉換器圓計算機隨機字母生成器隨機餐點產生器kg到lbs轉換器Instagram用戶ID查詢磅轉公斤轉換器校正鈣計算機圖片分割器年齡計算機隨機名稱生成器查找並替換文字坡度與傾斜度計算機ANC計算機加價計算機年份天數計算機 - 今天是今年的第幾天影片轉圖片擷取器反向文字百分比變化計算機調整影片速度最小公倍數計算機太陽、月亮與上升星座計算機 🌞🌙✨Bar to PSI 轉換器定期存款計算機凱薩密碼工具分數百分比轉換器CAGR計算機隨機日期生成器商和餘數計算機SRT時間偏移因子計算機AI內容檢測器簡單利息計算機樂透號碼生成器📅 日期計算機百分比誤差計算機隨機選擇器對數計算機合併影片剪刀石頭布產生器可整除測試計算機音訊分割器平均值計算機棒球打擊率計算機標準偏差計算機 - 高精度psi到kpa轉換器分數到小數計算機二次公式計算機影片壓縮器🌡️ 體感溫度計算機隨機數學題產生器賓果卡生成器YouTube頻道統計Facebook用戶ID查詢直角三角形計算機畢達哥拉斯定理計算機複利計算機輪胎尺寸計算機小字體生成器 ⁽ᶜᵒᵖʸ ⁿ ᵖᵃˢᵗᵉ⁾FPS 轉換器ERA計算機橢圓 周長計算機文字重複工具小數到分數計算機SRT轉換為TXT工具公因子計算機比率與百分比計算機最大公因子計算機壓力轉換器心算訓練器⏱️ 小時計算機組合計算機隨機英文單字產生器百分比到ppm轉換器跑步配速計算機迷宮產生器考拉茲猜想計算機HEX轉換器純利潤計算機我的幸運數字是什麼速度計算機OPS計算機填字遊戲製作器愛情兼容性計算機每局被上壘率 WHIP 計算機總和計算機MAC地址產生器斜率計算機摩爾斯電碼生成器數獨產生器與求解器刪除線文字產生器質數檢查器保齡球計分計算機常規時間到小數時間轉換器比率計算機隨機顏色生成器GUID / UUID 生成器cpm計算機乘積符號計算機 (Pi 記號)積分計算機發音音標轉換器圖片打碼工具汽車貸款計算機條碼產生器兩點間距離計算機文件大小轉換器科學記數法計算機上壘率計算機土星回歸計算機自行車齒輪比計算機血壓解讀器投球命中率計算機行數統計工具影片分割器線規計算機工時計算機汽車折舊計算機星期幾計算機演講時間計算機鋼筋計算機克到磅轉換器階乘計算機隨機錦標賽對陣生成器樓梯計算機騎行速度計算機SRT合併工具速度轉換器隨機字符串生成器AI標點符號添加器GIF 轉 MP4 轉換器分數簡化磅到克轉換器步數距離計算機弧長計算機百分比計算機PSI 轉 Bar 轉換器燃油費用計算機獲取字符串長度日曆隨機化數字預期壽命計算機高斯分布產生器厘米到英尺和英寸轉換器線性迴歸計算機💧 露點計算機YouTube收益估算器體型計算機數字轉分數轉換器旋轉影片中位數計算機慣性矩計算機斜率截距式計算機校正鈉計算機棒球長打率計算機翻轉圖片體積轉換器ppm到百分比轉換器半衰期計算機圖片壓縮器圖片旋轉器橢圓面積計算機為影片新增浮水印角度轉換器YouTube Shorts 收益計算機姓名靈數計算機🌐 時區轉換器WHtR計算機排序數字進製轉換器鋼材重量計算機隨機生日生成器債務資產比率計算機ASCII碼表Base64解碼器MP3循環播放器刪除空行變異係數計算機隨機時間生成器隱形文本生成器YouTube縮圖下載器比較分數計算機凱利公式計算機自行車尺寸計算機速通分段計時器斯特伯福計分計算機飛鏢殘局計算機板球投球經濟率計算機擊球率計算機淨得分率計算機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角色生成器