• 締切済み

数値解析　微分方程式

数値解析における常微分方程式を解くために用いる手法についてです。 オイラー公式、ホルン公式、ルンゲクッタ、ルンゲクッタ4次、有限差分法の関係 違いがよくわからなのでどなたか教えてください また常微分方程式を有限差分で解くとなったとき、結局オイラー公式などを使うと言う認識で間違いないでしょうか？それとも有限差分だけで解けるのでしょうか

みんなの回答

  • f272
  • ベストアンサー率46% (8005/17110)
回答No.1

有限差分法って偏微分方程式の解法として有名ですが，常微分方程式を解くのに使うという説明は聞かない。

  1. カテゴリ
  2. 学問・教育
  3. 数学・算数

関連するQ&A

  • 高階連立常微分方程式の数値計算

    4次のルンゲクッタ法を用いた数値計算を勉強しています． 1階連立常微分方程式と高階常微分方程式は理解でき，プログラムも作成することができました． 次に高階の連立常微分方程式を解こうと思ったら，頭が混乱してしまいました． 4次のルンゲクッタ法を用いて高解連立常微分方程式を解く考え方を教えて頂ければ嬉しいです． また何か良い参考書があれば教えて頂きたいと思います． よろしくお願いします．

  • 常微分方程式の数値計算

    実験で常微分方程式の数値計算をしました y'(x)=2xy,y(0)=1.0(解はy=ex2）←eのx二乗の微分方程式できざみ幅ｈ＝0.1、区間〔0，2〕の条件で オイラー法と改良オイラー法とルンゲクッタ法を使って 表計算ソフトを使って求めよとあるのですが、 表計算で近似値のy座標とか　Δｔ後のxの変化分K1とか求めるのですが、全然意味がわからなくて、 何かいい方法はないでしょうか？

  • 常微分方程式の問題

    微分方程式x'=x^2-t-1で初期値がx(0)=1の問題で、 オイラー法とルンゲクッタ法は出せたのですが、 それらと比較するために出す解析解がわかりません。 どなたかよろしくお願いします。

  • 差分法とオイラー法の違いについて

    最近微分方程式の数値解析について学びだした者です。 微分方程式の数値解法として差分法とオイラー法があると思うのですが、この2つの違いや互いの位置づけはどうなっているのでしょうか？ また、差分法には風上法などがありますが、これらとオイラー法の位置づけについても教えていただきたいです。 できればこれら近辺の全体的な体系について教えていただけるとうれしいです。よろしくお願いします。

  • オイラー法、２次ルンゲクッタ法、４次ルンゲクッタ法のＣ言語プログラムに

    オイラー法、２次ルンゲクッタ法、４次ルンゲクッタ法のＣ言語プログラムについて教えてください！ 課題なのですが、まったくわからず困ってます＞＜ 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つの近似精度を比較せよ。 以上の３つです。 休んでいた自分が悪いのですが、ネットで調べてもよくわからなくて… わかる方、よろしくおねがいします…

  • 微分を数値的に積分して近似解を求める方法

    du/dt=f(u,t) というような感じの単純な常微分方程式を逐次的にコンピュータで数値積分していく場合、単純陽解法とかルンゲクッタ法とか予測子・修正子法、台形公式などいろいろあると思います。あるいは反復計算を含むものなどもあると思います。これらについてグラフを作図して解析解との差を評価したりするわけですが、偏微分方程式で使う場合も全く同じ評価がそのまま当てはまるでしょうか。常微分方程式での知見がそのまま偏微分方程式の計算に当てはまるのか、ということなのですが。反復計算を用いる方法で常微分方程式では芳しくない結果だったのですが、偏微分の場合は反復計算の価値が高くなるのではないかと思ったりするものですから。反復計算、すなわち繰り返しをやるうちに陰解法的になっていくからなのですが。 いかがでしょうか。細かい内容の質問で恐縮ですが、よろしくお願いします。

  • 偏微分方程式の数値解析

    お忙しいところすみません． 現在，FDTD法による数値解析を行っています． FDTD法は空間及び時間微分に対して中心差分を行っていますがこれを4次のルンゲクッタ法で行うことは出来るのでしょうか？ある参考書に以下のような「連立1階微分方程式」を4次のルンゲ-クッタ法で解くというのがありました． dy/dt=f(t,y,z) dz/dt=g(t,y,z) これを4次のルンゲ-クッタ法で解くと y(n+1)=y(n)+1/6*(k1+2*k2+2*k3+k4) z(n+1)=z(n)+1/6*(j1+2*j2+2*j3+j4) ここで k1=Δt*f(t(n),y(n),z(n)) j1=Δt*g(t(n),y(n),z(n)) k2=Δt*f(t(n)+Δt/2,y(n))+k1/2,z(n)+j1/2) j2=Δt*g(t(n)+Δt/2,y(n))+k1/2,z(n)+j1/2) k3=Δt*f(t(n)+Δt/2,y(n))+k2/2,z(n)+j2/2) j3=Δt*g(t(n)+Δt/2,y(n))+k2/2,z(n)+j2/2) k4=Δt*f(t(n)+Δt,y(n)+k3,z(n)+j3) j4=Δt*g(t(n)+Δt,y(n)+k3,z(n)+j3) これをマクスウェルの方程式に対応させ ε∂E(t,x,y,z)/∂t=f(t,H)=∇×H μ∂H/∂t=g(t,E)=-∇×E として出来ないこともない気がするのですが． なぜこの方法を試みようとした理由は現在，通常のFDTD定式化により数値解析を行っていますが ・プログラムは正しい ・始めの部分は正常に計算出来ている ・途中から凹凸が出始め，それが成長して発散している という事態になり，数値的不安定の可能性が強いと思ったからです．オイラー法やルンゲ-クッタ法ではこのような異常振動は進み幅を小さくすれば止まる(ということが数学的に証明されている)ようです． 何卒回答の方，よろしくお願い致します．

  • この微分方程式は解析的に解けるでしょうか?

    下記の非線形常微分方程式を解析的に解くための方法は在るのでしょうか? それとも、数値計算でなければ無理ですか? 2・y '' ・y ' + x^2 = 0 y(0) = 0 , y(1) = 1 y = y(x) (yはxの関数) , y 'はxでの一回微分、y ''は二回微分です。

  • 数値解析についての質問です。お願いします

    以下の問題をプログラム、あるいはエクセルを用いて解きたいのですがどのように解いたらよいのかわかりません。お願いします dx/dt=1/2(1-x),x(0)=0.1の解を以下の方法で求めよ 1、オイラー法を用いて求めよ。時間刻みhを0.5、時間範囲を0<=t<=20とする。解析値、解析解、相対誤差を表にまとめよ。数値解、解析解の時間的変化をグラフで示せ（解析解）= 1/1+9e^-0.5t 2、オイラー法をルンゲクッタ法に変え同条件で解け

  • 数値微分法についてです。

    数値的に微分を評価する時に、中心差分を使っているのですが、 どう考えても数値誤差としか思えない結果しかでません。 区切り幅は誤差が最も小さくなるように選んでいます。 中心差分よりも精度の良い数値微分法があれば教えていただけないでしょうか。 評価する関数は解析的に与えられておらず、補間して得られるようなものです。 (補間の精度にもよるのだと思いますが・・・)

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する