三维空间的波动问题——泊松公式（婉转曲折，其妙无穷）

凉月松心

已于 2023-04-24 21:57:04 修改

阅读量2.6k

点赞数 3

分类专栏：数学物理方程文章标签： django python 后端机器学习算法

于 2022-10-09 22:55:09 首次发布

本文链接：https://blog.csdn.net/m0_52555753/article/details/127234948

版权

数学物理方程专栏收录该内容

6 篇文章 3 订阅

订阅专栏

简单分析

运用了平均值的思想，三维有x,y,z，我们做了一个平均值，就变成了一维的波动问题

这个思想以及方法非常关键

我们列出方程

$\small \mu_{tt}-a^2 \nabla^2\mu=0$

$\small \mu(M,0)=\phi(M)$

$\small \mu_t(M,0)=\psi(M),M=(x,y,z)$

$\small -\infty < x,y,z <+\infty$

无界意味着我们不考虑反射波的影响，我们要求的是 $\small \mu(M,t)$

这个方向各向同性，x，y，z三个方向是等价的

平均值方法

场量在某一点，以 $\small M_0$ 为秋心， $\small r$ 为半径，把球面上的值加起来，然后再除以总的面积，这个就是球面上u的平均值，我们再让 $\small r->0$ ，就是球心的值

这种方法可以很好的处理各向同性的问题

$\small \overline \mu(r,t)=\frac{1}{4 \pi r^2}\iint_{S_l^{M_0}}\mu(M,t)ds$ ,我们让其趋于0，就得到球心的值

$\small \mu(M,t)=\lim_{r->0}\frac{1}{4 \pi}\iint \mu(M,t)(\frac{ds}{r^2})=\lim_{r->0}\frac{1}{4 \pi}\iint \mu(M,t)d\Omega$

这个球面平均实际就是对立体角的积分，最后实际上是对于一个立体角的一个平均，球极坐标有 $\small (r,\theta, \phi)$ ，对立体角的积分实际上是把 $\small \theta,\phi$ 给积分掉了，我们可以联想到 $\small \int_{0}^{\pi} \int_0^{2 \pi} sin(\theta) d\theta d\phi=4 \pi$

也就是总的立体角等于 $\small 4 \pi$ ，对角度的积分实际是把这一部分积分掉

所以最后只剩下r的变量

在球极坐标里面，关系式

$\small x-x_0=rsin(\theta)cos(\phi)$

$\small y-y_0=rsin(\theta)sin(\phi)$

$\small z-z_0=rcos(\theta)$

求解

在泛定方程两边做球面平均，我们得到

$\small \frac{1}{4 \pi}\iint_{S_l^{M_0}}\mu_{tt}d\Omega=\frac{a^2}{4 \pi}\iint_{S_l^{M_0}}\nabla^2 \mu d\Omega$ ，对时间的二阶导数和球面上的积分是可以交换次序的

$\small \frac{\partial^2 }{\partial t^2}[ \frac{1}{4 \pi}\iint_{S_l^{M_0}}\mu d\Omega]-a^2 \nabla^2[\frac{1}{4 \pi} \iint_{S_l^{M_0}}ud\Omega]=0$

我们得到 $\small \frac{\partial^2 }{\partial t^2} \overline \mu (r,t)=a^2\nabla \overline \mu (r,t)$

$\small \nabla^2 \overline \mu (r,t)=\frac{1}{r}\frac{\partial^2 }{\partial t^2}(r \overline \mu)$

进一步，我们得到

$\small \frac{\partial^2 }{\partial t^2} \overline \mu (r,t)=\frac{a^2}{r}\frac{\partial^2 r \overline \mu}{\partial r^2}$ ,这个方程就变成

$\small \frac{\partial^2 }{\partial t^2}(r \overline \mu)=a^2 \frac{\partial^2 }{\partial r^2}(r \overline \mu)$ ,我们可以做一个变量代换,令 $\small r\overline \mu=v(r,t)$

我们可以得到 $\small v_{tt}-a^2v_{rr}=0$

它的解:

$\small f(\zeta,\eta)=f_1(\zeta)+f_2(\eta)$

其中 $\small \zeta=r+at,\eta=r-at$ ，我们可以得到

$\small v (x,t) = f_1 (\zeta)+f_2 ( \eta)<-\frac{\partial^2 v}{\partial \zeta \partial \eta}=0$

$\small =f_1(r+at)+f_2(r-at)=r\overline \mu(r,t)$ ，前者是向球心汇聚的球面波，后者是向球心发散的球面波

但是这只是一个通解，我们还需要做一个变换

我们有一个边界条件，当 $\small r=0$ 的时候， $\small v|_{r=0}=0=> f_1(at)+f_2(-at)=0........1$

$\small \overline \mu=\frac{f_1(r+at)+f_2(r-at)-[f_1(at)+f_2(at)]]}{r-0}$

令 $\small r->0,t=t_0 => \mu(M_0,t_0)=\lim_{r->0}\frac{\Delta (f_1+f_2)}{\Delta r}={f_1}'(at)+{f_2}'(-at)$

我们发现 $\small {f_1}'(at)={f_2}'(-at)$ （从1号等式发现）

$\small \mu(M_0,t_0)=2{f_1}'(at_0)$

泊松公式

我们对于 $\small f_1(r+at)+f_2(r-at)=r\overline \mu(r,t)$ 两边求导，列出两个等式

$\small \frac{\partial }{\partial r} r\overline \mu={f_1}'(r+at)+{f_2}'(r-at)=>r=at_0,t=0$

$\small \frac{1}{a}\frac{\partial }{\partial t} r\overline \mu={f_1}'(r+at)-{f_2}'(r-at)=>t=0,r=at_0$

我们需要利用初始条件两者一相加得到

$\small =\frac{\partial }{\partial r} r\overline \mu+\frac{1}{a} \frac{\partial }{\partial r} r\overline \mu |_{r=at_0}$

我们再进行计算，就可以得到泊松公式

$\small \mu(M,t)=\frac{1}{4 \pi a}[\frac{\partial }{\partial t} \iint_{S_{at}^{M}} \frac{\phi({M'})}{at}]d{S}'+\iint_{S_{at}^{M}} \frac{\psi({M'})}{at}]d{S}'$