二维、三维 Laplace 算子的极坐标表示
(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 算子的极坐标表示的更多相关文章
- 使用C语言实现二维,三维绘图算法(1)-透视投影
使用C语言实现二维,三维绘图算法(1)-透视投影 ---- 引言---- 每次使用OpenGL或DirectX写三维程序的时候, 都有一种隔靴搔痒的感觉, 对于内部的三维算法的实现不甚了解. 其实想想 ...
- 使用C语言实现二维,三维绘图算法(3)-简单的二维分形
使用C语言实现二维,三维绘图算法(3)-简单的二维分形 ---- 引言---- 每次使用OpenGL或DirectX写三维程序的时候, 都有一种隔靴搔痒的感觉, 对于内部的三维算法的实现不甚了解. 其 ...
- 使用C语言实现二维,三维绘图算法(2)-解析曲面的显示
使用C语言实现二维,三维绘图算法(2)-解析曲面的显示 ---- 引言---- 每次使用OpenGL或DirectX写三维程序的时候, 都有一种隔靴搔痒的感觉, 对于内部的三维算法的实现不甚了解. 其 ...
- ARCGIS二维三维导航
在使用代码前需要先安装arcgis10.0 或者10.1都可以 不过本人建议初学者安装10.0比较容易安装.. 安装方式和二维三维地图的加载网上都有,就不在此一一赘述了. 先从基本的功能开 ...
- ARCGIS二维三维互动
当对三维模型进行操作时(如导航.平移)二维地图自动跟进. private void Synckron() { m_pGlobe = this._GlobeControl.Globe; m_pMap = ...
- ARCGIS二维三维放大缩小
private void ULZoomPan() { ESRI.ArcGIS.SystemUI.ICommand com = new ControlsGlobeFixedZoomOutCommand( ...
- ARCGIS二维三维平移
private void glZoomPan() { ESRI.ArcGIS.SystemUI.ICommand com = new ControlsGlobePanTool(); com.OnCre ...
- HDU 3404&POJ 3533 Nim积(二维&三维)
(Nim积相关资料来自论文曹钦翔<从"k倍动态减法游戏"出发探究一类组合游戏问题>) 关于Nim积计算的两个函数流程: 代码实现如下: ][]={,,,}; int N ...
- VC、OpenGL、ArcGIS Engine开发的二维三维结合的GIS系统
一.前言 众所周知,二维GIS技术发展了近四十年,伴随着计算机软硬件以及关系型数据库的飞速发展,二维GIS技术已日臻完善.在对地理信息的分析功能上有着无可比拟的优势.一些宏观的地理信息,一维的地理信息 ...
随机推荐
- Windows 版 SourceTree 免登录跳过初始设置的方法
首先去官网下载最新的sourcetree安装包,点击打开下载地址. 下载完成后安装,等到他自启动开始提示你登录的时候,打开“我的电脑(此电脑)”,在最上边的输入栏输入%LocalAppData%\At ...
- solr 7.7 搭建和搜索
最近做了个solr搜索的demo, 用的是solr7.7,之前网上好多文章搭建solr都是5.5+tomcat.. 自己摆弄了下,发现solr7.7和5.5 的部分文件夹有些改动,没有深究原理也就没有 ...
- day14- 面向对象基础(一)
今天开始写关于面向对象的内容,面向对象是编程思想中一大思想,由于日常使用中经常会用到,本次主要是对于我个人认为重点的基础知识进行整理,不会特别详细的书写. 1.面向过程与面向对象的区别 2.类 3.类 ...
- 转://诊断 Grid Infrastructure 启动问题 (文档 ID 1623340.1) .
文档内容 用途 适用范围 详细信息 启动顺序: 集群状态 问题 1: OHASD 无法启动 问题 2: OHASD Agents 未启动 问题 3: OCSSD.BI ...
- SQL学习 DECODE
from 百度百科: DECODE有什么用途呢? 先构造一个例子,假设我们想给这些职员加工资,其标准是:工资在8000元以下的加20%:工资在8000元或以上的加15%,通常的做法是,先选出记录中的工 ...
- 分布式 cephfs
参考链接: http://docs.ceph.com/docs/mimic/cephfs/
- 将List按照指定大小等分的几种实现方式和效率对比及优化
今天碰到一个需求,定时任务,批量从表里取数据并做一些其他操作然后再存表,每次取1000条,由于计算过程比较耗时所以要起多个线程同时跑,需要将List按照指定大小等分,如每100条数据起一个线程,若最后 ...
- 【Swift 4.0】iOS 11 UICollectionView 长按拖拽删除崩溃的问题
正文 功能 用 UICollectionView 实现两个 cell 之间的位置交互或者拖拽某个位置删除 问题 iOS 11 以上拖拽删除会崩溃,在 iOS 9.10 都没有问题 错误 017-10- ...
- python 项目自动生成requirements.txt文件
主要使用目的: 任何应用程序通常需要设置安装所需并依赖一组类库来满足工作要求.通过requirements.txt可以一次性安装程序所需要和依赖的包. 为工程生成requirements.txt的两种 ...
- Sql Server2008如何让外网访问自己的数据库
1.打开 Sql Server配置管理器 2.修改Sql Server网络配置,将Tcp/IP协议修改为“已启用” 3.重启Sql Server服务即可(如果不知道怎么重启Sql Server服务,重 ...