方向場・傾き場プロッター

方向場・傾き場プロッター

方向場・傾き場プロッターは、任意の1階常微分方程式 y' = f(x, y) の幾何学的構造を、解析的に解くことなく視覚化する電卓です。カスタマイズ可能なグリッド上のすべての点で、その点での傾きが f(x, y) に等しい小さな接線分を描画し、解曲線の族全体を一目で明らかにします。インタラクティブな SVG キャンバスにより、クリックして RK4 積分された解曲線を生成したり、場に沿って流れる粒子をアニメーション化したり、結果を出版品質の画像としてエクスポートしたりできます。

方向場とは何ですか？

1階 ODE y' = f(x, y) が与えられたとき、方向場（傾き場とも呼ばれる）は、規則的な間隔で配置された点 (x_i, y_j) に置かれた短い線分のグリッドです。各線分は f(x_i, y_j) の傾きを持ち、これはその点を通過する解曲線の接線の傾きとなります。解はどこを通っても常に場に接していなければならないため、全体の図形は、明示的な公式を書き下ろす前に、アトラクタ、リペラ、平衡線、振動といった ODE の質的な振る舞いを示してくれます。

この手法は、微分方程式の定性的理論の一部として20世紀初頭に普及し、現在ではすべての初等常微分方程式コースの標準的な教育ツールとなっています。

このプロッターが優れている理由

解曲線の計算方法

指定された各初期条件 (x₀, y₀) に対して、ツールは古典的な4次ルンゲ＝クッタ（RK4）法を使用して ODE を積分します。RK4 は1ステップにつき、開始時に1回、中間に2回、最後に1回、計4回傾きをサンプリングし、それらを加重平均で組み合わせます。

RK4 の局所切り捨て誤差は O(h⁵)、大域誤差は O(h⁴) であるため、ステップサイズが小さくなるにつれて、オイラー法よりも4倍速く真の解に収束します。プロッターは (x₀, y₀) から前方および後方の両方向に積分するため、曲線は初期点の両側に伸び、表示領域全体を満たします。

プロットの読み方

平衡線とヌルクライン

線分が水平になる場所はヌルクライン（f(x, y) = 0 となる曲線）の上にあります。自励系 ODE y' = g(y) では、定数ヌルクラインは平衡解となります。色分けにより、これらは青い水平方向の帯として簡単に見つけることができます。

安定平衡 vs 不安定平衡

安定平衡点では、周囲の解がそこに向かって巻き込みます。上の矢印は下を向き、下の矢印は上を向きます。不安定平衡点では、逆のことが起こります。y' = y(1 − y) の場合、y = 1 は安定、y = 0 は不安定です。これはロジスティックプリセットで即座に確認できます。

急峻な領域とスティフネス

赤い線分は |f(x, y)| が大きい場所を示しており、そこでは解が急速に変化します。プロットが赤で占められている場合、その領域で方程式はスティフ（剛い）であり、数値積分器が精度を維持するためには小さなステップサイズが必要になります。

受け入れ可能な入力形式

1. 微分方程式

x と y を使用した有効な数学式として解析できるものであれば何でも可能です。一般的な例: y - x, x*y, sin(x) - y, exp(-x^2) + y, y*(1-y)。キャレット ^ は自動的に ** に変換されます。

2. ドメイン

x 範囲と y 範囲の4つの数字です。正方形のドメインが最も読みやすいプロットを提供します。一方の軸が極端に長い場合、傾きの値は正しくても接線分が歪んで見えることがあります。

3. 初期条件

セミコロンまたは改行で区切られた x, y のペアのリストです。各ペアが1つの RK4 解曲線になります。最大8つの初期条件を受け付けます。追加の曲線は、プロットをクリックしてインタラクティブに追加できます。

このプロッターの使い方

1. 発現式フィールドに y' = f(x, y) の右辺を入力するか、6つのプリセット例から1つを選択して古典的な振る舞いを確認します。
2. x 範囲と y 範囲を設定します。 興味深い振る舞いの近くを中心とした正方形の領域から始め、その後、範囲を絞って再送信することでズームインします。
3. 初期条件を セミコロンで区切られた x, y ペアとしてリストします。空白のままにして、プロット後に曲線を追加することもできます。
4. 「方向場をプロット」をクリックします。 傾き線分、色分けされた大きさ、および指定した解曲線を含む SVG が即座にレンダリングされます。
5. 操作: キャンバス上の任意の場所をクリックまたはタップして解曲線を追加し、ホバーして (x, y, 傾き) を読み取り、「フローをアニメーション化」を押して粒子が場に沿って流れるのを確認し、あるいは「SVG を保存」してエクスポートします。

計算例

古典的な方程式 y' = y − x を例に取ります。ヌルクラインは傾きがゼロになる直線 y = x です。この線より上では傾きが正（矢印が上を向く）、下では傾きが負（矢印が下を向く）であるため、すべての解曲線は垂直方向に y = x から漸近的に遠ざけられます。

プロッターはこの幾何学的構造を視覚的に裏付けます。特殊解 y = x + 1 を除くすべての軌跡は指数関数的に発散し、色分けによって直線 y = x が傾きが消失する鮮やかな青い筋として現れます。

一般的な用途

• ODE の概念の教育 — 平衡、安定性、吸引領域、サドルの挙動。
• 解析解の確認 — 手計算で導出した曲線を場に重ねて、接していることを確認。
• 個体群モデルの探索 — ロジスティック、アリー効果、収穫項はすべて特徴的な方向場のシグネチャを持ちます。
• 制御システムの視覚化 — 1階線形コントローラは y' = −k·y + u(x) に帰着し、その方向場は応答率を示します。
• 図表の作成 — 講義ノート、教科書、技術レポート用の図を作成（ロスレス出力には SVG 保存を使用）。

制限事項

このツールは1階の明示的 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) に等しい傾きを持ちます。これは1階 ODE y' = f(x, y) の右辺です。ODE の解曲線はすべての点で線分に接する必要があるため、方程式を解析的に解くことなく解の族全体を視覚化できます。

ツールはどのように解曲線を描画しますか？

提供されたすべての初期条件に対して、ツールは小さなステップサイズで古典的な4次ルンゲ＝クッタ（RK4）法を使用して ODE を数値的に積分します。RK4 は1ステップあたり4回傾きを評価し、それらを加重平均で組み合わせて、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 などがあります。

色の意味は何ですか？

「傾きの大きさで色分け」が選択されている場合、各傾き線分はその点での傾きの大きさに応じて対数スケールで色付けされます。青は小さな傾き（ほぼ水平な流れ）を示し、赤は大きな傾き（ほぼ垂直な流れ）を示します。これにより、平衡線、スティフな領域、アトラクタなどの特徴を一目で確認できます。

ヌルクラインとは何ですか？なぜ重要なのですか？

ヌルクラインは f(x, y) = 0 となる点の集合であり、方向場はヌルクラインに沿って水平になります。自励系の ODE では、ヌルクラインに平衡解が含まれることが多く、非自励系の方程式では解の転換点となります。ツールでは、色分けがオンのときに、これらの領域をほぼ水平な青い線分で強調します。

このツールはモバイルで使用できますか？

はい。レイアウトは小さな画面に適応し、SVG プロットはタッチイベントを使用するため、キャンバス上の任意の場所をタップして新しい解曲線を追加できます。すべての計算はサーバー側で実行されるため、ツールはスマートフォン、タブレット、デスクトップで同様に動作します。

参考文献

• 傾き場 — Wikipedia
• ルンゲ＝クッタ法 — Wikipedia
• Nullcline — Wikipedia (英語)
• 常微分方程式 — Wikipedia