利用概率密度函数生成随机变量过程中出现超几何分布函数问题
如题,已知概率密度函数,拟依据其生成随机变量,想法是先对其积分求累计概率分布,然后由(0,1)均匀分布随机数来生成随机变量,但是,其累计概率密度函数是一个含有hypergeom(超几何分布函数)的函数,怎么办呢???具体问题(符合贝塔分布的一个概率密度函数):
f=int((1/216815.2)*((x/196000)^-0.05)*(1-(x/196000))^-0.05);MATLAB运行结果是:
f =
2374999999999999993
-------------------
-5 2500000000000000000
0.8928931123 10 x hypergeom(
-5
[0.05000000000, 0.9500000000], [1.950000000], 0.5102040816 10 x)
当我继续运行
f=rand;
x=
x
这可怎么办呢???急死了,弄了一天了谢谢各位啊!
附上函数:负荷贝塔分布的概率密度函数,其中P为待求随机变量,R为已知数