Python 实现定积分与二重定积分的操作

更新时间：2021年05月26日 09:57:08 作者：qq_33203870

这篇文章主要介绍了Python 实现定积分与二重定积分的操作，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教

1.概述

最近项目需要使用程序实现数学微积分，最初想用java实现，后来发现可用文档太少，实现比较麻烦，后来尝试使用python实现，代码量较少，主要有sympy与scipy两种实现方式，本文主要记录scipy的实现方式。

2.内容

2.1 所求函数

在这里插入图片描述

2.2 python代码

# 引入需要的包
import scipy.integrate
from numpy import exp
from math import sqrt
import math

# 创建表达式
f = lambda x,y : exp(x**2-y**2)

# 计算二重积分：（p:积分值，err:误差）
# 这里注意积分区间的顺序
# 第二重积分的区间参数要以函数的形式传入
p,err= scipy.integrate.dblquad(f, 0, 2, lambda g : 0, lambda h : 1)	
print(p)

2.3 注意问题

1. exp尽量使用numpy的exp

2. 注意积分区间参数的顺序

3. 第二重积分的区间参数要以函数的形式传入

补充：python实现求解积分

例子 1：

假设有随机变量 x，定义域 X，其概率密度函数为 p(x)，f(x) 为定义在 X 上的函数，目标是求函数 f(x) 关于密度函数 p(x) 的数学期望。

蒙特卡洛法根据概率分布 p(x) 独立地抽样 n 个样本 x1,x2,…..xn，得到近似的 f(x) 期望为：

其实这个的理解就是要求一个拥有概率密度的函数期望值

期望=积分（每个点的密度函数*每个点的价值函数）

例子 2：

假设我们想要求解 h(x) 在 X 上的积分：

我们将 h(x) 分解成一个函数 f(x) 和一个概率密度函数 p(x) 的乘积，进而又将问题转换为求解函数 f(x) 关于密度函数 p(x) 的数学期望：

这里的Ep(x)是相当于把整个分布当时了概率分布，即总发生概率为1.

这里，f(x) 表示为，则有：

更一般的，假设我们想要求解，熟悉积分的同学肯定已经知道答案为，那么如何用采样的方法来得到这个值呢？

令，0<x<10，那么。

下面是代码：

'''import random
num=1000000
sum=0
for i in range(0,num):
    x=random.uniform(0,10)
    sum+=x*x*10
sum/=1000000
print(sum)'''
import random
numSamples=10000
samples=[random.uniform(0,10)for _ in range(numSamples)]
f_samples=[10*sample**2 for sample in samples]
result=1/10000.0*sum(f_samples)
print(result)

result=333.10527012455066

random.uniform(x,y)表示在[x,y)之间生成一个实数

对于复杂的 h(x)，这种方法计算起来显然就更加方便了（特别是忘记积分怎么算的同学）。

蒙特卡洛方法其实就是利用大数定理通过大量统计来算出最后的值。

到这里为止，我们简单的介绍了蒙特卡洛方法，但是依旧没有提到要怎么利用复杂的概率密度函数进行采样。

接下来我们来看一下接受-拒绝法（accept-reject sampling method），它也是蒙特卡洛法中的一种类型适用于不能直接抽样的情况。

以上为个人经验，希望能给大家一个参考，也希望大家多多支持脚本之家。

您可能感兴趣的文章:

Python使用定时调度任务的方式
Python 有几种方法可以定时调度一个任务，这就是我们将在本文中学习的内容。接下俩下边将给大家介绍5种Python使用定时调度任务方式,需要的朋友可以参考下面文章内容吧
2021-09-09
python中单双下划线的区别对比分析
Python中存在一些特殊的方法，有些方法以双下划线 “__” 开头和结尾，它们是Python的魔法函数，比如__init__()和__str__等等，不用要这种方式命名自己的变量或者函数，这篇文章主要介绍了python里面单双下划线的区别,需要的朋友可以参考下
2023-04-04
python基础入门之普通操作与函数（三）
这篇文章主要介绍了python基础入门之普通操作与函数
2021-06-06
详解Python是如何处理不同时区的
时区是指在地球上不同地方的时间差异，地球分为 24 个时区，每个时区都相对于格林威治标准时间或协调世界时（UTC）有所偏移。本文主要和大家来聊聊Python是如何处理不同时区的，希望对大家有所帮助
2023-02-02
python求绝对值的三种方法小结
今天小编就为大家分享一篇python求绝对值的三种方法小结，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2019-12-12
Pytorch可视化之Visdom使用实例
Visdom 是一个专门用于 PyTorch 的交互式可视化工具,可以对实时数据进行丰富的可视化,帮助我们实时监控在远程服务器上进行的科学实验,这篇文章主要给大家介绍了关于Pytorch可视化之Visdom使用的相关资料,需要的朋友可以参考下
2021-08-08
使用keras时input_shape的维度表示问题说明
这篇文章主要介绍了使用keras时input_shape的维度表示问题说明，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2020-06-06
配置python连接oracle读取excel数据写入数据库的操作流程
这篇文章主要介绍了配置python连接oracle，读取excel数据写入数据库,本文通过实例代码给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
2021-03-03
django框架CSRF防护原理与用法分析
这篇文章主要介绍了django框架CSRF防护原理与用法,结合实例形式分析了Django框架CSRF防护的概念、原理、使用方法及相关操作注意事项,需要的朋友可以参考下
2019-07-07
Hadoop中的Python框架的使用指南
这篇文章主要介绍了Hadoop中的Python框架的使用指南,Hadoop一般使用复杂的Java操作,但通过该框架使得Python脚本操作Hadoop成为了可能,需要的朋友可以参考下
2015-04-04