时间:2021-05-22
背景:
python函数库scipy的quad、dblquad实现一维二维积分的范例。需要注意dblquad的积分顺序问题。
代码:
import numpy as npfrom 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 = 10000x = np.linspace(0,2,num=N)#,endpoint=True)dh = (2-0)/Ny = 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实现一维二维积分的实例详解就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
声明:本页内容来源网络,仅供用户参考;我单位不保证亦不表示资料全面及准确无误,也不保证亦不表示这些资料为最新信息,如因任何原因,本网内容或者用户因倚赖本网内容造成任何损失或损害,我单位将不会负任何法律责任。如涉及版权问题,请提交至online#300.cn邮箱联系删除。
python生成二维码的实例详解版本相关操作系统:MacOSXEICaptionPython版本:2.7IDE:SublimeText3依赖库Python生成二
本文实例讲述了Python实现将文本生成二维码的方法。分享给大家供大家参考,具体如下:#coding:utf-8'''Python生成二维码v1.0主要将文本生
python调用二维列表中的一维列表的方法:访问二维列表中的一维列表可以用下标法“列表名(数字)”的方式获取到一维列表所有元素x=[[23,25,15,69,4
本文实例讲述了python实现把二维列表变为一维列表的方法。分享给大家供大家参考,具体如下:c=[[1,2,3],[4,5,6],[7,8,9]]1.用列表推导
本文实例讲述了C#实现的二维数组排序算法。分享给大家供大家参考,具体如下:classOrder{//////对二维数组排序//////排序的二维数组///排序根