我试图用scipy计算一个二重积分。被积函数有点复杂,因为它包含一些概率分布,以加权x和y的每个值的可能性(就像一个混合模型)。以下代码的计算结果为负数,但它应被[0,1]绑定。另外,计算大约需要半个小时。在

我有两个问题。在

1)有没有更好的方法来计算这个积分?在

2)这个负值从何而来?对我来说,最大的问题是如何加快计算速度,因为我可以在我的代码中找到导致我以后自己负的错误。在from scipy import stats

from scipy.integrate import dblquad

import itertools

p= [list whose entries are each different stats.beta(a,b) distributions]

def integrand(x,y):

delta=x-y

marg=0

for distA,distB in itertools.permutations(p,2):

first=distA.pdf(x)

second=distB.pdf(y)

weight1=0

weight2=0

for distC in p:

if distC == distA:

continue

w1=distC.cdf(x)-distC.cdf(y)

if weight1 == 0:

weight1=w1

else:

weight1=weight1*w1

marg+=(first*weight1*second)

I=delta*marg

return I

expect=dblquad(integrand,0,1,lambda x: 0, lambda x: x)

这本质上是问两点之间的最大距离在分布向量中的期望值是多少。积分的极限是y∊[0,x]和x∊[0,1]。这给了我-.49,估计的积分误差在10e-10级,所以它不应该是由于积分方法。在

我已经和这件事斗争了一段时间了,感谢你的帮助。谢谢。在

编辑:更正了打字错误

python可以计算复杂积分吗_python中二重积分的正确计算相关推荐

  1. Python基础_第3章_Python中的循环结构

    Python基础_第3章_Python中的循环结构 文章目录 Python基础_第3章_Python中的循环结构 Python中的循环结构 一.回顾分支练习题 1.判断是否为一个合法三角形 2.求世界 ...

  2. Python基础_第5章_Python中的数据序列

    Python基础_第5章_Python中的数据序列 文章目录 Python基础_第5章_Python中的数据序列 Python中的数据序列 一.字典--Python中的==查询==神器 1.为什么需要 ...

  3. python基本统计量_Python中简单统计量的计算

    本篇文章给大家带来的内容是关于Python中简单统计量的计算,有一定的参考价值,有需要的朋友可以参考一下,希望对你有所帮助. 1.这些操作都要确保已经在电脑中安装好了Anaconda集成库,如果安装好 ...

  4. python datetime计算时间差_Python中关于日期的计算总结

    1.获取当前时间的两种方法: 代码如下: 2.获取上个月最后一天的日期(本月的第一天减去1天) 代码如下: 3.获取时间差(时间差单位为秒,常用于计算程序运行的时间) 代码如下: 4.计算当前时间向后 ...

  5. python通过经纬度算出距离_Python中如何利用经纬度进行距离计算

    前些天由于工作中需求,要计算经纬度之间的距离,当我去网上搜索距离计算的方法时,发现很多文章中的方法乍一看都是很不同的,同是进行距离计算,为啥这么不一样呢?后来才发现问题在于很多文章没有进行相关的原理说 ...

  6. python里的shell是什么_python中的shell操作

    http://blog.csdn.net/meng_tianshi/article/details/6682317 首先介绍一个函数: os.system(command) 这个函数可以调用shell ...

  7. python里的log怎么表示_python中的对数log函数表示及用法

    在数学运算中,如果没有计算器,对于很大的数字相乘,我们花费大量的时间计算,而且一旦出错,就要重新计算,很是麻烦.其实对于数字相乘,不依靠靠计算器,想要准确简单的运算的方法不是没有,那就是对数和指数,他 ...

  8. python 获取用户的一个输入值_Python中,用于获取用户输入的命令为:

    [多选题]以下关于机器学习说法正确的是? [判断题]Python内置函数sum____用来返回数值型序列中所有元素之和. [单选题]关于自定义函数的下列说法不正确的是: [判断题]Python内置函数 ...

  9. python选取元音开头的单词_Python中的Regex,用于查找遵循以下模式的单词:元音、辅音、元音、辅音...

    如果将辅音有向图映射成单个辅音,则最长的单词是解剖病理学的10*VC字符串. 如果您正确映射y,那么您将得到完整的字符串,如乙酰丙酮作为8*VC,下胚轴作为8*CV. 如果不需要字符串是完整的,那么在 ...

最新文章

  1. 【自用】 sklearn 用 train_test_split 简单划分训练和测试集
  2. mysql双机热备份的实现步骤
  3. DrawerLayout + Toolbar + ViewPager
  4. php_os用法,PHP教程:PHPUnit学习笔记(二)PHPUnit基本用法
  5. 如何学好初中计算机,初中生怎么学习方法好 十大方法告诉你
  6. 教你如何使用redis分布式锁
  7. Python笔记-相关性分析(连续变量和分类变量)
  8. 给你的网站添加运行时间
  9. PageRank算法原理与实现
  10. 浪潮服务器2016销售额,Gartner:2016年前三季度浪潮服务器销售额中国第一
  11. 16. Use the same form in correspondng uses of new and delete
  12. 一个故事告诉你什么才是好的程序员
  13. Vb6 Modbus TCP通讯示例源代码
  14. 六道腾讯、百度、美团常爱问的面试智力题和答案
  15. 用Python弹奏周董的《菊花台》是什么效果
  16. PCB天线设计及射频布局设计指南
  17. 【解决方案】scrapy报错KeyError: ‘Spider not found‘
  18. 中国男人到底配不配得上中国…
  19. 2020美容师(中级)复审模拟考试及美容师(中级)模拟考试系统
  20. 使用 psycopg2 时遇到的两个坑

热门文章

  1. python比较两个xml_python比较2个xml内容的方法
  2. r语言ggplot2一夜多图_ggplot2简介
  3. flask uwsgi nginx 部署在 ubuntu 上
  4. flask传值给js
  5. 自然语言处理-Word2Vec
  6. Python-PyCharm 报错解决:ImportError: cannot import name 'InteractiveConsole' from 'code'
  7. LeetCode 91. Decode Ways--动态规划DP的Python和Java解法
  8. linux如何添加默认路由表_Linux-如何添加路由表
  9. 关于int main( int argc, char* argv[] ) 中arg和argv参数理解
  10. proftpd mysql_Proftpd(MySQL验证+配额)