2

这是一个简单的问题，我找不到更好的答案，也许其他人可以！这是代码：

示例变量

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矩阵。

欢迎任何想法。谢谢。

4

1 回答 1

3

你真的不需要使用它repmat。您可以像这样直接输入该“子矩阵” data-

increaseFactors = bsxfun(@power, basis.increaseRate, data(:,1:4,:));

bsxfun内部负责维度的扩展，在这种情况下singleton，这是第一个维度（行） 。data由于basis.increaseRate有nSim行并且data(:,1,4,:)只有一行，它将被扩展为具有相同数量的行，即nSim，因此在内部完成了repmat-ing/expanding的工作。

其余代码保持不变。

于 2014-10-23T10:00:44.380 回答