数学 旧 .com 迁移

偏微分方程 / 习题/考试 / 作业 / 5

从旧 .com 全量搬运的历史内容,来源路径:/math/课程/偏微分方程/exams/作业/5/

迁移来源

20251126 #

姓名: 刘欣楠

序号: 12

班级: 数学强基 2301

学号: 2233310237

\begin{tikzpicture}[remember picture,overlay] \node[anchor=north east] at ($(current page.north east)+(-25mm,-50mm)$) { ![图](/media/figures/pde/12.png) }; \end{tikzpicture}

3-16 #

Ql={0<x<l,0<tT}Q^l=\{0<x<l,0<t\leqslant T\}, 设 ulC2,1(Ql)C(Ql)u_l\in C^{2,1}(Q^l)\cap C(\overline{Q}^l) 是定解问题

{ult2ulx2=0,(x,t)Ql,ult=0=0,0xl,ul(0,t)=g1(t), ul(l,t)=0,0tT\begin{cases} \dfrac{\partial u_l}{\partial t}-\dfrac{\partial^2 u_l}{\partial x^2}=0, & (x,t)\in Q^l,\\ u_l|_{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(x)0g_1(x)\geqslant 0 . 证明 ul(x,t)u_l(x,t) 关于 ll 是递增的, 即对于 l1<l2l_1<l_2,

ul1(x,t)Ul2(x,t),(x,t)Ql1,u_{l_1}(x,t)\leqslant U_{l_2}(x,t),\quad (x,t)\in Q^{l_1},

并给出物理解释.

note

对每个 l>0l>0, ulu_l 满足

\begin{cases} u_t-u_{xx}=0,& 0<x<l,\ 0<t\leqslant T,\

2mm] u_l(x,0)=0,& 0\leqslant x\leqslant l,\

2mm] u_l(0,t)=g_1(t)\ge0,\quad u_l(l,t)=0,& 0\leqslant t\leqslant T. \end{cases}

v=ulv=-u_l, 则 vtvxx=(utuxx)=0v_t-v_{xx}=-(u_t-u_{xx})=0, 且在抛物边界上

v(x,0)=ul(x,0)=0,v(0,t)=g1(t)0,v(l,t)=ul(l,t)=0. v(x,0)=-u_l(x,0)=0,\quad v(0,t)=-g_1(t)\le0,\quad v(l,t)=-u_l(l,t)=0.

因此 vv 在边界上最大值为 00, 由最大值原理知 v0v\le0, 即 ul0u_l\ge0Ql\overline{Q^l}.

0<l1<l20<l_1<l_2, 设对应解为 ul1,ul2u_{l_1},u_{l_2}. 在区域 Ql1={0<x<l1,0<tT}Q^{l_1}=\{0<x<l_1,0<t\leqslant T\} 上考虑 w(x,t):=ul2(x,t)ul1(x,t).w(x,t):=u_{l_2}(x,t)-u_{l_1}(x,t). 两者都满足 utuxx=0u_t-u_{xx}=0, 故 wtwxx=0w_t-w_{xx}=0 ww 在抛物边界上的值:

w(x,0)=ul2(x,0)ul1(x,0)=0,0xl1, w(x,0)=u_{l_2}(x,0)-u_{l_1}(x,0)=0,\quad 0\leqslant x\leqslant l_1, w(0,t)=ul2(0,t)ul1(0,t)=g1(t)g1(t)=0, w(0,t)=u_{l_2}(0,t)-u_{l_1}(0,t)=g_1(t)-g_1(t)=0, w(l1,t)=ul2(l1,t)ul1(l1,t)=ul2(l1,t)0=ul2(l1,t)0, w(l_1,t)=u_{l_2}(l_1,t)-u_{l_1}(l_1,t)=u_{l_2}(l_1,t)-0=u_{l_2}(l_1,t)\ge0,

于是 ww 在抛物边界上满足 w0w\ge0, 且在内部满足 wtwxx=0w_t-w_{xx}=0. 再由最大值原理得 w(x,t)0,(x,t)Ql1,w(x,t)\ge0,\quad (x,t)\in \overline{Q^{l_1}},ul1(x,t)ul2(x,t),(x,t)Ql1.u_{l_1}(x,t)\le u_{l_2}(x,t),\quad (x,t)\in Q^{l_1}. 这就证明了 ulu_l 关于 ll 的递增性.

