2

我有一个 1xm 单元阵列A{},阵列的每个元素都是 NxN 矩阵和一个矩阵W(N1,m)
我需要计算

Sum(j) = W(j,1)*A{1,1} + W(j,2)*A{1,2}  

我正在做以下事情:

for j=1:N1
  sum=false(N);  
  for k=1:m  
    sum = sum + W(j,k)*A{1,k};  
  end  
  Sum(j)=sum  
end

或者更直观地: 有没有办法在不使用循环的情况下做到这一点?
Matrix W(let's say N1=2)
|W11 W12||A{1,1}| = |W11*A{1,1} + W12*A{1,2}|
|W21 W22||A{1,2}| = |W21*A{1,1} + W22*A{1,2}|

4

3 回答 3

0

要在没有 for 循环的情况下做到这一点,您可以强奸(请原谅表达式)arrayfun命令:

w_func = @(j)arrayfun(@(k)(W(j, k) * A{k}), 1:m, 'Un', 0)
sum_func = @(x)sum(cat(3, x{:}), 3)
S = arrayfun(@(j)sum_func(w_func(j)), 1:N1, 'Un', 0);

这将生成一个包含所有总和的元胞数组S,从S{1}S{N1}

于 2012-07-25T18:20:33.563 回答
0

我对你想要做什么感到困惑,但如果我理解正确,这段代码应该可以工作:

temp = cell2mat(A);
a_sum = temp*repmat(eye(n),m,1); % this reduces A by performing sum like operation so [1 1 1 3;0 1 0 2]
                                 % becomes [2 4; 0 3]

Sum = W * a_sum
于 2012-07-25T13:56:08.600 回答
0

我也不确定我是否理解了这个问题,但这里有一些代码需要考虑:

%# create some data resembling what you described
N = 2;
m = 4;
N1 = 5;
W = rand(N1,m);
A = cell(1,m); for i=1:m, A{i} = rand(N); end

%# do the multiplications
s = cell(N1,1);
for j=1:N1
    AA = cellfun(@times, A, num2cell(W(j,:)), 'UniformOutput',false);
    s{j} = sum(cat(3,AA{:}), 3);
end

元胞数组s现在包含这样的结果:

s{j} = W(j,1)*A{1} + W(j,2)*A{2} + ... + W(j,m)*A{m}

因此s是一个大小为 N1×1 的元胞数组,其中每个元胞包含一个 N×N 矩阵

于 2012-07-25T16:18:18.367 回答