- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我在 python 中创建连续分布时遇到问题,它真的开始让我烦恼。我已经阅读并重读了this python guide (scipy guide)这对我的问题没有帮助。
我的代码是这样的:
import sys
import scipy.stats
import numpy
def CDF_Random(N,NE,E,SE,S,SW,W,NW,Iterations):
WindDir = [0,45,90,135,180,225,270,315]
Freq = N,NE,E,SE,S,SW,W,NW
mydist = scipy.stats.rv_continuous(#My problem is what to write here)
cdf_rand=mydist.rvs(size=Iterations)
return (cdf_rand)
if __name__ == '__main__':
N = float(sys.argv[1])
NE = float(sys.argv[2])
E = float(sys.argv[3])
SE = float(sys.argv[4])
S = float(sys.argv[5])
SW = float(sys.argv[6])
W = float(sys.argv[7])
NW = float(sys.argv[8])
Iterations = float(sys.argv[9])
numpy.set_printoptions(threshold=Iterations)
sys.stdout.write(str(CDF_Random(N,NE,E,SE,S,SW,W,NW,Iterations)))
如您阅读代码时所见,我的问题是知道将什么放在括号中以创建连续分布。 scipy.stats.rv_continuous(#这里放什么)
。
我尝试了很多不同的东西,主要是 this document(scipy guide) 中建议的那些,就像设置我的上限和下限值 a=,b=
将其设置为 pdf
或 ppf
。我已经尝试使用在命令行中输入的那些或我写入代码本身的那些[arrays]
。
我从命令行运行这个命令python C:\Users\...\CDF.py 0.01 0.01 0.01 0.01 0.01 0.93 0.01 0.01 10
每次我得到;RuntimeError: maximum recursion depth exceeded
我曾尝试将递归深度重置为不同的值,但这不起作用或使 python 崩溃。sys.setrecursionlimit(10000)
所以基本上应该在 scipy.stats.rv_continuous()
之后的括号中输入什么来创建 [array]
的连续分布,称为 WindDir
对于给定的分布 freq
?老实说,我已经通过谷歌和 stackoverflow 网站仔细查看,使用关键字、带有标签的关键字和单独的标签进行搜索,但找不到解决方案。
编辑 1-期望的结果我希望输出是 0,360
或 0,2pi
最佳答案
好的,所以为了使用 rv_continuous
,您需要提供某种概率密度函数。在下面的示例中,我为给定的风向间隔 [0,360) 实现了累积密度函数。我通过在输入中指定的最近两个风向之间插入概率密度函数来实现这一点。请注意 rv_continuous
基类构造函数中指定的参数 a
和 b
...这些参数指定了所考虑区间的最小值和最大值。试用代码,如果您有任何问题,请提出,我会尽力帮助澄清。
编辑 我修改了 python 3 的代码,并更新了 cdf 以更准确地在基本方向给定的频率之间进行插值。
import scipy.stats
class rvc(scipy.stats.rv_continuous):
def __init__(self, freqs):
super().__init__(a=0,b=359.9999)
self.WindDir = [0.,45.,90.,135.,180.,225.,270.,315.,360.]
self.Freqs = freqs
def _cdf(self, x):
return [self.do_cdf(i) for i in x]
def do_cdf(self, x):
if x < 0: return 0.0
if x >= 360: return 1.0
v = 0.0
for i in range(9):
if x >= self.WindDir[i]:
v += self.Freqs[i]
else:
v += (self.Freqs[i]-self.Freqs[i-1])*(x-self.WindDir[i-1])/45.
break
return v
rv = rvc([0.01,0.01, 0.01, 0.01, 0.01, 0.01, 0.92, 0.01, 0.01])
print(rv.rvs(size=100))
运行它会产生类似的输出
[ 305.76400244 271.49390379 285.5514374 280.4422151 302.52158889
273.85068415 278.8377575 287.05260894 270.66219523 301.97131911
277.77725392 293.29633122 273.60832876 278.90703722 271.12590324
288.38285244 299.89362412 276.68015144 278.46429959 304.33532312
313.42248626 272.10226608 311.00385667 284.48822492 298.134523
300.08266465 272.17850056 270.07995773 359.9999 289.13032076
206.58066169 308.95821915 288.998036 280.93451276 294.33703562
278.02828894 307.96335447 292.50172701 313.67335738 291.99562061
275.98551449 307.72620259 283.11307243 309.6052904 283.65619152
305.80365684 303.43088801 314.74811251 302.8390404 288.16580807
299.84487396 302.34636138 291.09465231 310.74393544 279.23446355
309.83353391 292.62761642 273.22324646 291.16193395 298.51520679
299.87754111 280.29779055 276.76741796 285.5186257 301.7074023
274.9771402 280.6619726 276.53202603 289.50757382 313.61213159
98.86469637 298.41091812 83.62474126 290.42538277 306.70218844
270.81939255 290.25089647 305.3277742 278.03965968 280.5979916
307.9492377 284.30184233 307.14788891 283.33779011 270.86398644
297.70610336 277.48672772 300.35147777 271.07308885 304.41450287
312.5640489 273.54296504 311.62527023 298.95246144 275.199183
302.6063864 359.9999 311.25861396 306.83491033 313.52858514]
关于python - 在 python 中创建连续分布,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/17676710/
前言: 有时候,一个数据库有多个帐号,包括数据库管理员,开发人员,运维支撑人员等,可能有很多帐号都有比较大的权限,例如DDL操作权限(创建,修改,删除存储过程,创建,修改,删除表等),账户多了,管理
所以我用 Create React App 创建并设置了一个大型 React 应用程序。最近我们开始使用 Storybook 来处理和创建组件。它很棒。但是,当我们尝试运行或构建应用程序时,我们不断遇
遵循我正在创建的控件的代码片段。这个控件用在不同的地方,变量也不同。 我正在尝试编写指令来清理代码,但在 {{}} 附近插入值时出现解析错误。 刚接触 Angular ,无法确定我错过了什么。请帮忙。
我正在尝试创建一个 image/jpeg jax-rs 提供程序类,它为我的基于 post rest 的 Web 服务创建一个图像。我无法制定请求来测试以下内容,最简单的测试方法是什么? @POST
我一直在 Windows 10 的模拟器中练习 c。后来我改用dev C++ IDE。当我在 C 中使用 FILE 时。创建的文件的名称为 test.txt ,而我给出了其他名称。请帮助解决它。 下面
当我们创建自定义 View 时,我们将 View 文件的所有者设置为自定义类,并使用 initWithFrame 或 initWithCode 对其进行实例化。 当我们创建 customUITable
我正在尝试为函数 * Producer 创建一个线程,但用于创建线程的行显示错误。我为这句话加了星标,但我无法弄清楚它出了什么问题...... #include #include #include
今天在做项目时,遇到了需要创建JavaScript对象的情况。所以Bing了一篇老外写的关于3种创建JavaScript对象的文章,看后跟着打了一遍代码。感觉方法挺好的,在这里与大家分享一下。 &
我正在阅读将查询字符串传递给 Amazon 的 S3 以进行身份验证的文档,但似乎无法理解 StringToSign 的创建和使用方式。我正在寻找一个具体示例来说明 (1) 如何构造 String
前言:我对 C# 中任务的底层实现不太了解,只了解它们的用法。为我在下面屠宰的任何东西道歉: 对于“我怎样才能开始一项任务但不等待它?”这个问题,我找不到一个好的答案。在 C# 中。更具体地说,即使任
我有一个由一些复杂的表达式生成的 ILookup。假设这是按姓氏查找人。 (在我们简单的世界模型中,姓氏在家庭中是唯一的) ILookup families; 现在我有两个对如何构建感兴趣的查询。 首
我试图创建一个 MSI,其中包含 和 exe。在 WIX 中使用了捆绑选项。这样做时出错。有人可以帮我解决这个问题。下面是代码: 错误 error LGH
在 Yii 中,Create 和 Update 通常使用相同的形式。因此,如果我在创建期间有电子邮件、密码、...other_fields...等字段,但我不想在更新期间专门显示电子邮件和密码字段,但
上周我一直在努力创建一个给定一行和一列的 QModelIndex。 或者,我会满足于在已经存在的 QModelIndex 中更改 row() 的值。 任何帮助,将不胜感激。 编辑: QModelInd
出于某种原因,这不起作用: const char * str_reset_command = "\r\nReset"; const char * str_config_command = "\r\nC
现在,我有以下由 original.df %.% group_by(Category) %.% tally() %.% arrange(desc(n)) 创建的 data.frame。 DF 5),
在今天之前,我使用/etc/vim/vimrc来配置我的vim设置。今天,我想到了创建.vimrc文件。所以,我用 touch .vimrc cat /etc/vim/vimrc > .vimrc 所
我可以创建一个 MKAnnotation,还是只读的?我有坐标,但我发现使用 setCooperative 手动创建 MKAnnotation 并不容易。 想法? 最佳答案 MKAnnotation
在以下代码中,第一个日志语句按预期显示小数,但第二个日志语句记录 NULL。我做错了什么? NSDictionary *entry = [[NSDictionary alloc] initWithOb
我正在使用与此类似的代码动态添加到数组; $arrayF[$f+1][$y][$x+1] = $value+1; 但是我在错误报告中收到了这个: undefined offset :1 问题:尝试创
我是一名优秀的程序员,十分优秀!