网站建设资讯

NEWS

网站建设资讯

球谐函数c语言,球谐函数数学物理方法

球谐波函数怎么计算

向量球谐函数(Vector spherical harmonics)是应用于球坐标系的拉普拉斯方程式的向量解,是球谐函数的向量衍伸形式。在必须计算向量场的电动力学等领域中被广泛应用。

番禺网站建设公司创新互联,番禺网站设计制作,有大型网站制作公司丰富经验。已为番禺成百上千家提供企业网站建设服务。企业网站搭建\成都外贸网站建设要多少钱,请找那个售后服务好的番禺做网站的公司定做!

定义

在球坐标系下,拉普拉斯算符作用在一三维向量场上可以写为

利用分离变数法可以将此一方程式的解分解为一系列本征函数的线性组合

其中的径向解与标量球谐函数相同,而为一与角度相关的向量解,也就是向量球谐函数。

向量球谐函数依用途有很多定义方式。这边我们依照 Barrera 等人的定义,以对球谐函数Yℓm(θ, φ)为基础,将三个向量球谐函数表示为

这边是对应球座标 (r, θ, φ) 的向量,而则为其单位向量。

主要特性

依照上述 Barrera 的定义,向量球谐函数有以下特性:

对称性

与球谐函数相同,向量球谐函数有对称性

星号 * 代表共轭函数。

正交性

三种向量球谐函数彼此两两正交

另外同种类的球谐函数的内积为:

标量场的梯度

对一个标量场,若其多极展开可表示为:

则其梯度可以向量球谐函数表示为:

散度

三种向量球谐函数之散度分别为:

其中为球谐函数之径向分布,为球谐函数。

球谐函数方程

“球谐函数”,可能大家非数理专业说这是什么鬼东西,没错它正是特殊函数,也是拉普拉斯方程的球坐标系形式解的角度部分,而球谐函数在量子力学与动力学等方面有着广泛应用,可以解决实际生活问题与数理模型。

接下来八一从它的定义讲解开始,它有实数形式与复数形式的表达,并列出相关性质,这里我列出几项球谐函数值,在百度百科上有列出详细球谐函数表,球谐函数的宇称两种关系,对相反的两个方向证明以及相关三个球谐积分对称性的证明,然后给出了对于归一化后的球谐函数的正交形式表达,以及广义傅里叶级数的展开和非轴对称拉普拉斯方程的通解,最后给你们练习的三道习题。

球谐函数的基础数学理论

球谐函数在图形学光照计算等领域有着重要应用,因为目前在实际工作中接触较少,所以对其的理解仅仅停留在表面,本着越是基础的东西,其重要性越高的想法,特此开篇文章对其背后的数学理论进行拆解,拆解过程参考了大量其他同学的工作,相应链接在文末的参考文献中有列出,引用过程中如有表述不清晰的内容,可以通过原文辅助阅读。

调和函数指的是一种特殊的二阶连续可导函数(简称C2,在某个定义域存在二阶导数,且二阶导数连续),数学符号用 表达,其中 是 (表示n维实数域)的一个开子集(相当于一维数据空间中的开区间),其特殊在于需要满足拉普拉斯方程(下面有介绍),用(笛卡尔坐标系下)数学表达式来描述的话,就是对于任意 ,需要满足如下的二阶偏微分方程:

这里来回顾一下微分方程的相关知识,单个变量下,也就是一元变量情况下,函数与函数各阶导数组成的微分方程叫做常微分方程:

多元函数而言,函数以及函数对各个自变量的各阶偏导数组成的微分方程叫做偏微分方程:

这个公式也经常以如下的形式出现( 称为拉普拉斯算子, 称为向量微分算子,也就是nabla算子):

其中 叫做拉普拉斯算子,光看定义太抽象,我们来举个例子吧,下面两个函数都是二元的调和函数:

拉普拉斯方程也被称为调和方程、位势方程,这是一种偏微分方程,因为其可以用势函数的形式来描述电磁场、引力场、流场(统称为保守场或者有势场)的性质而被广泛应用。

笛卡尔坐标系下的表述形式前面已经写过了,下面给出球面坐标系下的拉普拉斯方程表述形式:

这个方程也常用如下的简化形式来代替:

或者

其中div指的是向量场(指的是空间中的每一点都有一个对应的带长度的向量)的散度(divergence),grad表示的是标量场的梯度(gradient)。

