(1) 设 $(r,\theta)$ 是 $\bbR^2$ 的极坐标, 即 $$\bex x=r\cos\theta,\quad y=r\sin \theta. \eex$$ 证明 Laplace 算子 $\dps{\lap=\frac{\p^2}{\p x^2}+ \frac{\p^2}{\p y^2}}$ 可以表示为 $$\bex \lap u=u_{rr}+\frac{1}{r}u_r+\frac{1}{r^2}u_{\theta\theta}. \eex$$ (2) 设 $(r,\theta,\phi)$ 是 $\bbR^3$ 的极坐标, 即 $$\bex x=r\sin \theta\cos \phi,\quad y=r\sin \theta\sin \phi,\quad z=r\cos \theta. \eex$$ 证明 Laplace 算子 $\dps{\lap=\frac{\p^2}{\p x^2}+ \frac{\p^2}{\p y^2}+\frac{\p^2}{\p z^2}}$ 可以表示为 $$\bex \lap u=\frac{1}{r^2}\frac{\p}{\p r}\sex{r^2\frac{\p u}{\p r}} +\frac{1}{r^2\sin \theta} \frac{\p}{\p \theta}\sex{\sin \theta\frac{\p u}{\p \theta}} +\frac{1}{r^2\sin^2\theta}\frac{\p^2u}{\p \phi^2}. \eex$$

证明: (1) 由 $$\bex x_r=\cos \theta,\quad x_\theta=-r\sin \theta=-y,\quad y_r=\sin \theta,\quad y_\theta=r\cos \theta=x \eex$$ 知 $$\beex \bea u_r&=u_x\cos \theta+u_y\sin \theta,\\ u_{rr}&=u_{xx}\cos^2\theta +2u_{xy}\sin \theta\cos \theta +u_{yy}\sin^2\theta,\\ u_\theta&=-yu_x+xu_y,\\ u_{\theta\theta}&= -xu_x-y(-yu_{xx}+xu_{xy})\\ &\quad-yu_y+x(-yu_{xy}+xu_{yy})\\ &=y^2u_{xx}-2xyu_{xy} +x^2u_{yy}-(xu_x+yu_y). \eea \eeex$$ 而 $$\bex u_{rr}+\frac{1}{r}u_r+\frac{1}{r^2}u_{\theta\theta} =u_{xx}+u_{yy}=\lap u. \eex$$ (2) 设 $\rho=r\sin \theta$, 则 $$\bex \ba{ll} x=\rho \cos \phi,&y=\rho\sin \phi,\\ z=r\cos \theta,&\rho=r\sin \theta. \ea \eex$$ 而由 (1), $$\beex \bea u_{xx}+u_{yy}&=u_{\rho\rho} +\frac{1}{\rho}u_\rho+\frac{1}{\rho^2}u_{\phi\phi},\\ u_{zz}+u_{\rho\rho} &=u_{rr}+\frac{1}{r}u_r+\frac{1}{r^2}u_{\theta\theta}. \eea \eeex$$ 两式相加得 $$\bee\label{3.2:1} \lap u=u_{rr}+\frac{1}{r}u_r+\frac{1}{\rho}u_\rho +\frac{1}{r^2}u_{\theta\theta}+\frac{1}{\rho^2}u_{\phi\phi}. \eee$$ 我们再计算 $u_\rho$ 如下: $$\bee\label{3.2:2} \bea u_\rho&=u_rr_\rho+u_\theta\theta_\rho\quad\sex{z=r\cos\theta,\ \rho=r\sin \theta,\quad u(z,\rho)=u(r,\theta)}\\ &=\frac{\rho}{r}u_r+\frac{z}{r^2}u_\theta\\ &=u_r\sin \theta+u_\theta\frac{\cos\theta}{r}. \eea \eee$$ 把 \eqref{3.2:2} 代入 \eqref{3.2:1}, 得 $$\beex \bea \lap u&=u_{rr}+\frac{1}{r}u_r+\frac{1}{r\sin\theta} \sex{u_r\sin \theta+u_\theta\frac{\cos \theta}{r}} +\frac{1}{r^2}u_{\theta\theta} +\frac{1}{r^2\sin^2\theta}u_{\phi\phi}\\ &=u_{rr}+\frac{2}{r}u_r+ \frac{1}{r^2\sin\theta}\cdot u_\theta\cos \theta +\frac{1}{r^2}u_{\theta\theta} +\frac{1}{r^2\sin^2\theta}u_{\phi\phi}\\ &=\frac{1}{r^2}(r^2u_{rr}+2ru_r) +\frac{1}{r^2\sin\theta} (u_\theta\cos \theta+u_{\theta\theta}\sin \theta) +\frac{1}{r^2\sin^2\theta}u_{\phi\phi}\\ &=\frac{1}{r^2}(r^2u_r)_r +\frac{1}{r^2\sin\theta}(u_\theta\sin \theta)_\theta +\frac{1}{r^2\sin^2\theta}u_{\phi\phi}. \eea \eeex$$

