跳到主要内容

東京工業大学 工学院 システム制御系 2021年度 2020年8月実施 問題4

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

2つの独立変数 x,tx,t(ともに実数)の関数 u(x,t)u(x,t) は、xx に関するフーリエ変換が可能であるとする。ここで、u(x,t)u(x,t) の xx に関するフーリエ変換とその逆変換は次式で表される。ただし、i=−1i=\sqrt{-1} とする。

F[u(x,t)]=U(ω,t)=12π∫−∞∞u(x,t)e−iωx dx.(1)\mathcal F[u(x,t)]=U(\omega,t)=\frac1{\sqrt{2\pi}} \int_{-\infty}^{\infty}u(x,t)e^{-i\omega x}\,dx.\tag{1}
F−1[U(ω,t)]=u(x,t)=12π∫−∞∞U(ω,t)eiωx dω.(2)\mathcal F^{-1}[U(\omega,t)]=u(x,t)=\frac1{\sqrt{2\pi}} \int_{-\infty}^{\infty}U(\omega,t)e^{i\omega x}\,d\omega.\tag{2}

u(x,t)u(x,t) が次の偏微分方程式と初期条件を満たすとき、以下の問いに答えよ。

∂u∂t=∂2u∂x2(−∞<x<∞,0≤t<∞),(3)\frac{\partial u}{\partial t}=\frac{\partial^2u}{\partial x^2} \qquad(-\infty<x<\infty,\quad0\le t<\infty),\tag{3}
u(x,0)=e−x2(−∞<x<∞).(4)u(x,0)=e^{-x^2}\qquad(-\infty<x<\infty).\tag{4}
  1. 式(3)の両辺を xx に関してフーリエ変換し、U(ω,t)U(\omega,t) の tt に関する微分方程式を導け。

  2. U(ω,0)U(\omega,0) を求めよ。ただし、以下の式を用いてもよい。ここで、α\alpha は正の実数、β\beta は複素数である。

    ∫−∞∞e−α(x+β)2 dx=πα.(5)\int_{-\infty}^\infty e^{-\alpha(x+\beta)^2}\,dx =\sqrt{\frac\pi\alpha}.\tag{5}
  3. 問1と問2で得られた結果を用いて、U(ω,t)U(\omega,t) を求めよ。

  4. 問3で求めた U(ω,t)U(\omega,t) を逆フーリエ変換し、u(x,t)u(x,t) を求めよ。

题目描述​

采用上述对称归一化的傅里叶变换与逆变换。函数 u(x,t)u(x,t) 满足全实轴上的热方程 ut=uxxu_t=u_{xx},t≥0t\ge0,初始条件为 u(x,0)=e−x2u(x,0)=e^{-x^2}。

  1. 对空间变量 xx 作傅里叶变换,导出 U(ω,t)U(\omega,t) 关于 tt 的微分方程。
  2. 求 U(ω,0)U(\omega,0),可使用式(5)的高斯积分,其中 α>0\alpha>0、β∈C\beta\in\mathbb C。
  3. 利用前两问求 U(ω,t)U(\omega,t)。
  4. 作逆傅里叶变换,求 u(x,t)u(x,t)。

Kai​

問1​

部分積分を2回行い、無限遠の境界項が消えることを用いると

F[uxx]=(iω)2F[u]=−ω2U.\mathcal F[u_{xx}]=(i\omega)^2\mathcal F[u]=-\omega^2U.

従って、求める微分方程式は

∂U∂t=−ω2U.\boxed{\frac{\partial U}{\partial t}=-\omega^2U}.

問2​

指数を平方完成し、式(5)を用いると

U(ω,0)=12π∫−∞∞e−x2−iωx dx=e−ω2/42π∫−∞∞e−(x+iω/2)2 dx=12e−ω2/4.\begin{aligned} U(\omega,0) &=\frac1{\sqrt{2\pi}}\int_{-\infty}^{\infty}e^{-x^2-i\omega x}\,dx\\ &=\frac{e^{-\omega^2/4}}{\sqrt{2\pi}} \int_{-\infty}^{\infty}e^{-(x+i\omega/2)^2}\,dx =\boxed{\frac1{\sqrt2}e^{-\omega^2/4}}. \end{aligned}

問3​

問1の解は U(ω,t)=C(ω)e−ω2tU(\omega,t)=C(\omega)e^{-\omega^2t} であり、初期値を代入して

U(ω,t)=12exp⁡[−(t+14)ω2].\boxed{U(\omega,t)=\frac1{\sqrt2}\exp\left[-\left(t+\frac14\right)\omega^2\right]}.

問4​

a=t+1/4>0a=t+1/4>0 とおく。逆変換の指数を平方完成すると

−aω2+iωx=−a(ω−ix2a)2−x24a.-a\omega^2+i\omega x =-a\left(\omega-\frac{ix}{2a}\right)^2-\frac{x^2}{4a}.

従って

u(x,t)=12π∫−∞∞e−aω2+iωx dω=e−x2/(4a)2ππa=11+4texp⁡(−x21+4t).\begin{aligned} u(x,t)&=\frac1{2\sqrt\pi}\int_{-\infty}^{\infty} e^{-a\omega^2+i\omega x}\,d\omega\\ &=\frac{e^{-x^2/(4a)}}{2\sqrt\pi}\sqrt{\frac\pi a} =\boxed{\frac1{\sqrt{1+4t}}\exp\left(-\frac{x^2}{1+4t}\right)}. \end{aligned}