• 締切済み

常微分方程式、4次のルンゲクッタ法

(d^2x/dt^2)-2(dy/dt)=f(x) (d^2y/dt^2)+2(dx/dt)=g(y) この連立常微分方程式を4次のルンゲクッタ法で解くためにはどうすればいいのでしょうか?

みんなの回答

  • Tacosan
  • ベストアンサー率23% (3656/15482)
回答No.2

左辺に 1階微分だけを残すと dz/dt = f(x) + 2w, dx/dt = z, dw/dt = g(y) - 2z, dy/dt = w という連立 1階微分方程式ができて, これにふつ~に Runge-Kutta法を適用するだけです.

全文を見る
すると、全ての回答が全文表示されます。
  • Tacosan
  • ベストアンサー率23% (3656/15482)
回答No.1

1. dx/dt = z, dy/dt = w とおく 2. Runge-Kutta-Nystrom法 (o にはウムラウトが付く) を使う

teketeke_09_1
質問者

補足

早速回答ありがとうございます。こう置くと与えられた方程式は dz/dt-2w=f(x) dw/dt+2z=g(y) となりますよね。 ここからrunge-kutta法へいくまでの手順をもう少し詳しく教えていただけないでしょうか?よろしくお願いします。

全文を見る
すると、全ての回答が全文表示されます。

関連するQ&A

  • オイラー法、2次ルンゲクッタ法、4次ルンゲクッタ法のC言語プログラムに

    オイラー法、2次ルンゲクッタ法、4次ルンゲクッタ法のC言語プログラムについて教えてください! 課題なのですが、まったくわからず困ってます>< 1 常微分方程式 dy/dx=f(x,y),y(0)=1 の数値解をオイラー法を用いて計算するプログラムを作為せよ。ただし、f(x,y)=3-6x^2-4x+2xyとする。 2 α=1,β=1,γ=1/2,σ=1/2 の場合の2次ルンゲクッタ法を考える。1と同じ常微分方程式(f(x,y)も同じ)を考え、その数値解を求めるプログラムを作成せよ。また、オイラー法と2次ルンゲクッタ法の実行結果を示して、2つの近似精度を比較せよ。 3 1と同じ常微分方程式(f(x,y)も同じ)を考え、その数値解を4次ルンゲクッタ法を使って求めるプログラムを作成せよ。また、オイラー法、2次ルンゲクッタ法、4次ルンゲクッタ法の実行結果を示して、3つの近似精度を比較せよ。 以上の3つです。 休んでいた自分が悪いのですが、ネットで調べてもよくわからなくて… わかる方、よろしくおねがいします…

  • 連立常微分方程式の問題。。。

    手元の参考書などを調べても、連立された微分方程式について書いていなくて困っています。 以下の問題なのですが、どのように進めていけばよいのでしょうか?? ------------------------------------------------------- 問) dx/dt + 2x - 3y = exp(t) dy/dt - 3x + 2y = exp(2t) について、以下の問に答えなさい。 (1)x に関する2階の非同次常微分方程式を求めなさい。 (2)(1)を解き、x の一般解を求めなさい。 (3)(2)を用い、y の一般解を求めなさい。 -------------------------------------------------------- 基本的なものなのかもしれませんが、連立微分方程式について、 一般的にどのように取り組んだらよいのかわからず困っています。 お手数ですが、よろしくお願いします。

  • 連立微分方程式

    連立微分方程式の問題です。 dx/dt = -3x+y dy/dt = 5x+y 回答お願いします。

  • 連立微分方程式

    連立微分方程式 dx/dt=5x+4y dy/dt=-x+y の一般解を求めて下さい。

  • 連立微分方程式の解き方を教えてください.

    連立微分方程式の解き方を教えてください. 2d(^2)y/dt^2-dx/dt-4y = t 4dx/dt+2dy/dt-3x = 0 ヒントとしてtで一回微分するとよいとありました. まだ勉強を初めて間もないので,解法が本当にわかりません. お手数ですが,御教授よろしくお願いいたします.

  • 常微分方程式の解き方

    以下の常微分方程式をどのように解けばいいか教えてください。 (dy/dx)^2 +(x −3y −1)dy/dx + 2y^2 −2xy −x + y=0

  • 2階微分方程式の解法

    d^2y/dx^2+2*x*dy/dx=0 境界条件 x=0: y=1、x→∞: y→0 上記の微分方程式をルンゲクッタ法を使って数値的に解きたいのですが、どのように x→∞ : y→0 の境界条件をいれればいいかわかりません。どなたか教えて頂けないでしょうか。回答よろしくお願いします。

  • 微分方程式についてです

    dy/dx = -y の微分方程式で ,x=10tと置いた場合. dx/dt = 10 ∴dt/dx = 1/10 を使って, dy/dx = (dy/dt)(dt/dx) から, -y =(dy/dt)(1/10) ∴dy/dt = -10y とするのはいいのでしょうか. 1次近似なので,x=10t は分割が小さくなっただけのような気がするんですけど.

  • 常微分方程式の問題

    常微分方程式の問題でいくつか解けなかったところがあるので教えていただきたいです。 この章で扱っているのは 変数分離系・同時系・線形1階微分方程式・完全微分形・線形2階微分方程式(同次形)・線形2階微分方程式(非同次形) を扱っていました。 その内、一般解を求める以下の問題 (1)dy/dx=xe^-y (2)x(dy/dx)-y=1 (3)(2y-x^2)dx+(2x-y^2)dy=0 と 与えられた条件をそれぞれ満たす微分方程式の解を求める以下の問題 (1)dy/dx=y/x (x=1のときY-2) (5)y''+5y'+6y=0 (x=0のときy=0、y'=1) の問題が解くことができませんでした。 どなたか解法をわかりやすく教えていただけないでしょうか?

  • 連立微分方程式

    点P(x,y)は連立微分方程式 dx/dt=y dy/dt=-x を満たすものとする。t=0で原点以外の点から出発した点P(x,y)は、tが増加するにつれてどのようにふるまうか述べよ。図を用いてもよい。 この問題の解き方がよく分かりません。 連立微分方程式について、色々な文献を見てみたのですが、どうもいまいちです。 上の連立方程式を2つともdt=のかたちにして、dx/y=dy/-xという式にし、変数を分離して両辺を積分して・・・すると、x^2+y^2=Cという式に なりました。 円の方程式っぽいです。 でも、tは消えてしまい・・・ よく分からなくなってきました。 そもそもここまでの解き方も自分は間違っているのでしょうか?? ご意見やヒント、解答ヨロシクお願いしますm(_ _)m

このQ&Aのポイント
  • Windows初期セットアップの「富士通パソコンのご利用登録(無料)」画面にて「次へ」ボタンが表示されずセットアップが完了できない問題が発生しています。チェックボックスは全てチェック済みです。Windows11の場合も同じ現象が起きる可能性があります。
  • そこで、解決方法としては以下の手順を試してみることができます。 1. パソコンを再起動してみる 2. インターネット接続を確認する 3. パソコンのドライバを最新版に更新する 4. パソコンのセットアップをスキップして後から設定する このように対処することで、セットアップが正常に完了する可能性があります。
  • Windows初期セットアップが完了できない問題について、富士通FMVのご利用登録画面で「次へ」ボタンが表示されない場合は、再起動やドライバの更新などの対処方法を試してみましょう。Windows11でも同じ現象が起きる可能性がありますので、注意が必要です。
回答を見る