二阶偏微分方程的化简思路

本文主要是对顾樵老师 数物方法 一书对应章节的内容的梳理(主要为了抛砖引玉),有一些自己的理解,如有不妥,还请慷慨指出。

化简的理论

这里所说的二阶偏微分方程主要是指二阶线性双变量偏微分方程,它的一般形式如下所示:
A2ux2+2B2uxy+C2uy2+Duy+Fu=G
其中大写的英文字母都是x,y的函数,这是一个非齐次方程,如果G=0,那么它就变成一个齐次的方程。

为了化简上述方程,我们进行变量的替换: ξ=ξ(x,y); η=η(x,y). 需要注意的是这里进行变量替换的目的是为了化简上述方程,它可以用图表示为:

而在多元函数的微积分中我们也有过进行变量替换的操作(主要是为了进行积分的方便): x=x(ξ,η); y=y(ξ,η), 它的图为:

需要注意这两种变量替换是不同的,在多元函数的微积分中的变量替换一般是线性的,而这里为了化简方程所进行的变量替换一般不是线性的。

进行变量替换后的方程可以写为:
a2uξ2+2b2uξη+c2uη2+duη+fu=g(1)
其中的小写字母如下所示:
小写项
双曲型的: 2uxy= 或者 2ux22uy2=
抛物线型的: 2ux2= 或者 2uy2=
椭圆型的:2ux2+2uy2=
根据以上的想法,我们可以看到我们定义的变换ξ,η只需要使得(1)式中的二阶导数项中有一项或者两项的系数为0就可以达到目的。
这里a,c有相同的形式,我们先考虑它,进而考虑齐次方程
A(Wx)2+2BWxWy+C(Wy)2=0(2)
W不是常数时,两边同时除以(Wy)2, 可以得到下面的方程:
A(WxWy)2+2BWxWy+C=0
此时,如果y=y(x),那么根据隐函数定理有dydx=WxWy
那么有:
A(dydx)22B(dydx)+C=0
这是一个关于y的二次方程,它的解为:
y=B+B2ACA 以及 y=BB2ACA
那么根据Δ=B2AC的情况我们有:
a. Δ>0时,有两个不同的dydx,从而解出两个不同的y1=y1(x)y2=y2(x), 那么通过整理可以得到ξ(x,y1)=γ1以及η(x,y2)=γ2
其中γ1γ2是任意的两个常数,它们的存在保证了ξ(x,y1)η(x,y2)使得a=c=0, 从而达到将(1)式化简为双曲型的目的。
b. Δ=0时,有两个相同的dydx, 从而解出一种 y=y(x), 整理得到ξ(x,y)=γ1这里 γ1的存在使得a=0(当然也可以是c=0, 但是它俩只有一个等于0)。另外,将B2=AC带入上面求b的式子中我们可以得到b=0, 结合a=0, 我们可以将(1)式化简为抛物线型的式子。
c. Δ<0时,dydx有一对共轭复函数解,从而解出一对共轭的复函数y1=y1(x)以及y2=y2(x)。这里假设得到的其中一个整理过后的复函数解为ψ(x,y)=ψ1(x,y)+iψ2(x,y), 这里ψ(x,y)满足方程(2),即A(ψx)2+2Bψxψy+C(ψy)2=0 , 根据这个等式,把ψ换为ψ1+iψ2然后进行化简,最后根据实部和虚部都为0,可以得到下列方程:
实部虚部为0
这个方程结合前面ab,c的方程可以得出,如果取ξ=ψ1(x,y), η=ψ2(x,y), 那么有a=cb=0,那么方程(1)就可以化简为椭圆型的方程。

化简的步骤

根据上面介绍的化简的理论,下面我们来总结一下给定一个二阶线性双变量的偏微分方程,我们如何对其进行化简:
①. 首先识别出给定的二阶偏微分方程中对应上面的A,B,C,然后计算出Δ
②. 根据计算出的Δ求解出y=y(x), 然后整理成为ξ(x,y)=γ1η(x,y)=γ2 (Δ=0 时,η(x,y)可以任取一个与ξ(x,y)不相关的函数)
③. 利用②中的结果求解出其他系数(b,d,e),最后带入到变换后的方程(1)中得到结果。

posted @   SiranLee  阅读(1743)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· winform 绘制太阳,地球,月球 运作规律
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· AI与.NET技术实操系列(五):向量存储与相似性搜索在 .NET 中的实现
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
点击右上角即可分享
微信分享提示