概率论实验 01 - | 使用Matlab产生随机数

一、实验目的

     1.了解随机数的产生方法；

     2.了解常用随机数的概率分布函数、分布律和概率密度函数。

二、实验原理

随机数的产生有好多方法，可以利用乘积法和同余法产生[0,1]之间的均匀分布，然后利用函数变换法产生所需不同分布的随机数。可以按照所产生的随机数，对落在不同区间的数据进行统计，从而画出所产生的随机数的统计特性。MATLAB工具提供了非常丰富的函数来产生经常使用的分布的随机数，比如rand,randn就是用来产生均匀分布随机数和高斯分布随机数的。

本实验充分利用MATLAB提供的工具来产生随机数，验证和观察其统计特性。

1. disttool ：分布函数和密度函数的可视化工具

分布函数和密度函数的工具能够产生22种常用分布的概率分布曲线和概率密度曲线，并通过图形方式显示。我们还可以通过修改参数产生同一种分布不同参数的概率分布曲线和概率密度曲线。

    在MATLAB工作环境下，在命令窗口键入disttool (distribution tool的缩写)，弹出图形窗口，在函数类型中选择ＰＤＦ（概率密度曲线）和ＣＤＦ（概率分布曲线）后，我们可以通过修改参数观察不同参数情况下的曲线。这些分布包括正态（Ｎormal）、二项（Binomial）、指数（Exponential）、均匀（Uniform）、泊松（Poisson）等。也可以选择cdf选项对分布函数进行观察。

2.randtool（随机变量模拟工具）

随机变量模拟工具能够模拟产生22种常用分布的随机数，并可以通过修改它们的参数产生同一种分布不同参数的随机数，并通过图形方式显示它们的概率密度统计。在MATLAB工作环境下，在命令窗口键入randtool并回车，将弹出随机变量模拟工具窗口。Distribution(分布)： 在此可选择多种分布。Samples （样本）：模拟的样点数，可选择。Resample（重新抽样）：刷新结果，重新模拟。Export（输出）: 输出结果到工作空间，输出结果后可直接在工作空间中查看数值。

3.dfittool：数据的分布拟合的可视化工具

    dfittool （distribution fit tool的缩写），其主要步骤如下：（1）生成数据集：点击Data弹出子菜单，子菜单中点击Data à 选择数据变量（这里选前面已输出的normrv）à Data set name 中输入数据集名à 点击Create Data Set 建立数据集 à close关闭子菜单。 （2）选择分布密度拟合：点击New Fit 弹出子菜单 à 在fit name输入拟合名（此处输入了“正态拟合”）à Distribution提供了多个分布的菜单选项，此处选中normal à 点击Apply进行拟合并观察效果 à 满意后点击close关闭子菜单。

三、实验内容与步骤

1.绘制正态分布密度函数曲线。

clear all
close all
x=-10:0.1:10;
u=0;c12=4;
c11=sqrt(c12);
pf1=1/(sqrt(2*pi)*c11)*exp((-(x-u).^2)/(2*c11.^2));
u=0;c22=1;
c21=sqrt(c22);
pf2=1/(sqrt(2*pi)*c21)*exp((-(x-u).^2)/(2*c21.^2));
u=0;c32=1/4;
c31=sqrt(c32);
pf3=1/(sqrt(2*pi)*c31)*exp((-(x-u).^2)/(2*c31.^2));
plot(x,pf1,'r-',x,pf2,'b-',x,pf3,'g-');
legend('C^2=4','C^2=1','C^2=0.25');
title('正态分布概率密度函数曲线');
xlabel('x');
ylabel('f(x)');
grid on;

2. disttool

 

3. randtool

 

4. dfittool