这是一个简单的问题,我找不到更好的答案,也许其他人可以!这是代码:
示例变量
nSim = 3000;
nRow = 10000;
data = zeros(1, 5, nRow);
data (:, 1:4, :) = rand(4, nRow)*0.5; % 4 columns of duration values
data (:, 5, :) = 1000; % 1 column of actual value
basis.increaseRate = 1 + (rand(nSim, 4)*0.1);
示例计算
dataWithSim = repmat(data(:, 1:4,:),nSim, 1, 1);
increaseFactors = bsxfun(@power, basis.increaseRate, dataWithSim);
Values = bsxfun(@times, data(:,5,:), prod(increaseFactors,2));
需要 repmat 感觉不对,但我找不到避免它的方法。
实际上我正在增加^数据,我真的不想遍历两个维度(sims
或data
行)。可以以我选择的任何方式对虚拟数据进行排序,但输出的值必须是nSim
按nRow
矩阵。
欢迎任何想法。谢谢。