简化您的工作流程:搜索 miniwebtool。
添加插件
主页 > 数学 > 微积分 > 常微分方程组求解器
 

常微分方程组求解器

通过符号和数值方法求解常微分方程组 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计算器磅转千克转换器厘米到英尺和英寸转换器随机选择器毛利率计算器样本量计算器VAT计算器cpm计算器定期存款计算器百分比折扣计算器📅 日期计算器百分比增长率计算器kg到lbs转换器斜边计算器音频提取器图片打码工具音频分割器比例计算器血糖转换器质数检查器FPS 转换器名人名言搜索 (英文)英寸到厘米转换器罗马数字转换器MAC地址生成器圆计算器卡方检验计算器🎮 游戏灵敏度转换器厘米到英寸转换器MAC 地址分析工具平方根计算器随机名字选择器随机字符串生成器合并视频SRT转为TXT工具HEX计算器DOY日历百分比计算器股票平均成本计算器年度天数计算器 - 今天是今年的第几天t检验计算器体脂百分比计算器视频转图片提取器对数计算器百分比增加计算器月亮星座计算器跑步配速计算器调整视频速度闰年清单年龄计算器英尺到米转换器两个日期之间复利计算机数字提取器图片压缩器条形码生成器真心话大冒险生成器线性回归计算器分数百分比转换器利润计算器鞋码转换器卧推计算器删除空格分数计算器最简分数计算器Facebook用户ID查询PSI 转 Bar 转换器HEX转换器srt时间偏移凯利公式计算器盎司到克转换器二进制计算器互补误差函数计算器图片分割器圆形面积计算器每个月的天数AI Token 计数器One Rep Max (1RM) 计算器椭圆周长计算器随机IMEI生成器为图片添加文字半衰期计算器年金现值计算器kpa到psi转换器组合计算器不可见字符移除器百分比到ppm转换器分贝 (dB) 计算器AI内容检测器质数列表随机分组生成器随机化数字日本算盘模拟器俄罗斯农民乘法吠陀数学技巧计算器古埃及乘法计算器罗马数字数学求解器心算训练器乘法口诀表测验进位与借位可视化工具数的分合生成器硬币应用题求解器距离速度时间三角形计算器工作效率问题求解器混合问题求解器年龄问题求解器火车相遇问题求解器补水计算器配速卡路里计算器药物剂量计算器酒精卡路里计算器身体重塑计算器随机辩论话题生成器随机猫狗名字生成器随机圣经经文生成器随机数学题生成器随机段落生成器随机英文句子生成器砾石、砂和表土计算器钢材重量计算器螺栓扭矩计算器管道流量计算器梁荷载计算器美元换黄金转换器期权概率计算器股票拆分计算器员工持股计划计算器发票滞纳金计算器自由职业者时薪计算器租赁与购买对比计算器高级小费分摊计算器装箱清单生成器时差反应计算器旅行预算计算器飞行距离计算器热损失计算器发电成本计算器用水量计算器家电用电成本计算器家庭能源审计计算器太阳能投资回报率计算器太阳能板计算器堆肥CN比计算器草坪肥料计算器霜冻日期计算器高床种植箱土壤计算器NPK肥料计算器种子发芽率计算器视频比特率计算器音乐调性转换器音乐BPM节拍点击器照片文件大小估算计算器百万像素到打印尺寸计算器裁切系数计算器曝光三角计算器车辆牵引能力计算器汽车租赁计算器0–60与四分之一英里计算器电动车充电时间计算器电动汽车续航计算器汽车油耗计算器服装尺码转换器纸张尺寸参考表戒指尺寸转换器天文单位转换器燃油效率转换器数据传输速率转换器扭矩转换器 (Nm, ft-lb, kgf-cm)删除线文字生成器空白字符可视化工具阅读时间计算器演讲时间计算器段落计数器句子计数器音节计数器文本转二进制/十六进制/ASCII转换器Lorem Picsum / 占位符图片生成器.env 文件生成器Git 命令生成器颜色代码转换器全格式Bcrypt 哈希生成器和校验器JWT生成器CSS Grid 生成器数值积分计算器z变换计算器快速傅里叶变换FFT计算器张量积计算器矩阵指数计算器约当标准形计算器环与域计算器群论阶数计算器常微分方程组求解器伯努利微分方程求解器欧拉方法计算器方向场斜率场绘图器二阶常微分方程求解器一阶常微分方程求解器稳定婚姻问题求解器网络最大流计算器平面图检查器哈密顿路径检查器旅行商问题求解器 TSP线性规划求解器容斥原理计算器递推关系求解器邻接矩阵计算器拓扑排序计算器图着色计算器逻辑门模拟器卡诺图 (K-Map) 求解器布尔代数化简器分拆函数计算器数字根计算器斐波那契数检查器埃及分数计算器莫比乌斯函数计算器哥德巴赫猜想验证器梅森素数检查器孪生素数查找器亲和数检查器完全数检查器模幂运算计算器重复排列计算器效果量计算器相对风险计算器优势比计算器列联表计算器费舍尔精确检验计算器斯皮尔曼等级相关系数计算器贝塔分布计算器威布尔分布计算器指数分布计算器几何分布计算器负二项分布计算器超几何分布计算器F检验/F分布计算器贝叶斯定理计算器特征多项式计算器矩阵幂计算器乔列斯基分解计算器QR分解计算器矩阵对角化计算器克莱姆法则计算器列空间计算器零空间计算器向量夹角计算器单位向量计算器向量模计算器向量叉积计算器向量点积计算器矩阵乘法计算器逆矩阵计算器RREF计算器行最简阶梯形牛顿迭代法计算器雅可比矩阵计算器曲面积分计算器线积分计算器旋度计算器散度计算器梯度计算器多变量优化计算器微积分相关变化率求解器瞬时变化率计算器平均变化率计算器无限级数求和计算器级数收敛判定计算器幂级数计算器麦克劳林级数计算器洛必达法则计算器广义积分计算器辛普森法则计算器梯形法则计算器黎曼和计算器参数曲线绘图器旋转体表面积计算器旋转体体积计算器坐标几何距离计算器海伦公式计算器圆的切线计算器角平分线计算器内切圆计算器三角形外接圆计算器大圆距离计算器3D距离计算器环面计算器圆台计算器不规则多边形面积计算器正多边形计算器圆锥曲线识别器双曲线计算器抛物线计算器二项式定理展开计算器帕斯卡三角形生成器乘积符号计算器 (Pi记号)西格玛求和计算器有理根定理计算器笛卡尔符号法则计算器平行线和垂直线计算器直线方程计算器标准形式转斜截式转换器点斜式计算器非线性方程组求解器有理方程求解器字母方程求解器三角方程求解器指数方程求解器对数方程求解器四次方程求解器三次方程求解器估算计算器数字转分数转换器跳数生成器单位费率计算器上取整和下取整计算器绝对值计算器数列模式查找器位值图生成器运算顺序计算器PEMDAS竖式加减法计算器长乘法计算器乘法表生成器🎮 游戏货币换算器🎲 掉落概率计算器🎰 抽卡保底计算器⚔️ DPS计算器❄️ 雪天计算器🚚 搬家费用估算器🔍 抄袭检测器📷 OCR / 图片文字识别📈 折线图制作工具🥧 饼图制作工具📊 柱状图制作工具🔊 音调发生器🖱️ 点击计数器在线记事本⬛ 宽高比计算器🌍 碳足迹计算器向 文胸尺码计算器轮胎尺寸计算器燃油费用计算器💧 露点计算器🌡️ 体感温度计算器🌬️ 风寒指数计算器⏰ 在线闹钟⏰ 考勤卡计算器📅 日期差计算器🕐 军事时间转换器⏱️ 小时计算器⏱️ 在线秒表⏱️ 倒计时器🌐 时区转换器地毯计算器挡土墙计算器HVAC容量计算器隔热材料计算器铺路石计算器钢筋计算器木材计算器平方英尺计算器交叉相乘计算器五数概括计算器百分位数计算器正态分布计算器p值计算器比率计算器配方法计算器四舍五入计算器长除法计算器Twitter/X 字符计数器YouTube评论抽选器YouTube标签提取器YouTube缩略图下载器youtube收益估算器随机RPG角色生成器