\noindent**物理解释: **

可以把问题看作一根长度为 ll 的均匀细杆, 其初温为 00, 左端 x=0x=0 的温度被控制为非负的给定函数 g1(t)g_1(t), 右端 x=lx=l 一直保持在 00^\circ. 当把杆从长度 l1l_1 延长到 l2>l1l_2>l_1 时, 对同一内部点 x(0,l1)x\in(0,l_1) 来说, 右端这个“冷端”被移得更远, 热量更不容易通过右端散失掉, 因此该点的温度在任何时刻都不会比短杆情况更低, 即

ul1(x,t)ul2(x,t). u_{l_1}(x,t)\le u_{l_2}(x,t).

3-18 #

uC(Q)C2,1(Q)u\in C(\overline{Q})\cap C^{2,1}(Q) 且满足

Luuta2uxx+c(x,t)u0,(x,t)Q,Lu\equiv u_t-a^2u_{xx}+c(x,t)u\leqslant0,\quad (x,t)\in Q,

其中 c(x,t)c(x,t) 有界, 且 c(x,t)0c(x,t)\geqslant 0. 试证明: 如果 uuQ\overline{Q} 上存在非负最大值, 则 uu 必在抛物边界 Γ\Gamma 上达到它在 Q\overline{Q} 上的非负最大值.

note

Q={(x,t):0<x<l, 0<tT}Q=\{(x,t):0<x<l,\ 0<t\leqslant T\}, 其抛物边界为 Γ={(x,t):0xl, t=0}{(x,t):x=0 或 x=l, 0tT}.\Gamma=\{(x,t):0\leqslant x\leqslant l,\ t=0\}\cup\{(x,t):x=0\text{ 或 }x=l,\ 0\leqslant t\leqslant T\}. 已知 uC(Q)C2,1(Q)u\in C(\overline Q)\cap C^{2,1}(Q)Lu:=uta2uxx+c(x,t)u0,c(x,t)0, c 有界.Lu:=u_t-a^2u_{xx}+c(x,t)u\le0,\qquad c(x,t)\ge0,\ c\ \text{有界}.

取常数 λ>0\lambda>0, 令 v(x,t)=eλtu(x,t).v(x,t)=e^{\lambda t}u(x,t).vC(Q)C2,1(Q)v\in C(\overline Q)\cap C^{2,1}(Q), 且计算得 vt=eλt(ut+λu),vxx=eλtuxx,v_t=e^{\lambda t}(u_t+\lambda u),\quad v_{xx}=e^{\lambda t}u_{xx}, 于是

vta2vxx+(c+λ)v=eλt(uta2uxx+cu)=eλtLu0. v_t-a^2v_{xx}+(c+\lambda)v =e^{\lambda t}\big(u_t-a^2u_{xx}+c\,u\big) =e^{\lambda t}Lu\le0.

Lλv:=vta2vxx+(c+λ)v0L_\lambda v:=v_t-a^2v_{xx}+(c+\lambda)v\le0, 其中 c+λλ>0c+\lambda\geqslant\lambda>0.

vv 在某个内部点 (x0,t0)Q(x_0,t_0)\in Q 取到正的最大值 M>0M>0, 则因为 (x0,t0)(x_0,t_0)vv 的局部极大点, 有 vt(x0,t0)=0,vx(x0,t0)=0,vxx(x0,t0)0.v_t(x_0,t_0)=0, v_x(x_0,t_0)=0, v_{xx}(x_0,t_0)\le0. 于是

Lλv(x0,t0)=vt(x0,t0)a2vxx(x0,t0)+(c+λ)v(x0,t0)0a2vxx(x0,t0)+(c+λ)M>0, L_\lambda v(x_0,t_0) = v_t(x_0,t_0)-a^2v_{xx}(x_0,t_0)+(c+\lambda)v(x_0,t_0) \ge 0- a^2\cdot v_{xx}(x_0,t_0)+(c+\lambda)M>0,

这与 Lλv0L_\lambda v\le0 矛盾. 故 vvQQ 内部不能有正的局部最大值.

M:=maxQu0M:=\max_{\overline Q}u\ge0uuQ\overline Q 上的非负最大值.

对任意 τ(0,T)\tau\in(0,T), 考虑截断柱体 Qτ={(x,t):0<x<l, 0<tτ},Q_\tau=\{(x,t):0<x<l,\ 0<t\leqslant\tau\}, 其抛物边界为

