偏微分方程式の問題と境界条件についての質問

このQ&Aのポイント
  • 偏微分方程式の問題についての質問です。境界条件や解法について教えてください。
  • 求めるべき偏微分方程式の途中式がわかりません。具体的な計算手順を教えてください。
  • 以前の質問で得た解答に関して、特定の途中式が理解できません。詳細な説明をお願いします。
回答を見る
  • ベストアンサー

偏微分方程式の問題です。

・問題 u=u(t,x) u_t=ku_{xx} k>0 (0<x<1,t>0) 初期条件 u(0,x) sin(πx)+1/2sin(3πx) 境界条件 u(t,0)=u(t,1)=0 ・答え 境界条件より正弦級数展開をする. (☆)u(t,x)=Σ_{n=1}^∞b_n(t)sin(nπx) (★)b_n(t)=2∫_0^1u(t,x)sin(nπx)dx すると, u_t=Σ_{n=1}^∞{db_n(t)/dt}sin(nπx) ku_{xx}=Σ_{n=1}^∞b_n(t)d^2{sin(nπx)}/dx^2 =Σ_{n=1}^∞b_n(t)k(-nπ)^2sin(nπx) u_t=ku_{xx}より db_n(t)/dt=b_n(t)k(-nπ)^2=-n^2π^2kb_n(t) ∴b_n(t)=b_n(0)e^{-n^2π^2kt} ☆に代入して (☆☆)u(t,x)=Σ_{n=1}^∞b_n(0)e^{-n^2π^2kt}sin(nπx) ∴u(0,x)=Σ_{n=1}^∞b_n(0)sin(nπx) これと初期条件u(0,x)=sin(πx)+(1/2)sin(3πx)を係数比較して b_1(0)=1,b_3(0)=1/2,b_n(0)=0(n≠1,3) ☆☆に代入して ∴u(t,x)=e^{-π^2kt}sin(πx)+(1/2)e^{-9π^2kt}sin(3πx)(答) ーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーーー 以前、上記の問題をここで質問し、解答を頂いたのですが、この部分の途中式がどうしてもわかりません。 ku_{xx}=Σ_{n=1}^∞b_n(t)d^2{sin(nπx)}/dx^2 =Σ_{n=1}^∞b_n(t)k(-nπ)^2sin(nπx) お手数ですが宜しくお願い致します!

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

  • ベストアンサー
回答No.1

以前回答したものです. >ku_{xx}=Σ_{n=1}^∞b_n(t)d^2{sin(nπx)}/dx^2 >=Σ_{n=1}^∞b_n(t)k(-nπ)^2sin(nπx) この1行目の右辺にkがかかるところを抜かしていました. ku_{xx}=Σ_{n=1}^∞b_n(t)kd^2{sin(nπx)}/dx^2 となります.失礼しました.

longgoodbye
質問者

お礼

何度もすみません。ありがとうございました!

