1.符号积分:#通过integrate()功能(facility),SymPy对基本和特殊函数定与不定积分有卓越的支持。

#该功能使用有力的扩展Risch-Norman算法,启发算法和模式匹配:

from sympy import integrate, exp, sin, log, oo, pi,symbols

x, y = symbols('x,y')#定义符号变量x,y

#求不定积分

integrate(6*x**5, x)

integrate(log(x), x)

#求定积分

integrate(x**3, (x, -1, 1))

#反常积分

integrate(exp(-x), (x, 0, oo))

#运行结果:x**6,x*log(x) - x,0,1

2.数值积分:

#符号积分与数值积分各有优劣,下面介绍数值积分

from scipy import integrate#导入数值积分模块

import numpy as np

#读者可以dir(integrate)来查看里面有哪些数值积分函数并help他们的用法

#quad命令(该命令使用了Fortran library QUADPACK)

f=lambda x:x**3#注意这里匿名函数的用法,不懂请百度

integrate.quad(f,-1,1) #(0.0, 5.551208455924673e-15)后面那个超小的是浮点运算产生的误差

#fixed_quad命令(使用fixed-order Gaussian quadrature)

f=lambda x:np.exp(-x)

integrate.fixed_quad(f,0,1) #(0.63212055882831653, None)

#trapz(复化梯形求积分)

f=lambda x:np.exp(-x)

x=np.linspace(0,1,10)

integrate.trapz(f(x),x) #0.63277075484819789

#simps(抛物线求积公式)

def f(x):

return

np.exp(-x)

x=np.linspace(0,1,10)

integrate.simps(f(x),x) #0.63219734897921753

#romberg(龙贝格求积)

f=lambda x:np.exp(-x)

integrate.romberg(f,0,1) #0.63212055882856977

python怎么算积分_python求积分相关推荐

  1. python怎么算积分_python计算积分

    python有多个方法计算积分,下面介绍其中三个,以下式为例: 方法一:直接用numpy计算 start = 1 stop = 2 length = 101 x = np.linspace(start ...

  2. matlab中对于xf(x)的积分,[matlab 积分]MATLAB求积分?

    MATLAB求积分? 问题补充:我想用MATLAB求一下这个式子的积分,谁能帮我一下?g(f)=(e^4kxf)/sinh(kctf)未知量是f,最后求出来的式子里保留其他几个字母. ●matcom ...

  3. python怎么算阶乘_python计算阶乘

    广告关闭 腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! 问题本身很简单,主要是通过这个小问题来演示python的一些用法,例如测试代码运 ...

  4. python求积分_python 求定积分和不定积分示例

    求f(x) = sin(x)/x 的不定积分和负无穷到正无穷的定积分 sin(x)/x 的不定积分是信号函数sig ,负无穷到正无穷的定积分为pi import math import numpy a ...

  5. python怎么计算积分_python 求定积分和不定积分示例

    求f(x) = sin(x)/x 的不定积分和负无穷到正无穷的定积分 sin(x)/x 的不定积分是信号函数sig ,负无穷到正无穷的定积分为pi import math import numpy a ...

  6. python怎么计算积分_Python 计算积分实例

    import numpy as np from scipy import integrate import matplotlib.pyplot as plt plt.figure(figsize=(9 ...

  7. python 列表算平均分_python平均列表

    广告关闭 腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! 我有一个列表:类似: data = , , , ] 我想把这个算出来... 本节 ...

  8. python计算多边形面积_Python求凸包及多边形面积教程

    一般有两种算法来计算平面上给定n个点的凸包:Graham扫描法(Graham's scan),时间复杂度为O(nlgn):Jarvis步进法(Jarvis march),时间复杂度为O(nh),其中h ...

  9. python中最小公倍数函数_python求最大公约数和最小公倍数的简单方法

    python怎么求最大公约数和最小公倍数 一.求最大公约数 用辗转相除法求最大公约数的算法如下: 两个正整数a和b(a>b),它们的最大公约数等于a除以b的余数c和b之间的最大公约数.比如10和 ...

最新文章

  1. android4被淘汰,这五类安卓手机已经被淘汰,购买时别贪便宜,需慎重考虑!
  2. hdu 5340(manacher+枚举)
  3. Django模板语言中的自定义方法filter过滤器实现web网页的瀑布流
  4. 深入react技术栈(10):受控组件和非受控组件
  5. ioc spring 上机案例_Spring的IoC入门案例
  6. 修改了xml要不要重新起服务器,关于设置:Eclipse每次运行项目时都会修改server.xml(运行-在服务器上运行)...
  7. 中国移动公布2019年智能家庭网关集采结果:华为、中兴中标
  8. php和xml区别,html与xhtml和xml有什么区别
  9. (web前端每日一练)用JS制作一个能够在桌面运行的倒计时插件
  10. 中达优控触摸屏编程视频教程_中达优控触摸屏软件下载 YKBuilder(中达优控触摸屏编程软件) v5.0.200 官方安装免费版 下载-脚本之家...
  11. redis获取哈希类型的值,hget、hmget、hgetall
  12. 【新手必看】全网最全平面设计理论基础知识,平面设计重要理论
  13. 商汤科技创业领导_从技术专家到领导者的创业之旅
  14. ecshop 添加php标签,ECSHOP模板标签【ecshop标签大全】ecshop标签手册
  15. UGUI_03_补充之_Image的属性(image type这个属性simple、Sliced、tiled、filled样式详解)
  16. Linux企业运维篇——git+gitlab+jenkins+docker构成持续集成环境
  17. 项目实训(十六)——总结
  18. TensorFlow.js简介
  19. Python-- error: [WinError 3] 系统找不到指定的路径。: 'C:\\Program Files (x86)\\Microsoft Visual Studio 14.0\\VC
  20. 常见html的标题含义(1)

热门文章

  1. c语言 北京时间转换utc时间_c语言实战 BJT时间转化位UTC时间
  2. autojs批量删除好友源码,安卓免root引流脚本
  3. 【操作系统原理】1.操作系统概述
  4. 大学课程 | 《计算机操作系统》详细知识点总结
  5. 企业库存过大应怎样处理?
  6. 将DIB复制到剪贴板
  7. html 素材表(模板、日期控件、对话框样式、 翻页 、框架)
  8. json_tuple,get_json_object提取日志埋点信息
  9. js时间戳转换为日期字符串
  10. Hibernate使用原生SQL语句(left join左连接查询)