• 締切済み

数学 微分方程式

微分方程式 d^2u/dx^2=1 0<x<π/2 境界条件u(0)=0、u(π/2)=0 を有限差分を用いて解く問題なのですがどなたか教えてください

みんなの回答

  • f272
  • ベストアンサー率46% (8010/17118)
回答No.3

> 手動でやる場合解である関数をある程度予想してやるしかないのでしょうか。 どのような関数になるのかわかっているのなら、数値計算で答えを出すのは馬鹿らしいです。分からないから数値計算で出そうとするのです。 図はエクセルで刻み幅π/20で求めた例です。高々10分割でも正解に一致しますね。

  • f272
  • ベストアンサー率46% (8010/17118)
回答No.2

d^2u/dx^2=1 を見て,2回微分したら定数1になるのだから2次の係数が1である2次関数で, 境界条件u(0)=0、u(π/2)=0を見て2点が決まるのだから,2次関数を特定できます。

19194545
質問者

お礼

ありがとうございます。一つ質問なのですが、uに関数差分方程式を立てた後各点における値を出せ、プロットまで行ったのですが、手動でやる場合解である関数をある程度予想してやるしかないのでしょうか。 そもそもコンピュータで細かくプロットすることで成立するものだとは思いますが

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

ぱっと見てu=(1/2)x^2-(π/4)xと思うが,いかが。

19194545
質問者

補足

どのように考えたのでしょうか

関連するQ&A

  • 微分方程式について

    L(U)=d^2U/dx^2 - U =0(0<x<1) 境界条件U(0)=0 U(1)=1 この微分方程式を1次元2次要素を用いて解くという問題がわかりません。 お願いします。教えてください。

  • 微分方程式について

    次のような微分方程式があります d^2 x/dx^2 - (dy/dx)(4+x)/x +y*(6+2x)/x^2 =0 問題は以下です y=ux^2(uはxの関数)がこの微分方程式の解となるために uの満たすべき微分方程式を求めなさい。 要は u''=u'=u になればいいということじゃないのでしょうか ですがこれだと微分方程式になりません もしくはこれが解答でいいのでしょうか? ヒントのみでもいいので教えてください。

  • 微分方程式の解法

    d^2y/dx^2+2*x*dy/dx=0 境界条件 x=0: y=1、x→∞: y→0 この2階の微分方程式を解けという問題ができません。 dy/dx=z と置いて、1階の微分方程式にして解こうとしたのですが、exp(-x^2)が出てきてしまいました。これは確率積分みたいに積分できるのでしょうか。 回答よろしくお願いします。

  • 微分方程式

    微分方程式 d^2y/dx^2=2y^3+2y をみたす関数y(x)を求めよ。ただし,境界条件は,y(0)=0,dy/dx(x=0)=1

  • 微分方程式

    こんにちは。微分方程式の問題が解けなくて困っています。 次のx(t)に関する微分方程式 d^2x/dt^2=-1/x^2 ただし初期条件はt=0でx=X0(x0>0),dx/dt=√2であるとする。 問題 与式の両辺にdx/dtを乗じて積分することにより、初期条件を満たすxについての1階微分方程式をもとめよ。 必要ならば、公式d/dt(dx/dt)^2=2*(dx/dt)*(d^2x/dt^2) 少し問題の書き方がおかしいかもしれませんが(微分の書き方)どなたかお願いします。 自分なりにといたのですが 与式の両辺にdx/dtをかけて dx/dt(d^2x/dt^2)=-1/x^2*(dx/dt) 与えられた公式をつかい (1/2)*d/dt*(dx/dt)^2=-dx/dt*(1/x^2) ∫(1/2)*d/dt*(dx/dt)^2=-∫dx/dt*(1/x^2) ????? と与えられたヒント通りにしてそこからどうしたらいいのかわからなくなってしまいました・・・

  • 偏微分方程式の解き方を教えていただけないでしょうか

    偏微分方程式の解き方を教えていただけないでしょうか。 u_t (tの一階微分) = u_xx (xの二階微分) x∈[0,1]のとき、 境界条件 u_x(0,t)=0 、u(1,t)=5t (↑xの一階微分) 初期条件が、 u(x,0)=0 自分で _____________________ du/dt = d^u/dx^2 x∈[0,1] du/dx(0,t)=0 、u(1,t)=5t u(x,0)=0 のとき、変数を分離して、 u=(X,Y) X''=-λXとしました。 X=c1 cos(√(λ) x) +c2 sin(√(λ) x) として、 X’=√(λ) *(ーc1 sin(√(λ) x) +c2 cos(√(λ) x) ) 境界条件をいれると、 X’(0)=√(λ) *(ーc1 sin(√(λ) 0) +c2 cos(√(λ) 0) ) より c2=0 X(1)=c1 cos(√(λ)*1) +c2 sin(√(λ)*1) =5t c1*cos(√(λ)*1) =5t ____________________________ と計算をしてみたのですが、5tの扱い方がわからず、躓いてしまいました。 どのように計算をすればよいか、教えていただけないでしょうか。

  • 微分方程式

    (d^2)φ/dx^2=eNa/ε0εs 境界条件 x=-xp:dφ(x)/dx=0  (pはxの添え字です) x=0:φ(x)=0   これをφについて微分方程式を解きたいのですが… ちなみに解答は φ(x)=eNa/ε0εs{(x^2/2)+xpx} (pはxの添え字です) です。お願いします。

  • 微分方程式

    微分方程式 u=f(x) d^4 u/dx^4=u、f(0)=0 uを求める

  • 固有値問題 微分方程式

    固有値問題です。 (d^2/dx^2)*u(x)=λ*u(x) 境界条件 du(0)/dx-u(0)=0 du(1)/dx-u(1)=0 (d^2u(x)/dx^2)→u(x)をxで二階微分 u(x)→xを変数とする関数(固有関数) λ→固有値 という問題をλが正の時、0の時、負の時にわけて解きたいのですが解き方がわかりません。 よろしければ教えて下さい。

  • 2階微分方程式の解法

    d^2y/dx^2+2*x*dy/dx=0 境界条件 x=0: y=1、x→∞: y→0 上記の微分方程式をルンゲクッタ法を使って数値的に解きたいのですが、どのように x→∞ : y→0 の境界条件をいれればいいかわかりません。どなたか教えて頂けないでしょうか。回答よろしくお願いします。