python数值积分_Python-Scipy进行数值积分
Python的Scipy模块中拥有大量的数值计算函数,方便我们快速进行数值计算。
Scipy中的integrate模块提供了几种数值积分算法,导入方式为:
from scipy import integrate
使用integrate时,需要先将要进行积分的方程定义为函数。求取一至三重积分的函数分别为:
integrate.quad(func,a,b,args,full_output)
integrate.dblquad(func,a,b,gfun,hfun,args,epsabs,epsrel)
integrate.tplquad(func,a,b,gfun,hfun,qfun,rfun,args,epsabs,epsrel)
以三重积分为例。func为运算对象函数,形式为func(z,y,x)。a,b对应变量x的积分区域,gfun,hfun对应变量y的积分区域,依次类推。
注意gfun,hfun等的形式应为函数,其中gfun,hfun是自变量为x的函数,qfun,rfun是自变量为x,y的函数。这些函数可以使用lambda函数进行定义,形式通常为:
lambda x,y:x*y
如果是常函数,则定义为:
lambda x:0
lambda x,y:1
args可选,为传递给func的格外参数;
full_output可选,非零则返回积分信息的dictionary。如果非零,则还会禁止显示警告消息,并将消息附加到输出元组。
epsabs可选,绝对容差直接传递到内部1-D正交积分。默认值为1.49e-8。
epsrel可选,内部1-D积分的相对容差。默认值为1.49e-8。
运算结果输出一个元组,第一项为积分结果,第二项为绝对误差,此外还有收敛情况等信息。
以下是一个用integrate.tplquad()计算立方体晶体的LJ势能的例子:
#To calculate the LJ potential energy between a cubic and an atom
#Author: Lewisbase
#Date: 2019.01.01
from __future__ import division
import numpy as np
from scipy import integrate
#Cubic length width hight A,B,C
A=10
B=10
C=10
#Average sigma(s,A) epsilon(e,K) rho(r,A^-3) of the cubic
s=np.array([4.24,3.84]); e=np.array([208.4541197,69.0609]); r=0.061540355
def LJ_cubic(x,y,z):
return 4*eps*rho*((sig**2/((xb-x)**2+(yb-y)**2+(zb-z)**2))**6-(sig**2/((xb-x)**2+(yb-y)**2+(zb-z)**2))**3)
lx=50
ly=50
lz=100
dx=lx/101
dy=ly/101
dz=lz/201
x0=lx/2-A/2
x1=lx/2+A/2
y0=ly/2-B/2
y1=ly/2+B/2
z0=lz/2-C/2
z1=lz/2+C/2
for m in np.linspace(0,1,2):
for ix in np.linspace(0,100,101):
for iy in np.linspace(0,100,101):
for iz in np.linspace(0,200,201):
xb=ix*dx
yb=iy*dy
zb=iz*dz
sig=s[m]
eps=e[m]
rho=r
if abs(xb-lx/2)
LJ = 1E100
else:
LJ=integrate.tplquad(LJ_cubic,
z0,
z1,
lambda y: y0,
lambda y: y1,
lambda y,x: x0,
lambda y,x: x1)
LJ=LJ[0]/298.15
print "%f %f %f %f\n"%(xb,yb,zb,LJ)
有关数值积分的方法还有高斯积分,龙贝格积分等,自己暂时还没有用到,待到以后再做详细学习。
参考资料
python数值积分_Python-Scipy进行数值积分相关推荐
- gamma分布python代码_python scipy stats.gamma用法及代码示例
伽玛连续随机变量. 作为一个实例rv_continuous类,gamma对象从中继承了通用方法的集合(完整列表请参见下文),并使用特定于此特定发行版的详细信息来完善它们. 注意: 的概率密度函数gam ...
- ks检验python代码_python scipy stats.kstest用法及代码示例
进行Kolmogorov-Smirnov测试是否合身. 这将针对给定的分布G(x)对观察到的随机变量的分布F(x)进行测试.在原假设下,两个分布相同,F(x)= G(x).替代假设可以是" ...
- python奇异值分解_python scipy 奇异值分解 SVD
svd函数返回三个分解部分u,s,vh,由于s是一个对角矩阵,因此s是一个一维数组,保存对角线上的非零值 s中的每个值与vh的行向量以及u中的列向量对应,默认按照从大到小的顺序排列,他表示与其对应向量 ...
- python求数值积分_Python大数据处理-Scipy基础入门,数值积分计算
温馨提示:阅读本文只需要1分钟,您就可以掌握Scipy进行定积分计算.二重.三重积分.多重积分的计算.继续承接上文学习Scipy科学数据处理,为我们后面Python大数据处理开发打基础.今天主要学习分 ...
- python数值积分_python与计算物理:实现数值积分的Simpson方法
1.[代码]python与计算物理:实现数值积分的Simpson方法. #!/usr/bin/env python #coding=utf-8 ''' Created on 2012年10月27日 @ ...
- python中的scipy基础知识_python中SciPy是什么?
python中Numpy常用于计算二维数组计算,而python的另一个库SciPy库与Numpy有着密切的关系,是需要通过Numpy为基础,同时也是通过Numpy数据来操控科学计算.常见的是插值运算. ...
- 【Python图像处理】SciPy库
目录 1.ndimage模块 1.1 gaussian_filter()函数 1.2 convolve()函数 2.signal模块 2.1 convolve2d()函数 2.2 gaussian() ...
- Python库之Scipy库的简介、安装详细
目录 Scipy库的简介 Scipy库的安装 Scipy库的简介 Scipy高级科学计算库:和Numpy联系很密切,Scipy一般都是操控Numpy数组来进行科学计算.统计分析,所以可以说是基于Num ...
- pandas是基于什么的python库_Python的五大常用库——numpy,pandas,matplotlib等
Numpy Numpy 是python科学计算的基础包,本书大部分内容都基于numpy以及构建于其上的库.其功能有:快速高效的多维数组对象ndarray 用于对数组执行元素级计算以及直接对数组执行数字 ...
- python画代码-Python教程_Python画Mandelbrot集 代码
Python教程_Python画Mandelbrot集 代码 作者:Comet 来源: 课课家 www.kokojia.com点击数:278发布时间:2015-06-19 11:17:19 曼德勃罗集 ...
最新文章
- 河科大c语言上机实验答案,2016年河南科技学院信息工程学院C语言上机编程考研复试题库...
- 223.主成分分析PCA
- 服务器系统重装后要配置吗,服务器系统重装后要配置吗
- angular入门-ngOptions
- 2003DC 升级到2008R2 DC实验
- powerBi odbc 连接impala 实现自助分析
- RedHat el5.0 搭建 Postfix 邮件服务器系统一
- Html.RadioButtonFor和Html.DropDownListFor 用法--备忘
- SQLite: 关于日期的标准化问题
- 0基础Java自学之路
- metro3000 OSN2500 OSN3500 华为传输 华为光端机
- python根据经纬度计算两点之间距离
- 中国 vs 卡塔尔 一场幸运的比赛
- 活动回顾 | 大咖云集“开源安全治理模型和工具”线上研讨会
- 民营股份制企业是什么意思
- Tableau-仪表板
- Python-直线检测-提取邮票编码
- 共享汽车管理系统设计软便件研究
- OpenGL学习脚印:伽马校正(Gamma Correction)
- linux ./ 执行run文件,如何在Ubuntu中执行.bin和.run文件