How to do failure probability estimation using the direct Monte Carlo simulation. my code does not run it gives me "untitled'' as ans

Question

Ashwini Subhash 2022-11-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1854218-how-to-do-failure-probability-estimation-using-the-direct-monte-carlo-simulation-my-code-does-not

评论： Ashwini Subhash 2022-11-17

%USE MONTE CARLO METHOD TO CALCULATE PROBABILITY OF FAILURE OF WIND TURBINE

%The given limit state function for the given beam is

%g = -3.689 + 0.0184*D^2 + 0.2938*V ;

%Given :

D = 90.0219343 ;

V = 6.672408424 ;

mean_D = 90 ;

mean_V = 6.6724 ;

sigma_D = 10.29023 ;

sigma_V = 0.67708 ;

%Now use normrnd function

D = normrnd(mean_D, sigma_D, [n,1]);

V = normrnd(mean_V, sigma_V, [n,1]);

%

g = zeros(n,1);

K = 0;

%Create for loop

for i=1:n

g(i,1) = -3.689 + 0.0184*(D(i,1)^2) + 0.2938*(V(i,1)) ;

if g(i,1)<0

K = K + 1;

end

P_fail = K/n;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2022-11-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1854218-how-to-do-failure-probability-estimation-using-the-direct-monte-carlo-simulation-my-code-does-not#answer_1103013

在 MATLAB Online 中打开

%USE MONTE CARLO METHOD TO CALCULATE PROBABILITY OF FAILURE OF WIND TURBINE
%The given limit state function for the given beam is 
%g = -3.689 + 0.0184*D^2 + 0.2938*V ;
%Given :
D = 90.0219343 ;
V = 6.672408424 ;
mean_D = 90 ;
mean_V = 6.6724 ;
sigma_D = 10.29023 ;
sigma_V = 0.67708 ;
%Now use normrnd function
n = 1e6;
D = normrnd(mean_D, sigma_D, [n,1]);
V = normrnd(mean_V, sigma_V, [n,1]);
%
%Create for loop
K = nnz(-3.689 + 0.0184*D.^2 + 0.2938*V < 0);
P_fail = K/n
P_fail = 0