levy distribution 是我们经常学习使用的一个密度函数,那么对于这样的密度函数一般定义都是由一些超越函数定义的,但是超越函数又是由一些无穷级数定义的。所以。。。操作起来比较麻烦。下面介绍一种略微简单的方法。
方法/步骤
1
首先确定定义,选择一个与你定义向一致的特征函数定义。
2
然后设定相关参数 alpha beta mu 以及 sigma 的具体参数。
3
如果首先确定该密度函数在某一处的取值,可以通过如下:PDF[StableDistribution[1, 1.5, 1, 0, 2], x]
4
对于密度函数,我们可以使用plot函数进行画出来。Plot[Table[ PDF[StableDistribution[0, \[Alpha], 1.5/2, 0, 2], x], {\[Alpha], {1/3, 2/3, 2}}] // Evaluate, {x, -6, 6}, Filling -> Axis]
5
另一个我们比较关注的均值和方差可以通过mean和Variance 进行求得。
6
如果需要生成对应的随机变量可以通过a = RandomVariate[StableDistribution[0, 1.45, -0.5, 1, 2], 10^2];