monte carlo 仿真及其对应编程

看了一篇论文，后面用到蒙特卡洛仿真，但是文中没给出具体方法和相关程序，只有最终结果。谁能帮我看一下这个到底怎么解决，万分感激。

生成Z的模拟值的基本步骤：                              

步骤1：随机生成随机变量序列{W(r),r=0,1,2,…..55},其中W(r)~N(0,r);

步骤2：根据v(r)=exp(- )生成随机变量序列{v(r),r=0,1,2,…..55};

步骤3：生成10000个[0,1]上的均匀分布随机数，根据 =U,生成10000个成员的取整余命 的模拟样本{ ，i=1,2,…8; j=1,2,…  },上式中，T代表个体的余命，U代表[0,1]上的均匀分布随机数；

步骤4：根据式（1），式（2）得到随机变量序列{ }；

步骤5：根据式（3），生成Z的一个模拟值。

其中式（1），式（2），式（3）如下

    （ ）     （1）   ；                       （2）；      Z=          （3）

式中出现的数据如表

    

序号

    

  

     

  

  

  

     

  

  

  

     

  

    

1

    

  

0

  

  

  

0.6001

  

  

  

1100

  

    

2

    

  

0

  

  

  

0.5710

  

  

  

1250

  

    

3

    

  

0

  

  

  

0.5445

  

  

  

1300

  

    

4

    

  

0

  

  

  

0.5181

  

  

  

1400

  

    

5

    

  

25

  

  

  

0.2250

  

  

  

1000

  

    

6

    

  

20

  

  

  

0.2447

  

  

  

1100

  

    

7

    

  

10

  

  

  

0.4501

  

  

  

1500

  

    

8

    

  

10

  

  

  

0.4079

  

  

  

1350

  

步骤1~步骤5重复10000次，即可生成10000个Z值。根据10000个模拟值得到Z的经验分布函数图和经验密度函数图。

怎么回事呀？怎么乱码了呀？

不好意思，重新传一下，在这个回复的附件中。

Thanks for sharing!