(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. gulp 自动ftp至服务器时,处理开发 测试服务器地址问题

    var gulp=require('gulp'), babel = require('gulp-babel'), gulpSequence = require('gulp-sequence'), ht ...

  2. 【转】localStorage使用总结

    原文地址:https://www.cnblogs.com/st-leslie/p/5617130.html 一.什么是localStorage.sessionStorage 在HTML5中,新加入了一 ...

  3. (hdu)4858 项目管理 (vector)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4858 Problem Description 我们建造了一个大项目!这个项目有n个节点,用很多边连接起 ...

  4. PS制作漂亮紫色霓虹灯光文字

    一.新建画布,大小1500 * 950像素,分辨率为300,置入墙壁图像,大小适合. 二.调整图层的色阶,色相/饱和度. 三.新建文字图层,颜色为#a33e88,大小为103,字体为Beon Medi ...

  5. 证明与计算(4): 完美散列函数(Perfect Hash function)

    原文:wiki: 完美散列函数 假设,写一个SQL语句解析器,词法分析对SQL语句解析,把语句分成了多个token,一般这个时候会需要查询这个token是否是一个关键字token. 例如keyword ...

  6. Django rest framework 使用haystack对接Elasticsearch

    Elasticsearch 介绍 ElasticSearch是一个基于Lucene的搜索服务器.它提供了一个分布式多用户能力的全文搜索引擎,基于RESTful web接口.Elasticsearch是 ...

  7. Ubuntu 系统安装详解 19.04最新版本

    Ubuntu 19.04版本系统安装详解 1 .镜像的下载 推荐 阿里云镜像下载 2.安装 1.1.新建虚拟机 注意硬件的兼容性问题 当前只有5.x可以用,其他兼容各位可以尝试下,我也都试过,但只有5 ...

  8. Spring-framework

    1.spring注解驱动开发 官方文档 @Configuration 告诉spring这是一个配置类,配置类=配置文件 @Bean 给容器中注入一个bean,类型为返回值类型,id默认用方法名作为id ...

  9. Activiti6作业执行器Job Executor配置(学习笔记)

    内容概况: 异步执行配置相关: asyncExecutorActivate:这个属性是激活作业执行器,它的默认参数是false,只有设为true,activiti启动的时候才会开启线程池去扫描定时操作 ...

  10. 关于mysql主从架构master宕机后,请求转移问题解决办法

    mysql架构:一主一从 问题一:有两台mysql数据库,已做好主从.如果运行某一天master服务器mysql故障导致前端请求无法处理怎么办? 答:将前端需要数据库处理的请求转移到slave机上. ...