2009年4月17日星期五

problem

最近想编一个用蒙特卡罗法计算土石坝漫顶风险的程序,琢磨一阵子后发现问题很多。problem 1:在产生正态分布的随机数时,竟发现书上的算例不知道是怎么算出来的,复核了一下貌似不太对啊,寒,于是给作者DOC.WU发了个短信,约个时间探讨下,嘿嘿,其实是有很多东西不知道怎么下手,想请这位可靠度的专家给指导下。。。
problem 2 更是离谱了,按照这种算法,计算次数应该在100/pf 以上,如果取pf为 E-6 量级的话,就得计算一亿次,o,mygod。一亿次是什么概念啊,师兄说没事,慢慢算就是了,寒寒寒寒。。。
看来得找找方法改进下,否则我的电脑非废了不可。。。