简化您的工作流程:搜索 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)\)。支持函数如 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}\) 退化节点

迹-行列式平面

对于一个迹为 \(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 \text{ 时为中心点} \\ T^2 - 4D = 0 \Rightarrow \text{退化节点 (重特征值)} \end{array} \]

这就是为什么结果面板显著显示 \(T\)、\(D\) 和 \(\Delta = T^2 - 4D\)——这三个数字足以命名平衡点。

非线性系统与相图

大多数现实世界的 ODEs 是非线性的,没有闭式解。本工具通过四阶龙格-库塔 (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,快速衰减的稳定系统为 3–6。
  4. 点击求解。将出现包含分类、特征值、特征向量、闭式解(线性模式)、动画相图和时间序列图的完整结果页面。
  5. 使用相图下方的按钮重播轨迹,如果您想再次观察游标穿过轨迹的过程。

计算示例 — 阻尼谐振子

阻尼振子 \(\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} \)。

该工具可以符号化求解非线性系统吗?

非线性模式使用四阶龙格-库塔 (RK4) 积分器对系统进行数值求解并绘制相图。大多数非线性系统没有闭式解,因此这是标准方法。您仍然可以通过线性化来读取平衡点附近的局部行为,这可以通过 2×2 线性模式处理 — 在不动点处计算雅可比矩阵并将其作为 \(A\) 输入即可。

什么是相图?

相图是 2D 系统在 \(x\)–\(y\) 平面上的解的几何图形。每个解描绘出一条称为轨迹的曲线,轨迹集合与向量场箭头一起揭示了定性行为:解是向内螺旋、向外发散、振荡还是趋于平衡点。相图使系统的全局结构一目了然。

延伸阅读

引用此内容、页面或工具为:

"常微分方程组求解器" 于 https://MiniWebtool.com/zh-cn/常微分方程组求解器/,来自 MiniWebtool,https://MiniWebtool.com/

由 miniwebtool 团队开发。更新日期:2026年4月23日

您还可以尝试我们的 AI数学解题器 GPT,通过自然语言问答解决您的数学问题。

微积分:

常用工具:

随机信用卡生成器MAC地址查找彩票号码生成器凯利公式计算器网址提取器英尺英寸转换为厘米音频提取器太阳、月亮与上升星座计算器 🌞🌙✨CAGR计算器SRT转为TXT工具🌡️ 体感温度计算器相对标准偏差计算器随机扑克牌生成器MAC地址生成器样本量计算器厘米到英尺和英寸转换器cpm计算器图片打码工具随机IMEI生成器百分比折扣计算器考拉兹猜想计算器FPS 转换器🎮 游戏灵敏度转换器调整视频速度定期存款计算器音频分割器kg到lbs转换器HEX计算器VAT计算器对数计算器SHA256 哈希生成器毛利率计算器质数检查器t检验计算器📅 日期计算器Facebook用户ID查询Log Base 10 计算器血糖转换器农历转换器比例计算器斜边计算器随机字符串生成器One Rep Max (1RM) 计算器太阳位置计算器职位查找器随机选择器为图片添加文字真心话大冒险生成器随机字母生成器密码强度测试器百分比变化计算器卡方检验计算器图片压缩器随机化数字DOY日历根式化简器椭圆周长计算器英寸到厘米转换器随机装备生成器标准偏差计算器 - 高精度磅转千克转换器颜料调色计算器因子计算器年度天数计算器 - 今天是今年的第几天罗马数字转换器分数计算器六西格玛过程能力计算器半衰期计算器🖱️ 点击计数器三角函数绘图器位数计算器余弦定理计算器变异系数计算器泰勒级数计算器百分比到ppm转换器花样字体生成器视频转图片提取器💧 露点计算器螺栓扭矩计算器srt时间偏移两点间距离计算器十六进制到十进制转换器⏱️ 小时计算器新月和满月日历PSI 转 Bar 转换器合并视频抛硬币年龄计算器最简分数计算器积分计算器MD5哈希生成器分贝 (dB) 计算器弧长计算器条形码生成器米到英尺转换器利率计算器配色方案生成器移除标点符号在线工具多项式展开计算器日出日落计算器利润计算器图片分割器圆柱体体积计算器 高精度填字游戏制作器极坐标方程绘图器翻转视频视频分割器📷 OCR / 图片文字识别双重积分计算器高斯分布生成器Instagram用户ID查询卧推计算器原根计算器发音音标转换器圆计算器复数计算器模计算器死链检查器AI Token 计数器SRT合并工具伊斯兰历转换器体脂百分比计算器复合增长率计算器年金现值计算器百分比增长率计算器百分比计算器相关系数计算器随机虚假地址生成器雷诺数计算器kpa到psi转换器⏱️ 倒计时器总和计算器每个月的天数获取字符串长度视频压缩器🔊 音调发生器🥧 饼图制作工具MAC 地址分析工具十六进制转CMYK转换器摩尔斯电码生成器百分比增加计算器磅到克转换器圆形面积计算器循环播放视频📈 折线图制作工具数织生成器-picross随机坐标生成器隐形文本生成器黄金时刻和蓝调时刻计算器GIF 转 MP4 转换器凯撒密码工具图片添加线条快速傅里叶变换FFT计算器指数计算器-高精度数字提取器月亮星座计算器邻接矩阵计算器Argon2哈希生成器Webhook测试器YouTube频道统计中点计算器圆锥展开图模板生成器杀手数独生成器阅读时间计算器atan2计算器⚔️ DPS计算器Mann-Whitney U 检验计算器误差函数计算器质数列表AI标点符号添加器Zalgo文本生成器删除线文字生成器anova计算器Bar to PSI 转换器JWT生成器商和余数计算器极限计算器矩形计算器迷宫生成器随机数学题生成器黄金分割计算器WHOIS查询圆台计算器康奈尔笔记生成器排序数字样本标准差计算器特征值和特征向量计算器百分比减少计算器随机RPG角色生成器随机域名生成器p值计算器RC时间常数计算器不可见字符移除器两个日期之间中位数计算器二进制计算器图片Base64转换器度分秒到十进制转换器惯性矩计算器按位计算器板球击球率计算器净得分率计算器KD比率计算器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计算器每条线索成本计算器营运资本计算器youtube收益估算器