方向场斜率场绘图器

方向场斜率场绘图器

方向场斜率场绘图器无需通过解析法解方程，即可直观展示任何一阶常微分方程 y' = f(x, y) 的几何特征。它在自定义网格的每个点上绘制一段微小的切线段，其斜率等于 f(x, y)，从而让您一眼洞察整个解曲线族的性状。交互式 SVG 画布允许您通过点击生成 RK4 积分的解曲线、开启动画观察粒子沿场流动，并将结果导出为出版质量的图像。

什么是方向场？

给定一阶常微分方程 y' = f(x, y)，方向场（也称为斜率场）是放置在规则间隔点 (x_i, y_j) 处的一组短线段。每条线段的斜率为 f(x_i, y_j)，即通过该点的任何解曲线的切线斜率。由于解曲线在经过的任何地方都必须与向量场相切，因此在您写出显式公式之前，整体图像就能向您展示常微分方程的定性行为 —— 吸引子、排斥子、平衡线和振荡。

这种技术在 20 世纪初作为微分方程定性理论的一部分而流行，现在已成为每一门常微分方程入门课程的标准教学工具。

为什么这个绘图器与众不同

解曲线是如何计算的

对于您提供的每个初始条件 (x₀, y₀)，该工具会使用经典的四阶龙格-库塔法 (RK4) 对常微分方程进行积分。RK4 在每步采样四次斜率 —— 步初一次、步中两次、步末一次 —— 并以加权平均值组合它们：

RK4 的局部截断误差为 O(h⁵)，全局误差为 O(h⁴)，因此随着步长的缩小，它收敛到真实解的速度比欧拉法快四倍。绘图器从 (x₀, y₀) 开始进行向前和向后积分，因此曲线会延伸到初始点的两侧并填满整个可见区域。

如何阅读图表

平衡线和零斜率线 (Nullclines)

凡是线段变为水平的地方，您就处于零斜率线上 —— 即 f(x, y) = 0 的曲线。在自治常微分方程 y' = g(y) 中，常数零斜率线即为平衡解；着色功能使它们很容易被识别为蓝色水平带。

稳定平衡与不稳定平衡

在稳定平衡点处，相邻的解会向其靠拢：上方的箭头向下指，下方的箭头向上指。在不稳定平衡点处则相反。对于 y' = y(1 − y)，y = 1 是稳定的，而 y = 0 是不稳定的 —— 您可以在逻辑回归预设中立即看到这一点。

陡峭区域与刚性 (Stiffness)

红色线段标记了 |f(x, y)| 较大的地方，解在那里变化迅速。如果您的图表大部分是红色的，说明该方程在对应区域是刚性的，任何数值积分器都需要极小的步长才能保持精度。

接受的输入格式

1. 微分方程

任何可以解析为包含 x 和 y 的有效数学表达式。常见示例：y - x, x*y, sin(x) - y, exp(-x^2) + y, y*(1-y)。脱字符 ^ 会自动转换为 **。

2. 范围 (Domain)

x 和 y 范围的四个数值。正方形定义域能提供最易读的图表；如果一个轴比另一个轴长得多，即使斜率值正确，切线段看起来也会有失真感。

3. 初始条件

以分号或换行符分隔的 x, y 对列表。每一对都会生成一条 RK4 解曲线。最多接受 8 个初始条件；额外的曲线可以通过点击图表进行交互式添加。

如何使用本计算器

1. 在表达式字段中输入 y' = f(x, y) 的右侧部分，或从六个预设示例中任选其一以查看经典行为。
2. 设置 x 和 y 范围。 建议从以感兴趣行为为中心的方形区域开始，然后通过提交更窄的范围来进行缩放。
3. 列出初始条件，以分号分隔的 x, y 对形式。您也可以留空，在绘图后手动添加曲线。
4. 点击“绘制方向场”。 SVG 将立即渲染斜率段、颜色编码的大小以及您指定的任何解曲线。
5. 交互： 点击或触摸画布上的任何位置以添加更多解曲线，悬停查看 (x, y, 斜率) 值，点击“动画流”观察粒子流向，或点击“保存 SVG”导出。