Γτ={t=0, 0xl}{x=0 或 x=l, 0tτ}Γ. \Gamma_\tau=\{t=0,\ 0\leqslant x\leqslant l\}\cup\{x=0\text{ 或 }x=l,\ 0\leqslant t\leqslant\tau\}\subset\Gamma.

由第一步的结论, uuQτ\overline Q_\tau 上的正最大值不能在内部取得, 只能在 Γτ\Gamma_\tau 上取得, 即 maxQτu=maxΓτu.\max_{\overline Q_\tau}u=\max_{\Gamma_\tau}u.τT\tau\uparrow T, 由 uu 的连续性及单调性可得

maxQu=limτTmaxQτu=limτTmaxΓτumaxΓu. \max_{\overline Q}u=\lim_{\tau\uparrow T}\max_{\overline Q_\tau}u =\lim_{\tau\uparrow T}\max_{\Gamma_\tau}u \le\max_{\Gamma}u.

另一方面显然有 maxQumaxΓu\max_{\overline Q}u\geqslant\max_{\Gamma}u, 故 maxQu=maxΓu.\max_{\overline Q}u=\max_{\Gamma}u.

因此, 只要 uuQ\overline Q 上存在非负最大值, 它必定在抛物边界 Γ\Gamma 上达到该最大值.

3-21 #

证明半无界问题

