小编dev*_*raf的帖子

我正在尝试生成 100 个场景 (X, Y) 的样本，其中 X 和 Y 均呈正态分布 X=N(50,5^2)、Y=N(30,2^2) 并且 X 和 Y 相关Cov(X,Y)=0.4。

我已经能够使用 Cholesky 分解生成 100 个场景：

# We do a Cholesky decomposition to generate correlated scenarios\nnScenarios = 10\n\n\xce\xa3 = [25 0.4; 0.4 4]\n\xce\xbc = [50, 30]\nL = cholesky(\xce\xa3)\nv = [rand(Normal(0, 1), nScenarios), rand(Normal(0, 1), nScenarios)]\n\nX = reshape(zeros(nScenarios),1,nScenarios)\nY = reshape(zeros(nScenarios),1,nScenarios)\n\nfor i = 1:nScenarios\n    X[1, i] = sum(L.U[1, j] *v[j][i] for j = 1:nBreadTypes) + \xce\xbc[1]\n    Y[1, i] = sum(L.U[2, j] *v[j][i] for j = …