作業フローを簡素化：miniwebtoolを検索。
追加
こちらもどうぞ
一階常微分方程式ソルバー二階常微分方程式ソルバーベルヌーイ方程式ソルバールンゲクッタRK4法電卓オイラー法電卓誘導リアクタンス計算機キネマティクス計算するする機
ホームページ > 数学 > 微積分 > 常微分方程式系ソルバー
 

常微分方程式系ソルバー

常微分方程式系 x' = Ax を記号的および数値的に解きます。平衡点（サドル、ノード、スパイラル、センター）の自動分類、固有値と固有ベクトルの段階的な導出、一般解と特殊解の表示、そしてアニメーション付きのインタラクティブな相図描画に対応しています（2×2、3×3の線形系および2D非線形系）。

常微分方程式系ソルバー
クイック例：
安定ノード 中心（振動子） 安定渦状点 鞍点 重解（退化） 振り子 ファン・デル・ポール 3×3 3階微分
\(\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 相図を生成します。

常微分方程式系とは何ですか？

常微分方程式（ODE）系は、単一の変数（通常は時間 \(t\)）に対する複数の未知関数を、その導関数を通じて関連付けます。最も簡潔な形式では、

\[ \mathbf{x}'(t) = \mathbf{F}(t, \mathbf{x}(t)) , \qquad \mathbf{x} \in \mathbb{R}^n \]

定数行列 \(A\) に対して \(\mathbf{F}(t, \mathbf{x}) = A\mathbf{x}\) である場合、その系は線形かつ自励的です。ここが理論の最も美しい部分であり、長期的な挙動全体が \(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系の3つのケース

固有値 一般解 相図
異なる実数 \(\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行列の場合、分類全体が1つの図に収まります：

\[ \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\) が目立つように表示されるのはこのためです。これら3つの数字だけで平衡点の名前を特定できます。

非線形系と相図

現実世界の多くの常微分方程式は非線形であり、閉じた形式の解を持ちません。このツールでは、局所切断誤差が \(O(h^5)\) であり滑らかなベクトル場のデフォルトの主力手法である4次ルンゲ＝クッタ（RK4）法を用いて、方程式を数値的に積分することでこれらを処理します。

\[ \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. 係数または方程式を入力。クイック例をクリックして、標準的な系（安定ノード、中心、鞍点、振り子、ファン・デル・ポールなど）を事前入力できます。
  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}\) と置くことで2次元系として書き換えることができます：

\[ \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\) を持つ安定渦状点が得られます。軌道は原点に向かって渦を巻き、時系列は指数関数的に減衰する正弦波を示します。

応用例

ヒントと注意点

よくある質問

常微分方程式系とは何ですか？

常微分方程式（ODE）系は、単一の独立変数（通常は時間）に対する複数の未知関数の導関数を関連付ける連立方程式のセットです。古典的な形式は \( \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\) が1つしかない場合、 \( (A - \lambda I) w = v \) を解く一般化固有ベクトル \(w\) も必要になります。このとき一般解は \( \mathbf{x}(t) = c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t} \) という形式をとります。固有空間が2次元である場合、行列はその不変部分空間上で単位行列のスカラー倍となり、解は \( \mathbf{x}(t) = (c_1 v1 + c_2 v2) e^{\lambda t} \) に簡略化されます。

このツールは非線形系を記号的に解くことができますか？

非線形モードでは、4次ルンゲ＝クッタ（RK4）積分法を使用してシステムを数値的に解き、相図をプロットします。ほとんどの非線形系には閉じた形式の解が存在しないため、これが標準的なアプローチです。平衡点付近の局所的な挙動は、2×2線形モードで処理可能な線形化によって読み取ることができます。固定点でのヤコビ行列を計算し、それを \(A\) として入力してください。

相図とは何ですか？

相図（フェーズポートレート）は、\(x\)–\(y\) 平面における2次元系の解の幾何学的な図です。各解は「軌道」と呼ばれる曲線を描き、軌道の集合とベクトル場の矢印を合わせることで、解が渦を巻いて吸い込まれるのか、鞍点のように離れていくのか、振動するのか、あるいは平衡点に落ち着くのかといった定性的な挙動を明らかにします。相図により、系の全体構造を一目で把握できるようになります。

参考文献

このコンテンツ、ページ、またはツールを引用する場合は、次のようにしてください：

"常微分方程式系ソルバー"（https://MiniWebtool.com/ja/常微分方程式系ソルバー/） MiniWebtool からの引用、https://MiniWebtool.com/

miniwebtool チーム作成。更新日：2026年4月23日

また、AI 数学ソルバー GPT を使って、自然言語による質問と回答で数学の問題を解決することもできます。

微積分:

おすすめ:

InstagramユーザーID検索パーセンテージ減少電卓パーセント増加電卓弧長電卓ランダムカラージェネレーターwar電卓動画を結合MACアドレス検索画像分割ツールフィートとインチからセンチメートルへのコンバーター標準偏差電卓 - 高精度円錐展開図テンプレートジェネレーターエンジェルナンバー電卓ランダム誕生日ジェネレーター番号を並べ替える引用検索 (英語)動画を逆再生ランダム日付ジェネレーターランダムポーカーハンドジェネレーター合計電卓シグマ記法電卓 総和FPSコンバーターランダム超能力ジェネレーターパーセント誤差電卓平方完成電卓売上総利益率電卓相対標準偏差電卓YouTubeチャンネル統計モジュロ電卓英単語ランダム生成ツールランダムトーナメント表作成ツール外れ値電卓熱膨張計算機HEX電卓ASCIIコード表血糖値コンバーター動画を回転オンライン句読点削除ツールMP3ルーパーデシベル (dB) 電卓逆テキスト中央値電卓マスターナンバー電卓ai句読点追加kva計算機楕円円周電卓手数料電卓センチメートルからフィートとインチへのコンバーターボウリングスコア計算機表面積電卓CAGR電卓太陽位置計算機空の行を削除する労働時間計算ツール圧力電卓🖱️ クリックカウンターマン・ホイットニーのU検定計算機ビデオ速度を調整小数時間から普通の時間へのコンバーター筆算足し算・引き算計算機ボルト締付トルク計算機ビンゴカードジェネレーターt検定電卓比率電卓コラージュメーカー桁数電卓筆算割り算電卓四分位範囲電卓私のIPアドレスは何ですかHEXコンバーター中間日計算機加速度電卓土星回帰電卓ポアソン分布電卓配管流量電卓パーセント成長率電卓変化率電卓求人検索ランダムクレジットカードジェネレーター不可視文字除去ツール熱伝達計算機パレットジェネレーターランダム絵文字ジェネレーターオーディオ スプリッター画像回転ツール関数電卓TikTok収益計算ツール愛の相性電卓SRT 時間シフト 電卓指数電卓（高精度）斜面計算機身長パーセンタイル電卓動画から画像抽出ツール配当利回り電卓散布図作成ツールSHA256 ハッシュジェネレーター年の日電卓 - 今日は今年の何日目Twitch収益計算ツール正多角形電卓車両重量配分計算機週番号計算機底2の対数電卓パーセントからppmへのコンバーター分数電卓アンテナ長計算機階段電卓XMLバリデーター世界時計atan2電卓相関係数計算機ANC電卓CRC32チェックサム電卓fena電卓BUN対クレアチニン比電卓GIFメーカーヒストグラムメーカーサイコロローラータンジェント電卓飛行距離電卓ランダムトランプカードジェネレーター角度変換ツールジニ係数電卓平方根電卓CMYKからHEXへの変換ツール動画分割ツール小文字生成器 ⁽ᶜᵒᵖʸ ⁿ ᵖᵃˢᵗᵉ⁾沸点計算ツールAIテキストヒューマナイザーアナグラム生成器ランダム名前ジェネレーターFacebookユーザーID検索カラースキームジェネレータークロスワードパズルメーカーワイヤーゲージ電卓対数電卓💧 露点電卓PSIからbarへの変換器クーパー12分間走計算ツールサッカーxg期待ゴール電卓バーコードジェネレーターローマ数字のコンバーター三相電力計算機取り消し線テキスト生成ツールトルク電卓ランニングペース電卓三角関数グラフ作成ツール割り切れるテスト電卓周波数波長変換ツール直線の方程式電卓角速度計算機じゃんけんジェネレータービデオをループ再生ランダム時刻ジェネレーター服のサイズ変換猫カロリー電卓ベーカーズパーセント電卓SRTからTXTへの変換ツール動画圧縮積分電卓迷路ジェネレーターBCD?????psiからkPaへのコンバーターそろばんシミュレーター文字数による改行日割り家賃計算魔方陣ジェネレーター🔊 トーンジェネレーター点つなぎジェネレーターhba1c電卓タイル目地計算機バイナリ電卓平均寿命電卓ポンドからキログラム変換分圧器計算電卓素数ですかカレンダー電卓hcg倍加時間計算ツール変動係数電卓水星逆行カレンダー自転車サイズ電卓IPアドレスから16進数への変換慣性モーメント計算機斜辺電卓AIお礼状ジェネレーターダイスロール確率電卓ペース カロリー電卓分散電卓 高精度直角三角形電卓10進数から16進数へのコンバーターHTMLからテキストコンバータmcgからmgへの変換ツールVTTからtxtへのコンバーターボロノイ図ジェネレーターヴィジュネル暗号ツール外接円電卓平方数リスト直列抵抗計算機グロスアップ計算機ジェスチャーゲームジェネレーターレシピ栄養計算ツールTwitter (X)用動画変換ゲーマータグジェネレーターチーム名ジェネレーターニックネームジェネレーターラップネームジェネレーターバンド名ジェネレーターファンタジー名前ジェネレーターフォーチュンクッキージェネレーターお絵描きアイデアジェネレータージャーナルプロンプトジェネレーターデイリーアファメーションジェネレーター褒め言葉ジェネレーター口説き文句ジェネレーター親父ギャグジェネレーターランダムジョークジェネレーターランダム雑学ジェネレータートリビアクイズジェネレーターハングマン単語ジェネレーターピクショナリー単語ジェネレーターアイスブレイク質問ジェネレーターNever Have I Ever ジェネレーターどっちを選ぶジェネレーターランダム質問ジェネレーターAI夢占いAI結婚式スピーチジェネレーターAIプレスリリースジェネレーターAI商品説明ジェネレーターAI動画スクリプトジェネレーターAI YouTubeタイトル・概要欄ジェネレーターInstagramキャプションジェネレーターAI LinkedInヘッドライン・自己紹介ジェネレーターAI職務経歴書ブレットジェネレーターAIカバーレター作成ツールAIラップ歌詞ジェネレーターPDFフラット化ツールPDFロック解除（オーナーパスワード削除）PDFテキスト抽出PDFページ抽出PDF単語カウンターPDFにパスワードを設定PDFにページ番号を追加PDFページ並べ替えPDFページ削除PDF回転JPGをPDFに変換PDFをJPGに変換PDF圧縮PDF分割PDF結合PTからPXへの変換器PXからREMへの変換器トラからグラムへの変換器インドの土地単位変換器AI歌詞ジェネレーターAI詩ジェネレーターAIストーリージェネレーターAI翻訳AIテキスト要約ツールmgからmlへの変換器MPHからKMHへの変換器ノットからマイル毎時変換器立方メートルから立方フィート立方フィートから立方ヤード変換器平方フィートから平方メートル換算平方メートルから平方フィート変換器オンスからmlへの変換器mlからオンスへの変換器ガロンからリットル変換器リットルからガロン変換器キログラムからストーンへの変換器ストーンからキログラム変換器インチからミリメートルへの変換器ミリメートルからインチへの変換器マイルからキロメートル変換器キロメートルからマイル変換器華氏から摂氏への変換器摂氏から華氏への変換器バックパックサイズ電卓サーフボード容量電卓テニスグリップサイズ電卓ヘルメットサイズ電卓手袋サイズ電卓帽子サイズ変換スノーボードサイズ電卓スキーサイズ電卓スピードランスプリットタイマーステーブルフォード計算機ダーツチェックアウト電卓ボウリング・エコノミーレート電卓打撃ストライクレート電卓ネットランレート計算機KDレシオ計算機Eloレーティング計算機心拍数回復電卓ハイキング時間計算機ローイングペース計算機ランニング年齢グレード計算機FTP・パワーゾーン計算機Beepテスト電卓ACFTスコア電卓Wilks & DOTS 電卓ホイールオフセット電卓タイヤ負荷指数・速度レーティング検索マイルあたりコスト電卓リース買取電卓オクタン価ブレンド電卓2ストロークオイル混合電卓エンジン排気量計算ソファ搬入電卓薪コード計算機空気清浄機CADR電卓除湿機サイズ計算シーリングファンサイズ電卓カーテンサイズ計算機ラグサイズ電卓額縁の掛け高さ計算機テレビ取り付け高さ計算機テレビサイズ電卓池の容量・ライナー電卓プール塩量計算機プール容量計算給湯器サイズ電卓エポキシ樹脂計算機手すり子間隔電卓巾木とトリムの電卓サイディング計算機デッキ用ステインの計算芝生の種計算機芝生電卓アスファルト電卓立方ヤード電卓電線管充填率電卓直列並列コンデンサ電卓誘導リアクタンス計算機部屋の照明計算機ルクスからルーメン電卓ルーメンからワット変換器発電機サイズ計算機mAh Wh 変換器 電卓アンペアからワット電卓ワットからアンペア電卓摩擦計算機機械的倍率計算機音速計算機波の速度計算機浮力電卓終端速度計算機ド・ブロイ波長計算機光子エネルギー計算機emc2電卓時間の遅れ計算機ケプラーの第三法則計算機脱出速度計算機万有引力計算機ベール・ランベルトの法則電卓ネルンストの式計算機浸透圧計算機沸点上昇電卓凝固点降下計算パーセント組成計算機規定度計算機質量モル濃度計算機pKa→Ka変換器ヘンダーソン・ハッセルバルヒ電卓理論収量計算機制限反応物計算機電子配置計算機インタラクティブ周期表AI授業計画ジェネレーターAIクイズ作成ツール引用ジェネレーター (APA/MLA/Chicago)出席率計算APスコア計算機ACTスコア計算機SATスコア電卓youtube収益見積もりツールランダムRPGキャラクタージェネレーター