看了一篇论文,后面用到蒙特卡洛仿真,但是文中没给出具体方法和相关程序,只有最终结果。谁能帮我看一下这个到底怎么解决,万分感激。
生成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的经验分布函数图和经验密度函数图。