在python中创建连续分布

2023-11-21

我在 python 中创建连续分布时遇到了麻烦,它真的开始让我烦恼。我已经阅读并重新阅读这个Python指南(scipy指南)这对我的问题没有帮助。

我的代码如下:

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(#what to put here).

我尝试了很多不同的方法,主要是中建议的方法本文档(scipy 指南),就像设置我的范围值上限和下限一样a=,b=将其设置为pdf or a 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] called WindDir对于给定的分布freq? 老实说,我已经仔细浏览了 Google 和 stackoverflow 网站,使用关键字、带标签的关键字和单独的标签进行搜索,但找不到解决方案。

编辑1-期望的结果我希望输出是介于0,360 or 0,2pi


好吧,为了使用rv_continuous您需要提供某种概率密度函数。在下面的示例中,我针对给定的风向间隔 [0,360) 实现了累积密度函数。我通过在输入中指定的最近两个风向之间插值概率密度函数来实现此目的。注意参数a and b中指定的rv_continuous基类构造函数...这些指定所考虑的间隔的最小值和最大值。尝试一下代码,如果您有任何疑问,请询问,我会尽力帮助澄清。

Edit我修改了 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]
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

在python中创建连续分布 的相关文章

随机推荐