题六·方形区域拉普拉斯(B 类模板)
抄题(0≤x≤a, 0≤y≤b):
uxx+uyy=0,ux∣x=0=0, u∣x=a=2y,uy∣y=0=0, uy∣y=b=0.
四条边里只有 x=a 一条非齐次,其余齐次 → 不必叠加,直接分离。
挑齐次方向的诀窍:选那条”两端都齐次”的轴做本征值。 这里 y 轴两端都是齐次 Neumann。
1. y 方向本征值问题
Y′′+μY=0,Y′(0)=Y′(b)=0(N–N)→ 查 本征值问题速查表(注意含零本征值):
μn=(bnπ)2,Yn=cosbnπy,n=0,1,2,…(Y0=1).
2. x 方向方程 + 齐次端
令 u=∑n≥0Xn(x)cosbnπy,代回得 Xn′′−μnXn=0,配 Xn′(0)=0:
- n=0:X0′′=0⇒X0=A0+B0x,X0′(0)=0⇒X0=A0。
- n≥1:Xn=Ancoshbnπx+Bnsinhbnπx,Xn′(0)=0⇒Bn=0。
故 u=A0+n≥1∑Ancoshbnπxcosbnπy。
3. 用 x=a 边定系数(2y 的余弦展开)
u∣x=a=2y。零模 A0=b1∫0b2ydy=b。对 n≥1,
Ancoshbnπa=b2∫0b2ycosbnπydy=n2π24b[(−1)n−1]=⎩⎨⎧0,−n2π28b,n 偶n 奇.
(用 ∫0bycosbnπydy=n2π2b2[(−1)n−1],模方 ∫0bcos2=2b。)
4. 答案
u(x,y)=b−n 奇∑n2π2coshbnπa8bcoshbnπxcosbnπy
题五同型但四边非全齐次 → 先按 齐次化三套路 套路 B 拆成两个”单边非齐次”子问题各做一遍再相加。