我正在尝试生成以下模型的蒙特卡罗样本:
Y = 20 / SQRT(X1)* 10 ^((X2-X3)/ 20)* X4 * 10 ^((X 5 -X 6)/ 20)。
在Matlab中,我使用以下代码为k变量生成N大小的随机样本。假设输入变量均匀分布在(0; 1)
中在(0; 1)
中生成均匀样品的%for i=1:N
X=[];
T = rand(1,2*k);
%PREPARATION OF THE SAMPLE MATRIX X
A=T(1:k);
B=T(k+1:2*k);
X=[A;B];
for j=1:k;
Ab=A;
Ab(j)=B(j);
X=[X;Ab];
end
end
我有两个问题: 1)如何定义模型y = f(x1,..,x6)? 2)如果我想将参数的pdf从uniform更改为例如,该怎么办?日志均匀?
感谢您的帮助。