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 感觉不对,但我找不到避免它的方法。

实际上我正在增加^数据,我真的不想遍历两个维度(simsdata行)。可以以我选择的任何方式对虚拟数据进行排序,但输出的值必须是nSimnRow矩阵。

欢迎任何想法。谢谢。

4

1 回答 1

3

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

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

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

其余代码保持不变。

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