関連するQ&A

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

    偏微分方程式の解き方を教えていただけないでしょうか。 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の扱い方がわからず、躓いてしまいました。 どのように計算をすればよいか、教えていただけないでしょうか。

  • 偏微分方程式の問題です。

    下記の問題の解き方をどうか教えて下さい。 u=u(t,x) ut=kuxx k>0 (0<x<1,t>0) 初期条件 u(0,x) sin(πx)+1/2sin(3πx) 境界条件 u(t,0)=u(t,1)=0 よろしくお願いします!!!

  • 偏微分方程式の問題について質問です。

    偏微分方程式の問題について質問です。 http://www.mech.tohoku.ac.jp/examination/grad/pdf/Problem2007.pdf の数学Bの1.の問題で、 (1)は、g(x)=-1/3x^3+1/3xとなった(多分、合っていると思う)のですが、 (2)で、 u(x,t)=v(x,t)+g(x)において、g(x)は求まっているのでv(x,t)を求めるために、 ∂v/∂t=∂^2v/∂x^2 の式を XT'=X''T X(0,t)=X(1,t)=0 として、X=Asin(nπ)x T=Cexp(-nπt) となり、重ね合わせから v(x,t)=?An*sin(nπ)x*exp(-nπt) よって、u(x,t)=?An*sin(nπ)x*exp(-nπt)-1/3x^3+1/3x u(x,0)=x^2(1-x)/3を使って、Anを求めればいいのか、とも思ったのですが、 ここで止まってしまいました。 どうすればいいのでしょうか。 よろしくお願いします。

  • 偏微分方程式

    閲覧ありがとうございます。 学校の課題で出された問題です。 【問題】 ∂^2U/∂x^2 = ∂^2U/∂t^2 を解け。(Uはx、tで定義されている関数です) なお、初期条件は ∂U(0,t)/∂x = ∂U(π,t)/∂x = 0 U(x,0) = 2+cosx ∂U(x,0)/∂t = x とする。 スキャナが使えない上、自分の字が汚い事で読みにくい事この上ないのですが、画像の通り一応解らしきものは求まりました。 【解】 U = 2+cosx*cost+Σ(Bn*cos(nx)*sin(nt)) (n : 0~∞) ただし n=0   → n*Bn = π/2 n=2m  → Bn = 0 n=2m-1→ Bn = -4/(π*n^3) (m:自然数) ここで、『n=0   → n*Bn = π/2』と、自分で書いたものの 実際にn=0を代入すると『0*B0 = π/2』となり、おかしいですよね? おかしいとしたらどこの計算が間違っているのかを指摘してほしいです。 根本的に間違っているようでしたら模範解答をいただけるとありがたいです。 長文失礼しました。

  • 応用数学、熱伝導方程式について

    大学でこのような問題がでたのですが、解き方や考え方を教えて下さい。 問、次の方程式の解を求めよ。 ∂u/∂t = 2 ((∂^2 u)/ (∂x^2))・・・(1) (0<x<2, t>0) u(0,t)=u(2,t)=0・・・(2), u(x,0)=5 sin 2πx - 3 sin 5πx・・・(3) 以下は自分でといた過程です。 (1)より X(x)*T'(t)=2*X"(x)*T(t) T'(t)/T(t) = 2X"(x)/X(x)=Κと置く。 これより T' - ΚT=0・・・(4) 2X" - ΚX = 0・・・(5) (Κ<0) (5)の解はΚ<0より X = A cos(√(-Κ/2))x +B sin(√(-Κ/2))x ・・・(6)(A,Bは任意定数) 境界条件(2)より X(0)=X(2)=0・・・(7) (T(t)は0ではない) よって(6)、(7)より A =0, Κ= -(n^2*π^2)/2 X=B sin(nπ/√2)x (4)よりT' = ΚT    ↓ T = D e ^(-(n^2 * π^2)/2)t (Dは任意定数) 次に偏微分方程式の解を求める。 →ここから先が全くわかりません。 正しい答えは5e^(-8(π^2)t)2πx - 3e^(-50(π^2)t)sin5πxになります。

  • 偏微分方程式の初期値問題

    x(ux)+(uy)=0 u(ξ,0)=sinξ (ux),(uy)はuのx,yそれぞれについての偏微分です。 という問題なんですが、 dx/dt=x dy/dt=1 (x,y,u)lt=0=(ξ,0,sinξ) x=ξe^t y=t となりますよね? そこから初期値の無い問題と同様に定数を作ると、 x*e^-y=ξ でu=f(x*e^-y)という答えになってしまったんですが、 これなら初期値の無い時と一緒の答えですよね? どこで初期値を使ってやればいいのでしょうか?

  • 1次元熱伝導方程式

    以下の問題で参考書に解法が載ってなかったので悩んでいます。 両端が熱的に絶縁されている場合 u_t=u_xx (0<x<1, 0<t<∞) u_x(0, t)=0 (0<t<∞) u_x(1, t)=0 (0<t<∞) u(x, 0)=x (0≦x≦1) 基本解まで求めたところ、 u_n(x, t)=a_n*exp(-(n*pi)^2*t)*cos(n*pi*x) となりました。 参考書には右辺のcosではなくsinの例しか載っておらず、 固有関数展開の係数a_nを求めるのに安直にt=0として両辺にcos(m*pi*x)をかけてxについて0から1まで積分したところ、 u(x, t)=Σ(-4/((2n+1)^2*pi^2))*exp(-((2n+1)*pi)^2*t)*cos((2n+1)*pi*x) nは0から∞まで となりました。 しかし、この解が直感と異なる(定常解が0なのはおかしい)ので困っている次第です。 実際、u(0, 0)=-0.5 u(1, 0)=0.5 なのでちょうど初期条件から0.5ずれています。 原因はやはり固有関数展開の係数を求めるところにあるのでしょうか。 どうか解説をお願いします。

  • 微分方程式の固有値問題

    僕は物理系の学生なので、なんとなくそれっぽい雰囲気の質問になってます。 ふと気になったのですが、例えば固有値問題 -(d^2/dx^2)ψ(x) = kψ(x) ただし、k>0 0<x<L (複素)内積<ψ(x),ψ(x)> = 1 の解は ψ(x) = (1/L)sin(kx) k = nπx/L nは(0を含まない)自然数 となります。 一方、 {-(d^2/dx^2)+x^2}ψ(x) = εψ(x) ただし 0<x<L (複素)内積<ψ(x),ψ(x)> = 1 の解は ψ(x) = cHn(x)exp{-(x^2)/2} ε = 2n + 1 ただし、cは規格化定数、Hn(x)はエルミート多項式、nはゼロ以上の整数です。 という解になることはわかるのですが、 ψ(x) = (1/L)sin(kx) k = nπx/L + x^2 nは(0を含まない)自然数 としても、解になるのではないでしょうか。 長くなりましたが、質問としては 『固有値に変数xが入ってはいけないのか』ということです。 よろしくお願いしますm(__)m

  • 定数係数の2階線形微分方程式の問題です。

      x'' + 2x' + 5x = 2cos(3t), x(0) = 1, x'(0) = 0.・・・・・・・(1)   k^2 + 2k + 5 = 0.  k = -1±2i.  x'' + 2x' + 5x = 0 の一般解 x0 は   x0 = e^(-t)(Acos2t + B2sin2t).  x を複素数と見なして(1)を書き直すと   (D^2 + 2D + 5)x = 2e^(3it).   1/(D^2 + 2D + 5)[2e^(3it)]  = 1/( (3i)^+ 2(3i) + 5 )[2e^(3it)] = 1/(-9 + 6i + 5)[2e^(3it)]  = 1/( -4 + 6i)[2e^(3it)] = 1/(-2 + 3i)[e^(3it)]  = -(2 + 3i)/13[e^(3it)] = -(2+3i)/13(cos3t + isin3t)  = (-2/13)(cos3t + isin3t) - (3i/13)(cos3t + isin3t)  = (-2/13)cos3t + (-2/13)isin3t - (3i/13)cos3t + (-3i/13)isin3t  = (-2/13)cos3t + (-2/13)isin3t - (3i/13)cos3t + (3/13)sin3t  = (-2/13)cos3t + (3/13)sin3t) + (-2i/13)sin3t - (3i/13)cos3t  = (-1/13)(2cos3t - 3sin3t) - (i/13)(2sin3t + 3cos3t ).  (1)の特解を x1 とすると   x1 = (1/13)(3sin3t - 2cos3t).  よって(1)の一般解 x は   x = e^(-t)(Acos2t + Bsin2t) + (1/13)(3sin3t - 2cos3t).   x' = -e^(-t)(Acos2t + Bsin2t) + e^(-t)(-2Asin2t + 2Bcos2t)                  + (1/13)(9cos3t + 6sin3t).   x(0) = A - 2/13 = 1.  A = 15/13.   x'(0) = 15/13 + 2B + 9/13 = 0.  2B = - 24/13. B = -12/13.  求める特解を改めて x とおくと   x = (e^(-t)/13)(15cos2t + 12sin2t) + (1/13)(3sin3t - 2cos3t).  ラプラス変換による解法では   x = (e^(-t)/13)(15cos2t + 3sin2t) + (1/13)(3sin3t - 2cos3t). となりました。どちらかが計算ミスしていると思いますのでご指摘下さい。  ラプラス変換での解法は分数表示が煩雑なので以下のテキストをご参照下さい。 http://ichigo-up.com/Sn2/download/1416878709.txt pass は 1234

  • 偏微分方程式の問題

    次の問題に解答が付いていませんでした。どうか答えが正しいか教えてください。よろしくお願いします。 問題 以下の条件を満たす偏微分方程式を解け     ∂^2f/∂t^2=∂^2f/∂x^2 f(0,t)=f(1,t)=0 (t>=0) f(x,0)=sin(πx) ∂f/∂t=1 (0<=x<=1)          答えはf(x、t)=sin(πx)cos(πx)+∑{2(1-(-1)^k)/(k^2π^2) }sin(kπx)sin(kπt)          となりました。Σはk=1から∞です