散度是向量分析中常用的向量算子,用于实现向量场到标量场的转换映射,也就是说,经过散度算子处理后,得到的是一个标量场(每一点有一个不带方向的数值)。以静电场为例,空间中的电场强度是一个向量场,电场线正出负归,在正电荷附近,对应的散度为正值,且电荷带电量越大,散度越大,负电荷附近则反之,其散度为负值,且电荷带电量越大,散度绝对值越大。更为通用的概括是,散度可以看成是向量场在某一点的通量密度,当散度大于0的时候,就表示该点有流量留出,此时这一点可以被称为源点,当散度小于0的时候,表示此点有流量流入,此时此点被称为汇点,散度为0,表示该点无流入也无流出,如果整个向量场的散度都是0,那么这个向量场可以称为无源场。

对于某个向量场 而言,其散度可以通过如下公式求得:

梯度是对多元函数的导数的一种描述,单元函数(只有一个自变量)的导数是标量值函数,而多元(多个自变量)函数的导数则是一个向量值函数,这里多元函数的导数,我们也称为多元函数的梯度,多元函数f在点P处的梯度指的是以f在P处的偏微分作为分量的向量,如一个三维空间函数 ,其梯度函数可以用如下的形式来表述:

单元函数的导数对应的是函数在某一点切线的斜率,对应到梯度上,如果多元函数在某点P的梯度不为0的话,那么计算出来的梯度方向指的是这个函数在P点处增长最快的方向(超平面的切线),而梯度的长度则是函数在此点处的增长率(超平面的斜率)。

举个例子,如果某个房间内的温度用一个函数来表示,那么这个函数在三维空间中的梯度就对应于房间中某点处温度上升最快的方向,而其长度则对应于温度增长率。

可以看到,一个多元函数的标量场,经过梯度转化后,得到的是一个向量场。

从调和函数的定义我们可以看到,所谓的调和函数,实际上就是拉普拉斯方程的解,而我们日常所说的球谐函数(Spherical Harmonics Function)实际上就是拉普拉斯方程在球坐标系空间下的解。

拉普拉斯方程是一个偏微分方程,而解偏微分方程常用的策略是分离变量法,即将偏微分方程分解成几个常微分方程进行求解,下面我们通过将半径跟角度进行分离来进行求解。

设 ,将之代入前面的拉普拉斯方程,可以得到:

上面公式乘上 之后可以得到:

对于上面公式中后面的等式,我们继续使用分离变量法,令(这里是假设Y具有可以分离的形式,当然这个假设不是必然成立的,只是为了简化计算而给出的,只有一些特殊的函数才具有这种假设的可分离的形式) ,代入前面公式可以得到:

简化后,令左右两边均等于 ,可以得到:

一个先验知识是m是一个复数常量(怎么得到的?),且由于 是一个周期函数,其周期可以整除 ,因此m就会是一个整数,而 则是复数指数 的线性组合,Y的常规解出现在极点,也就是 的时候,而在上面的第二个方程中求解 时的常规状态出现在Sturm-Liouville problem的边界点上,在这个边界点中会将 ,其中l是非负整数,且 ,此外,将上面公式中的 用t来替代,就能够得到勒让德公式(Legendre equation),而勒让德公式的解就是伴随勒让德多项式 的倍数。

对于满足前面假设的Y,对于给定的 ,我们总共有 个独立解,这些角度上的解可以表示为三角函数的乘积,这里可以用复数指数与伴随勒让德多项式来表示:

其中这个解需要满足:

上述公式中的 就被称为一个m阶(order)l度(degree)的球谐函数, 就是一个伴随勒让德多项式,N是一个归一化的常量, 则代表着球上的经纬度

所有的球谐函数组成了一组正交基,所谓的正交基指的是,两两基函数相乘的积分只有当两个基函数是同一个基函数的情况下结果为1,否则为0。

上图给出了不同的SH基函数的几何形状展示,这个图是通过以方向为自变量,到球心的距离作为因变量绘制的。

而其他函数都可以通过使用不同系数来对SH基函数进行线性组合来实现近似模拟,这个过程有点像是周期函数的傅里叶展开。

未完待续……

[1] Rendering-球谐光照推导及应用

[2] 调和函数

[3] 拉普拉斯方程

[4] 散度

[5] 梯度

[6] Spherical harmonics

[7]. Laplace's equation

什么是球谐函数

;url=http%3A//physics%2Ebnu%2Eedu%2Ecn/include/elearn/physics/slff/archive22%2Ehtmlb=0a=12user=baidu 看看吧 好运!!


当前文章:球谐函数c语言,球谐函数数学物理方法
标题路径:http://cdweb.net/article/dscedci.html