计算实例

以经典方程 y' = y − x 为例。零斜率线是直线 y = x，该处斜率为零。在此线上方，斜率为正（箭头向上）；在此线下方，斜率为负（箭头向下）。因此，除了特解 y = x + 1 之外，每条解曲线在垂直方向上都会渐进式地被推离 y = x。

绘图器直观地确认了这一几何特征：除 y = x + 1 外，所有轨迹都呈指数级爆炸，着色功能将直线 y = x 变成了一条明显的蓝色条纹，表示斜率消失。

常见应用场景

• 教学 ODE 概念 —— 平衡、稳定性、吸引盆、鞍点行为。
• 验证解析解 —— 将您手算的解曲线叠加在斜率场上，确认是否相切。
• 探索人口模型 —— 逻辑回归、阿利效应 (Allee effect)、收获项在斜率场中都有独特的特征。
• 可视化控制系统 —— 一阶线性控制器可简化为 y' = −k·y + u(x)，其斜率场展示了响应速率。
• 准备图表 —— 适用于讲义、教科书和技术报告（使用“保存 SVG”获得无损输出）。

局限性

本工具仅处理显式一阶 ODE —— 像 dy/dx = f(x, y), dz/dx = g(x, y, z) 这样的系统需要相图工具。隐式方程 F(x, y, y') = 0 在绘图前必须改写为 y' = f(x, y) 的形式。在奇点附近（即 f(x, y) 为无穷大或未定义的点），网格会变稀疏，RK4 追踪会直接停止而不是进行推测。

常见问题

什么是方向场（斜率场）？

方向场或斜率场是放置在 x-y 平面上规则间隔点处的一组短线段。在每个点 (x, y) 处，线段的斜率等于 f(x, y)，即一阶常微分方程 y' = f(x, y) 的右侧部分。常微分方程的解曲线在每个点必须与这些线段相切，这让您无需通过解析法解方程即可直观地看到整个解族。

该工具如何绘制解曲线？

对于您提供的每个初始条件，该工具都会使用经典的四阶龙格-库塔法（RK4）以较小的步长对常微分方程进行数值积分。RK4 在每步评估四次斜率，并使用加权平均值将它们组合，从而产生精度达到 O(h^4) 的轨迹。曲线从起始点向后和向前追踪，直到离开绘图区域或斜率变为无穷大。

我可以在表达式中使用哪些函数？

您可以使用算术运算符 + - * / ^ 以及变量 x 和 y，加上三角函数 (sin, cos, tan, asin, acos, atan)、双曲函数 (sinh, cosh, tanh)、指数和对数函数 (exp, ln, log, log10)、平方根 (sqrt)、绝对值 (abs) 以及常数 pi 和 e。有效的表达式示例包括 y - x, x*y, sin(x)*cos(y) 和 exp(-x^2) + y。

颜色代表什么意思？

当选择“按 |斜率| 着色”时，每个斜率段根据该点斜率的大小使用对数刻度进行着色。蓝色表示斜率较小（接近水平流动），红色表示斜率较大（接近垂直流动）。这可以一目了然地揭示平衡线、刚性区域和吸引子等特征。

什么是零斜率线（nullcline），它为什么重要？

零斜率线是 f(x, y) = 0 的点集，因此斜率场沿零斜率线是水平的。在自治常微分方程中，零斜率线通常包含平衡解；在非自治方程中，它们标记了解的转折点。开启“按斜率着色”时，该工具会用接近水平的蓝色段突出显示这些区域。

我可以在手机上使用这个工具吗？

可以。布局适应小屏幕，且 SVG 绘图支持触摸事件，因此您可以点击画布上的任何位置来添加新的解曲线。所有计算均在服务器端执行，因此该工具在手机、平板电脑和桌面上运行效果完全一致。

延伸阅读

• 斜率场 — 维基百科 (英文)
• 龙格-库塔法 — 维基百科 (英文)
• 零斜率线 — 维基百科 (英文)
• 常微分方程 — 维基百科 (英文)