python怎么算积分_Python求解数值积分-定积分求解
本小节求解下述定积分:
数值积分-integrate
integrate模块提供了好几种数值积分的方法,包括常微分方程组(ODE)的数值积分。相关函数列表如下:quad() - 一元定积分
dblquad() - 二元定积分
triquad() - 三元定积分
odeint() - 计算常微分方程组的数值解
本文节选自作者的《Python编程基础及应用》视频教程。Python编程基础及应用_哔哩哔哩 (゜-゜)つロ 干杯~-bilibiliwww.bilibili.com
微实践 - 定积分求解
本小节求解下述定积分:
, 为了方便说明,我们先使用下述代码画出示意图:
import numpy as np
from matplotlib import pyplot as plt
x = np.linspace(0,6,1000)
y = np.cos(2*np.pi*x)*np.exp(-x)+1.2
plt.axis([np.min(x),np.max(x),0,np.max(y)]) #坐标范围
plt.plot(x,y,label="$cos(2πx)e^{-x}+1.2$") #画曲线,带图示
plt.fill_between(x,y1=y,y2=0,where=(x>=0.7)&(x<=4), #填充积分区域
facecolor='blue',alpha=0.2)
plt.text(0.5*(0.7+4),0.4,r"$\int_{0.7}^4(cos(2πx)e^{-x}+1.2)\mathrm{d}x$",
horizontalalignment='center',fontsize=14) #增加说明文本
plt.legend() #显示图示
plt.show()
执行结果:
plt.axis()函数设定了图的坐标范围。fill_between(x,y1=y,y2=0,where=(x>=0.7)&(x<=4)...)则用于填充积分区域,其中,x和y1构成曲线1; x和y2=0构成曲线2(也就是横坐标线);该函数填充两条曲线之间x值域为[0.7,4]的部分,where参数指明了这个值域。facecolor指定填充颜色,alpha参数指定透明度。
plt.text()则在图上添加文本,前两个参数指定了文本的坐标位置,horizontalalignment='center'要求文本在指定的位置水平居中摆放(指定位置位于文本的水平中心)。r"$...$"为文本内容:字符串前加表示放弃对字符串内的内容进行\转义;两个"为文本内容:字符串前加r表示放弃对字符串内的内容进行\转义;两个$包含起来说明其中的内容为LaTeX格式的公式。
显然,上述定积分就是上图中阴影部分的面积。
方法1:分成小矩形,计算面积和
import numpy as np
x = np.linspace(0.7,4.0,1000)
y = np.cos(2*np.pi*x)*np.exp(-x)+1.2
dx = x[1] - x[0] #每个矩形的宽度
fArea = np.sum(y*dx) #矩形宽*高,再求和
print("Integral area:",fArea)
执行结果:
Integral area: 4.032803310221616
上述代码中,把曲线的阴影部分分成1000个矩形,每个矩形的宽都是dx,第i个矩形的高则是yi。每个矩形的长乘宽,再求和,得积分面积。
方法2:使用quad()函数进行积分
import math
from scipy import integrate
def func(x):
print("x=",x) #用于展示quad()函数对func的多次调用
return math.cos(2*math.pi*x)*math.exp(-x)+1.2
fArea,err = integrate.quad(func,0.7,4)
print("Integral area:",fArea)
执行结果:
x= 2.35
x= 0.7430542279466668
x= 3.9569457720533334
x= 2.4613227224815875
...
x= 3.4178741117287044
Integral area: 4.029065401143393
首先,我们定义了一个函数func(),它根据x计算y值。当对单个数值进行计算时,numpy的ufunc并不具备速度优势,所以我们使用了math模块。
integrate.quad()专门用于计算一元定积分,fArea,err = integrate.quad(func,0.7,4)取x值域[0.7,4]进行数值积分,在积分过程中,会反复调用func()函数计算y值。其返回一个元组,包括积分结果及误差。
integrate.quad()计算的积分会比方法1的矩形面积求和方法更加精确。
本文节选自作者的B站MOOC及同名教材:Python编程基础及应用 — 重庆大学 高等教育出版社,作者亲授_哔哩哔哩 (゜-゜)つロ 干杯~-bilibiliwww.bilibili.com
版权声明 本文可以在互联网上自由转载,但必须:注明出处(作者:海洋饼干叔叔)并包含指向本页面的链接。 本文不可以以纸质出版为目的进行改编、摘抄。
python怎么算积分_Python求解数值积分-定积分求解相关推荐
- python怎么算积分_python使用sympy不定积分入门及求解
本文章首发于古客 如果想要获得更好的阅读体验,请前往我的博客! 概述 为了解决数学中一些不定积分的求解方便自己做作业,让我们一起学习如何使用python求不定积分吧! <matlab版点这里&g ...
- python怎么算积分_python求积分
1.符号积分:#通过integrate()功能(facility),SymPy对基本和特殊函数定与不定积分有卓越的支持. #该功能使用有力的扩展Risch-Norman算法,启发算法和模式匹配: fr ...
- python怎么算积分_python计算积分
python有多个方法计算积分,下面介绍其中三个,以下式为例: 方法一:直接用numpy计算 start = 1 stop = 2 length = 101 x = np.linspace(start ...
- python怎么计算积分_python 求定积分和不定积分示例
求f(x) = sin(x)/x 的不定积分和负无穷到正无穷的定积分 sin(x)/x 的不定积分是信号函数sig ,负无穷到正无穷的定积分为pi import math import numpy a ...
- python怎么算积分_利用python求积分的实例
利用python求积分的实例 python的numpy库集成了很多的函数.利用其中的函数可以很方便的解决一些数学问题.本篇介绍如何使用python的numpy来求解积分. 代码如下: # -*- co ...
- python 列表算平均分_python平均列表
广告关闭 腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! 我有一个列表:类似: data = , , , ] 我想把这个算出来... 本节 ...
- python怎么算阶乘_python专家写阶乘 ()用python计算阶乘
请编写一个程序实现分享n的阶乘(即n.),并打印输出结果.要分享包括两个函def little_than_50(x): if x 语言很多时候是假的,一起经历过的事情才是真的. 这个用python写的 ...
- python怎样算入门_python初学者怎么入门
据百度大数据统计每月有4260705的用户在百度检索Python相关介绍,在大众周围,也总是充斥着各种如何学习Python的声音,许多人已经对于是否"应该学习Python"也都发表 ...
- python怎么算阶乘_Python 计算阶乘的算法
希望你的坚持是因为热爱,而不是不甘 by:缘分落地 阶乘在高中的数学排列组合曾经出现过,相信各位都不陌生.举个栗子: N!= N * (N-1) * (N-2)......1 带入数字形象的看一遍,我 ...
- python怎么算阶乘_python计算阶乘
广告关闭 腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! 问题本身很简单,主要是通过这个小问题来演示python的一些用法,例如测试代码运 ...
最新文章
- 网络配置——Linux入门到精通的难关之一
- 移植fatfs上电复位前两次不能正确打开文件解决方法
- ActiveMQ整合spring
- centos软件软件包
- 手机闪存速度排行_2020年双十二3000-4000元高性价比手机推荐!
- javacc案例之统计字符
- flask及扩展源码解读
- 二元偏导数存在的条件_偏导数连续怎么证明
- shell脚本头,#!/bin/sh与#!/bin/bash的区别.
- OutofMemory异常
- 图像 理想低通滤波_高通滤波、低通滤波、带通滤波 python实现
- 【论文笔记】From Facial Parts Responses to Face Detection: A Deep Learning Approach
- MySQL可视化工具使用
- Java生态技术体系科普
- 邮件中的抄送和密送的区别
- 移动云迁移工具:Hyper-V虚拟化迁移到移动云
- property中ref、value、name的区别
- 一道疯狂bypass的题目
- 华为2285v2服务器修改SN,华为RH2285 V2服务器升级bios
- ts-学习类型-class
热门文章
- 微软自带dns服务器,win10微软设置哪个DNS服务器地址最快
- Flink1.13.0 + Hudi 0.11.1 + Hive2.1.1 + presto0.273.3 + yanagishima 18.0
- 毕业生见习期考核鉴定表(单位考核意见)
- php判断是否连续出现数字,php通过gbk编码判断 含有连续数字 可用于判断QQ号,手机号等。...
- urllib爬虫封装(可设置代理,记录日志)方法
- 实体消歧(链接到实体库)
- IExplore.exe应用程序错误解决方法
- 【项目管理】项目管理四要素
- 8.	Zigbee应用程序框架开发指南 - 事件
- Ubuntu常用命令汇总