Latin hypercube sampling coupled with Monte Carlo simulation

3 次查看（过去 30 天）

Barry Akenga 2017-4-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/335935-latin-hypercube-sampling-coupled-with-monte-carlo-simulation

评论： Duc Thanh Tran 2022-7-20

Hi I am trying to formulate a Monte carlo simulation for corrosion initiation time and at the same time use Latin Hypercube sampling to come up with a more precise description of data ie to form a curve fitting pdf(probability density function). The relative variables I have are coefficient of variation and mean, however I am struggling to find the code for using Cov as the option provided is for matrix options. Below is some of the code I have managed to formulate: X = normrnd(40, 0.1, 1, 1000); Co = lognrnd(3.5, 0.5, 1, 1000); D = 5*(10^-5); Cth = unifrnd(0.6, 1.2, 1, 1000); E = 1/[erfc(Cth / Co)]; T = X .^ 2 / [4 * D .* E .^ 2]; a=mean(log(T)); b=std(log(T)); To=lognpdf(T,a,b); plot(T,To); grid; xlabel('T'); ylabel('To')

I am still a long way to go but I believe if I can find a code template for the simulation I can be able to progress further into the corrosion modelling.