二维、三维 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$$

时间: 2024-10-01 10:31:54

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

改写函数实现PHP二维/三维数组转字符串_php技巧

由于工作需要,自己在手册给定的示例函数基础上改写出了这样一个函数,代码如下: 复制代码 代码如下: //将多维数组中所有的数值转换成字符串---->最多支持三维数组 function implodex( $glue, $array, $separator='' ) { if ( ! is_array( $array ) ) return $array; $string = array(); $count = 0; foreach ( $array as $key => $val ) { if

海量二维图纸 与三维模型齐头并进

&http://www.aliyun.com/zixun/aggregation/37954.html">nbsp;   [IT168 资讯]当下三维时代正以迅猛态势前进着,与此同时,是否就代表着二维时代的结束呢? 其实不然,日前,3DSource云应用中心平台重磅推出特色资源:二维图纸http://apps.3dsource.cn/manufacturing/drawing.html.在已有海量的三维模型资源的基础上,进一步拓展出海量的二维图纸工程资源.据悉,二维图纸能在一定程度

二维CAD市场访谈之UGS篇

前言:e-workshttp://www.aliyun.com/zixun/aggregation/32268.html">调查显示,近两年二维CAD市场的营业额有了较大增长,为了对这种市场现象进行深入的分析,e-works正在策划<二维CAD是复兴还是回光返照?>专题.UGS以问卷的形式发表了自己对二维CAD市场的看法,以下是问答内容. 1.调查显示,二维CAD市场近两年有较大增长,而UGS在这期间反而将Solid Edge 2D Drafting免费,这是基于什么原因? U

(七十八)嵌套和二维数组

例如int a[5]是一维数组,他是由a[0].a[1].a[2].a[3].a[4]这五个int类型的数组组成的.   如果放到excel里面,那么就是: a[0] a[1] a[2] a[3] a[4]   这样一个样子,一维表示线. 一维数组的表示方法可以理解为: 变量名[横坐标x] x初始为0,表示原点   假如表格从线到面,就升格为二维. 二维数组的的特点是既有横坐标,又有纵坐标.声明格式为: int 变量名[纵坐标y][横坐标x] 表示每行有x个元素,每列有y行. 例如int a[3

怎么使用vtk将二维图片合成三维图片

问题描述 怎么使用vtk将二维图片合成三维图片 30C 我的毕业设计是做图片处理的 ,现在图片处理完了,老师说要将二维的图片合成三维的(之前的二维图片也是用matlab将三维图片分割来的)vtk工具没有用过,网上搜的教程也看不太懂 求会vtk的兄弟姐妹帮忙用vc代码写一下 顺便说说vc里面需要怎么配置 不胜感激!

用AutoCAD将三维弯管模型转成二维图纸

以前我们介绍了<AutoCAD转二维工程图为三维弯管模型>,今天我们介绍反方向操作:用AutoCAD将三维的弯管模型转成二维图纸.

ZBrush二维图片转变三维浮雕的方法

  制作原理 ZBrush软件有个很不错的工具,叫Alpha通道.它是用来控制笔刷的形状的.更不错的是,它还能以图片作为Alpha通道.其实,它是将图片的颜色信息作为强度的控制,黑色表示笔刷强度为0%,白色表示笔刷强度为100%.这样的处理方式正好吻合了浮雕灰度图1制作的本意,故利用此工具可以直接将二维浮雕灰度图转化为三维浮雕的效果. 操作演示 这里我们将霸气侧漏的牛浮雕贴在一个3D扫描人像的身上. 步骤一 导入模型:在右侧的托盘中点击Tool>Import[工具>导入],选择模型,在画布中拖

solidworks三维图怎么转二维图?

  solidworks三维图怎么转二维图?            1.打开需要转换的SolidWorks三维图形 2.点击标题栏"新建"图标,选择"从零件/装配体制作工程图" 3.选择模板点击确认 4.进入"工程图"界面,选择右侧的"等轴测",再选择左侧确认. 5.点击标题栏"文件(F)"下拉菜单,选择"另存为(A)" 6.在"另存为"对话框的"保存类型(

cad-纬地6.8怎么把二维地形图转三维地形图

问题描述 纬地6.8怎么把二维地形图转三维地形图 纬地6.8怎么把二维地形图转三维地形图,求大神!具体步骤是什么? 解决方案 等高地形图 解决方案二: 本质上是没办法的,没有高程数据 或许找到该地区的等高线图,可以进行转换 经纬度,加上高程 才是个完整的坐标 才好转换成三维的地心坐标,不过要考虑旋转椭球面的问题,地球表面不是正圆球体,