可以采用泰勒展开式的方法:
http://en.wikipedia.org/wiki/Taylor_expansions_for_the_moments_of_functions_of_random_variables
编辑:
取,V = log (Y )。U=log(X)V=log(Y)
使用多元泰勒展开来计算的近似值(类似于链接中“ First Moment”末尾的示例,该示例处理E(X .1 / Y )的情况更简单),并使用单变量展开来计算E(U )和E(V )的近似值(在同一部分的第一部分中给出),以达到类似的精度。根据这些内容,计算(近似)协方差。E(UV)E(X.1/Y))E(U)E(V)
扩展到与链接中的示例相似的近似程度,我认为您最终得到了每个(未转换)变量的均值和方差以及它们的协方差的术语。
编辑2:
但是这里有一些技巧可以节省一些精力:
注意和X = exp (U )和Y = exp (V )。E(XY)=Cov(X,Y)+E(X)E(Y)X=exp(U)Y=exp(V)
ë(EXP(Û))≈EXP(μÚ)+EXP(μÛ)
E[f(X)]≈f(μX)+f′′(μX)2σ2X
E(exp(U))≈exp(μU)+exp(μU)2σ2U≈exp(μU+12σ2U)
编辑:最后一步是从泰勒逼近出发,这对小很好(采用)。exp(b)≈1+bbb=12σ2U
(该近似值对于,法线是精确的:))UVE(exp(U))=exp(μU+12σ2U)
令W=U+V
E(XY)=E(exp(U).exp(V))=E(exp(W))
≈exp(μW)+exp(μW)2σ2W≈exp(μW+12σ2W)
并给定,然后Var(W)=Var(U)+Var(V)+2Cov(U,V)
(编辑:)
1+Cov(X,Y)E(X)E(Y)=E(XY)E(X)E(Y)
≈exp(μW+12σ2W)exp(μU+12σ2U).exp(μV+12σ2V)
≈exp(μU+μV+12(σ2U+σ2V+2Cov(U,V)))exp(μU+12σ2U).exp(μV+12σ2V)
≈exp[Cov(U,V)]
因此。对于双变量高斯这应该是精确的。U,VCov(U,V)≈log(1+Cov(X,Y)E(X)E(Y))U,V
如果您使用第一个近似值而不是第二个近似值,则此处将获得不同的近似值。