二维、三维 Laplace 算子的极坐标表示的更多相关文章

  1. 使用C语言实现二维,三维绘图算法(1)-透视投影

    使用C语言实现二维,三维绘图算法(1)-透视投影 ---- 引言---- 每次使用OpenGL或DirectX写三维程序的时候, 都有一种隔靴搔痒的感觉, 对于内部的三维算法的实现不甚了解. 其实想想 ...

  2. 使用C语言实现二维,三维绘图算法(3)-简单的二维分形

    使用C语言实现二维,三维绘图算法(3)-简单的二维分形 ---- 引言---- 每次使用OpenGL或DirectX写三维程序的时候, 都有一种隔靴搔痒的感觉, 对于内部的三维算法的实现不甚了解. 其 ...

  3. 使用C语言实现二维,三维绘图算法(2)-解析曲面的显示

    使用C语言实现二维,三维绘图算法(2)-解析曲面的显示 ---- 引言---- 每次使用OpenGL或DirectX写三维程序的时候, 都有一种隔靴搔痒的感觉, 对于内部的三维算法的实现不甚了解. 其 ...

  4. ARCGIS二维三维导航

    在使用代码前需要先安装arcgis10.0    或者10.1都可以    不过本人建议初学者安装10.0比较容易安装.. 安装方式和二维三维地图的加载网上都有,就不在此一一赘述了. 先从基本的功能开 ...

  5. ARCGIS二维三维互动

    当对三维模型进行操作时(如导航.平移)二维地图自动跟进. private void Synckron() { m_pGlobe = this._GlobeControl.Globe; m_pMap = ...

  6. ARCGIS二维三维放大缩小

    private void ULZoomPan() { ESRI.ArcGIS.SystemUI.ICommand com = new ControlsGlobeFixedZoomOutCommand( ...

  7. ARCGIS二维三维平移

    private void glZoomPan() { ESRI.ArcGIS.SystemUI.ICommand com = new ControlsGlobePanTool(); com.OnCre ...

  8. HDU 3404&POJ 3533 Nim积(二维&三维)

    (Nim积相关资料来自论文曹钦翔<从"k倍动态减法游戏"出发探究一类组合游戏问题>) 关于Nim积计算的两个函数流程: 代码实现如下: ][]={,,,}; int N ...

  9. VC、OpenGL、ArcGIS Engine开发的二维三维结合的GIS系统

    一.前言 众所周知,二维GIS技术发展了近四十年,伴随着计算机软硬件以及关系型数据库的飞速发展,二维GIS技术已日臻完善.在对地理信息的分析功能上有着无可比拟的优势.一些宏观的地理信息,一维的地理信息 ...

随机推荐

  1. applicationSettings设置和appsttings

    applicationSettings 可以和sttings一样在配置文件中,设置参数.支持定义参数的类型“serializaAs=string”,并可以使用  . 语法. 可以使用.语法

  2. 安装Gradle(Windows & Linux)

    Gradle 是以 Groovy 语言为基础,面向Java应用为主.基于DSL(领域特定语言)语法的自动化构建工具.在github上,gradle项目很多,有的是gradel跟maven构建一块儿使用 ...

  3. 一次CMS GC问题排查过程(理解原理+读懂GC日志)

    这个是之前处理过的一个线上问题,处理过程断断续续,经历了两周多的时间,中间各种尝试,总结如下.这篇文章分三部分: 1.问题的场景和处理过程:2.GC的一些理论东西:3.看懂GC的日志 先说一下问题吧 ...

  4. 实现element-ui中table点击一行展开

    转:https://www.jianshu.com/p/e51ba4cb11d6 先上效果   效果图 三要素 1.row-click 点击行 2.ref 自行了解vue 3.toggleRowExp ...

  5. docker面试整理

    为什么要使用docker  https://www.cnblogs.com/AshOfTime/p/10755479.html docker的使用场景 docker和虚拟机比较的优势   https: ...

  6. hotspot目录结构

    Hotspot的目录结构 ├─agent Serviceability Agent的客户端实现 ├─make 用来build出HotSpot的各种配置文件 ├─src HotSpot VM的源代码 │ ...

  7. C#自定义应用程序上下文对象+IOC自己实现依赖注入

    以前的好多代码都丢失了,加上最近时间空一些,于是想起整理一下以前的个人半拉子项目,试试让它们重生.自从养成了架构师视觉 搭建框架之后,越来 越看不上以前搭的框架了.先撸个上下文对象加上实现依赖注入.由 ...

  8. js中的枚举

    在JavaScript中,对象的属性分为可枚举和不可枚举之分,它们是由属性的enumerable值决定的.可枚举性决定了这个属性能否被for…in查找遍历到. js中基本包装类型的原型属性是不可枚举的 ...

  9. vue.js实战——$event

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  10. js 实现数据结构 -- 栈

    原文: 在 Javascript 中学习数据结构与算法. 概念: 栈是一种遵从先进后出 (LIFO) 原则的有序集合:新添加的或待删除的元素都保存在栈的末尾,称作栈顶,另一端为栈底.在栈里,新元素都靠 ...