数値計算の演習問題で以下の二階微分方程式をルンゲ・クッタで解けという問題があります。
-y"+x^2・y=e・y(eは定数、”・”は単なる掛け算)
y(0)=1, y'(0)=0, 0<=x<=2までを計算せよ。
これは
y’=z・・・(1)
z'=(x^2-e)y・・・(2)
この2つの連立方程式を解けばよいところまではわかります。
まず(2)を解くときにルンゲ・クッタの場合
(k1+2k2+2k3+k4)/6の項(←公式の右辺第二項)のk(1~4)を求めなければいけません。
質問はkの求め方です。
本にはy'=f(x,y,z) , z'=g(x,y,z)とおけば
(2)の場合だと例えばk1は
k1=g(xn,yn,zn)dxで計算する。と書いてあります
しかしz'=(x^2-e)y(←zが入ってない) なので、計算すると
k1=g(xn,yn)dxとなってしまうんですがどうなんでしょう?
おそらくどこかで勘違いしてると思うんです。
長い質問になってしまいましたがどうかご教授のほどよろしくお願いします。
No.1
- 回答日時:
良いと思いますよ。
ちなみに、私は、dx=hとして、
k1=f(x,y,z)・h
l1=g(x,y,z)・h
k2=f(x+h/2,y+k1/2,z+l1/2)・h
l2=g(x+h/2,y+k1/2,z+l1/2)・h
k3=f(x+h/2,y+k2/2,z+l2/2)・h
l3=g(x+h/2,y+k2/2,z+l2/2)・h
k4=f(x+h,y+k3,z+l3)・h
l4=g(x+h,y+k3,z+l3)・h
としてプログラムを書きました。
どうもありがとうございます。
同じ問題やられてたんですか(笑)
お答えの中で
k1=f(x,y,z)・h
l1=g(x,y,z)・h
とありますが
実際プログラムで計算するとき
f(x,y,z)の中身は
f(x,y,z)=z
g(x,y,z)の中身は
g(x,y,z)=(x^2-e)・y
ということでいいんですか?
って何がいいたいかと言うと
例えばg(x,y,z)・hを実際に計算する場合、
g(x,y,z)=(x^2-e)・y
に変数(xとy)に数を代入して計算しますよね。
その場合g(x,y,z)(←zが入っている)にはどこにもzを代入する
ところがないんですがg(x,y,z)と記述するのは単に形式的なものと
考えていいんですか?
No.2ベストアンサー
- 回答日時:
卒論で微分方程式の数値計算をやったので多少の覚えはあります。
結局僕の場合はy'の項が入ってるのでリープフロッグの方が収束が早いという事でそっちにしましたが。> k1=g(xn,yn)dxとなってしまうんですがどうなんでしょう?
問題ありません。
> 本にはy'=f(x,y,z) , z'=g(x,y,z)とおけば
というのはy', z'がx,y,zの関数という一般の場合であり、仮に
f(x, y, z) = x + y
という関数が存在しても何の問題もないのと同じです。
k1を計算する際にも
k1 = (x^2 - e)y * dx
と値が一意に決まり、問題ないでしょう。
ちなみにこの程度の微分方程式なら刻み幅dxにも寄りますが、プログラムなど組まなくてもExcelで一発で出ますよ。グラフも。
どうもありがとうございます。
もしかしたら私は勘違いしてるんじゃないかと思ってたので・・
しつこく質問してしまいすいませんでした。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 連立微分方程式の解き方について 7 2022/12/16 13:39
- 数学 x=r・cosθの2回微分 θ=ωtとすると? 5 2022/05/10 23:53
- 数学 数学の問題の解説お願いします! 4 2022/08/28 05:22
- 物理学 ポテンシャルが有限で不連続の時、右側の波動関数をφ1(x)、左側をφ2(x)とする。境界条件の「波動 2 2023/06/04 13:53
- 電気工事士 6.6kVケーブル単芯325sq-1.5kmの遮蔽銅テープ抵抗値は何Ω? 1 2023/05/02 21:06
- 化学 化学のエンタルピ変化を求め方について ある例題では各物質のモール数を換算して計算することもあり、ある 1 2022/06/20 23:22
- 数学 数2Bの数列の問題です。 自分は、 まず数列 an=ar^(n-1)と置き こちらの問題の、y= の 1 2022/07/07 16:26
- 数学 京都大学教授が証明。 「ABC予想・宇宙際タイヒミューラー予想」を、ザックリで説明お願致出来ますか? 1 2022/04/11 20:52
- 数学 情報処理詳しい人!! A4縦のレポート文書に4:3の大きさの横向きの写真画像を貼り付けることにした。 2 2022/12/18 02:30
- 物理学 電磁気学、TEMモードでマクスウェル方程式を満たさなくなる。 1 2022/07/19 23:58
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
公共工事の現場管理費率(%)...
-
eのマイナス無限大乗
-
「割る」と「割りかえす」の違い
-
30パーセントオフで371円だった...
-
ラプラス変換に関して
-
一個当たり15秒の製品を1時間で...
-
中学生の数学を習う順番に並べ...
-
面積から辺の長さを出す計算式
-
2割負担の計算。
-
楕円の円周の長さの計算の仕方...
-
袋のサイズから容量を計算する方法
-
10進法で時間の計算で30分が0.5...
-
(x^2-x+1)^10の展開式におけるx...
-
積分のエクセル計算式を教えて...
-
分数の計算で分子が0になったら...
-
展開について
-
プール計算って何ですか?
-
計算手順について
-
クリストッフェル記号
-
a=9 b=8 c=7の三角形ABCの角Aを...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
30パーセントオフで371円だった...
-
公共工事の現場管理費率(%)...
-
「割る」と「割りかえす」の違い
-
中学生の数学を習う順番に並べ...
-
楕円の円周の長さの計算の仕方...
-
2割負担の計算。
-
10進法で時間の計算で30分が0.5...
-
プール計算って何ですか?
-
体積凹 理解しない 頭掻きむしる
-
面積から辺の長さを出す計算式
-
一個当たり15秒の製品を1時間で...
-
袋のサイズから容量を計算する方法
-
eのマイナス無限大乗
-
積分のエクセル計算式を教えて...
-
赤で囲んだ式でもできますか? ...
-
積分での計算ミス直す方法。
-
映画を1.3倍速で見た時の時間計...
-
半径の計算方法を教えてください。
-
最後の指針がわかりません
-
量子力学の運動量pの微分演算子...
おすすめ情報