数学 偏微分方程

押题

偏微分方程复习题与证明提纲,涵盖变分、波动方程、热方程和位势问题。

本页目录19 节

复习汇总 #

题目 #

叙述热方程基本解的定义.

题目 #

叙述位势方程基本解的定义.

题目 #

叙述广义函数及其弱导数的定义.

题目 #

调和函数的定义.

题目 #

设

J(v)=12∫Ω(∣∇v∣2+v2)dx+12∫∂Ωα(x)v2ds−∫Ωfvdx−∫∂ΩgvdsJ(v)=\frac 1 2\int_\Omega (|\nabla v|^2+v^2)\text{d} x+\frac 12 \int_{\partial \Omega} \alpha(x)v^2\text{d} s-\int_\Omega fv\text{d} x-\int_{\partial\Omega}gv\text{d} s

其中 α(x)⩾0\alpha(x)\geqslant 0 考虑以下三个问题:

  • I 变分问题: 求 u∈M=C1(Ω‾)u\in M=C^1(\overline{\Omega}), 使得
J(u)=min⁡v∈MJ(v).J(u)=\min\limits_{v\in M}J(v).
  • II 求 u∈M=C1(Ω‾)u\in M=C^1(\overline{\Omega}), 使得它对于任意 v∈Mv\in M, 都满足
