2012-07-25 70 views
2

的元素我有一个1xm单元阵列A{},与所述阵列是NxN矩阵和矩阵W(N1,m)的每个元素。
我需要计算乘法细胞与基质Matlab的

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}|

是否有这样做不使用循环的方法吗?

+2

什么'M'?是不是只是'2'?是'N1'和'N'的意思是相同的变量?我认为你正在使用'N'在此表示两回事 – Drodbar 2012-07-25 11:31:16

+1

不知道我理解的数据结构,你能发布一些代码,定义了两个元素? – Thor 2012-07-25 11:43:56

回答

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{1}S{N1}单元阵列S

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

+0

你还有一个'for'环... – 2012-07-25 18:19:27

+0

@EitanT:双回路真的,如果算上CELLFUN :)我只是想,我提出了一个可行的解决方案。这样OP可以确认它是否给出预期的输出。 – Amro 2012-07-25 18:22:05

+0

我认为for循环比'arrayfun'或'cellfun'慢,所以它不是相同的循环。但是是啊:) – 2012-07-25 18:24:19