美文网首页我的大学简友广场如切如磋,如琢如磨
极坐标下二维拉普拉斯方程和泊松方程的求解

极坐标下二维拉普拉斯方程和泊松方程的求解

作者: Raow1 | 来源:发表于2021-04-01 23:22 被阅读0次
1. 一般形式

对于方程\begin{cases} \Delta u =f(r,\varphi) \\ (\alpha u+ \beta u_r)|_{r=R_0}=u_0(\varphi) \end{cases}如果f(r,\varphi) \equiv 0,则方程为拉普拉斯方程,否则为泊松方程

如果\beta =0,则称其为狄利克雷问题;如果\alpha=0,则称其为诺伊曼问题;否则为混合问题。

此外,根据r的范围,还可分为圆域内,圆环,圆外三种。

2. 拉普拉斯方程

使用分离变量的方法,我们能够得出如下结论:

  1. 对于圆域内的问题,其解的形式为u(r,\varphi)=C_0+\displaystyle\sum_{n=1}^{\infty}r^n(A_n\cos n\varphi +B_n\sin n\varphi)
  2. 对于圆外的问题,其解的形式为u(r,\varphi)=C_0+\displaystyle\sum_{n=1}^{\infty}\frac{1}{r^n}(A_n\cos n\varphi +B_n\sin n\varphi)
  3. 对于圆环上的问题,其解的形式为u(r,\varphi)=C_0+D_0\ln r+\displaystyle\sum_{n=1}^{\infty}(C_nr^n+\frac{D_n}{r^n})\cos n\varphi + \displaystyle\sum_{n=1}^{\infty}(A_nr^n+\frac{B_n}{r^n})\sin n\varphi

所以,对于拉普拉斯方程,我们按形式设解,代入边值条件即可。(可能需要对边值条件中的u_0(\varphi)进行傅里叶展开)

关于拉普拉斯方程三个不同形式的解:显然圆环上的解的形式是最全的。
对于圆域内的问题,因为r可以取到0,如果留下\ln rr^{-n}项,则此时当r \to 0时,u\to \infty,这与实际的物理含义不符。
对于圆域外的问题,同理留下\ln rr^n会使得在无穷远处u \to \infty,也是与实际不符。

3. 泊松方程

显然,泊松方程是拉普拉斯方程的非齐次形式。我们的求解步骤为:找特解,转化为拉普拉斯方程,求拉普拉斯方程,综合。

对于f(r,\varphi)=ar^n\cos m\varphi或是f(r,\varphi)=ar^n\sin m\varphi

1.如果(n+2)^2\neq m^2

则存在特解u^{*}=br^{n+2}\cos m\varphiu^{*}=br^{n+2}\sin m\varphi

其中,b=\frac{a}{(n+2)^2-m^2}

2.如果(n+2)^2 = m^2

则设特解u^{*}=\varphi(r)\cos m\varphiu^{*}=\varphi(r)\sin m\varphi

求出特解后,作变换v=u-u^{*}即可将泊松方程转换为拉普拉斯方程,然后求解拉普拉斯方程得到v,从而u=v+u^{*}

相关文章

  • 极坐标下二维拉普拉斯方程和泊松方程的求解

    1. 一般形式 对于方程如果,则方程为拉普拉斯方程,否则为泊松方程。 如果,则称其为狄利克雷问题;如果,则称其为诺...

  • 泊松图像编辑

    理论基础:拉普拉斯算子泊松方程泊松编辑、泊松融合泊松图像编辑matlab注:白马负金羁的博客很有意思 论文Pois...

  • 拉普拉斯和泊松方程题

    2016-1-10. 求解下列拉普拉斯方程的狄利克雷问题 由题, 由积化和差公式,易得 同理, 所以, 所以, 2...

  • 使用scipy库的root和fsolve函数求解方程

    求解非线性方程 使用scipy库的 root, fsolve 函数求解非线性方程。 求解传热方程(热辐射+热传导)...

  • 高中数学基础20190307

    圆的参数方程,圆的极坐标方程,圆的普通方程,能否分辨清楚。你看看上面的方程分别是圆的什么方程?一定要搞清楚。

  • 线性方程组(四)- 矩阵方程

    小结 矩阵方程的定义 矩阵方程的求解 矩阵方程、向量方程和线性方程组拥有相同的解集 的计算、行-向量规则和性质 若...

  • 2020-01-16

    机器学习概念统计学样本 —— 数据集待求解方程(待求解方程中的未知数为参数) —— 预训练模型求解参数后的N元方程...

  • [现代控制理论]2_state-space状态空间方程

    state-space状态空间方程 质量-弹簧-阻尼经典系统的微分方程,拉普拉斯变换 介绍了状态空间方程的表示 推...

  • 二分法

    本次使用二分法求解一个简单的方程的根 方程如下 方程的定义域为[0,1] 编程求解 ``` def bisect(...

  • 弦振动方程的解法(分离变量法+行波法)

    弦振动方程的解法 微分方程基础 本文所言的方程基本上是微分方程,而非中学阶段所言的代数方程,微分方程的求解目标一般...

网友评论

    本文标题:极坐标下二维拉普拉斯方程和泊松方程的求解

    本文链接:https://www.haomeiwen.com/subject/bjuxkltx.html