欢迎光临散文网 会员登陆 & 注册

6.1求解笛卡尔坐标系下的拉普拉斯方程

2023-07-16 14:45 作者:JJGYLevel2  | 我要投稿

6 拉普拉斯方程求解

在本节,我们将解开拉普拉斯方程在笛卡尔坐标系和球面坐标系下。前者解对平面或者区域应用是十分有用的。后者的解是全球方案。两者都将发展为以正交基函数表示(分别为傅里叶序列和球谐序列)

拉普拉斯方程的解是最重要的一步去解决边值问题。第二步将边界方程用同一组级数表示并确定这些级数的系数。

6.1 笛卡尔坐标系

考虑到以x和y为水平坐标,z为垂直坐标情况, z<0意味着位置在地球内部,z>0在地球外部(这里可能会有疑惑?)。我们的任务就是找出在z>0的解。对于这组微分方程,有个很重要与实用的策略就是分离变量法。这是一个很好的方法,常用于求解微分方程。分离变量法本身的意义就是建立互相正交的基。这是我一个很朴素但是又贴合实际的解释。首先,我们把这个拉普拉斯方程的解三个函数的乘积,这三个函数分别为三个独立变量的函数。

拉普拉斯

这里短暂回顾一下拉普拉斯算子的形式:

拉普拉斯算子

由此可以看出,拉普拉斯算子是个对标量才起作用的,输出标量的算符。

将算子用在公式6.1上,我们得到

拉普拉斯方程分离变量后

为了方便表示,我们采用了简单的形式 (), 并对公式进行调整得到以下结果。

简化形式

变量可以很清楚地分离出被求导的变量。我们不加解释的直接先定各个子方程的常数。

子方程

6.3的方程都是二阶常微分方程(ODE),因此直接可以求得它们的基础解

基础解系

当然,每一个解都可以跟一个常数或者说振幅相结合。这个6.2方程的解将会表示为f,g,h的乘积形式。但是对于每一个n和m,我们都会有一组新的解,因此,我们不得不叠加每个n和m的所有可能组合的所有解。

组合

用复数指数去替换这些sin和cos基函数,我们可以得到更紧凑的表达式:

紧凑形式

由欧拉公式 

欧拉公式

三角恒等变换

三角恒等变换

如果看不明白,可以尝试把n,m组成的四个象限的情况分开求和即可。


6.1求解笛卡尔坐标系下的拉普拉斯方程的评论 (共 条)

分享到微博请遵守国家法律