python中quad_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中quad_python用quad、dblquad实现一维二维积分的实例详解相关推荐

python生成二维码_python生成二维码的实例详解
python生成二维码的实例详解版本相关操作系统:Mac OS X EI Caption Python版本:2.7 IDE:Sublime Text 3 依赖库 Python生成二维码需要的依赖库 ...
python中代理模式分为几种_Python设计模式之代理模式实例详解
本文实例讲述了Python设计模式之代理模式.分享给大家供大家参考,具体如下: 代理模式(Proxy Pattern):为其他对象提供一种代理以控制对这个对象的访问 #!/usr/bin/env py ...
java前端长连接框架_Java中Spring Boot+Socket实现与html页面的长连接实例详解
Spring Boot+Socket实现与html页面的长连接,客户端给服务器端发消息,服务器给客户端轮询发送消息,附案例源码功能介绍客户端给所有在线用户发送消息客户端给指定在线用户发送消息服务器 ...
python从date目录导入数据集_PyTorch加载自己的数据集实例详解
数据预处理在解决深度学习问题的过程中,往往需要花费大量的时间和精力. 数据处理的质量对训练神经网络来说十分重要,良好的数据处理不仅会加速模型训练, 更会提高模型性能.为解决这一问题,PyTorch提供 ...
python 录制web视频_Python django框架 web端视频加密的实例详解
视频加密流程图: 后端获取保利威的视频播放授权token,提供接口api给前端参考文档:http://dev.polyv.net/2019/videoproduct/v-api/v-api-play ...
python 经典脚本文件_Python3.5文件读与写操作经典实例详解
本文实例讲述了Python3.5文件读与写操作.分享给大家供大家参考,具体如下: 1.文件操作的基本流程: (1)打开文件,得到文件句柄并赋值给一个变量 (2)通过句柄对文件进行操作 (3)关闭文件 ...
python如何调用文件进行换位加密_python 换位密码算法的实例详解
python 换位密码算法的实例详解一前言: 换位密码基本原理:先把明文按照固定长度进行分组,然后对每一组的字符进行换位操作,从而实现加密.例如,字符串"Error should neve ...
python二维元素向量_详解python Numpy中求向量和矩阵的范数
在python Numpy中求向量和矩阵的范数实例 np.linalg.norm(求范数):linalg=linear(线性)+algebra(代数),norm则表示范数. 函数参数 x_norm=n ...
在python中使用关键字define定义函数_python自定义函数def的应用详解
这里是三岁,来和大家唠唠自定义函数,这一个神奇的东西,带大家白话玩转自定义函数自定义函数,编程里面的精髓! def 自定义函数的必要函数:def 使用方法:def 函数名(参数1,参数2,参数-): ...

python中quad_python用quad、dblquad实现一维二维积分的实例详解

python中quad_python用quad、dblquad实现一维二维积分的实例详解相关推荐

最新文章

热门文章