monte.m

来自「应用与系统开发设计中的统计运算 蒙特卡洛算法」· M 代码 · 共 13 行

clear;clc
N=input('输人模拟次数N=');
t=cputime; %计时起点
p=rand(1,3*N);%产生随机数
x1=norminv(p(1:N),551.4,44.1);
x2=norminv(p(N+1:2*N),0.70,0.05);
x3=norminv(p(2*N+1:3*N),0.85,0.09);
S=x1.*x2.*x3 ;
E=sum(S)/N %求均值
sigma=std(S) %求标准差
sort(S) ; %排序
hist(S,30) %画直方图
t=cputime-t %计时终点