简化您的工作流程:搜索 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计算器太阳、月亮与上升星座计算器 🌞🌙✨随机选择器cpm计算器样本量计算器百分比折扣计算器VAT计算器HEX计算器比例计算器随机扑克牌生成器二进制计算器毛利率计算器音频提取器FPS 转换器月亮星座计算器定期存款计算器随机名字选择器时间持续计算器SRT转为TXT工具音频分割器kg到lbs转换器斜边计算器图片打码工具年龄计算器血糖转换器股票平均成本计算器厘米到英寸转换器MAC地址生成器厘米到英尺和英寸转换器📅 日期计算器随机IMEI生成器两个日期之间磅转千克转换器视频转图片提取器随机字符串生成器合并视频分数计算器英寸到厘米转换器MAC 地址分析工具卡方检验计算器💧 露点计算器螺栓扭矩计算器凯利公式计算器圆计算器百分比变化计算器One Rep Max (1RM) 计算器🎮 游戏灵敏度转换器随机分组生成器英尺到米转换器调整视频速度srt时间偏移质数检查器椭圆周长计算器三角函数绘图器百分比增长率计算器罗马数字转换器图片分割器年度天数计算器 - 今天是今年的第几天百分比计算器石头剪刀布生成器随机用户代理生成器对数计算器随机化数字SHA256 哈希生成器t检验计算器最简分数计算器商和余数计算器条形码生成器利润计算器MD5哈希生成器百分比增加计算器线性回归计算器每个月的天数随机超能力生成器日历DOY日历半衰期计算器圆形面积计算器职位查找器Facebook用户ID查询随机坐标生成器快速傅里叶变换FFT计算器电池续航计算器分数百分比转换器移除标点符号在线工具随机端口号生成器位数计算器为图片添加文字平方根计算器随机装备生成器六西格玛过程能力计算器逻辑门模拟器随机颜色生成器体脂百分比计算器数字提取器两点间距离计算器分贝 (dB) 计算器视频压缩器SRT合并工具复合增长率计算器股息收益率计算器行数统计工具AI Token 计数器自酿啤酒酒精度计算器随机字母生成器黄金分割计算器kpa到psi转换器命运数字计算器圆台计算器方向场斜率场绘图器质数分解计算器闰年清单Log Base 10 计算器图片添加线条圆柱体体积计算器 高精度填字游戏制作器查找并替换文本年金现值计算器组合计算器随机数字选择器AI改写工具AI标点符号添加器卧推计算器根式化简器箱线图生成器PSI 转 Bar 转换器方差计算器 高精度获取字符串长度迷宫生成器AI内容检测器Instagram用户ID查询🎲 掉落概率计算器楼梯计算器百分比减少计算器英亩到平方米转换器随机域名生成器YouTube频道统计克到磅转换器工资转换计算器IPv4/IPv6到十六进制转换器为视频添加水印十进制到十六进制转换器平均值计算器盎司到克转换器字符计数器相关系数计算器磅到克转换器视频分割器图片压缩器模计算器死链检查器黄金时刻和蓝调时刻计算器AI语言检测器标准偏差计算器 - 高精度科学计算器atan2计算器克到盎司转换器🎰 抽卡保底计算器🌍 碳足迹计算器线性方程组求解器误差函数计算器平均偏差计算器指数计算器-高精度旋转视频沸点计算器百分比到ppm转换器隐形文本生成器YouTube缩略图下载器⏱️ 倒计时器农历转换器双重积分计算器真心话大冒险生成器翻转视频随机圣经经文生成器鞋码转换器文本去重工具十六进制到十进制转换器坡度与倾斜度计算器幻方生成器折扣计算器斜率截距式计算器步数距离计算器随机生日生成器anova计算器原根计算器⏱️ 小时计算器工作日计算器文本格式化工具标准误差计算器直方图生成器竖式加减法计算器配色方案生成器随机物品生成器标准杯计算器葡萄酒搭配建议器攀岩难度等级转换器自行车齿轮比计算器钓鱼结强度计算器瑜伽体式保持计时器游泳SWOLF计算器跑步成绩预测计算器拳击出拳力量计算器橄榄球得分计算器板球得分率计算器足球 xG预期进球计算器网球计分器Wells评分计算器 (DVT/PE)格拉斯哥昏迷评分计算器阿普加评分计算器FFMI计算器库珀12分钟跑步计算器一英里步行测试Rockport计算器瘦体重力量计算器碳水化合物胰岛素比例计算器胰岛素敏感系数计算器希伯来历转换器伊斯兰历转换器跨文化年龄计算器多久以前计算器还有多久倒计时计算器日期模式生成器中间日期计算器日期添加工作日词频分析器句子长度方差分析器海明威风格可读性编辑器发音音标转换器维吉尼亚密码工具埃特巴什密码工具ROT13编码解码器EXIF数据查看与移除工具猪拉丁文翻译器倒推首字母缩写生成器首字母缩写生成器全字母句检查器漏字文检测器图像转SVG描摹器图片转 ASCII 艺术转换器JSON Schema 生成器TypeScript 在线演练场Less 到 CSS 编译器SCSS转CSS编译器SVG 转 React/JSX 转换器查询字符串生成器URL解析器UUID验证和解码器HTTP状态码参考cURL命令构建器谢尔宾斯基三角形生成器3D曲面绘图器极坐标方程绘图器朱利亚集合生成器曼德博集合探索器L-System分形生成器Delaunay 三角剖分生成器Voronoi 图生成器万花尺图案生成器镶嵌图案生成器帕累托图生成器NPS净推荐值计算器留存率同期群计算器客户流失率计算器客户获取成本CAC计算器客户终身价值 CLV 计算器转化率计算器A/B测试样本量计算器A/B测试显著性计算器透镜方程计算器导线磁场计算器电场计算器库仑定律计算器斯涅尔定律计算器惯性矩计算器角速度计算器向心力计算器单摆周期计算器弹簧劲度系数计算器多普勒效应计算器索提诺比率计算器特雷诺比率计算器股票贝塔系数计算器通胀保值美国国债 (TIPS) 计算器房贷重新摊还计算器远期利率计算器债券久期计算器 (麦考利和修正)债券凸性计算器固定指数年金计算器变额年金计算器反向抵押贷款计算器年金支付计算器日本算盘模拟器俄罗斯农民乘法吠陀数学技巧计算器古埃及乘法计算器罗马数字数学求解器心算训练器乘法口诀表测验进位与借位可视化工具数的分合生成器硬币应用题求解器距离速度时间三角形计算器工作效率问题求解器混合问题求解器年龄问题求解器火车相遇问题求解器补水计算器配速卡路里计算器药物剂量计算器酒精卡路里计算器身体重塑计算器随机辩论话题生成器随机猫狗名字生成器随机数学题生成器随机段落生成器随机英文句子生成器砾石、砂和表土计算器钢材重量计算器管道流量计算器梁荷载计算器美元换黄金转换器期权概率计算器股票拆分计算器员工持股计划计算器发票滞纳金计算器自由职业者时薪计算器租赁与购买对比计算器高级小费分摊计算器装箱清单生成器时差反应计算器旅行预算计算器飞行距离计算器热损失计算器发电成本计算器用水量计算器家电用电成本计算器家庭能源审计计算器太阳能投资回报率计算器太阳能板计算器堆肥CN比计算器草坪肥料计算器霜冻日期计算器高床种植箱土壤计算器NPK肥料计算器种子发芽率计算器视频比特率计算器音乐调性转换器音乐BPM节拍点击器照片文件大小估算计算器百万像素到打印尺寸计算器裁切系数计算器曝光三角计算器车辆牵引能力计算器汽车租赁计算器0–60与四分之一英里计算器电动车充电时间计算器电动汽车续航计算器3D距离计算器环面计算器不规则多边形面积计算器正多边形计算器圆锥曲线识别器双曲线计算器长除法计算器Twitter/X 字符计数器YouTube评论抽选器YouTube标签提取器youtube收益估算器随机RPG角色生成器