python用quad、dblquad实现一维二维积分的实例详解
更新时间：2019年11月20日 09:54:54 作者：潜水的飞鱼baby
今天小编大家分享一篇python用quad、dblquad实现一维二维积分的实例详解，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
背景：
python函数库scipy的quad、dblquad实现一维二维积分的范例。需要注意dblquad的积分顺序问题。
代码：
import numpy as np
from scipy import integrate
def half_circle(x):
"""
原心：(1,0),半径为1
半圆函数：(x-1)^2+y^2 = 1
"""
return (1-(x-1)**2)**0.5
"""
梯形法求积分：半圆线和x轴包围的面积
"""
N = 10000
x = np.linspace(0,2,num=N)#,endpoint=True)
dh = (2-0)/N
y = half_circle(x)
"""
梯形法求积分：（上底+ 下底）*高/2
"""
S = sum((y[1:]+y[:-1])*dh/2)
print("=========%s=========="%"梯形法")
print("面积：%f"%S)
"""
直接调用intergrate的积分函数quad
"""
S2,err = integrate.quad(half_circle,0,2)
print("=========%s=========="%"quad")
print("面积：%f"%S2)
"""
多重定积分:注意积分顺序
"""
def half_sphere(y,x):
"""
球心：（1，0，0）
半径：1
半球：(x-1)^2+y^2+z^2=1
"""
return (1-(x-1)**2-y**2)**0.5
"""
积分顺序：
v = V x in [0,2] :V y in [-g(x),h(x)]
"""
V3,err = integrate.dblquad(half_sphere,0,2,lambda x:-half_circle(x),lambda x:half_circle(x))
print("========%s==========="%"dblquad")
print("体积：%f"%V3)
结果：
========
=========梯形法==========
面积：1.570638
=========quad==========
面积：1.570796
========dblquad===========
体积：2.094395
以上这篇python用quad、dblquad实现一维二维积分的实例详解就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持脚本之家。
相关文章
作为常用的 python 自带数据库管理模块,shelve 和 dbm 都是非常方便的对象持久化存储和检索工具,本文将从用法、优势以及不同点等方面进行介绍,希望对大家有所帮助2023-10-10
这篇文章主要为大家详细介绍了python实现ftp文件传输功能，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2020-03-03
这篇文章主要为大家详细介绍了TensorFlow如何实现反向传播，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2018-02-02
python基于paramiko库远程执行 SSH 命令,实现 sftp 下载文件
这篇文章主要介绍了python基于paramiko库远程执行 SSH 命令,实现 sftp 下载文件的方法，帮助大家更好的理解和学习使用python，感兴趣的朋友可以了解下2021-03-03
这篇文章主要介绍了Pycharm如何对python文件进行打包，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教2023-02-02
这篇文章主要介绍了Python利用代码计算2个坐标之间的距离,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教2023-08-08
Python unittest如何生成HTMLTestRunner模块
这篇文章主要介绍了Python unittest如何生成HTMLTestRunner模块,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下2020-09-09
Python压缩包处理模块zipfile和py7zr操作代码
目前对文件的压缩和解压缩比较常用的格式就是zip格式和7z格式，这篇文章主要介绍了Python压缩包处理模块zipfile和py7zr,需要的朋友可以参考下2022-06-06
这篇文章主要介绍了Python中else的用法总结，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教2022-05-05
这篇文章主要介绍了python批量处理PDF文档，输出自定义关键词的出现次数，文中有详细的代码示例，需要的朋友可以参考阅读2023-04-04
最新评论