python用quad、dblquad实现一维二维积分范例
2017-09-25 13:43
387 查看
背景:
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实现二维曲线拟合
- python实现数值积分的Simpson方法实例分析
- python实现数值积分的Simpson方法实例分析
- PHP实现判断数组是一维、二维或几维的方法
- Python实现App自动签到领取积分功能
- 【python】一维二维插值
- 实现动态分配一维,二维,三维数组
- python实现二维函数高次拟合
- python3实现一维数组求最大值返回参数需要注意的地方
- Python实现二维有序数组查找的方法
- EM算法逼近GMM参数针对二维数据点的python实现
- Python契约式设计范例 - AOP风格的实现
- 【python】双向二维PCA(2D-2D PCA)算法实现
- python OpenCV学习笔记实现二维直方图
- python二维列表一维列表的互相转换
- Python 实现「食行生鲜」签到领积分
- Python的“二维”字典 (two-dimension dictionary)定义与实现方法
- python利用elaphe制作二维条形码实现代码
- 0-1背包的二维实现及一维实现
- python 将一维列表转二维列表