非線形連立方程式ソルバー

ニュートン・ラフソン法を使用して、2つ以上の未知数を含む非線形連立方程式を解きます。ステップごとの反復計算、ヤコビ行列、収束分析、およびインタラクティブな等高線グラフで、すべての解を見つけます。

非線形連立方程式ソルバー

非線形連立方程式ソルバーは、ニュートン・ラフソン法を使用して、2つ以上の非線形方程式からなるシステムのすべての解を見つけます。方程式を入力すると、ソルバーは詳細なステップごとの反復、ヤコビ行列分析、収束の視覚化、および2変数システム用のインタラクティブな等高線グラフを使用して、すべての解を自動的に検索します。

非線形連立方程式ソルバーの使い方

方程式を入力する: 変数 x、y（および3変数システムの場合は z）を使用して各方程式を入力します。方程式は x^2 + y^2 - 25（= 0 と見なされます）または x^2 + y^2 = 25 のように入力できます。累乗には ^、乗算には * を使用し、sin、cos、exp、log、sqrt などの標準的な関数を使用してください。
方程式の数を選択する: ドロップダウンから 2 または 3 を選択します。適切に決定されたシステムの場合、方程式の数は変数の数と等しくなければなりません。
初期推定値を設定する (オプション): x₀、y₀（および z₀）の開始値を入力します。ソルバーはこれらをニュートン・ラフソン反復の開始点として使用します。空白のままにすると、デフォルトで 1 になります。
「連立方程式を解く」をクリックする: ソルバーは、初期推定値からニュートン・ラフソン法を実行するほか、[-5, 5] の範囲でマルチスタート検索を実行してすべての解を見つけます。
結果を確認する: 見つかったすべての解、収束を示す反復テーブル、解におけるヤコビ行列、およびインタラクティブな等高線グラフ（2変数システムの場合）を調べます。

非線形連立方程式とは何ですか？

非線形連立方程式は、少なくとも1つの方程式に非線形項（$x^2$、$\sin(x)$、$e^x$、$xy$など）が含まれる2つ以上の方程式で構成されます。一般形式は以下の通りです：

$$\begin{cases} f_1(x_1, x_2, \ldots, x_n) = 0 \\ f_2(x_1, x_2, \ldots, x_n) = 0 \\ \vdots \\ f_n(x_1, x_2, \ldots, x_n) = 0 \end{cases}$$

（最大1つの解しか持たない）線形システムとは異なり、非線形システムは解がゼロ、1つ、または複数存在することがあり、解くのが大幅に難しくなります。

連立方程式におけるニュートン・ラフソン法

ニュートン・ラフソン法（ニュートン法とも呼ばれます）は、よく知られた単変数根探索アルゴリズムを連立方程式に拡張したものです。反復公式は次のとおりです：

$$\mathbf{x}_{k+1} = \mathbf{x}_k - J(\mathbf{x}_k)^{-1} \mathbf{F}(\mathbf{x}_k)$$

ここで、$\mathbf{F}$ は方程式のベクトル、$J$ はヤコビ行列です。実際には、逆行列を計算するのではなく、各ステップで線形システム $J \cdot \Delta\mathbf{x} = -\mathbf{F}$ を解きます。

ヤコビ行列

ヤコビ行列は、導関数を多変数ベクトル関数に一般化したものです。$n$ 個の未知数を持つ $n$ 個の方程式のシステムの場合、次のようになります：

$$J = \begin{pmatrix} \frac{\partial f_1}{\partial x_1} & \frac{\partial f_1}{\partial x_2} & \cdots & \frac{\partial f_1}{\partial x_n} \\ \frac{\partial f_2}{\partial x_1} & \frac{\partial f_2}{\partial x_2} & \cdots & \frac{\partial f_2}{\partial x_n} \\ \vdots & \vdots & \ddots & \vdots \\ \frac{\partial f_n}{\partial x_1} & \frac{\partial f_n}{\partial x_2} & \cdots & \frac{\partial f_n}{\partial x_n} \end{pmatrix}$$

このソルバーは、記号微分を必要とせず、良好な精度が得られる中心差分を用いて数値的にヤコビ行列を計算します。

収束特性

ニュートン・ラフソン法は、ヤコビ行列が非特異である解の近くで二次収束を示します。これは、各反復で正しい桁数がおよそ2倍になることを意味します。ただし、収束は以下に依存します：

初期推定値が解に十分に近いこと
解の近くでヤコビ行列が非特異であること (det(J) ≠ 0)
関数が滑らかであること（連続微分可能）

ヤコビ行列が特異または特異に近い場合、収束は線形に低下するか、メソッドが完全に失敗することがあります。

複数の解とマルチスタート戦略

ニュートン・ラフソン法は開始点に最も近い解に収束するため、このソルバーはマルチスタート戦略を使用しています。各変数について [-5, 5] の範囲のグリッド上で多くの異なる初期推定値を試行します。（異なる開始点から）複数回見つかった解は重複排除されます。このアプローチにより、検索範囲内のほとんどの解が見つかりますが、すべての解を見つけることを保証するものではありません。

等高線グラフの理解

2変数システムの場合、ソルバーはインタラクティブな等高線グラフを表示します。各方程式 $f_i(x,y) = 0$ は、xy平面上の曲線（ゼロレベルセット）を定義します。解はこれらの曲線の交点です。グラフには、初期推定値からのニュートン・ラフソン反復経路も表示され、アルゴリズムがどのように収束するかを示しています。