概率论实验 02 - | 基于Matlab随机数的产生和特性统计

最新推荐文章于 2023-10-07 10:26:55 发布

嵌入式逍遥

最新推荐文章于 2023-10-07 10:26:55 发布

阅读量1k

点赞数 2

分类专栏：电子科学与技术文章标签： matlab 概率论随机数数据可视化

本文链接：https://blog.csdn.net/Neutionwei/article/details/117932027

版权

电子科学与技术专栏收录该内容

65 篇文章 189 订阅

订阅专栏

该实验旨在通过MATLAB编程生成不同分布的随机数，包括均匀分布和正态分布，并通过统计分析验证其概率特性。实验内容包括绘制随机数分布图、计算均值与方差、建立直方图以及对比理论概率密度函数。通过实验，学生能够理解并掌握随机数生成方法及其概率特性验证过程。

摘要由CSDN通过智能技术生成

一、实验目的

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

2．了解随机数的概率特性；

3．产生实际随机数进行验证。

二、实验原理

随机数的产生有好多方法，可以利用乘积法和同余法产生[0,1]之间的均匀分布，然后利用函数变换法产生所需不同分布的随机数。可以按照所产生的随机数，对落在不同区间的数据进行统计，从而画出所产生随机数的统计特性。所有这些工作我们可以自己动手用MATLAB，VC或VB等语言进行编程实现。

三、实验内容与步骤

close all
clear all
n=5000;
m=100;
xmin=0,xmax=1;
x=rand(1,n);
subplot(3,1,1);
plot(x,'.');
cc=sprintf('Uniform Distribution,n=%d',n);
title(cc);
subplot(3,1,2);
m1=mean(x);
c2=var(x);
x1=xmin:(xmax-xmin)/m:xmax;
h=hist(x,x1);
bar(x1,h);
grid on
cc=sprintf('Uniform Distribution Histogram n=%d',n);
title(cc);
cc=sprintf('mean=%5.2f',m1);
gtext(cc);
cc=sprintf('C^2=%5.3f',c2);
gtext(cc);
subplot(3,1,3);
pdf1=h/n/((xmax-xmin)/m);
plot(x1,pdf1);
grid on

close all
clear all
n=5000;
m=100;
xmin=-5,xmax=5;
x=randn(1,n);
subplot(3,1,1);
plot(x,'.');
cc=sprintf('Uniform Distribution,n=%d',n);
title(cc);
subplot(3,1,2);
m1=mean(x);
c2=var(x);
x1=xmin:(xmax-xmin)/m:xmax;
h=hist(x,x1);
bar(x1,h);
grid on
cc=sprintf('Uniform Distribution Histogram n=%d',n);
title(cc);
cc=sprintf('mean=%5.2f',m1);
gtext(cc);
cc=sprintf('C^2=%5.3f',c2);
gtext(cc);
subplot(3,1,3);
pdf1=h/n/((xmax-xmin)/m);
plot(x1,pdf1);
grid on

close all
clear all
n=5000;
m=100;
xmin=-5,xmax=5;
x=rand(12,n);
x=sum(x)-6;
subplot(3,1,1);
plot(x,'.');
cc=sprintf('Uniform Distribution,n=%d',n);
title(cc);
subplot(3,1,2);
m1=mean(x);
c2=var(x);
x1=xmin:(xmax-xmin)/m:xmax;
h=hist(x,x1);
bar(x1,h);
grid on
cc=sprintf('Uniform Distribution Histogram n=%d',n);
title(cc);
cc=sprintf('mean=%5.2f',m1);
gtext(cc);
cc=sprintf('C^2=%5.3f',c2);
gtext(cc);
subplot(3,1,3);
pdf1=h/n/((xmax-xmin)/m);
plot(x1,pdf1);
grid on

close all
clear all
n=5000;
m=100;
xmin=-5,xmax=5;
u1=rand(1,n);
u2=rand(1,n)
x=sqrt(-log(u1)).*cos(2*pi*u2);
subplot(3,1,1);
plot(x,'.');
cc=sprintf('Uniform Distribution,n=%d',n);
title(cc);
subplot(3,1,2);
m1=mean(x);
c2=var(x);
x1=xmin:(xmax-xmin)/m:xmax;
h=hist(x,x1);
bar(x1,h);
grid on
cc=sprintf('Uniform Distribution Histogram n=%d',n);
title(cc);
cc=sprintf('mean=%5.2f',m1);
gtext(cc);
cc=sprintf('C^2=%5.3f',c2);
gtext(cc);
subplot(3,1,3);
pdf1=h/n/((xmax-xmin)/m);
plot(x1,pdf1);
grid on

close all
clear all
n=5000;
m=100;
xmin=-3,xmax=8;
y=randn(1,n);
x=y*1+2;
subplot(3,1,1);
plot(x,'.');
cc=sprintf('Uniform Distribution,n=%d',n);
title(cc);
subplot(3,1,2);
m1=mean(x);
c2=var(x);
x1=xmin:(xmax-xmin)/m:xmax;
h=hist(x,x1);
bar(x1,h);
grid on
cc=sprintf('Uniform Distribution Histogram n=%d',n);
title(cc);
cc=sprintf('mean=%5.2f',m1);
gtext(cc);
cc=sprintf('C^2=%5.3f',c2);
gtext(cc);
subplot(3,1,3);
pdf1=h/n/((xmax-xmin)/m);
plot(x1,pdf1);
grid on