How to do failure probability estimation using the direct Monte Carlo simulation. my code does not run it gives me "untitled'' as ans
显示 更早的评论
%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
end
P_fail = K/n;
回答（1 个）
%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
4 个评论
Ashwini Subhash
2022-11-17
Ashwini Subhash
2022-11-17
Torsten
2022-11-17
first order reliability method ?
compare beta function and pf of FORM ?
Ashwini Subhash
2022-11-17
类别
在 帮助中心 和 File Exchange 中查找有关 Wind Power 的更多信息
2022-11-17
2022-11-17
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
