How can I do a Monte Carlo simulation for 1000 runs on this operation?
3 ビュー (過去 30 日間)
古いコメントを表示
Please can someone help to verify that this is the proper way to run a monte carlo simulation on this code?
% A = T(:,1); %Column vector
A = [71.213; 74.499; 79.175; 54.163; 83.008; 52.615];
B = length(A); %size of column (number of elements)
t = 0; %Variable to store harvest time
number_of_runs = 1000;
for n = 1:number_of_runs
for i = 1:B
H = A(i);
while H > 0
H = H - randi([5 10]);
t = t + randi([2 4]);
end
end
end
t_mean = t/number_of_runs;
disp(t_mean);
2 件のコメント
John D'Errico
2023 年 3 月 8 日
編集済み: John D'Errico
2023 年 3 月 9 日
Hard to say. That is the proper way to simulate what you simulated. Is it the proper way to simulate what you think you were simulating? That is impossible to tell, since we see only the code you wrote. So perhaps you need to explain what is the goal.
採用された回答
Torsten
2023 年 3 月 9 日
Maybe you mean:
% A = T(:,1); %Column vector
A = [71.213; 74.499; 79.175; 54.163; 83.008; 52.615];
B = length(A); %size of column (number of elements)
number_of_runs = 1000;
T = zeros(number_of_runs,1); %Variable to store harvest time
for n = 1:number_of_runs
t = 0;
for i = 1:B
H = A(i);
while H > 0
H = H - randi([5 10]);
t = t + randi([2 4]);
end
end
T(n) = t;
end
plot(1:number_of_runs,T)
t_mean = sum(T)/number_of_runs;
disp(t_mean);
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Monte-Carlo についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!