傅里叶方法求解偏微分方程2
题目
使用傅里叶方法求解正方形区域 0<x<π 0 < x < \pi 0<x<π,0<y<π 0 < y < \pi 0<y<π 中以下问题的所有解:
4uxx−9uyy=0,
4u_{xx} - 9u_{yy} = 0,
4uxx−9uyy=0,
边界条件:
- 在 y=0 y = 0 y=0 和 y=π y = \pi y=π 时,u=0 u = 0 u=0,
- 在 x=0 x = 0 x=0 和 x=π x = \pi x=π 时,u=0 u = 0 u=0.
问题解决
给定偏微分方程(PDE):
4∂2u∂x2−9∂2u∂y2=0,0<x<π,0<y<π,
4 \frac{\partial^2 u}{\partial x^2} - 9 \frac{\partial^2 u}{\partial y^2} = 0, \quad 0 < x < \pi, \quad 0 < y < \pi,
4∂x2∂2u−9∂y2∂2u=0,0<x<π,0<y<π,
边界条件为齐次 Dirichlet 边界条件:
- u(x,0)=u(x,π)=0 u(x, 0) = u(x, \pi) = 0 u(x,0)=u(x,π)=0 对所有 x∈[0,π] x \in [0, \pi] x∈[0,π],
- u(0,y)=u(π,y)=0 u(0, y) = u(\pi, y) = 0 u(0,y)=u(π,y)=0 对所有 y∈[0,π] y \in [0, \pi] y∈[0,π].
这是一个线性齐次 PDE,且区域为正方形,因此使用分离变量法(Fourier 方法)求解。假设解的形式为:
u(x,y)=X(x)Y(y).
u(x, y) = X(x) Y(y).
u(x,y)=X(x)Y(y).
代入 PDE 得:
4X′′(x)Y(y)−9X(x)Y′′(y)=0.
4 X''(x) Y(y) - 9 X(x) Y''(y) = 0.
4X′′(x)Y(y)−9X(x)Y′′(y)=0.
整理为:
4X′′X=9Y′′Y.
4 \frac{X''}{X} = 9 \frac{Y''}{Y}.
4XX′′=9YY′′.
令该常数等于 −λ -\lambda −λ(为满足边界条件,通常取负常数),则:
4X′′X=−λ,9Y′′Y=−λ.
4 \frac{X''}{X} = -\lambda, \quad 9 \frac{Y''}{Y} = -\lambda.
4XX′′=−λ,9YY′′=−λ.
于是得到两个常微分方程(ODE):
- X′′+λ4X=0 X'' + \frac{\lambda}{4} X = 0 X′′+4λX=0, 边界条件 X(0)=X(π)=0 X(0) = X(\pi) = 0 X(0)=X(π)=0,
- Y′′+λ9Y=0 Y'' + \frac{\lambda}{9} Y = 0 Y′′+9λY=0, 边界条件 Y(0)=Y(π)=0 Y(0) = Y(\pi) = 0 Y(0)=Y(π)=0.
求解特征值问题
对于 X(x) X(x) X(x) 的 ODE:
X′′+μX=0,μ=λ4,X(0)=X(π)=0.
X'' + \mu X = 0, \quad \mu = \frac{\lambda}{4}, \quad X(0) = X(\pi) = 0.
X′′+μX=0,μ=4λ,X(0)=X(π)=0.
这是标准的 Sturm-Liouville 问题。特征值和特征函数为:
μn=n2,Xn(x)=sin(nx),n=1,2,3,…
\mu_n = n^2, \quad X_n(x) = \sin(n x), \quad n = 1, 2, 3, \dots
μn=n2,Xn(x)=sin(nx),n=1,2,3,…
因此:
λ4=n2 ⟹ λ=4n2.
\frac{\lambda}{4} = n^2 \implies \lambda = 4n^2.
4λ=n2⟹λ=4n2.
对于 Y(y) Y(y) Y(y) 的 ODE:
Y′′+νY=0,ν=λ9,Y(0)=Y(π)=0.
Y'' + \nu Y = 0, \quad \nu = \frac{\lambda}{9}, \quad Y(0) = Y(\pi) = 0.
Y′′+νY=0,ν=9λ,Y(0)=Y(π)=0.
特征值和特征函数为:
νm=m2,Ym(y)=sin(my),m=1,2,3,…
\nu_m = m^2, \quad Y_m(y) = \sin(m y), \quad m = 1, 2, 3, \dots
νm=m2,Ym(y)=sin(my),m=1,2,3,…
因此:
λ9=m2 ⟹ λ=9m2.
\frac{\lambda}{9} = m^2 \implies \lambda = 9m^2.
9λ=m2⟹λ=9m2.
匹配特征值
由于 λ \lambda λ 必须同时满足两个方程,即:
4n2=9m2 ⟹ (2n)2=(3m)2 ⟹ 2n=3m(取正值,因为 n,m>0).
4n^2 = 9m^2 \implies (2n)^2 = (3m)^2 \implies 2n = 3m \quad (\text{取正值,因为 } n, m > 0).
4n2=9m2⟹(2n)2=(3m)2⟹2n=3m(取正值,因为 n,m>0).
令 m=2k m = 2k m=2k, n=3k n = 3k n=3k,其中 k=1,2,3,… k = 1, 2, 3, \dots k=1,2,3,…(因为 m m m 和 n n n 必须是正整数)。则:
λ=4(3k)2=9(2k)2=36k2.
\lambda = 4(3k)^2 = 9(2k)^2 = 36k^2.
λ=4(3k)2=9(2k)2=36k2.
对应的特征函数为:
X3k(x)=sin(3kx),Y2k(y)=sin(2ky).
X_{3k}(x) = \sin(3k x), \quad Y_{2k}(y) = \sin(2k y).
X3k(x)=sin(3kx),Y2k(y)=sin(2ky).
因此,基本解为:
uk(x,y)=sin(3kx)sin(2ky),k=1,2,3,…
u_k(x, y) = \sin(3k x) \sin(2k y), \quad k = 1, 2, 3, \dots
uk(x,y)=sin(3kx)sin(2ky),k=1,2,3,…
通解
所有解为这些基本解的线性组合:
u(x,y)=∑k=1∞cksin(3kx)sin(2ky),
u(x, y) = \sum_{k=1}^{\infty} c_k \sin(3k x) \sin(2k y),
u(x,y)=k=1∑∞cksin(3kx)sin(2ky),
其中 ck c_k ck 是任意常数。
验证解
- PDE 满足:计算二阶导数:
uxx=−(3k)2sin(3kx)sin(2ky)=−9k2uk, u_{xx} = - (3k)^2 \sin(3k x) \sin(2k y) = -9k^2 u_k, uxx=−(3k)2sin(3kx)sin(2ky)=−9k2uk,
uyy=−(2k)2sin(3kx)sin(2ky)=−4k2uk. u_{yy} = - (2k)^2 \sin(3k x) \sin(2k y) = -4k^2 u_k. uyy=−(2k)2sin(3kx)sin(2ky)=−4k2uk.
代入 PDE:
4(−9k2uk)−9(−4k2uk)=−36k2uk+36k2uk=0. 4(-9k^2 u_k) - 9(-4k^2 u_k) = -36k^2 u_k + 36k^2 u_k = 0. 4(−9k2uk)−9(−4k2uk)=−36k2uk+36k2uk=0. - 边界条件满足:
- 当 y=0 y = 0 y=0 或 y=π y = \pi y=π 时,sin(2ky)=0 \sin(2k y) = 0 sin(2ky)=0,故 u=0 u = 0 u=0。
- 当 x=0 x = 0 x=0 或 x=π x = \pi x=π 时,sin(3kx)=0 \sin(3k x) = 0 sin(3kx)=0(因为 3kπ 3k\pi 3kπ 是 π \pi π 的整数倍),故 u=0 u = 0 u=0.
解的唯一性
当 m m m 和 n n n 不满足 2n=3m 2n = 3m 2n=3m 时,λ \lambda λ 无法同时满足两个特征值问题,导致只有零解。例如:
- 若 n=1 n = 1 n=1,则 λ=4 \lambda = 4 λ=4,代入 Y Y Y 的方程:Y′′+49Y=0 Y'' + \frac{4}{9} Y = 0 Y′′+94Y=0,边界条件 Y(0)=Y(π)=0 Y(0) = Y(\pi) = 0 Y(0)=Y(π)=0,仅零解(因为 49 \frac{4}{9} 94 不是特征值)。
- 类似地,其他不满足 2n=3m 2n = 3m 2n=3m 的组合均无解。
零特征值 (λ=0 \lambda = 0 λ=0) 也仅得零解。因此,所有非平凡解均由上述级数给出。
答案
所有解为:
u(x,y)=∑k=1∞cksin(3kx)sin(2ky)
\boxed{u(x, y) = \sum_{k=1}^{\infty} c_k \sin(3k x) \sin(2k y)}
u(x,y)=k=1∑∞cksin(3kx)sin(2ky)
其中 ck c_k ck 是任意常数,k=1,2,3,… k = 1, 2, 3, \dots k=1,2,3,…。1