- ベストアンサー
計算機シュミレーション
tintagelの回答
- tintagel
- ベストアンサー率77% (214/277)
専門家ではありませんがこんなＨＰを見つけたので参考にして下さい。 「Runge-Kutta 法による微分方程式解法」 上記ＨＰのリンクに下記の項目がありました。 「Runge-Kutta法のコードと計算のページ」 「上記のclassを少し変更して利用した差分プログラミングによるアプレット」
関連するQ&A
- 高階連立常微分方程式の数値計算
4次のルンゲクッタ法を用いた数値計算を勉強しています． 1階連立常微分方程式と高階常微分方程式は理解でき，プログラムも作成することができました． 次に高階の連立常微分方程式を解こうと思ったら，頭が混乱してしまいました． 4次のルンゲクッタ法を用いて高解連立常微分方程式を解く考え方を教えて頂ければ嬉しいです． また何か良い参考書があれば教えて頂きたいと思います． よろしくお願いします．
- ベストアンサー
- 数学・算数
- この方程式、stiffなんでしょうか・・・
僕は以下のような連立微分方程式を数値計算で解こうと試みています。 現在微分方程式の解を級数展開と変数変換とルンゲクッタ法で求めようとしています。 しかし、どのように初期値をいじっても必ず解が発散してしまいます。 色々と自分で調べてみる過程で微分方程式を数値計算で解くときにstiff問題というものが生じることがあると知りました。 解の挙動が急に変わる場合や境界条件に大きく依存する場合などで途中からルンゲクッタ法などで数値計算できなくなり、解が発散してしまうという条件がぴったり自分の場合と一致しました。 以下の連立微分方程式はstiffなのでしょうか？ (x-1)^4*f ''(x)+(2x-1)(x-1)^3/x*f '(x)-{ ((x-1)/x+A(x))^2+f^2-1 }f=0, (x-1)^4*A''(x)+(2x-1)(x-1)^3/x*A'(x)-(x-1)^2/x^2A(x)-f^2(A+(x-1)/x)=0 stiffだとルンゲクッタ法は使えなくなってしまうんですよね・・・・？
- ベストアンサー
- 数学・算数
- 常微分方程式、４次のルンゲクッタ法
(d^2x/dt^2)-2(dy/dt)=f(x) (d^2y/dt^2)+2(dx/dt)=g(y) この連立常微分方程式を４次のルンゲクッタ法で解くためにはどうすればいいのでしょうか？
- 締切済み
- 数学・算数
- 解き方を教えてください。
2粒子があってｒ１粒子がｒ２粒子から受ける力を考えて 連立2階微分方程式を数値的に解く。 解く方法は「ルンゲ・クッタ・ニストローム法」です。 どうやってｒのベクトルをルンゲクッタの公式を利用して 解いたらいいのかわかりません。 誰か助けてください。
- ベストアンサー
- 物理学
- ルンゲクッタ法について
ルンゲクッタ法についての問題なのですがよくわかりません。 常微分方程式 du/dt=u について、ルンゲクッタ法（3次.4次）による誤差分析を行えという問題です。
- ベストアンサー
- 数学・算数
- オイラー法、２次ルンゲクッタ法、４次ルンゲクッタ法のＣ言語プログラムに
オイラー法、２次ルンゲクッタ法、４次ルンゲクッタ法のＣ言語プログラムについて教えてください！ 課題なのですが、まったくわからず困ってます＞＜ 1 常微分方程式 dy/dx=f(x,y),y(0)=1 の数値解をオイラー法を用いて計算するプログラムを作為せよ。ただし、f(x,y)=3-6x^2-4x+2xyとする。 2 α=1,β=1,γ=1/2,σ=1/2 の場合の２次ルンゲクッタ法を考える。1と同じ常微分方程式(f(x,y)も同じ)を考え、その数値解を求めるプログラムを作成せよ。また、オイラー法と２次ルンゲクッタ法の実行結果を示して、2つの近似精度を比較せよ。 3 1と同じ常微分方程式(f(x,y)も同じ)を考え、その数値解を４次ルンゲクッタ法を使って求めるプログラムを作成せよ。また、オイラー法、２次ルンゲクッタ法、４次ルンゲクッタ法の実行結果を示して、3つの近似精度を比較せよ。 以上の３つです。 休んでいた自分が悪いのですが、ネットで調べてもよくわからなくて… わかる方、よろしくおねがいします…
- 締切済み
- C・C++・C#
- エクセルを使ってのシミュレーション
連立方程式をエクセルで解き，シミュレーションをしたいと思ってます．最初の状態としては制約条件をすべてゼロに等しくなるようにしています（そのため，制約条件を書き込んだセルは０になります）．その状態から，パラメータを変化させて，解きたいのですが，パラメータの数値を初期値と変えてしまうと制約条件のセルが０にどうしてもならないことになってしまいます．どのようにシミュレーションをしたら宜しいでしょうか？
- ベストアンサー
- オフィス系ソフト
- 数値解析 微分方程式
数値解析における常微分方程式を解くために用いる手法についてです。 オイラー公式、ホルン公式、ルンゲクッタ、ルンゲクッタ4次、有限差分法の関係 違いがよくわからなのでどなたか教えてください また常微分方程式を有限差分で解くとなったとき、結局オイラー公式などを使うと言う認識で間違いないでしょうか？それとも有限差分だけで解けるのでしょうか
- 締切済み
- 数学・算数
- 数値シミュレーションの勉強に役立つ文献
数値シミュレーションの勉強に役立つ文献を教えてください。 数値シミュレーションを勉強したいと考えています。 ※ここでいう数値シミュレーションは ・Flow3D ・CAE といった、既存のソフトを使用するものではなく、 C言語などを用いて実際に自分でプログラムをくむことを言っています。 伝わりますかね。 実際にシミュレーションを組む流れは １：支配方程式を導き ２：数学的な計算方法を盛り込む(ex ルンゲクッタ) といったものという認識でいます。 上記らを体系的に学びたいと考えています。 それに必要な文献や資料を、ご教授願います。ⵁ
- ベストアンサー
- 物理学
- 常微分方程式の数値計算
実験で常微分方程式の数値計算をしました y'(x)=2xy,y(0)=1.0(解はy=ex2）←eのx二乗の微分方程式できざみ幅ｈ＝0.1、区間〔0，2〕の条件で オイラー法と改良オイラー法とルンゲクッタ法を使って 表計算ソフトを使って求めよとあるのですが、 表計算で近似値のy座標とか Δｔ後のxの変化分K1とか求めるのですが、全然意味がわからなくて、 何かいい方法はないでしょうか？
- ベストアンサー
- 数学・算数
補足
ありがとうございました。 僕が今、扱ってる方程式は拡散方程式です。 その時間微分がどうにもなりません。 大変参考になりました。 でももう少し突っ込んだものがありましたらよろしくお願い致します。