[交流] 【转帖】使用matlab对 某正态函数进行 蒙特卡罗模拟

function func_baidu_monte_carlo

randn('state', sum(100*clock));
%利用时钟设置随机种子，这样每次产生的随机数就不同了

N = 1000000;
% 设置样本个数

m = 0.43 + sqrt(0.008) * randn(N,1);
n = 0.19 + sqrt(0.008) * randn(N,1);
L = 0.5 + sqrt(0.008) * randn(N,1);
%随机生成m, n, L，根据他们各自的分布

g = 16.4445 .* (m.*m.*L./(m.*m+n.*n))-8.2;
%计算g

fprintf('the probability of g>0 is %.2f\n', sum(g>0)/length(g));
%统计 g>0的概率，结果为0。18

已经测试通过，可以运行，结果为
the probability of g>0 is 0.18