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!