• 締切済み

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

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

みんなの回答

  • alice_44
  • ベストアンサー率44% (2109/4759)
回答No.2

数学の問題としては、  1、オイラー法を使う。  2、ルンゲ・クッタ法を使う。 で終わりです。 解法を指定されているのだから、 その手順に従って作業するだけで、 頭を使う箇所は全くありません。 馬鹿でもコンピュータにでもできることなので、 プログラムに書きあらわすことが可能なのです。 貴方の言う「解く」が、プログラムを書くこと を指しているのであれば、 質問するカテゴリーが違います。

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

「どのように」もなにも, そこに書いてあることが全てですが. 当然「オイラー法」とか「ルンゲクッタ法」は完全に理解できているんでしょ?

関連するQ&A

  • 数値解析の問題

    dx/dt=1-x^2をオイラー法とルンゲクッタ法で0<t<2.8の区間におけるX(t)がわからないです。わかるかた是非教えてください。

  • 常微分方程式の問題

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

  • 単振動の微分方程式を刻みhについてルンゲクッタで求める。

    m*d^2x/dt^2=-kx x(0)=1 dx/dt=0 というのが与えられて二階微分だから一階微分にするために dx(t)/dt=v(t) dv(t)/dt=-k*x(t)/mという式を立てました。オイラー法ではできたのですが2次、4次のルンゲクッタだとできません。どなたか回答お願いします。

  • 数値解析の刻み幅

    偏微分方程式を数値解析にて解く場合, 差分法を使って解こうとすると, dtやdxを設定してやる必要があると思います. たとえば, δa/δt = -δa/δx + a*cos(x) + (a^2+b^2)*a δb/δt = -δb/δx + b*cos(x) + (a^2+b^2)*b っていう方程式を解こうと考えた場合, (仮に初期条件を孤立波としておく) とりあえず線形項を考慮して,dt/dx<<1 は満たす必要があると思いますが,非線形項を含む場合は どのようにdtやdxを決定してやればいいのでしょうか? dxは小さければ小さいほど良いというのは聞いたことがあるのですが・・・ テキトーに決定してやるのでいいのでしょうか?

  • 数値解析 微分方程式

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

  • オイラー法、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つです。 休んでいた自分が悪いのですが、ネットで調べてもよくわからなくて… わかる方、よろしくおねがいします…

  • オイラー法の誤差

    「オイラー法の誤差を最小にするtの刻み幅」とは、どうやって求めるのでしょうか?丸め誤差の求め方も知らず、どうしたらいいかわかりません。 与えられている式は dx/dt=-2t^3+12t^2-20t+8.5 です。

  • 微分方程式・・・。

    問題はdt/dx=x^2 , x(0)=a この微分方程式の解x(t)が任意の時刻t>0までに存在するまでの初期値aの満たすべき条件を求めよ。 という問題です。オイラー法を使ってやるのかなってとこぐらいまでしかわかりません・・。dx/dt=f(t,x)とx(to)=xo から先に進みません。 アドバイスいただけませんか。

  • 解析学の問題

    次の一般解を求めよ。 d**2x/dt**2 - 4dx/dt + 6x = 0 こういった問題が宿題に出されました。 そこで積分をして ∫(d**2x/dt**2)dt - ∫(4dx/dt)dt + ∫(6x)dt = C ↓ dx/dt - 4x + 6tx = Ct ↓さらに積分して x - 4tx + 3(t**2)x = 1/2*(Ct**2) Cは全ての定数なので、0でもいい。 よって x - 4tx + 3(t**2)x = 0 ここまで考えてみたのですが、あっている自信がほとんどありません。 よろしければどうすればいいか教えていただきたいです。 よろしくお願いします。

  • オイラー法

    ①dx(t)/dt=1-x²(t),Δt=1/10,x(0)=0であるときx(0.3)の値をオイラー法とエクセルを用いてグラフにせよ オイラー法を使いエクセルで上記の微分方程式を解きたいのですがエクセルの使い方がわからず悩んでいます 教えてください、お願いします