∫Ω(∇u⋅∇v+u⋅v−fv)dx+∫∂Ω(α(x)uv−gv)ds=0.\int_\Omega (\nabla u\cdot\nabla v+u\cdot v-fv)\text{d} x+\int_{\partial\Omega}(\alpha(x)uv-gv)\text{d} s=0.
  • III 第三边值问题: 求 u∈C2(Ω)∩C1(Ω‾)u\in C^2(\Omega)\cap C^1(\overline{\Omega}), 满足以下边值问题
{−Δu+u=f,x∈Ω,∂u∂n+α(x)u=g,x∈∂Ω.\begin{cases} -\Delta u+u=f, & x\in\Omega,\\ \dfrac{\partial u}{\partial n}+\alpha(x)u=g, & x\in\partial \Omega. \end{cases}
  • (1) 证明问题 I 与问题 II 等价.
  • (2) 当 u∈C2(Ω)∩C1(Ω‾)u\in C^2(\Omega)\cap C^1(\overline{\Omega}) 时, 证明问题 I,II,III 等价.
证明
  • (1) ⇒\Rightarrow: 令 j(ε)=J(u+εv)j(\varepsilon)=J(u+\varepsilon v), u∈M,v∈Mu\in M, v\in M. 那么有
j′(0)=∫Ω∇v⋅∇u+vudx⃗+∫∂Ωαuvdx⃗−∫Ωfvdx⃗−∫∂Ωgvdx⃗=0j'(0)=\int_\Omega \nabla v\cdot \nabla u+vu\text{d}\vec x+\int_{\partial \Omega} \alpha u v\text{d} \vec x-\int_\Omega fv\text{d} \vec x-\int_{\partial\Omega} gv\text{d} \vec x=0

即

∫Ω(∇u⋅∇v+u⋅v−fv)dx+∫∂Ω(α(x)uv−gv)ds=0.\int_\Omega (\nabla u\cdot\nabla v+u\cdot v-fv)\text{d} x+\int_{\partial\Omega}(\alpha(x)uv-gv)\text{d} s=0.

再验证 j′′(ε)=∫Ω(∣∇v∣2+v2)dx⃗+∫∂Ωαv2dx⃗>0j''(\varepsilon)=\int_\Omega(|\nabla v|^2+v^2)\text{d}\vec x+\int_{\partial\Omega}\alpha v^2\text{d} \vec x> 0.

⇐\Leftarrow: 令 w=v−u∈Mw=v-u\in M, 考虑

J(u+w)−J(u)=∫Ω∇u⋅∇w+uw+12(∣∇w∣2+w2)dx⃗+∫∂Ωα(uw+12w2)dx⃗−∫Ωfwdx⃗−∫∂Ωgwdx⃗.J(u+w)-J(u)=\int_\Omega \nabla u\cdot\nabla w+uw+\frac{1}{2}(|\nabla w|^2+w^2)\text{d}\vec x+\int_{\partial\Omega}\alpha (uw+\frac12 w^2)\text{d}\vec x-\int_{\Omega} fw\text{d}\vec x-\int_{\partial\Omega} gw\text{d}\vec x.

由条件带入 v=wv=w.

得到

J(u+w)−J(u)=12∫Ω(∣∇w∣2+w2)dx⃗+12∫∂Ωαw2dx⃗⩾0J(u+w)-J(u)=\frac12\int_\Omega(|\nabla w|^2+w^2)\text{d}\vec x+\frac12\int_{\partial \Omega}\alpha w^2\text{d}\vec x\geqslant 0

从而问题 II 的解 uu 是 J(u)J(u) 的最小值点.

综上问题 I, II 等价.

  • (2) II ⇒\Rightarrow III: 取 v∈C02v\in C_0^2, 由 Green 公式
∫Ω(∇u⋅∇v+uv−fv)dx⃗=∫Ω(−Δu+u−f)vdx=0\int_\Omega (\nabla u\cdot\nabla v+uv-fv)\text{d}\vec x=\int_\Omega (-\Delta u+u-f)v\text{d} x=0

那么由 vv 任意可得 −Δu+u=f-\Delta u+u=f. 再带回 II 得到

∫∂Ω(∂u∂n+αu−g)vds=0\int_{\partial \Omega}(\frac{\partial u}{\partial n}+\alpha u-g)v\text{d} s=0

从而 ∂u∂n+αu=g\dfrac{\partial u}{\partial n}+\alpha u=g.

III ⇐\Leftarrow II: 由 Green 公式, 可得 II 条件成立.

题目 #

用特征线法求解下述问题

  • (1) {ut+2ux=0,t>0,−∞<x<∞u∣t=0=x2,−∞<x∞\begin{cases} u_t+2u_x=0, & t>0,-\infty<x<\infty\\ u\big|_{t=0}=x^2, & -\infty<x\infty \end{cases}
  • (2) {ut+(1+x2)ux−u=0,t>0,−∞<x<∞u∣t=0=arctan⁡x,−∞<x<∞.\begin{cases} u_t+(1+x^2)u_x-u=0, & t>0,-\infty<x<\infty\\ u\big|_{t=0}=\arctan x, & -\infty<x<\infty. \end{cases}

题目 #

(!!!) 设 u∈C2(Rn×[0,∞))u\in C^2(\mathbb{R}^n\times [0,\infty)) 是方程 utt−a2Δu=0u_{tt}-a^2\Delta u=0 的解. 令

Mu(x,r,t)=1ωn∫∣ξ∣=1u(x+rξ,t)dSξM_u(x,r,t)=\frac{1}{\omega_n}\int_{|\xi|=1} u(x+r\xi,t)\text{d} S_\xi

其中 ωn\omega_n 是单位球面的面积. 证明

∂2∂t2Mu=a2(∂2Mu∂r2+n−1r∂Mu∂r)\frac{\partial^2}{\partial t^2}M_u=a^2\left(\frac{\partial^2 M_u}{\partial r^2}+\frac{n-1}{r}\frac{\partial M_u}{\partial r}\right)
证明

一方面,

∂2∂t2Mu=1ωn∫∣ξ∣=1utt(x+rξ,t)dSξ=1ωn∫∣ξ∣=1a2Δx+rξu(x+rξ,t)dSξ=a2ωn∫∂B(x,r)Δu(y,t)1rn−1dS\begin{aligned} \frac{\partial^2}{\partial t^2}M_u=\frac{1}{\omega_n}\int_{|\xi|=1}u_{tt}(x+r\xi,t)\text{d} S_\xi\\ =\frac{1}{\omega_n}\int_{|\xi|=1}a^2\Delta_{x+r\xi} u(x+r\xi,t)\text{d} S_\xi\\ =\frac{a^2}{\omega_n}\int_{\partial B(x,r)}\Delta u(y,t) \frac{1}{r^{n-1}} \text{d} S \end{aligned}

另一方面,

∂Mu∂r=1ω∫∣ξ∣=1∇u(y,t)⋅ξ⃗dSξ=\begin{aligned} \frac{\partial M_u}{\partial r}=\frac{1}{\omega}\int_{|\xi|=1}\nabla u(y,t)\cdot\vec\xi\text{d} S_\xi\\ = \end{aligned}

题目 #

计算下列函数的 Fourier 变换.

  • (1) f(x)=e−Ax2,(A>0,x∈R)f(x)=e^{-Ax^2},(A>0, x\in\mathbb{R}).
  • (2) f(x⃗)=e−A∣x⃗∣2,(A>0,x∈Rn)f(\vec x)=e^{-A|\vec x|^2}, (A>0,x\in\mathbb{R}^n).

题目 #

计算/证明广义函数.

  • (1) (∣x∣)(m),(m⩾1)(|x|)^{(m)}, (m\geqslant 1).
  • (2) (H(x)eax)′′(H(x)e^{ax})''.
  • (3) xmδ(m)(x)=(−1)mm!δ(x)x^m\delta^{(m)}(x)=(-1)^m m! \delta(x).

题目 #

用分离变量法求解

  • (1)
{ut=a2uxx,0<x<π, t>0,u∣t=0=sin⁡x,0⩽x⩽π,ux∣x=0=ux∣x=π=0,t>0.\begin{cases} u_t = a^2 u_{xx}, & 0 < x < \pi,\ t > 0, \\ u\big|_{t=0} = \sin x, & 0 \leqslant x \leqslant \pi, \\ u_x\big|_{x=0} = u_x\big|_{x=\pi} = 0, & t > 0. \end{cases}
  • (2) {ut=a2uxx,0<x<l, t>0,u∣t=0=x2(l−x),0⩽x⩽l,ux∣x=0=u∣x=l=0,t>0.\begin{cases} u_t = a^2 u_{xx}, & 0 < x < l,\ t > 0, \\ u\big|_{t=0} = x^2 (l - x), & 0 \leqslant x \leqslant l, \\ u_x\big|_{x=0} = u\big|_{x=l} = 0, & t > 0. \end{cases}
答案
u(x,t)=2π+4π∑m=1∞11−4m2e−a2(2m)2tcos⁡(2mx)u(x,t)=\frac{2}{\pi}+\frac{4}{\pi}\sum\limits_{m=1}^\infty \frac{1}{1-4m^2}e^{-a^2(2m)^2t}\cos(2mx)

注意此时是以 {cos⁡nx}\{\cos nx\} 作为正交基, 所以会有 n=0n=0 的项.

答案

该边值条件算出以 {cos⁡μnx}\{\cos\mu_n x\} 为基, 其中 μn=(2n+1)π2\mu_n=\frac{(2n+1)\pi}{2}, 也有 (cos⁡μnx,cos⁡μnx)=l2(\cos\mu_n x,\cos\mu_n x)=\frac{l}{2}, 故展开的系数还是 2l\frac 2l. 且没有 n=0n=0 的项.

un(x,t)=∑n=1∞φne−(a(2n+1)π/(2l))2tcos⁡(2n+1)π2lx,φn=(−1)n64(2n+1)πl3−192l3(2n+1)4π4u_n(x,t)=\sum\limits_{n=1}^\infty\varphi_n e^{-\left(a(2n+1)\pi/(2l)\right)^2t}\cos\frac{(2n+1)\pi}{2l}x,\quad \varphi_n=\frac{(-1)^n 64(2n+1)\pi l^3-192l^3}{(2n+1)^4\pi^4}

题目 #

(!!!) 设 φ(x)\varphi(x) 是连续函数且有界, 则

u(x,t)=∫−∞∞K(x−ξ,t)φ(ξ)dξu(x,t)=\int_{-\infty}^\infty K(x-\xi,t)\varphi(\xi)\text{d}\xi

给出初值问题

{ut−a2uxx=0,−∞<x<∞,t>0u(x,0)=φ(x),−∞<x<∞\begin{cases} u_t-a^2u_{xx}=0, & -\infty<x<\infty, t>0\\ u(x,0)=\varphi(x), & -\infty<x<\infty \end{cases}

的解.

证明

注意积分和求导可交换的条件, 即对 KK 求导后, 固定 tt 均对 xx 可积.

K(x,t)={12aπte−x2/(4a2t),t>00,t⩽0K(x,t)=\begin{cases} \frac{1}{2a\sqrt{\pi t}}e^{-x^2/(4a^2t)}, & t>0\\ 0, & t\leqslant 0 \end{cases}

题目 #

设 u,ux∈C2,1(Q)∩C1,0(Q‾)u,u_x\in C^{2,1}(Q)\cap C^{1,0}(\overline{Q}), uu 满足第三边值问题

{Lu=ut−uxx=f(x,t)(x,t)∈Qu∣t=0=φ(x),0⩽x⩽l[−∂u∂x+αu]x=0=g1(t),0⩽t⩽T[∂u∂x+βu]x=l=g2(t),0⩽t⩽T\begin{cases} Lu=u_t-u_{xx}=f(x,t) & (x,t)\in Q\\ u|_{t=0}=\varphi(x), & 0\leqslant x\leqslant l\\ [-\frac{\partial u}{\partial x}+\alpha u]_{x=0}=g_1(t), & 0\leqslant t\leqslant T\\ [\frac{\partial u}{\partial x}+\beta u]_{x=l}=g_2(t), & 0\leqslant t\leqslant T \end{cases}

其中 α,β⩾0\alpha,\beta\geqslant 0, 给出 max⁡∣∂u∂x∣\max|\frac{\partial u}{\partial x}| 的估计.

证明

设 v=uxv=u_x, 首先由 uu 的第三边值问题可以给出估计

∣u∣⩽C(F+B)|u|\leqslant C(F+B)

, F=sup⁡∣f∣,B=max⁡{sup⁡∣φ∣,sup⁡∣g1∣,sup⁡∣g2∣}F=\sup|f|, B=\max\{\sup|\varphi|,\sup|g_1|,\sup|g_2|\}, C=max⁡{T,1+2a2Tl+l4}C=\max\{T,1+\frac{2a^2T}{l}+\frac l4\}.

则考虑 vv 满足第一边值问题

{Lv=vt−vxx=fx,(x,t)∈Qv(x,0)=φ′(x),0⩽x⩽lv(0,t)=αu(0,t)−g1(t),0⩽t⩽Tv(l,t)=−βu(l,t)−g2(t),0⩽t⩽T\begin{cases} Lv=v_t-v_{xx}=f_x, & (x,t)\in Q\\ v(x,0)=\varphi'(x), & 0\leqslant x\leqslant l\\ v(0,t)=\alpha u(0,t)-g_1(t), & 0\leqslant t\leqslant T\\ v(l,t)=-\beta u(l,t)-g_2(t), & 0\leqslant t\leqslant T \end{cases}

从而有估计

∣v∣⩽FT+B|v|\leqslant FT+B

其中 F=sup⁡∣fx∣,B=max⁡{sup⁡∣φ′∣,sup⁡∣u(0,t)−g1(t)∣,sup⁡∣u(l,t)−g2(t)∣}F=\sup|f_x|, B=\max\{\sup|\varphi'|,\sup|u(0,t)-g_1(t)|,\sup|u(l,t)-g_2(t)|\}, 由 uu 的估计可得 BB 存在.

题目 #

记 Ql={0<x<l,0<t⩽T}Q^l=\{0<x<l,0<t\leqslant T\}, 设 ul∈C2,1(Ql)∩C(Ql‾)u_l\in C^{2,1}(Q^l)\cap C(\overline{Q^l}) 是定解问题

{Lu=ut−uxx=0,(x,t)∈Qlul∣t=0=0,0⩽x⩽lul(0,t)=g1(t),ul(l,t)=0,0⩽t⩽T\begin{cases} Lu=u_t-u_{xx}=0, & (x,t)\in Q^l\\ u_l\big|_{t=0}=0, & 0\leqslant x\leqslant l\\ u_l(0,t)=g_1(t), u_l(l,t)=0, & 0\leqslant t\leqslant T \end{cases}

的解, 其中 g1(t)⩾0g_1(t)\geqslant0 证明: l1<l2l_1<l_2 时 ul1(x,t)⩽ul2(x,t),(x,t)∈Ql1u_{l_1}(x,t)\leqslant u_{l_2}(x,t),\quad (x,t)\in Q^{l_1}.

证明

先证明 ulu_l 非负, 直接使用弱极值原理即可得到.

于是定义 w(x,t)=ul2(x,t)−ul1(x,t)w(x,t)=u_{l_2}(x,t)-u_{l_1}(x,t), 在 Ql1Q^{l_1} 上.

那么再用一次弱极值原理, 对于 w(x,t)w(x,t) 的边界 w(0,t)=0,w(l1,t)=ul2(l1,t)⩾0w(0,t)=0, w(l_1,t)=u_{l_2}(l_1,t)\geqslant 0, 从而 w(x,t)w(x,t) 非负. 证毕.

题目 #

设 uu 满足

Lu=ut−a2uxx+c(x,t)u⩽0Lu=u_t-a^2u_{xx}+c(x,t)u\leqslant 0

其中 c(x,t)c(x,t) 有界, 且 c(x,t)⩾0c(x,t)\geqslant 0. 证明如果 uu 在 Q‾\overline{Q} 上存在非负最大值, 则必在边界上达到非负最大值.

证明

按弱极值原理证明过程走一遍即可, 先考虑严格小于 0.

再考虑辅助函数 v(x,t)=u(x,t)+εe−Ntv(x,t)=u(x,t)+\varepsilon e^{-Nt}, 变为严格小于 0. 其中 N>∣c(x,t)∣N>|c(x,t)|, 由 c(x,t)c(x,t) 有界, 保证存在这样的 N.

题目 #

设 u(x,t)∈C2,1(Q‾)u(x,t) \in C^{2,1}(\overline{Q}) 是问题

{∂u∂t−∂2u∂x2=f,(x,t)∈Q,u(x,0)=φ(x),0⩽x⩽l,u(0,t)=u(l,t)=0,0⩽t⩽T \begin{cases} \frac{\partial u}{\partial t} - \frac{\partial^2 u}{\partial x^2} = f, & (x,t) \in Q, \\ u(x,0) = \varphi(x), & 0 \leqslant x \leqslant l, \\ u(0,t) = u(l,t) = 0, & 0 \leqslant t \leqslant T \end{cases}

的解, 证明 uu 满足以下估计

sup⁡0⩽t⩽T∫0l(∂u∂x)2dx+∫0T∫0l(∂u∂t)2dxdt⩽M[∫0l(φ′(x))2dx+∫0T∫0lf2(x,t)dxdt], \sup_{0 \leqslant t \leqslant T} \int_0^l \left(\frac{\partial u}{\partial x}\right)^2 \mathrm{d}x + \int_0^T \int_0^l \left(\frac{\partial u}{\partial t}\right)^2 \mathrm{d}x \mathrm{d}t \leqslant M \left[ \int_0^l (\varphi'(x))^2 \mathrm{d}x + \int_0^T \int_0^l f^2(x,t) \mathrm{d}x \mathrm{d}t \right],

其中 MM 只依赖于 T,lT,l.

题目 #

设 u(x,t)∈C1,0(Q‾)∩C2,1(Q)u(x,t) \in C^{1,0}(\overline{Q}) \cap C^{2,1}(Q) 且满足以下定解问题

{ut−a2uxx=f(x,t),(x,t)∈Q,u(x,0)=φ(x),0⩽x⩽l,[−∂u∂x+αu]x=0=[∂u∂x+βu]x=l=0,0⩽t⩽T, \begin{cases} u_t - a^2 u_{xx} = f(x,t), & (x,t) \in Q, \\ u(x,0) = \varphi(x), & 0 \leqslant x \leqslant l, \\ \left[-\frac{\partial u}{\partial x} + \alpha u\right]_{x=0} = \left[\frac{\partial u}{\partial x} + \beta u\right]_{x=l} = 0, & 0 \leqslant t \leqslant T, \end{cases}

其中 α⩾0,β⩾0\alpha \geqslant 0, \beta \geqslant 0, 证明

sup⁡0⩽t⩽T∫0lu2(x,t) dx+∫0T∫0l(∂u∂x)2dxdt⩽M(∫0lφ2(x)dx+∫0T∫0lf2(x,t)dxdt), \sup_{0 \leqslant t \leqslant T} \int_0^l u^2(x,t) \, \mathrm{d}x + \int_0^T \int_0^l \left(\frac{\partial u}{\partial x}\right)^2 \mathrm{d}x \mathrm{d}t \leqslant M \left( \int_0^l \varphi^2(x) \mathrm{d}x + \int_0^T \int_0^l f^2(x,t) \mathrm{d}x \mathrm{d}t \right),

其中 MM 只依赖于 T,aT, a.

题目 #

(!!!) 二维/三维 Poisson 公式推导 (讲义21 P3-4).

题目 #

(!!!) 设 u∈C2(Ω‾)u\in C^2(\overline{\Omega}) 是方程

{Δu=0,{∣x∣>1}≜Ωlim⁡∣x∣→∞u(x)=0\begin{cases} \Delta u=0, & \{|x|>1\}\triangleq \Omega\\ \lim\limits_{|x|\to \infty} u(x)=0 \end{cases}

证明: max⁡Ω‾∣u∣=max⁡∂Ω∣u∣\max\limits_{\overline{\Omega}}|u|=\max\limits_{\partial\Omega}|u|.

位势的极值和调和函数未整理!

讨论

评论

正在加载评论…

输入关键词开始搜索。