How to find beta and Pf values using Monte Carlo

1 次查看(过去 30 天)
%USE MONTE CARLO METHOD TO CALCULATE PROBABILITY OF FAILURE OF WIND TURBINE
%The given limit state function for the given 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 = 10e6;
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

回答(1 个)

Aditya
Aditya 2023-8-31
Hi Ashwini,
From the code that you have shared I can see that you have already calculated the Pf (probability of failure) value. To calculate the beta value, you can utilize the "norminv" function in MATLAB. 
For more detailed information on how to use the "norminv" function, please refer to the documentation available at: https://in.mathworks.com/help/stats/norminv.html
Hope this helps!

类别

Help CenterFile Exchange 中查找有关 Spline Postprocessing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by