- ベストアンサー
微分方程式の数値解法
chukanshiの回答
- chukanshi
- ベストアンサー率43% (186/425)
>やっぱりそういう方法しかないですよね。 >なんか超絶技巧な技があるのかと思いましたが。 あとで、自分で回答を読んで、わかりにくかったかな、と思いましたが、 わかっていただけてよかったです。 というよりも、誰でも真っ先に思いつく方法ですよね。 英語で、「shooting method」とちゃんと名前がついていること自体 始めは驚きでした。 でも、周囲の数値計算の専門の仲間に聞いても、これ以外の方法は みんな知らないって言うので、私は論文(理論物理です)の 微分方程式の数値計算をこの方法でやって、英文専門誌に投稿して受理 されましたから、認められた方法だと今では思っています。 tgbさんのように、なるほど、両側から解いていって、間で継ぎ合わせる という方法もありますよね。 収束の速さや、解の精度を挙げるには、いろいろと検討しなければ ならないかもしれません。 私もいまだに、超絶技巧があるのでは?と調べたり、研究していますが、 新しい方法は、見つかりません。 このページに画期的な方法が現れるのを期待していますし、 ibm_111さんが、新しい方法を手に入れられたら、こちらが 教えていただきたいほどです。
関連する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) の問題が解くことができませんでした。 どなたか解法をわかりやすく教えていただけないでしょうか?
- ベストアンサー
- 数学・算数