最初的问题是模拟一个24/7全天候使用的灯泡,并且通常会持续25天。一盒灯泡包含12.盒子可能持续一年以上的概率是多少?使用Matlab查找高斯分布的概率
我不得不使用MATLAB来建模一个基于指数变量的高斯曲线。 下面的代码生成一个平均值= 300和std = sqrt(12)* 25的高斯模型。 我不得不使用如此多的不同变量并加上它们的原因是因为我应该证明中心极限定理。高斯曲线表示一盒灯泡持续一天的概率,其中300是盒子将持续的平均天数。
我在使用生成的高斯函数和查找天> 365的概率时遇到了问题。声明1-normcdf(365,300,sqrt(12)* 25)是试图找出概率的预期值,我得到了.2265。 有关如何根据生成的高斯I来查找天> 365的概率的任何提示将不胜感激。
谢谢!
clear all
samp_num=10000000;
param=1/25;
a=-log(rand(1,samp_num))/param;
b=-log(rand(1,samp_num))/param;
c=-log(rand(1,samp_num))/param;
d=-log(rand(1,samp_num))/param;
e=-log(rand(1,samp_num))/param;
f=-log(rand(1,samp_num))/param;
g=-log(rand(1,samp_num))/param;
h=-log(rand(1,samp_num))/param;
i=-log(rand(1,samp_num))/param;
j=-log(rand(1,samp_num))/param;
k=-log(rand(1,samp_num))/param;
l=-log(rand(1,samp_num))/param;
x=a+b+c+d+e+f+g+h+i+j+k+l;
mean_x=mean(x);
std_x=std(x);
bin_sizex=.01*10/param;
binsx=[0:bin_sizex:800];
u=hist(x,binsx);
u1=u/samp_num;
1-normcdf(365,300, sqrt(12)*25)
bar(binsx,u1)
legend(['mean=',num2str(mean_x),'std=',num2str(std_x)]);
你的高斯代表什么?这是一个灯泡将持续多少小时?剩余在一个盒子里或其他东西的灯泡数量?请更新您的信息。 – kkuilla 2014-10-30 10:22:05
我相信高斯表示一盒灯泡持续某些天数的概率。 300通常会持续多长时间,平均值。 – azumakazuma 2014-10-30 10:31:11
为什么你的'std = sqrt(12)* 25'?你确定答案.2265是错误的吗? – kkuilla 2014-10-30 11:14:04