lnX~N(mu,sig2);
X=exp(lnX) 即服从对数正态分布;
所以先生成二元正态分布,之后指数化即得到二元对数正态分布;
生成多元正态分布的内容,我在这篇帖子里写过
https://bbs.pinggu.org/thread-2157211-1-1.html
这里需要注意的是前后均值与协方差的变化;
假设生成的二元正态分布:
{X1 X2} 的均值为[mu1 mu2]
协方差阵为[var1 cov
cov var2];
指数化之后,得到
Y1=exp(X1) Y2=exp(X2)
{Y1 Y2} 服从二元对数正态分布;
此时
{Y1 Y2}的均值为[exp(mu1+0.5*var1) exp(mu2+0.5*var2)]
协方差阵为
[(exp(var1)-1)*exp(2mu1+var1) (exp(cov)-1)*exp(mu1+mu2+0.5var1+0.5var2)
(exp(cov)-1)*exp(mu1+mu2+0.5var1+0.5var2) (exp(var2)-1)*exp(2mu2+var2)]