当前位置: 首页 > news >正文

球面坐标系下的三重积分

涉及知识点

  • 三重积分
  • 球面坐标系
  • 点火公式
  • 一些常见积分处理手法

球面坐标系定义

球面坐标系由方位角φ\varphiφ、仰角θ\thetaθ和距离rrr构成

直角坐标系(x,y,z)(x,y,z)(x,y,z)到球面坐标系的(r,φ,θ)(r,\varphi,\theta)(r,φ,θ)的转化规则如下:

{x=rsin⁡φcos⁡θy=rsin⁡θsin⁡φz=rcos⁡φ\left\{ \begin{aligned} x & = & r\sin φ\cosθ \\ y & = & r\sin θ\sin φ \\ z & = & r\cos φ \end{aligned} \right. xyz===rsinφcosθrsinθsinφrcosφ

在这里插入图片描述

适用

适用于积分区域为球或球的部分、锥或锥的部分。

处理方法

按规则直角坐标系的积分式转换成球面坐标系就行

∭Ωf(x,y,z)dxdydz=∭Ωf(rsin⁡φcos⁡θ,rsin⁡θsin⁡φ,rcos⁡φ)r2sin⁡φdθdφdr\iiint \limits_{\Omega} f(x,y,z)dxdydz=\iiint \limits_{\Omega}f(r\sin φ\cosθ,r\sin θ\sin φ,r\cos φ)r^2\sin \varphi d\theta d\varphi drΩf(x,y,z)dxdydz=Ωf(rsinφcosθ,rsinθsinφ,rcosφ)r2sinφdθdφdr

然后一般按如下顺序写出积分式:

∫dθ∫dφ∫f(r,θ,φ)dr\int d\theta \int d\varphi \int f(r,\theta,\varphi)drdθdφf(r,θ,φ)dr

由于“后积先定限”,所以先处理方位角,即下图中1的轨迹,随后处理仰角,即下图中2的轨迹,两个角取值范围都是[0,2π][0,2\pi][0,2π]

在这里插入图片描述

例题

计算三重积分∭Ω(x2+y2)dv\iiint \limits_{\Omega}(x^2+y^2)dvΩ(x2+y2)dv其中Ω\OmegaΩ是右半球面x2+y2+z2=a2(y≥0,a>0)x^2+y^2+z^2=a^2\text{ }(y\ge 0,a>0)x2+y2+z2=a2 (y0,a>0)xOzxOzxOz面所围成的区域


【解析】

Ω={0≤r≤a,0≤θ≤π,0≤φ≤π}\Omega = \{0\le r\le a,0\le \theta \le \pi,0\le \varphi \le \pi \}Ω={0ra,0θπ,0φπ}

本题即解如下积分

∭Ωr2sin⁡2φ⋅r2sin⁡φdrdθdφ\iiint \limits_{\Omega}r^2\sin^2\varphi ·r^2\sin \varphi drd\theta d\varphiΩr2sin2φr2sinφdrdθdφ

即:

∫0πdθ∫0πdφ∫0ar4sin⁡3φdr\int_0^\pi d\theta \int_0^\pi d\varphi \int_0^a r^4\sin^3\varphi dr0πdθ0πdφ0ar4sin3φdr

其中在drdrdrsin⁡3φ\sin^3\varphisin3φ是常量,可提出,剩下就是对r4r^4r4积分,即变为:

∫0πdθ∫0πsin⁡3φ⋅a55dφ\int_0^\pi d\theta \int_0^\pi \sin^3\varphi · \frac{a^5}{5} d\varphi0πdθ0πsin3φ5a5dφ

a55\frac{a^5}{5}5a5是常数可提出,并且这个对φ\varphiφ积分完要对θ\thetaθ积分,可以先变换顺序先对θ\thetaθ积分,则原式变为:

π5a3∫0πsin⁡3φdφ\frac{\pi}{5}a^3 \int_0^\pi \sin^3\varphi d\varphi5πa30πsin3φdφ

对于sin⁡3φ\sin^3\varphisin3φ的积分步骤中用到了点火公式,过程如下:

∫0πsin⁡3φdφ=∫0π2sin⁡3φdφ+∫π2πsin⁡3φdφ=23+∫π2πsin⁡3φdφ\int_0^\pi \sin^3\varphi d\varphi=\int_0^{\frac\pi 2} \sin^3\varphi d\varphi+\int_{\frac\pi 2}^\pi \sin^3\varphi d\varphi=\frac23+\int_{\frac\pi 2}^\pi \sin^3\varphi d\varphi0πsin3φdφ=02πsin3φdφ+2ππsin3φdφ=32+2ππsin3φdφ

对于右侧的积分继续进行处理,令φ=π−t\varphi = \pi - tφ=πt(好像算是区间再现公式)

∫π2πsin⁡3φdφ=∫π20sin⁡3(π−t)d(π−t)=∫π20sin⁡3(π−t)d(−t)\int_{\frac\pi 2}^\pi \sin^3\varphi d\varphi=\int_{\frac\pi 2}^0 \sin^3(\pi - t) d(\pi - t)= \int_{\frac\pi 2}^0 \sin^3(\pi - t) d(- t)2ππsin3φdφ=2π0sin3(πt)d(πt)=2π0sin3(πt)d(t)

提出负号,上下限颠倒,则右侧积分式等于:

∫0π2sin⁡3(π−t)dt\int_0^{\frac\pi 2} \sin^3(\pi - t) dt02πsin3(πt)dt

根据sin⁡3x\sin^3xsin3x的对称性,该式子又等于:

∫0π2sin⁡3tdt=23\int_0^{\frac\pi 2} \sin^3t dt=\frac2302πsin3tdt=32

故原式等于

π5a3∫0πsin⁡3φdφ=π5a3⋅(23+23)=415πa5\frac{\pi}{5}a^3 \int_0^\pi \sin^3\varphi d\varphi=\frac{\pi}{5}a^3·(\frac 23+\frac 23)=\frac4{15}\pi a^55πa30πsin3φdφ=5πa3(32+32)=154πa5

即最终结果

http://www.lryc.cn/news/205.html

相关文章:

  • 谷歌 Jason Wei | AI 研究的 4 项基本技能
  • excel数据整理:合并计算快速查看人员变动
  • vit-pytorch实现 MobileViT注意力可视化
  • Python将字典转换为csv
  • EasyX精准帧率控制打气球小游戏
  • 你知道 GO 中什么情况会变量逃逸吗?
  • 一篇文章学懂C++和指针与链表
  • TPGS-cisplatin顺铂修饰维生素E聚乙二醇1000琥珀酸酯
  • 【20230206-0209】哈希表小结
  • c++11 标准模板(STL)(std::multimap)(一)
  • python进阶——自动驾驶寻找车道
  • 男,26岁,做了一年多的自动化测试,最近在纠结要不要转行,求指点。?
  • 源码级别的讲解JAVA 中的CAS
  • JUC锁与AQS技术【我的Android开发技术】
  • 【问题代码】顺序点的深入理解(汇编剖析+手画图解)
  • BinaryAI全新代码匹配模型BAI-2.0上线,“大模型”时代的安全实践
  • nvidia设置wifi和接口
  • PostgreSQL 变化数据捕捉(CDC)
  • Spring 事务【隔离级别与传播机制】
  • HTTP和HTTPS协议
  • day3——有关java运算符的笔记
  • Git多人协同远程开发
  • Chapter4:机器人仿真
  • python(14)--集合
  • 【Spark分布式内存计算框架——Spark Core】4. RDD函数(中)Transformation函数、Action函数
  • Mysql 数据类型
  • 运行Whisper笔记(1)
  • 2023年最强大的12款数据可视化工具,值得收藏
  • LeetCode刷题系列 -- 523. 连续的子数组和
  • LeetCode刷题系列 -- 525. 连续数组