常见的一元分布随机数
)1) MATLAB统计工具箱中函数名称以rnd三个字符串结尾的函数用于生成一般分布的随机数。 如下表所示
函数名称的说明
binornd二元分布随机数
chi2rnd卡侧分布随机数
exprnd指数分布随机数
frnd F分布随机数
lognrnd对数分布随机数
mnrnd多项式分布随机数
normrnd正态分布随机数
………。
示例:
调用normrnd函数生成1000x3的正态分布随机函数,其中平均值u=75,标准差为8,绘制了各列的频数直方图
调用normrnd函数生成1000行3列的随机数矩阵x,其元素遵循平均值75、标准差8的正态分布
x=normrnd(75、8、1000、3 ); %调用-rnd函数需要传递给其分布函数的参数以及生成矩阵的维数
图形;
HIST(x ); %绘图矩阵每x列的度数直方图
xlabel ('正态分布随机数(mu=75,sigma=8) );
ylabel (度数);
legend (“第一列”、“第二列”和“第三列”);
)2) random函数
random函数是计算每个分布随机数的常用函数,格式如下:
y=random (‘name’,A1,A2,A3,m,n ) )。
name指定是这样的分布函数,取值的是上表MATLAB中常见分布函数名称的简称
A1、A2、A3是分布函数的参数
m,n用于指定生成随机数的行和列
例:用normrnd和random函数分布生成1000x1正态分布随机函数,其中均值u=75,标准差8绘制频数直方图
x=normrnd(75、8、1000、1 ); %normrnd函数生成平均u=75、标准偏差为8的正态随机数
xx=随机(norm )、75、8、1000、1 ); %random函数生成平均u=75、标准偏差为8的正态随机数
图形;
HIST(x )、title、“normrnd频数分布直方图”;
图形;
hist(xx )、title、“随机频数分布直方图”;
示例:调用random函数生成1000x1的二元分布随机数向量,制作频率直方图。 二元分布的参数为n=10,p=0.3
注:要计算频率直方图,首先调用ecdf函数来计算经验累积概率分布函数值,然后调用ecdfhist函数来创建频率直方图。 此参数是ecdf函数的返回值
%调用random函数生成1000行1列的随机数向量x,其要素遵循二元分布b (10,0.3 )
x=随机(bino )、10、0.3、1000、1 );
[fp,XP]=ecdf(x );
ecdfhist(FP,xp,20 ); p表示条形图的数量,默认值为10
xlabel ()二元分布(n=10,p=0.3 )随机数);
ylabel(f ) x ) );
明天继续发表MATLAB如何生成任意一元分布的随机数。