我想知道是否有办法提高矩阵的功率一个 作为一个数组?
假设我们有这个矩阵
A =
5 4
3 6
然后我们重复它的形状。
>> repmat(A, 5, 1)
ans =
5 4
3 6
5 4
3 6
5 4
3 6
5 4
3 6
5 4
3 6
现在我想提高功率,所以长时间重复的矩阵如下所示:
>> [A^1; A^2; A^3; A^4; A^5]
ans =
5 4
3 6
37 44
33 48
317 412
309 420
2821 3740
2805 3756
25325 33724
25293 33756
在 MATLAB/Octave 中没有 for 循环是否可以做到这一点?
最佳答案
使用 arrayfun
的另一个选项
B = cell2mat(arrayfun(@(x)A^x,1:5,'UniformOutput',0).')
结果:
B =
5 4
3 6
37 44
33 48
317 412
309 420
2821 3740
2805 3756
25325 33724
25293 33756
但在这种情况下,基本的 for 循环可能是最快的选择。
Octave 基准测试:
tic
iif = @(varargin) varargin{2*find([varargin{1:2:end}], 1, 'first')}();
recPower = @(A, B, n, f) iif(n > 1, @() [B; f(A, A * B, n - 1, f)], true, @() B);
nPower = @(A, n) recPower(A, A, n, recPower);
for ii = 1:1000
% Calculate for arbitrary n.
nPower(A, 5);
end
toc
耗时是 1.023 秒。
tic
for ii = 1:1000
B = cell2mat(arrayfun(@(x)A^x,1:5,'UniformOutput',0).');
end
toc
耗时是 4.8684 秒。
tic
for ii = 1:1000
B=[];
for jj = 1:5
B = [B;A^jj];
end
end
toc
耗时是 0.039371 秒
关于matlab - 有没有办法在矩阵 repmat 数组上执行 mpower?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/55347294/