使用特征法求解一阶线性偏微分方程时的推导错误排查

你好！你的特征法思路方向是对的，但问题出在特征方程积分时对“常数”的理解错误，以及后续联立方程的逻辑偏差，导致你只得到了一个不满足原PDE的错误特解形式。下面一步步帮你梳理问题并修正推导：

一、先明确特征方程的正确转换（这部分你是对的） #

原PDE：
$$(c_1t+c_2)\frac{\partial u}{\partial t}-\left(c_1x+c_3\right)\frac{\partial u}{\partial x}=c_1u$$
对应的特征方程为：
$$\frac{dt}{c_1t+c_2} = \frac{dx}{-(c_1x+c_3)} = \frac{du}{c_1u}$$
令$r=c_1t+c_2$，$s=c_1x+c_3$，则$dr=c_1dt$，$ds=c_1dx$，代入后特征方程可简化为：
$$\frac{dr}{r} = -\frac{ds}{s} = \frac{du}{u}$$
这一步的转换是完全正确的。

二、核心错误：积分常数不是全局实数，而是特征不变量的函数 #

特征法的关键是：从特征方程积分得到的“常数”并非固定实数，而是沿特征线保持不变的量（特征不变量）的任意可微函数，你错误地把它当成了全局固定常数，导致推导偏离正确方向。

正确的积分步骤： #

1. 先求特征不变量
   从$\frac{dr}{r} = -\frac{ds}{s}$积分：
   $$\int \frac{dr}{r} = -\int \frac{ds}{s} \implies \ln|r| + \ln|s| = \ln|\xi| \implies \xi = rs = (c_1t+c_2)(c_1x+c_3)$$
   这里$\xi$就是特征不变量——沿着每条特征线，$\xi$的值固定不变，不同特征线对应不同的$\xi$。

2. 再联立$\frac{dr}{r} = \frac{du}{u}$求通解
   积分$\frac{dr}{r} = \frac{du}{u}$时，注意这里的“积分常数”不是固定实数，而是依赖于特征不变量$\xi$的任意可微函数（因为不同特征线的“常数”可以不同）：
   $$\ln|u| = \ln|r| + \ln|f(\xi)| \implies u = r \cdot f(\xi)$$
   其中$f$是任意一元可微函数。

把$\xi=rs$代入，就得到原PDE的通解：
$$\boldsymbol{u = (c_1t+c_2)f\left((c_1t+c_2)(c_1x+c_3)\right)}$$

你也可以用另一种形式表示通解：从$-\frac{ds}{s} = \frac{du}{u}$积分得$u = \frac{g(\xi)}{s}$，这和上面的通解是等价的（只需令$g(\xi) = \xi f(\xi)$即可转换）。

三、验证通解的正确性 #

我们来验证这个通解满足原PDE：
令$u=(c_1t+c_2)f(\xi)$，$\xi=(c_1t+c_2)(c_1x+c_3)$，计算偏导数：

• $\frac{\partial u}{\partial t} = c_1f(\xi) + (c_1t+c_2)f'(\xi) \cdot c_1(c_1x+c_3)$
• $\frac{\partial u}{\partial x} = (c_1t+c_2)f'(\xi) \cdot c_1(c_1t+c_2)$

代入原PDE左边：
$$
\begin{align*}
&(c_1t+c_2)\frac{\partial u}{\partial t} - (c_1x+c_3)\frac{\partial u}{\partial x} \
=& (c_1t+c_2)\left[c_1f(\xi) + c_1(c_1t+c_2)(c_1x+c_3)f'(\xi)\right] - (c_1x+c_3)\left[c_1(c_1t+c_2)^2f'(\xi)\right] \
=& c_1(c_1t+c_2)f(\xi) + c_1(c_1t+c_2)^2(c_1x+c_3)f'(\xi) - c_1(c_1t+c_2)^2(c_1x+c_3)f'(\xi) \
=& c_1(c_1t+c_2)f(\xi) = c_1u
\end{align*}
$$
完全等于原PDE的右边，说明这个通解是正确的。

四、你的错误根源总结 #

1. 积分常数的理解偏差：你把特征方程积分得到的“常数”当成了全局固定实数（比如$A,B,k$），但特征法中这个“常数”必须是特征不变量的任意函数——因为不同特征线上的积分常数可以不同，由特征不变量决定。
2. 联立方程的逻辑错误：你从$u=Ar$和$u=B/s$联立得到$u=k\frac{r}{s}$，忽略了$A$和$B$并非独立的固定常数，而是由同一个特征不变量绑定的函数，不能直接用固定比例联立。

这样修正后，就能得到满足原PDE的通解啦！

备注：内容来源于stack exchange，提问作者Gomes93