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

(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$$

目录
相关文章
|
算法 前端开发
三维形体投影面积
三维形体投影面积
57 0
二维坐标系空间变换(详细解读,附MATLAB代码)
二维坐标系空间变换(详细解读,附MATLAB代码)
959 0
二维坐标系空间变换(详细解读,附MATLAB代码)
|
7月前
|
算法 图形学
【计算机图形学】实验四 二维图形的缩放、旋转,平移,组合变换
【计算机图形学】实验四 二维图形的缩放、旋转,平移,组合变换
333 2
|
7月前
halcon基础系列之二维量测相关算子
halcon基础系列之二维量测相关算子
258 0
|
7月前
[Halcon&几何] 直线的垂线与延长线的计算
[Halcon&几何] 直线的垂线与延长线的计算
377 1
|
数据采集 数据管理 大数据
真正用起来,二维比三维更重要
今天这篇文章起这个标题似乎有点逆潮流,但是这确实是开发的表达,同时也是我这两天和开发交流之后最直接的感受,当然这个开发的群体是有一定的限定的,就是时空大数据的开发团队主管,但是这些经验对于建设其他数字孪生类的应有也有着很好的借鉴和参考意义。(本人工作于交通行业,现在本人单位也在建设数字孪生相关内容,有感而发故此想谈谈感受)
|
Python
点云在任意平面上获取二维投影
点云在任意平面上获取二维投影
1177 0
点云在任意平面上获取二维投影
|
7月前
|
算法
[Halcon&几何] 矩形顶点和对角连线角度计算
[Halcon&几何] 矩形顶点和对角连线角度计算
147 0
|
人工智能 算法 计算机视觉
基于Fringe-Projection环形投影技术的人脸三维形状提取算法matlab仿真
基于Fringe-Projection环形投影技术的人脸三维形状提取算法matlab仿真
|
算法 图形学
【计算机图形学】实验三:二维图形变换
【计算机图形学】实验三:二维图形变换
276 0
【计算机图形学】实验三:二维图形变换