{uta2uxx=f(x,t),0<x<,t>0,ut=0=φ(x),0x<,ux=0=μ(t),t0.\begin{cases} u_t-a^2u_{xx}=f(x,t),& 0<x<\infty,t>0,\\ u|_{t=0}=\varphi(x), & 0\leqslant x<\infty,\\ u|_{x=0}=\mu(t), & t\geqslant 0. \end{cases}

的有界解是唯一的.

note

u1,u2u_1,u_2 是题中问题的两个有界解, 记 w(x,t)=u1(x,t)u2(x,t),w(x,t)=u_1(x,t)-u_2(x,t),ww 亦有界, 并满足齐次问题

{wta2wxx=0,0<x<, t>0,w(x,0)=0,x0,w(0,t)=0,t0. \begin{cases} w_t-a^2w_{xx}=0,& 0<x<\infty,\ t>0,\\ w(x,0)=0,& x\ge0,\\ w(0,t)=0,& t\ge0. \end{cases}

要证唯一性, 只需证明 w0w\equiv0.

取任意 λ>0\lambda>0, 定义 F(t)=0w(x,t)2eλxdx.F(t)=\int_0^\infty w(x,t)^2 e^{-\lambda x}\,dx. 由于 ww 有界, F(t)F(t) 对每个 tt 都有限;又由 w(x,0)=0w(x,0)=0, 得 F(0)=0F(0)=0.

wt=a2wxxw_t=a^2w_{xx} 计算:

F(t)=20wwteλxdx=2a20wwxxeλxdx. F'(t)=2\int_0^\infty w w_t e^{-\lambda x}\,dx =2a^2\int_0^\infty w w_{xx} e^{-\lambda x}\,dx.

xx 分部积分两次可得

0wwxxeλxdx=0wx2eλxdx+λ220w2eλxdx. \int_0^\infty w w_{xx} e^{-\lambda x}\,dx = -\int_0^\infty w_x^2 e^{-\lambda x}\,dx +\frac{\lambda^2}{2}\int_0^\infty w^2 e^{-\lambda x}\,dx.

于是

F(t)=2a2(0wx2eλxdx+λ22F(t))a2λ2F(t). F'(t) =2a^2\Big(-\int_0^\infty w_x^2 e^{-\lambda x}\,dx +\frac{\lambda^2}{2}F(t)\Big) \le a^2\lambda^2 F(t).

C=a2λ2C=a^2\lambda^2, 考虑 G(t)=eCtF(t),G(t)=e^{-Ct}F(t),G(t)=eCt(F(t)CF(t))0,G'(t)=e^{-Ct}\big(F'(t)-CF(t)\big)\le0, 说明 G(t)G(t) 非增. 由 F(0)=0F(0)=0G(0)=0G(0)=0, 再由 G(t)0G(t)\ge0, 得到 G(t)0F(t)0, t0.G(t)\equiv0\quad\Rightarrow\quad F(t)\equiv0,\ \forall t\ge0. 从而对任意 tt, 0w(x,t)2eλxdx=0,\int_0^\infty w(x,t)^2 e^{-\lambda x}\,dx=0, 由被积函数非负可知 w(x,t)0w(x,t)\equiv0. 因此 u1u2u_1\equiv u_2, 半无界问题的有界解唯一.

3-22 #

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

{utuxx=f,(x,t)Q,u(x,0)=φ(x),0xl,u(0,t)=u(l,t)=0,0tT\begin{cases} u_{t}-u_{xx}=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 满足以下估计

sup0tT0l(ux)2dx+0T0l(ut)2dxdtM[0l(φ(x))2dx+0T0lf2(x,t)dxdt],\begin{aligned} & \sup\limits_{0\leqslant t\leqslant T} \int_0^l (u_x)^2\text{d} x + \int_0^T\int_0^l (u_t)^2\text{d} x\text{d} t\\ \leqslant & M\left[\int_0^l (\varphi'(x))^2\text{d} x + \int_0^T\int_0^l f^2(x,t)\text{d} x\text{d} t\right], \end{aligned}

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

note

对方程两边乘以 utu_t, 在 (0,l)(0,l) 上积分, 得 0lut2dx0luxxutdx=0lfutdx.\int_0^l u_t^2\,dx-\int_0^l u_{xx}u_t\,dx=\int_0^l f\,u_t\,dx. 对第二项分部积分:

0luxxutdx=[uxut]0l+0luxutxdx. -\int_0^l u_{xx}u_t\,dx = -\big[u_x u_t\big]_{0}^{l}+\int_0^l u_x u_{tx}\,dx.

由于边界条件 u(0,t)=u(l,t)=0u(0,t)=u(l,t)=0 与时间无关, 对 tt 求导得 ut(0,t)=ut(l,t)=0,u_t(0,t)=u_t(l,t)=0, 故边界项 [uxut]0l=0[u_x u_t]_{0}^{l}=0, 于是

0luxxutdx=0luxuxtdx=12ddt0lux2dx. -\int_0^l u_{xx}u_t\,dx=\int_0^l u_x u_{xt}\,dx =\frac12\frac{d}{dt}\int_0^l u_x^2\,dx.

代回原式:

0lut2dx+12ddt0lux2dx=0lfutdx. \int_0^l u_t^2\,dx + \frac12\frac{d}{dt}\int_0^l u_x^2\,dx = \int_0^l f\,u_t\,dx.

应用 Cauchy—Schwarz 与 Young 不等式:

0lfutdx120lut2dx+120lf2dx, \left|\int_0^l f\,u_t\,dx\right| \le \frac12\int_0^l u_t^2\,dx + \frac12\int_0^l f^2\,dx,

从而

12ddt0lux2dx+0lut2dx120lut2dx+120lf2dx, \frac12\frac{d}{dt}\int_0^l u_x^2\,dx + \int_0^l u_t^2\,dx \le \frac12\int_0^l u_t^2\,dx + \frac12\int_0^l f^2\,dx,

ddt0lux2dx+0lut2dx0lf2dx. \frac{d}{dt}\int_0^l u_x^2\,dx + \int_0^l u_t^2\,dx \le \int_0^l f^2\,dx.

tt00 积分到任意 τ(0,T]\tau\in(0,T], 得到

0lux2(x,τ)dx0lux2(x,0)dx+0τ0lut2dxdt0τ0lf2dxdt. \int_0^l u_x^2(x,\tau)\,dx - \int_0^l u_x^2(x,0)\,dx + \int_0^\tau\int_0^l u_t^2\,dx\,dt \le \int_0^\tau\int_0^l f^2\,dx\,dt.

ux(x,0)=φ(x)u_x(x,0)=\varphi'(x), 并放宽上限 τT\tau\leqslant T, 得

0lux2(x,τ)dx+0τ0lut2dxdt0lφ(x)2dx+0T0lf2dxdt. \int_0^l u_x^2(x,\tau)\,dx + \int_0^\tau\int_0^l u_t^2\,dx\,dt \le \int_0^l \varphi'(x)^2\,dx + \int_0^T\int_0^l f^2\,dx\,dt.

τ[0,T]\tau\in[0,T] 取上确界即可得到

sup0tT0lux2(x,t)dx+0T0lut2(x,t)dxdtM[0lφ(x)2dx+0T0lf2dxdt], \sup_{0\leqslant t\leqslant T}\int_0^l u_x^2(x,t)\,dx +\int_0^T\int_0^l u_t^2(x,t)\,dx\,dt \le M\left[\int_0^l\varphi'(x)^2\,dx + \int_0^T\int_0^l f^2\,dx\,dt\right],

其中可取 M=1M=1.

4-17 #

求边值问题

{Δu=f(x,y),(x,y)Ω,uΩ=φ\begin{cases} -\Delta u =f(x,y), & (x,y)\in \Omega,\\ u|_{\partial\Omega}=\varphi & \end{cases}

的 Green 函数, 其中

  • (1) Ω\Omega 是上半平面;
  • (2) Ω\Omega 是第一象限;
note
  • (1)
G(x,y;ξ1,ξ2)=12πln(xξ1)2+(yξ2)2(xξ1)2+(y+ξ2)2,y>0. G(x,y;\xi_1,\xi_2) =-\frac{1}{2\pi} \ln\frac{\sqrt{(x-\xi_1)^2+(y-\xi_2)^2}} {\sqrt{(x-\xi_1)^2+(y+\xi_2)^2}}, \qquad y>0.
  • (2)
G(x,y;ξ1,ξ2)=12π[ln(xξ1)2+(yξ2)2ln(xξ1)2+(y+ξ2)2ln(x+ξ1)2+(yξ2)2+ln(x+ξ1)2+(y+ξ2)2]. \begin{aligned} G(x,y;\xi_1,\xi_2) &=-\frac{1}{2\pi}\big[ \ln\sqrt{(x-\xi_1)^2+(y-\xi_2)^2} -\ln\sqrt{(x-\xi_1)^2+(y+\xi_2)^2}\\ &-\ln\sqrt{(x+\xi_1)^2+(y-\xi_2)^2} +\ln\sqrt{(x+\xi_1)^2+(y+\xi_2)^2} \big]. \end{aligned}

4-19 #

B+(R)={(x,y)x2+y2<R2,y>0}B^+(R)=\{(x,y)|x^2+y^2<R^2,y>0\}, 求定解问题

{Δu=f(x,y),(x,y)B+(R),uB+(R){y>0}=φ(x,y),uyy=0=ψ(x,0),RxR\begin{cases} -\Delta u = f(x,y), & (x,y)\in B^+(R),\\ u|_{\partial B^+(R)\cap\{y>0\}}=\varphi(x,y), & \\ u_y|_{y=0}=\psi(x,0), & -R\leqslant x\leqslant R \end{cases}

的 Green 函数. 如果 uC1(B+(R))C2(B+(R))u\in C^1(\overline{B^+}(R))\cap C^2(B^+(R)) 是上述问题的解, 试给出解的表达式.

note

B(R)={(x,y):x2+y2<R2},B+(R)=B(R){y>0},B(R)=\{(x,y):x^2+y^2<R^2\},\quad B^+(R)=B(R)\cap\{y>0\}, 其边界分解为 B+(R)=ΓDΓN,ΓD=B(R){y>0},ΓN={(x,0):x<R},\partial B^+(R)=\Gamma_D\cup\Gamma_N,\quad \Gamma_D=\partial B(R)\cap\{y>0\},\quad \Gamma_N=\{(x,0):|x|<R\},ΓD\Gamma_D 上为 Dirichlet 条件, 在 ΓN\Gamma_N 上为 Neumann 条件 uy=ψu_y=\psi.

先给出圆盘 B(R)B(R) 的 Dirichlet Green 函数 GRG_R. 用复数记号 z=x+iy, ζ=ξ+iηz=x+iy,\ \zeta=\xi+i\eta, GR(z,ζ)=12πlnzζzR2ζ,z<R, ζ<R.G_R(z,\zeta) =-\dfrac{1}{2\pi}\ln\left| \frac{z-\zeta}{z-\frac{R^2}{\overline\zeta}} \right|, \quad |z|<R,\ |\zeta|<R. 它满足 ΔzGR(z,ζ)=δ(zζ),zB(R);GR(z,ζ)=0,z=R.-\Delta_z G_R(z,\zeta)=\delta(z-\zeta),\quad z\in B(R);\quad G_R(z,\zeta)=0,\quad |z|=R.

设极点 ζ=ξ+iη\zeta=\xi+i\eta 在上半圆内:ζ<R, η>0|\zeta|<R,\ \eta>0. 对其关于 xx 轴的镜像为 ζ=ζ=ξiη.\zeta^*=\overline\zeta=\xi-i\eta.G+(z,ζ):=GR(z,ζ)+GR(z,ζ),zB+(R).G^+(z,\zeta):=G_R(z,\zeta)+G_R(z,\zeta^*),\quad z\in B^+(R).

则有: 对 zB+(R)z\in B^+(R), ζ\zeta^* 在下半圆, 不在 B+(R)B^+(R) 内, 故

ΔzG+(z,ζ)=ΔzGR(z,ζ)=δ(zζ). -\Delta_z G^+(z,\zeta) =-\Delta_z G_R(z,\zeta) = \delta(z-\zeta).

zΓDz\in\Gamma_D, 即 z=R, z>0|z|=R,\ \Im z>0, 有 GR(z,)=0G_R(z,\cdot)=0, 故 G+(z,ζ)=0,zΓD.G^+(z,\zeta)=0,\quad z\in\Gamma_D.

z=xΓNz=x\in\Gamma_N, 可见 G+(x+iy,ζ)=GR(x+iy,ζ)+GR(x+iy,ζ)G^+(x+iy,\zeta)=G_R(x+iy,\zeta)+G_R(x+iy,\zeta^*)

yy 为偶函数, 从而在 y=0y=0G+y(x,0;ζ)=0,x<R,\frac{\partial G^+}{\partial y}(x,0;\zeta)=0,\quad |x|<R, 即在 ΓN\Gamma_N 上满足齐次 Neumann 条件.

因此 G+(z,ζ)G^+(z,\zeta) 即为半圆 B+(R)B^+(R) 上满足

ΔzG+(z,ζ)=δ(zζ),G+ΓD=0,G+nΓN=0 -\Delta_z G^+(z,\zeta)=\delta(z-\zeta),\quad G^+\big|_{\Gamma_D}=0,\quad \frac{\partial G^+}{\partial n}\Big|_{\Gamma_N}=0

的 Green 函数.

P=(ξ,η)B+(R)P=(\xi,\eta)\in B^+(R), Q=(x,y)Q=(x,y). Green 表示公式为

u(P)=B+(R)G+(P,Q)f(Q)dQ+B+(R)(u(Q)G+nQ(P,Q)G+(P,Q)un(Q))dsQ. u(P)=\int_{B^+(R)}G^+(P,Q)\,f(Q)\,dQ +\int_{\partial B^+(R)}\bigl(u(Q)\frac{\partial G^+}{\partial n_Q}(P,Q) -G^+(P,Q)\frac{\partial u}{\partial n}(Q)\bigr)\,ds_Q.

由于在 ΓD\Gamma_Du=φu=\varphi、在 ΓN\Gamma_Nuy=ψu_y=\psi, 且

G+ΓD=0,G+nΓN=0, G^+\big|_{\Gamma_D}=0,\qquad \frac{\partial G^+}{\partial n}\Big|_{\Gamma_N}=0,

上式化简为

u(P)=B+(R)G+(P,Q)f(Q)dQ+ΓDφ(Q)G+nQ(P,Q)dsQ+ΓNG+(P,Q)ψ(Q)dsQ. u(P)=\int_{B^+(R)}G^+(P,Q)\,f(Q)\,dQ +\int_{\Gamma_D}\varphi(Q)\frac{\partial G^+}{\partial n_Q}(P,Q)\,ds_Q +\int_{\Gamma_N}G^+(P,Q)\,\psi(Q)\,ds_Q.

综上, uC1(B+(R))C2(B+(R))u\in C^1(\overline{B^+}(R))\cap C^2(B^+(R)) 的解可表示为

u(P)=B+(R)G+(P,Q)f(Q)dQ+B+(R){y>0}φ(Q)G+nQ(P,Q)dsQ+RRG+((ξ,η);(x,0))ψ(x,0)dx. \begin{aligned} u(P)&=\int_{B^+(R)}G^+(P,Q)\,f(Q)\,dQ +\int_{\partial B^+(R)\cap\{y>0\}}\varphi(Q)\frac{\partial G^+}{\partial n_Q}(P,Q)\,ds_Q\\ &+\int_{-R}^R G^+((\xi,\eta);(x,0))\,\psi(x,0)\,dx. \end{aligned}

讨论

评论

正在加载评论...