我有一个函数Y(k)的矩阵Bi,我计算如下:
function Y = specm(A,k)
% A is a 3-D array, and k is a number from 1 to p( p is the number of matrices in A)
B = eye(size(A,2));
p = size(A,3);
if k > 0;
B = B' * A(:,:,k);
end
for l = 1:(p-k);
B = B + A(:,:,l)' * A(:,:,k+l);
end
if k > 0;
Y = 2*B;
end
end
问题是我希望从1到p的所有k的Y之和。有谁知道怎么做?
提前多多感谢。