带指定边界条件的二阶偏微分方程求解思路咨询

你选择的w = u + 3变量代换方向完全正确，出现4个无法确定的待定常数，是代换后的分离变量流程操作有误，和代换思路本身无关。

原问题梳理 #

待求解方程为：
$$\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = u + 3$$
对应第一类齐次Dirichlet边界条件：

• 当$0<y<1$时，$u(0,y)=0$、$u(1,y)=0$
• 当$0<x<1$时，$u(x,0)=0$、$u(x,1)=0$

代换后的形式验证 #

将w = u + 3代入原方程，常数项求导为0，可直接化简得到w满足的齐次Helmholtz方程：
$$\frac{\partial^2 w}{\partial x^2} + \frac{\partial^2 w}{\partial y^2} = w$$
注意代换后w的边界条件不再是齐次的：

• 当$0<y<1$时，$w(0,y)=3$、$w(1,y)=3$
• 当$0<x<1$时，$w(x,0)=3$、$w(x,1)=3$

待定常数无法确定的原因 #

你得到4个待定常数，本质是把分离变量得到的单组特解形式w=X(x)Y(y)直接当成了方程通解：

1. 对w的齐次方程分离变量后，会得到x方向、y方向各带2个待定常数的二阶常微分方程通解，合计4个常数
2. 齐次线性偏微分方程在矩形域上的通解是无穷多组分离特解的级数和，不是单组分离特解，仅靠单组特解无法匹配4段边界条件，自然会出现常数无法确定的问题
3. 这类问题的正确处理逻辑是先求解对应齐次边界的本征值问题，得到正交本征函数系，再将非齐次边界条件展开成本征函数的级数形式，逐阶匹配系数即可消去不确定的常数，最终得到收敛的级数解。

如果觉得代换后处理非齐次边界麻烦，也可以跳过代换，直接将齐次边界下的u展开为双重正弦级数代入原方程，直接求解级数系数，两种路径最终得到的结果完全一致。

内容的提问来源于stack exchange，提问作者spaldix