SciPy 积分的计算

2024-05-15 19:31:02

Scipy中的integrate模块提供了很多数值积分方法，例如，一重积分、二重积分、三重积分、多重积分、高斯积分等等。

一重积分

SciPy积分模块中，quad函数是一个重要函数，用于求一重积分。例如，在给定的a到b范围内，对函数f(x)求一重积分。

$\int_{a}^{b}f(x)dx$

quad的一般形式是scipy.integrate.quad(f, a, b)，其中f是求积分的函数名称，a和b分别是下限和上限。

求0到5范围内的积分,首先需要定义函数→ $f(x)=e^{-x^{2}}$

import scipy.integrate
from numpy import exp
f = lambda x:exp(-x**2)
i = scipy.integrate.quad(f, 0, 5)
print(i)

结果：

(0.8862269254513955, 2.3183115159980698e-14)

注意：quad函数返回两个值，第一个值是积分的值，第二个值是对积分值的绝对误差估计。

积分的函数f带系数参数：

$I(a,b)=\int_{0}^{1}(ax^{2}+b)dx$

那么a和b可以通过args传入quad函数：

from scipy.integrate import quaddef f(x, a, b):return a * (x ** 2) + bret = quad(f, 0, 1, args=(3, 1))
print (ret)

结果：

(2.0, 2.220446049250313e-14)

重积分

要计算二重积分、三重积分、多重积分，可使用dblquad、tplquad和nquad函数。

二重积分

dblquad的一般形式是scipy.integrate.dblquad(func, a, b, gfun, hfun)，其中，func是待积分函数的名称，a、b是x变量的上下限，gfun、hfun为定义y变量上下限的函数名称。

求二重积分：

$\int_{0}^{2}dy\int_{0}^{\sqrt{1-4y^{2}}}19xydx$

使用lambda表达式定义函数f、g和h。注意，在很多情况下g和h可能是常数，但是即使g和h是常数，也必须被定义为函数。

import scipy.integrate
from numpy import exp
from math import sqrt
f = lambda x, y : 19*x*y
g = lambda x : 0
h = lambda y : sqrt(1-4*y**2)
i = scipy.integrate.dblquad(f, 0, 0.5, g, h)
print (i)

结果：

(0.59375, 2.029716563995638e-14)

Scipy的integrate模块中还有许多其他积分方法，例如nquad，用于求多重积分。不过，大部分场景quad和dblquad就够用了。

SciPy 积分的计算相关推荐

python 二重积分_Python机器学习（五十七）SciPy 积分
Scipy中的integrate模块提供了很多数值积分方法,例如,一重积分.二重积分.三重积分.多重积分.高斯积分等等. 下面介绍几种常用的积分函数. 一重积分 SciPy积分模块中,quad函数是一 ...
python应用-scipy,numpy,sympy计算微积分
python应用-scipy,numpy,sympy计算微积分今天来讲一下使用python进行微积分运算,python有很多科学计算库都可以进行微积分运算,当然如果知晓微积分计算的原理也可以自己编程 ...
SciPy 积分 | Python技能树征题
SciPy 积分 | Python技能树征题 0. 前言 1. 第 1 题:具有函数表达式的被积函数求积分 2. 第 2 题:函数表达式未知的积分求解试题代码地址 0. 前言积分在科学和工程应用中 ...
运用数学软件matlab求无穷积分,matlab积分的计算及其简单应用论文.doc
积分的计算及其简单应用摘要:本文简要的概述了MATLAB 在高等数学中积分的计算及应用:利用MATLAB 中符号积分和数值积分的命令,计算定积分和不定积分.同时,也可以通过这些命令来解决一些实际问题 ...
曲面积分的投影法_在家学|第一类曲面积分与第二类曲面积分的计算
利用投影法计算第一类曲面积分设函数为定义在曲面上的连续函数.曲面的方程为.具有对和的连续偏导数,即此曲面是光滑的,且其在平面上的投影为可求面积的.则如果曲面由方程给出,在坐标面上的投影区域为,函数 ...
对圆柱面的曲面积分_积分曲面为圆柱面的曲面积分的计算
数理医药学杂志 2013年第 26卷第2期文章编号 :1004-4337(2013)02-0133-02 中图分类号:0172．2 文献标识码:C 积分曲面为圆柱面的曲面积分的计算贺勇 (武汉 ...
Simpson积分方法计算NURBS曲线弧长，详细原理+代码实现
Simpson积分方法计算NURBS曲线弧长,详细原理+代码实现 Simpson 积分方法是一种数值积分方法,可以用于计算曲线的弧长.它的基本思想是将曲线分成若干小段,对每一小段采用 Simpson ...
高等数学复习笔记（五）-一元函数积分的计算
本节为高等数学复习笔记的第五部分,一元函数积分的计算,主要包括: 四大积分基本法即凑微分法.换元法(第一类和第二类).分部积分法以及有理函数积分法. 1. 凑微分法 1.1 基本导数公式(反过来用) ...
Scipy中积分的计算
一重积分的运算形如f(x)={k(x**n)dx,首先积分运算需先引用scipy.intergrate,一重积分引入quad(), 用法:quad(func,a,b,args(n,k)),func是 ...

最新文章

热门文章