python中scipy.intergrate的部分积分用法

scipy.integrate中提供了几种积分算法,使用help命令来获取该模块的简述
在这里插入图片描述

一般积分(quad)

quad提供了对一个函数一段距离(边界可以为正无穷或负无穷)内积分的功能。如下:对J2.5(x)函数在[0,4.5]范围内进行积分运算。
在这里插入图片描述

import scipy.integrate as integrate
import scipy.special as special
result = integrate.quad(lambda x: special.jv(2.5,x), 0, 4.5)
result=(1.1178179380783249, 7.8663172481899801e-09)
from numpy import sqrt, sin, cos, pi
I = sqrt(2/pi)*(18.0/27*sqrt(2)*cos(4.5) - 4.0/27*sqrt(2)*sin(4.5) +
                sqrt(2*pi) * special.fresnel(3/sqrt(pi))[0])
>>>I=1.117817938088701
print(abs(result[0]-I))
1.03761443881e-11

quad的第一个参数是一个“可调用的”python对象(即函数,方法或类实例)。注意在这种情况下使用lambda函数作为参数。接下来的两个参数是集成的极限。返回值是一个元组,第一个元素保存积分的估计值,第二个元素保存误差的上限。注意,在这种情况下,该积分的真实值为
在这里插入图片描述
Si(x)是菲涅耳正弦积分。
如果被积分函数中有其他参数需要设定,我们可以通过以下方法来处理。例如:
在这里插入图片描述

from scipy.integrate import quad
def integr
  • 3
    点赞
  • 34
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值