python用quad、dblquad实现一维二维积分的实例详解
作者：潜水的飞鱼baby 发布时间：2022-02-17 05:32:51
标签：python,一维,二维,积分
背景：
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
来源：https://blog.csdn.net/u011702002/article/details/78084245
0
投稿猜你喜欢
- check.asp 代码如下:<％ '''''''''
本文实例讲述了Python实现读取及写入csv文件的方法。分享给大家供大家参考，具体如下：新建csvData.csv文件，数据如下：具体代码- 1、吃金币源码分享：import osimport cfgimport sysimport pygameimport randomfrom
- 一、制作思路由于注册的时候常常会用到注册码来防止机器恶意注册，这里我发表一个产生png图片验证码的基本图像，简单的思路分析：1、产生一张pn
- Data URIData URI是由RFC 2397定义的一种把小文件直接嵌入文档的方案。通过如下语法就可以把小文件变成指定编码直接嵌入到页
- 在my.cnf中设置了sql_mode='STRICT_TRANS_TABLES';却发现还是有非法数据的写入，比如无符号整
- 从Python3.2引入的concurrent.futures模块，Python2.5以上需要在pypi中安装futures包。future
- Pytorch 使用GPU训练使用 GPU 训练只需要在原来的代码中修改几处就可以了。我们有两种方式实现代码在 GPU 上进行训练方法一 .
- 变量插入字符串的方法Python中的format()函数是一种将变量插入字符串的方法，能够使字符串更易于阅读和理解。它支持许多不同的用法，以
- 本文中的示例主要是解决在函数间不能传递多个(32个以上)参数的问题，解题的具体思路就是采用记录类型作为函数的输入和返回值，所以我们需要先定义
- 一、介绍ADB（Android调试桥）是一个命令行工具（CLI），可用于控制Android设备并与之通信。您可以执行许多操作，例如安装应用程
- 仿豆瓣分页原型(Javascript版)写了个分页的样式。自我感觉，这样的分页前后兼顾，对于用户的体验是蛮好使的Javascript分页代码
- 本文实例总结了php随机取mysql记录方法。分享给大家供大家参考。具体分析如下：在php中要随机取mysql记录我们可以直接使用mysql
- 这篇文章主要介绍了python Opencv计算图像相似度过程解析,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价
- 在做我的友情链接批量检查工具过程中，碰到一些情况，就是对方网页会用gzip压缩。用gzip压缩的好处是，能压缩网页大小，加快网页的浏览速度，
- 本文实例讲述了Python基于正则表达式实现检查文件内容的方法分享给大家供大家参考，具体如下：这个是之前就在学python，欣赏python
- 本文实例为大家分享了python实现打砖块小游戏的具体代码，供大家参考，具体内容如下开发益智的打砖块小游戏，你可以试一下能打几块import
- 通常我们做网站，设计版面布局是第一步，如何做到版面布局具有创意又美观大方呢?这就需要一定的版面处理功底。让我们先来了解一下版面布局的步骤:一
- 在Python3中已经有很大一部分语句与Python2不互通了，运行暂停的方法也有所不同。1、input();这种方法不用包含模块，因此这也
- 当服务器必须提供与两个或更多个网络或网络子网的连接时，典型的方案是使用多宿主计算机。此计算机通常位于外围网络（也称为 DMZ、外围安全区域或
