Fourier变换解热方程

考虑

(∗){∂tu−Δu=f(x,t),x∈Rnu(x,0)=φ(x)

由于算子 −ΔRn 有绝对连续谱:

−ΔRn(e−ix⋅ξ)=|ξ|2e−ixξ,∀x⩾0⟹绝对连续谱为[0,+∞)

可以用 Fourier 变换处理。

注:分离变量法与 Fourier 变换
分离变量法和 Fourier 变换是一回事
离散谱:Fourier 级数
绝对连续谱:Fourier 变换

用 Fourier 变换研究初值问题

在 (∗) 两边关于 x 作 Fourier 变换,可得

{∂tu^(ξ,t)+4π2|ξ|2u^(ξ,t)=f^(ξ,t)u^(ξ,0)=φ^(ξ)

这是一阶线性常微分方程,解得

u^(ξ,t)=e−4π2|ξ|2tφ^(ξ)+∫0te−4π2|ξ|2(t−s)f^(ξ,s)ds

若 f≡0,则 u(x,t)=(e−4π2|ξ|2t)∨∗φ.,∨ 表示逆变换

对 Schwarz 函数的对偶空间 f∈(S(Rn))′=S′(Rn) )(缓增分布)
定义 Fourier 变换为 ⟨f^,φ⟩=⟨f,φ^⟩ ,∀φ∈S(Rn)

例

解方程

{ut−Δu=0,x∈Rnu(x,0)=φ(x)

解

两边作 Fourier 变换得

{∂tu^+4π2|ξ|2u^(ξ,t)=0u^(ξ,0)=φ^(ξ)

得 u^(ξ,t)=e−4π2|ξ|2tφ^(ξ)。

两边同时取 Fourier 逆变换,得

u(x,t)=(e−4π2|ξ|2tφ^)∨=(e−4π2|ξ|2t)∨∗φ=例41(2πt)nπn2e−π2|x2πt|2=e−|x|24t(4πt)n2

这就解出了方程:

u(x,t)−1(4πt)n2∫Rne−|x−y|24tφ(y)dy

性质探究:是不是解?

令 K(x)=1(4π)n2e−|x|24,Kt(x)=1(t)nK(xt),则

u(x,t)=(Kt∗φ)(x)

对于 {Kt(x)}t>0,有

∫RnKt(x)dx=∫RnK(x)dx=1
∫Rn|Kt(x)|dx=1(bounded)
  1. ∀η>0,
∫|x|>η|Kt(x)|dx=y=x/t∫|y|>η/tK(y)dy⟶t→0+0.

Claim .

若 φ 是连续、有界函数,则 limt→0+u(x,t)=φ(x).

说明解是解这个初值问题的。

证明

u(x,t)−φ(x)=∫RnKt(y)φ(x−y)dx−∫RnKt(y)φ(x)dy=z=y/t∫RnK(z)φ(x−tz)−∫K(z)φ(x)dz=∫RnK(z)(φ(x−t)−φ(x))dz.

∀ε>0,∃δ s.t. |y|<δ,则

|φ(x−y)−φ(x)|<ε

且由于 K(z) 在 Rn 上可积, ∀ε>0,∃R s.t.

∫|z|>RK(z)dz<ε

于是

|u(x,t)−φ(x)|=|∫RnK(z)(φ(x−tz)−φ(x))dz|≤∫|z|>R+∫|z|<R

后一项:由于 t→0+,|tz|<δ,因此

|φ(x−tz)−φ(x)|<ε

故 |后一项|≤ε∫|K(z)|dz⩽ε。

前一项:设 |φ|≤M,∀x∈Rn,则 第一项≤2M∫|z|>R|K(z)|dz≤2Mε.
故 |u−φ|<ε,Claim 得证。

解的性质

原方程的解是 φ 与一个光滑函数的卷积,因此

  1. ∀t>0 ,解 u(x,t)∈C∞ 是光滑的。(卷积磨光)
  2. supx|u(x,t)|⩽supx|φ(x)|;若 φ(x)>0 则 u>0
  3. 无限传播速度。
  4. 沿时间不能反向演化:u(x,−t) 不满足方程。(反观波方程:可以)
    (热方程的熵:∫Rnulog⁡udx)

非齐次

{∂tu−Δu=f(x,t),x∈Rnu(x,0)=φ(x)

作 Fourier 变换

{∂tu^+4π2|ξ|2u^(ξ,t)=f^(ξ,t)u^(ξ,0)=φ^(ξ)u^(ξ,t)=e−4π2|ξ|2tφ^(ξ)+∫0te−4π2|ξ|2(t−s)f^(ξ,s)ds

再作 Fourier 逆变换

u(x,t)=1(4πt)n2∫Rne−|x−y|24tφ(y)dy+∫0t1(4π(t−s))n2∫Rne−|x−y|24(t−s)f(y,s)dyds