跳到主要内容

東京工業大学 環境・社会理工学院 融合理工学系 2022年8月実施 微分方程

Author

思齐塾, 祭音Myyura

Description

  1. 0xπ,t00 \le x \le \pi, t \ge 0 で定義された次の偏微分方程式 (1) を,与えられた境界条件・初期条件 (2) のもと,変数分離法により解け.
u(x,t)t=42u(x,t)x2(1)\frac{\partial u(x, t)}{\partial t} = 4 \frac{\partial^2 u(x, t)}{\partial x^2} \tag{1}
u(0,t)=0,u(π,t)=0,u(x,0)=2sin3x4sin5x(2)u(0, t) = 0, \quad u(\pi, t) = 0, \quad u(x, 0) = 2 \sin 3x - 4 \sin 5x \tag{2}

题目描述

在区域

0xπ,t00\leq x\leq\pi,\qquad t\geq0

上,使用分离变量法求解热方程

u(x,t)t=42u(x,t)x2,(1)\frac{\partial u(x,t)}{\partial t} =4\frac{\partial^2u(x,t)}{\partial x^2}, \tag{1}

并满足齐次固定端边界条件与初始条件

u(0,t)=0,u(π,t)=0,u(x,0)=2sin3x4sin5x.(2)u(0,t)=0,\qquad u(\pi,t)=0,\qquad u(x,0)=2\sin3x-4\sin5x. \tag{2}

Kai

2.

境界条件は u(0,t)=u(π,t)=0u(0,t)=u(\pi,t)=0 と読む。 u(x,t)=X(x)T(t)u(x,t)=X(x)T(t) とおくと、

XT=4XTXT'=4X''T

より

T4T=XX=λ\frac{T'}{4T}=\frac{X''}{X}=-\lambda

である。したがって

X+λX=0,T+4λT=0.X''+\lambda X=0,\qquad T'+4\lambda T=0.

境界条件 X(0)=X(π)=0X(0)=X(\pi)=0 を満たす非自明解は

λ=n2,Xn(x)=sinnx(n=1,2,)\lambda=n^2,\qquad X_n(x)=\sin nx\qquad(n=1,2,\ldots)

であり、対応する時間因子は Tn(t)=e4n2tT_n(t)=e^{-4n^2t} である。よって一般に

u(x,t)=n=1cne4n2tsinnx.u(x,t)=\sum_{n=1}^{\infty}c_ne^{-4n^2t}\sin nx.

初期条件

u(x,0)=2sin3x4sin5xu(x,0)=2\sin3x-4\sin5x

から、直交性により c3=2, c5=4c_3=2,\ c_5=-4 、その他の cnc_n00 である。したがって

u(x,t)=2e36tsin3x4e100tsin5x\boxed{u(x,t)=2e^{-36t}\sin3x-4e^{-100t}\sin5x}

を得る。各項は ut=4uxxu_t=4u_{xx} を満たし、 x=0,πx=0,\pi で消え、 t=0t=0 で所定の初期関数になるので、すべての条件を満たしている。