• ベストアンサー

常微分方程式の問題

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

質問者が選んだベストアンサー

  • ベストアンサー
  • KappNets
  • ベストアンサー率27% (1557/5688)
回答No.1

難しそうな微分方程式ですね。mathematicaを使ってみましたが、 AiryAiPrime、AiryBi、AiryBiPrime、AiryAi などの関数が出てきます。Primeは微分形ですし、早い話解けないもようです。如何でしょうか。

deedeedee
質問者

お礼

mathematica使っても出ないなら、あきらめるしかないようですね。 助かりました、ありがとうございました。

関連するQ&A

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

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

  • 偏微分方程式と常微分方程式

    物質濃度をC、時間をt、座標をx、物質の分子拡散係数をνとすると分子拡散による物質濃度の時空間変化は以下の偏微分方程式によって記述される。これについて以下の問いに答えよ。 ∂C/∂t=ν((∂^2)C/∂x^2) (1)C(x,t)=X(x)T(t)と仮定することにより、X(x)およびT(t)に関する常微分方程式をそれぞれ導出せよ。 (2)(1)での2つの常微分方程式の一般解をそれぞれ求めよ。 (3)上記拡散方程式は一般に放物型と言われる偏微分方程式に分類される。これとは別の楕円型と言われる偏微分方程式を1つ、数式で記述せよ。 困っているのは(2)の問題です。 以下のようなwebサイトを見つけました。 http://www12.plala.or.jp/ksp/mathInPhys/partial/ これに沿って問題を解いていったとき、一般解をどのようにするべきか迷いが生じました。今回の問題では初期条件や境界条件はないため、一般解はλが正、ゼロ、負のとき全ての場合の一般解を求めなければならないということですか? 後もう1点、もしよければ、楕円型の微分方程式として有名な物理現象、あるいは式を教えていただけないでしょうか? ヨロシクお願いしますm(_ _)m 特に(2)の問題に関する質問、ヨロシクお願いします。。。

  • 常微分方程式の問題

    常微分方程式の問題です。 dx/dt =x^3 + e^-t^2 (初期条件x(0) =1)の解を求めよ。という問題なのですが、御回答よろしくおねがいします

  • 微分方程式・・・。

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

  • 微分方程式の問題2

    微分方程式x''-4x'+20x=0 の初期条件x(0)=1、x(0)=6を満たす解x=x(t)をあらわす式を求めよ

  • 数値解析 微分方程式

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

  • 1回の常微分方程式解法について

    常微分方程式df/dt=g(t)で、関数fとtの数値的な関係を求める場合、ルンゲクッタ法が用いられることが多いと思います。特にtが時間の場合、g(t)は未来には影響しますが、過去には影響できないので解法としてもルンゲクッタ形式の解法になることはわかります。計算が一方向に進むというイメージです。 df/dx=g(x)でxが空間の関数だった場合、g(x)がxの±領域全体に影響を与えると考える場合、反復計算によって全体に影響しながら解を求めていくことは可能でしょうか。ただし、微分は1回しかありません。 2回微分だとSOR法のような反復計算が可能です。全体が影響を与え合いながら式を満たすfを求めていくことができますが、1回微分でもそのようにできるでしょうか。 独立変数がx(空間)に限定ですから時間という概念はありませんが、時間が含まれないということは時間的に変化しなくなった解と見ることもできます。 例えば、シューティングメソッドとか言うのだろうと思いますが、ルンゲクッタ法で一方向に向かって計算し、計算が領域の端部になったらその結果と境界条件を比較して収束が足りないならば少し変化させて再度ルンゲッタで計算する、すなわちルンゲクッタ法の収束計算というやり方もあるようです。 ルンゲクッタ法で計算するような問題を反復計算で計算領域がお互いに影響を与えながら計算する方法があったら教えて頂きたいのですが。境界条件が計算領域の両端にある場合、2点境界値問題というのでしょうか。これを解く方法ということにもなろうかと思います。手前の境界で与えられたものから出発し、もう片方の境界で計算結果と境界値を見合わせることになるのかもしれませんが。 よろしくお願いします。

  • 常微分方程式の問題で。

    常微分方程式の問題で オイラー法、ルンゲ・クッタ法、ルンゲ・クッタ・ギル法、ミルン法、アダムス法より2つの方法を用いて解け。 y''+y'/x+y=0 初期条件y'/x=-1/2 0<=x<=4 h=0.4 という問題なのですが、基本的に解き方がわかりません。 二階微分なので連立ルンゲ・クッタで解くのかな?っと思ったのですが 結局それからどうすればいいのかわかりません。 2つの方法とありますが、これはどの方法でも解けるものなのでしょうか? 本当はC言語のプログラムを組まなければいけないのですが、 この理論がわからないのでおねがいします。

  • 2階常微分方程式

    x*y''-y'=0 この微分方程式を解きたいのですが解法がわかりません・・・ y"の係数がx^2であればオイラーの微分方程式としてy=x^αを代入して解を得ることができますが この形ですとα=0となってしまい解を求めることができません。 このような定数係数以外の時の2階常微分方程式の解の求め方について説明をお願いします。

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