python用quad、dblquad實(shí)現(xiàn)一維二維積分的實(shí)例詳解
更新時(shí)間：2019年11月20日 09:54:54 作者：潛水的飛魚(yú)baby
今天小編大家分享一篇python用quad、dblquad實(shí)現(xiàn)一維二維積分的實(shí)例詳解，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧
背景：
python函數(shù)庫(kù)scipy的quad、dblquad實(shí)現(xiàn)一維二維積分的范例。需要注意dblquad的積分順序問(wèn)題。
代碼：
import numpy as np
from scipy import integrate
def half_circle(x):
"""
原心：(1,0),半徑為1
半圓函數(shù)：(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)
"""
直接調(diào)用intergrate的積分函數(shù)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)
結(jié)果：
========
=========梯形法==========
面積：1.570638
=========quad==========
面積：1.570796
========dblquad===========
體積：2.094395
以上這篇python用quad、dblquad實(shí)現(xiàn)一維二維積分的實(shí)例詳解就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。
相關(guān)文章
本文主要介紹了Pandas中inf值替換的方法，文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值，需要的朋友們下面隨著小編來(lái)一起學(xué)習(xí)學(xué)習(xí)吧2022-08-08
PyCharm搭建一勞永逸的開(kāi)發(fā)環(huán)境
這篇文章主要介紹了PyCharm搭建一勞永逸的開(kāi)發(fā)環(huán)境,本文通過(guò)圖文并茂的形式給大家介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值，需要的朋友可以參考下2021-04-04
關(guān)于文件Permission denied解決方案(pip)
這篇文章主要介紹了文件Permission denied解決方案(pip),具有很好的參考價(jià)值,希望對(duì)大家有所幫助,如有錯(cuò)誤或未考慮完全的地方,望不吝賜教2023-08-08
這篇文章主要介紹了Python tkinter和exe打包的方法，本文給大家介紹的非常詳細(xì)，具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2020-02-02
Python的collections模塊中namedtuple結(jié)構(gòu)使用示例
namedtuple顧名思義,就是名字+元組的數(shù)據(jù)結(jié)構(gòu),下面就來(lái)看一下Python的collections模塊中namedtuple結(jié)構(gòu)使用示例2016-07-07
Python實(shí)現(xiàn)加解密,編碼解碼和進(jìn)制轉(zhuǎn)換(最全版)
這篇文章主要為大家詳細(xì)介紹了Python實(shí)現(xiàn)加解密、編碼解碼、進(jìn)制轉(zhuǎn)換、字符串轉(zhuǎn)換的最全版操作方法，文中的示例代碼講解詳細(xì)，大家可以收藏一下2023-01-01
python geopandas讀取、創(chuàng)建shapefile文件的方法
shapefile是GIS中非常重要的一種數(shù)據(jù)類(lèi)型，在ArcGIS中被稱為要素類(lèi)(Feature Class)，主要包括點(diǎn)(point)、線(polyline)和多邊形(polygon)，本文重點(diǎn)給大家介紹python geopandas讀取、創(chuàng)建shapefile文件的方法，需要的朋友參考下吧2021-06-06
這篇文章主要介紹了python如何提取xml指定內(nèi)容，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。如有錯(cuò)誤或未考慮完全的地方，望不吝賜教2023-01-01
Python實(shí)現(xiàn)把類(lèi)當(dāng)做字典來(lái)訪問(wèn)
今天小編就為大家分享一篇Python實(shí)現(xiàn)把類(lèi)當(dāng)做字典來(lái)訪問(wèn)，具有很好的參考價(jià)值，希望對(duì)大家有所幫助。一起跟隨小編過(guò)來(lái)看看吧2019-12-12
這篇文章主要介紹了Python中列表(List)的詳解操作方法,包含創(chuàng)建、訪問(wèn)、更新、刪除、其它操作等,需要的朋友可以參考下2021-10-10