蒙特卡洛方法模拟小例子.docxVIP

蒙特卡洛方法模拟小例子 蒙特卡洛方法模拟小例子 例 在我方某前沿防守地域，敌人以一个炮排（含两门火炮）为单位对我方进行干扰和破坏．为躲避我方打击，敌方对其阵地进行了伪装并经常变换射击地点． 经过长期观察发现，我方指挥所对敌方目标的指示有50％是准确的，而我方火力单位，在指示正确时，有1/3的射击效果能毁伤敌人一门火炮，有1/6的射击效果能全部毁伤敌人火炮． 现在希望能用某种方式把我方将要对敌人实施的20次打击结果显现出来，确定有效射击的比率及毁伤敌方火炮的平均值。 使用蒙特卡洛方法模拟50次打击结果： function [out1 out2 out3 out4]=Msc(N) % N开炮次数 % out1射中概率 % out2平均每次击中次数 % out3击中敌人一门火炮的射击总数 % out4击中敌人2门火炮的射击总数 for i=1:N x0=randperm(2)-1; y0=x0(1); if y0==1 fprintf(第%d次：指示正确||,i); x1=randperm(6); y1=x1(1); if y1==1|y1==2|y1==3 fprintf(第%d次：击中0炮||,i); k1=k1+1; elseif y1==4|y1==5 fprintf(第%d次：击中1炮||,i); k2=k2+1; fprintf(第%d次：击中2炮||,i); k3=k3+1; fprintf(第%d次：指示错误，击中0炮||,i); fprintf(\n); out1=(k2+k3)/N; out2=(0*k1+k2+2*k3)/20; out3=k2/N; out4=k3/N; 1. [out1 out2 out3 out4]=Msc(50) 结果： 1. 第1次：指示正确||第1次：击中2炮|| 2. 第2次：指示错误，击中0炮|| 3. 第3次：指示错误，击中0炮|| 4. 第4次：指示正确||第4次：击中0炮|| 5. 第5次：指示错误，击中0炮|| 6. 第6次：指示正确||第6次：击中1炮|| 7. 第7次：指示正确||第7次：击中0炮|| 8. 第8次：指示错误，击中0炮|| 9. 第9次：指示正确||第9次：击中2炮|| 10. 第10次：指示正确||第10次：击中1炮|| 11. 第11次：指示正确||第11次：击中1炮|| 12. 第12次：指示正确||第12次：击中2炮|| 13. 第13次：指示错误，击中0炮|| 14. 第14次：指示正确||第14次：击中1炮|| 15. 第15次：指示错误，击中0炮|| 16. 第16次：指示错误，击中0炮|| 17. 第17次：指示正确||第17次：击中0炮|| 18. 第18次：指示错误，击中0炮|| 19. 第19次：指示正确||第19次：击中1炮|| 20. 第20次：指示错误，击中0炮|| 21. 第21次：指示正确||第21次：击中0炮|| 22. 第22次：指示正确||第22次：击中1炮|| 23. 第23次：指示正确||第23次：击中0炮|| 24. 第24次：指示错误，击中0炮|| 25. 第25次：指示正确||第25次：击中1炮|| 26. 第26次：指示错误，击中0炮|| 27. 第27次：指示正确||第27次：击中1炮|| 28. 第28次：指示正确||第28次：击中0炮|| 29. 第29次：指示正确||第29次：击中0炮|| 30. 第30次：指示正确||第30次：击中0炮|| 31. 第31次：指示错误，击中0炮|| 32. 第32次：指示错误，击中0炮|| 33. 第33次：指示正确||第33次：击中0炮|| 34. 第34次：指示错误，击中0炮|| 35. 第35次：指示正确||第35次：击中0炮|| 36. 第36次：指示正确||第36次：击中0炮|| 37. 第37次：指示错误，击中0炮|| 38. 第38次：指示正确||第38次：击中0炮|| 39. 第39次：指示错误，击中0炮|| 40. 第40次：指示正确||第40次：击中0炮||