概率论与数理统计实验报告

概率论与数理统计

实验报告

概率论部分实验二

《正态分布综合实验》

实验名称：正态分布综合实验

实验目的：通过本次实验，了解Matlab在概率与数理统计领域的应用，学会用matlab做概率密度曲线，概率分布曲线，直方图，累计百分比曲线等简单应用；同时加深对正态分布的认识，以更好得应用之。

实验内容：

实验分析：

本次实验主要需要运用一些matlab函数，如正态分布随机数发生器normrnd函数、绘制直方图函数hist函数、正态分布密度函数图形绘制函数normpdf函数、正态分布分步函数图形绘制函数normcdf等；同时，考虑到本次实验重复性明显，如，分别生成100,1000,10000个服从正态分布的随机数，进行相同的实验操作，故通过数组和循环可以简化整个实验的操作流程，因此，本次实验程序中要设置数组和循环变量。

实验过程：

1.直方图与累计百分比曲线

1）实验程序

m=[100,1000,10000]; 产生随机数的个数

n=[2,1,0.5]; 组距

for j=1:3

for k=1:3

x=normrnd(6,1,m(j),1); 生成期望为6，方差为1的m(j)个

正态分布随机数

a=min(x); a为生成随机数的最小值

b=max(x); b为生成随机数的最大值

c=(b-a)/n(k); c为按n(k)组距应该分成的组数 subplot(1,2,1); 图形窗口分两份

hist(x,c);xlabel('频数分布图'); 在第一份里绘制频数直方图

yy=hist(x,c)/1000; yy为各个分组的频率

s=[];

s(1)=yy(1);

for i=2:length(yy)

s(i)=s(i-1)+yy(i);

end s[]数组存储累计百分比

x=linspace(a,b,c);

subplot(1,2,2); 在第二个图形位置绘制累计百分

比曲线

plot(x,s,x,s);xlabel('累积百分比曲线');

grid on; 加网格

figure; 另行开辟图形窗口，为下一个循

环做准备

end

end

2）实验结论及过程截图

实验结果以图像形式展示，以下分别为产生100,1000,10000个正态分布随机数，组距分别为2,1,0.5的频数分布直方图和累积百分比曲线，从实验结果看来，随着产生随机数的数目增多，组距减小，累计直方图逐渐逼近正态分布密度函数图像，累计百分比逐渐逼近正态分布分布函数图像。

N=100,组距为1的频数分布图与累计百分比曲线

N=100,组距为0.5的频数分布图与累计百分比曲线

N=1000,组距为1的频数分布图与累计百分比曲线

N=1000,组距为0.5的频数分布图与累计百分比曲线

N=10000,组距为2的频数分布图与累计百分比曲线

N=10000,组距为1的频数分布图与累计百分比曲线

N=10000,组距为0.5的频数分布图与累计百分比曲线

1）实验程序

clear all;

x=[-0.15:0.000002:0.2]';

y1=[];y2=[];

mul=[0.05 0.05 0.05];

sigmal=[0.01 0.02 0.03];

for i=1:length(mul)

y1=[y1,normpdf(x,mul(i),sigmal(i))]; y2=[y2,normcdf(x,mul(i),sigmal(i))]; end

subplot(1,2,1);plot(x,y1);

xlabel('(a) 概率密度函数');

subplot(1,2,2);plot(x,y2);

xlabel('(b) 分布函数');

2）实验结果截图

1）实验程序

clear all;

x=[-0.15:0.000001:0.2]';

y1=[];y2=[];

mul=[0.03 0.05 0.07];

sigmal=[0.02 0.02 0.02];

for i=1:length(mul)

y1=[y1,normpdf(x,mul(i),sigmal(i))]; y2=[y2,normcdf(x,mul(i),sigmal(i))]; end

subplot(1,2,1);plot(x,y1);

xlabel('(a) 概率密度函数');

subplot(1,2,2);plot(x,y2);

xlabel('(b) 分布函数');

2）实验结果截图