黑马程序员技术交流社区

标题: 【上海校区】scipy计算三重积分 [打印本页]

---

作者: 小影姐姐    时间: 2018-7-20 13:51
标题: 【上海校区】scipy计算三重积分

圆柱体直径4cm，长20cm，两端各有一个球，每个球直径20cm，密度为5600kg/m3计算rotational inertia matrix .坐标原点为圆柱的中心，沿长度的方向为z轴。

先在草稿纸上写出直角坐标的积分公式，将其写入程序即可。计算三重积分用到的函数是tplquad，此外还有二重积分，定积分，n重积分。f为被积函数，cyl为求圆柱的积分，sphereup求上球的积分，spheredown求下球的积分，返回一个元组，第一个元素是积分值，第二个元素是精度。

• 
  
  from scipy import integrate
  
  
• 
  
  f = lambda x, y, z: x**2+y**2
  
  
• 
  
  
  
  
• 
  
  cyl=integrate.tplquad(f, -0.1, 0.1,
  
  
• 
  
               lambda y: -0.02, lambda y: 0.02,
  
  
• 
  
       lambda x, y: -np.sqrt(0.0004-y**2), lambda x, y: np.sqrt(0.0004-y**2))
  
  
• 
  
  sphereup=integrate.tplquad(f, 0.1, 0.3,
  
  
• 
  
           lambda z: -np.sqrt(0.40*z - z**2 - 0.03),
  
  
• 
  
           lambda z: np.sqrt(0.40*z - z**2 - 0.03),
  
  
• 
  
       lambda z,y: -np.sqrt(0.40*z - z**2 - 0.03 - y**2),
  
  
• 
  
       lambda z,y: np.sqrt(0.40*z - z**2 - 0.03 - y**2))
  
  
• 
  
  spheredown=integrate.tplquad(f, -0.3, -0.1,
  
  
• 
  
           lambda z: -np.sqrt(-0.40*z - z**2 - 0.03),
  
  
• 
  
           lambda z: np.sqrt(-0.40*z - z**2 - 0.03),
  
  
• 
  
       lambda z,y: -np.sqrt(-0.40*z - z**2 -0.03 - y**2),
  
  
• 
  
       lambda z,y: np.sqrt(-0.40*z - z**2 - 0.03 - y**2))
  
  
• 
  
  
  
  
• 
  
  print((cyl[0]+sphereup[0]+spheredown[0])*5600)
  
  
  

---

作者: 不二晨    时间: 2018-7-20 17:45
优秀，奈斯

---

作者: wuqiong    时间: 2018-7-23 11:34

---

作者: 摩西摩西OvO    时间: 2018-7-23 13:44

---

作者: 不二晨    时间: 2018-7-23 13:52
奈斯，优秀

---

作者: 摩西摩西OvO    时间: 2018-7-26 09:40

---

作者: 吴琼老师    时间: 2018-7-26 16:05

---

欢迎光临 黑马程序员技术交流社区 (http://bbs.itheima.com/)

黑马程序员IT技术论坛 X3.2