任意一元分布随机数
(1)离散分布随机数
给出只取有限个值的离散总体x的分布列
x x1 x2 . xn
PPP1P2.pn
通过函数randsample函数，可以生成遵从分布的随机数。
例：离散母集团x的分布如下
X -2 -1 0 1 2
P 0.05 0.2 0.5 0.2 0.05
调用randsample函数生成100个遵从该分布的随机数，调用tabulate函数统计各数字的出现频率和频率
xvalue=[-2 -1 0 1 2]； %定义取值的向量
xp=[0.05 0.2 0.5 0.2 0.05]； %定义概率向量，与取值向量一一对应
调用randsample函数生成100个遵循指定离散分布的随机数
x=randsample(xvalue，100，true，xp )； %为true表示存在返回样本，false表示没有返回样本
reshape(x，[ 10，10 ] )将向量x置换为10行10列的矩阵，便于观察
ans=
2 0 1 2 1 0 1 0 1 0
0-10-22-10
1 0 0 0 2 0 0 -1 -2 1
1 1 1 1 0 0 0 0 0 0
0 -1 1 -2 1 0 -1 0 2 0
1 0 1 -1 0 -1 -1 0 0 1
0 2 -1 1 -1 -1 -1 0 0 0
0 0 1 0 0 0 1 0 -1 -1
0 0 0 0 -2 0 0 -1 -2 -1
1 0 -1 0 0 -2 0 0 2 1
tabulate(x ) %调用tabulate函数统计各数字出现的频率和频率
值计数百分比
-2 6 6.00%
-1 17 17.00%
0 49 49.00%
1 21 21.00%
2 7 7.00%
从统计结果来看，随机数的频率分布与真实分布相差不大，但不完全相同，生成足够多的随机数后，就会逐渐接近真实分布。
函数randsrc和无情的睫毛膏也可以实现同样的功能。
)2)连续分布随机数
MATLAB统计工具箱提供了slicesample函数，用于生成任何指定的分布随机数
rand=slecesample(initial，nasample，‘pdf’，pdffun () ) ) ) ) ) ) ) ) ) )。
初始：指定初始值
nasample :产生的随机数的个数
“pdf”:概率分布函数
pdffun :概率分布函数的句柄
例：整体x服从抛物线分布，其概率密度函数如下
然后调用slicesample函数生成1000个遵循其分布的随机数，创建频率直方图并将其与实际密度函数曲线进行比较。
pdffun=@(x )6*x ) )1-x )； %用匿名函数定义密度函数，取pdffun=@(x )6*x ) (1-x ) ) *(x0x1) )
对于%x1或x0，密度函数不能小于0，因为pdffun0。 MATLAB直接
%为pdffun=0，因此不需要写为pdffun=@(x )6*x ) (1-x ) ) *(x0x1)，而
%否则，必须指定x的范围，如下例所示
调用slicesample函数生成1000个遵循指定密度函数的随机数
x=slice sample (0.5，1000，' pdf '，pdffun )； %指定的初始值为0.5，需要{pdffun(0.5 )0}
[fp，XP]=ecdf(x )； %计算经验累积概率分布函数值
ecdfhist(FP，xp，20 )； %频率创建直方图
保持接通；
FPlot(pdffun，[0 1]，' r ' )； %fplot用于绘制数学函数，plot为向量或矩阵列绘制曲线
xlabel(x )；
ylabel(f ) x ) )；
legend (“频率直方图”、“密度函数曲线”)；
例：整体x服从三角分布，其概率密度函数如下
调用slicesample函数生成1000个遵循其分布的随机数，创建频率直方图，并将其与真正的密度函数进行比较。
pdffun=@(x ) x * (x=0x1)2-x ) ) *(x=1x2 )； %密度函数
调用slicesample函数生成1000个遵循分布的随机数
x=slice sample (1.5，1000，' pdf '，pdffun )；
[fp，XP]=ecdf(x )； %计算经验累积概率分布函数值
ecdfhist(FP，xp，20 )； %频率创建直方图
保持接通；
FPlot(pdffun，[0 2]，' r ' )； %绘制真正的密度函数曲线
xlabel(x )；
ylabel(f ) x ) )；
legend (“频率直方图”、“密度函数曲线”)；