%E I服从正态分布，P服从极值I型分布
clear all;clc;
N=2000000;
n=0;
E=normrnd(2*10^7,0.5*10^7,N,1);
I=normrnd(10^(-4),0.2*10^(-4),N,1);
uP=4;dP=1;
aEv=sqrt(6)*dP/pi;uEv=0.5772*aEv-uP;
y=rand(N,1);%产生[0,1]区间均匀分布随机数
for i=1:N
P=-evinv(1-y(i,1),uEv,aEv);
g=E(i,1)*I(i,1)-78.12*P;
if g<0
n=n+1;
end
end
pf=n/N;
beta=-norminv(pf,0,1);
fprintf('pf=%f,beta=%f',pf,beta)
评论2