- ベストアンサー
微分方程式の数値解法
motsuanの回答
- motsuan
- ベストアンサー率40% (54/135)
基本的に境界値問題のことをおっしゃっているのだとおもいますが 線形であれば 線型方程式を解く方法もあると思います。 Σa_n(x) (∂ /∂x)^n f(x) = C というようなものであれば、結局、各場所の値 f(x_j) (j=1,2,...,N)を未知数として、 方程式から境界条件を除いて適当な数の連立方程式がえられますよね(適切な問題であれば)。 たとえば、端の点の値が与えられているのであれば、1階の微分方程式であれば完全にきまって、 2個の境界条件と微分方程式を差分式になおしたN-2個の連立方程式が得られて 結局、線形連立方程式に帰着できることがわかります。 でも、線形じゃないばあいはやっぱり「いろいろやってみる」方式のようで、 Mathematica などでは(マニュアルによると) 非剛性アダムズ(Adams)法と剛性ギア(Gear)法とを交互に切り替えるのだそうです。 http://www.wolfram.com/products/calculationcenter/features/solvers.ja.html ちなみに 線型境界値問題の場合はゲルファンド・ロクチェフスキー(Gel'fand-Lokutsiyevskii)追跡法 だそうです。 確認できなかったのですが、方程式には硬い軟らかいのがあるそうで http://grape.c.u-tokyo.ac.jp/~makino/kougi/system_suuri4_1999/note8/node1.html 局所的にそれを切り替えて収束をはやくしているようです。 非線形数値計画法とかの(準)ニュートン法の微分方程式版みたいなものなのでしょうか? あと、L(f(x) )=0を微分方程式として 発展方程式 (∂/∂t ) f(x) = L(f(x) ) +「境界では境界を満たすようにする」 を考えて、(良く収束するとして)、 発展方程式の時間を進めて 定常解 (∂/∂t ) f(x) =0を解とするという考え方もありますよね。
関連するQ&A
- 偏微分方程式の数値解法
偏微分方程式の込み入った質問です。 2次元(x,y)の空間で2つの関数f(x,y),g(x,y)を考えます。 そこで、それぞれにラプラス方程式を立てました。 fxx+fyy = 0 (1) gxx+gyy = 0 (2) です。これは境界値問題で、差分式からSOR法を使って収束計算によって数値解を求めることができます。f, gはそれぞれ独立という形にはなります。 そこにもう1つ式が出てきました。 fxfy + gxgy = 0 (3) というものです。f,gをx,yで1回微分してできる式です。 都合3つの式が出てきました。 この数値解を求めるにはどのような方法があるでしょうか。 数値解ですから近似解です。 3つ目の拘束条件の下でのラプラス方程式とみると、ペナルティ関数とかラグランジュの未定係数法とかいろいろあるかもなと思いますが。 3つ目の式は完全に満たすというより、できるだけ満足するようにしたいというものです。 よろしくお願いします。
- 締切済み
- 数学・算数
- 微分方程式の解法について
微分方程式の解法について X^3-3xy+y^3=0 ならば、dy/dxはいくつになるか どうやったら、いいんでしょうか。 解法を教えてください。
- ベストアンサー
- 数学・算数
- 微分方程式の解法について
以下の微分方程式の解法が分かりません。y '=pと置いて変数分離で解くことを試みたのですが、行き詰まってしまいました…。 2xy"=√{1+(y ')^2}、y(1)=y '(1)=0 回答よろしくお願いします。
- ベストアンサー
- 数学・算数
- 微分方程式の数値解法 1段解法収束性の証明
微分方程式の1段解法の収束性を示す過程で、 h:刻み幅 L:リプシッツ条件を満たす定数 としたとき、(1+hL)<e^(hL)という式が出てきていました。 どの本を見てもこの式が使われているのですが、どうしてこの関係が言えるのかがわかりません。 まだ数値解法を勉強し始めたばかりで知識が足りないだけかもしれませんが、分かる方、ぜひ教えてください! お願いします(><)
- ベストアンサー
- 数学・算数
- 常微分方程式の問題
常微分方程式の問題でいくつか解けなかったところがあるので教えていただきたいです。 この章で扱っているのは 変数分離系・同時系・線形1階微分方程式・完全微分形・線形2階微分方程式(同次形)・線形2階微分方程式(非同次形) を扱っていました。 その内、一般解を求める以下の問題 (1)dy/dx=xe^-y (2)x(dy/dx)-y=1 (3)(2y-x^2)dx+(2x-y^2)dy=0 と 与えられた条件をそれぞれ満たす微分方程式の解を求める以下の問題 (1)dy/dx=y/x (x=1のときY-2) (5)y''+5y'+6y=0 (x=0のときy=0、y'=1) の問題が解くことができませんでした。 どなたか解法をわかりやすく教えていただけないでしょうか?
- ベストアンサー
- 数学・算数