Thread Subject:

Matix maths

Does anyone know how I can sum all the diagonals of a matrix and get a vector of the results?

For example if a had the following matrix

18 36 24 12
12 24 16 8
 6 12 8 4

I would get a vector of:

18 48 54 40 16 4

All help appreciated!
Dan

On 11/15/2012 2:09 PM, Daniel Robbins wrote:
> Does anyone know how I can sum all the diagonals of a matrix and get a
> vector of the results?
> For example if a had the following matrix
> 18 36 24 12
> 12 24 16 8
> 6 12 8 4
>
> I would get a vector of:
>
> 18 48 54 40 16 4
>

Not terribly clever, but...given your array as 'm'--

 >> n=fliplr(m)
n =
     12 24 36 18
      8 16 24 12
      4 8 12 6
 >> [nr,nc]=size(n);d=nc-1;for i=1:nr+nc-1,sum(diag(n,d)),d=d-1;end
ans =
     18
ans =
     48
ans =
     54
ans =
     40
ans =
     16
ans =
      4
 >>

--

On 11/15/2012 2:09 PM, Daniel Robbins wrote:
> Does anyone know how I can sum all the diagonals of a matrix and get a vector of the results?
>
> For example if a had the following matrix
>
> 18 36 24 12
> 12 24 16 8
> 6 12 8 4
>
> I would get a vector of:
>
> 18 48 54 40 16 4
>
> All help appreciated!
> Dan
>

completley not clear question. Which diagonals? positive slope or
negative slope? which sub matrices are you summing its diagonals?

I can see where the '18' and the '48' and the '54' came from. But
the rest is not clear now how you got it.

"Daniel Robbins" <d.w.e.robbins@gmail.com> wrote in message <k83i5d$rqg$1@newscl01ah.mathworks.com>...
> Does anyone know how I can sum all the diagonals of a matrix and get a vector of the results?
> For example if a had the following matrix
>
> 18 36 24 12
> 12 24 16 8
> 6 12 8 4
>
> I would get a vector of:
>
> 18 48 54 40 16 4
- - - - - - - - -
 Let A be your matrix. Then v will be your desired vector:

 [m,n] = size(A);
 B = [zeros(m,m-1),A,zeros(m,m-1)];
 v = sum(B(bsxfun(@plus,(0:m-1:(m-1)^2)',m:m:m*(m+n-1))),1);

Roger Stafford

A=[18 36 24 12
12 24 16 8
6 12 8 4]

[m n]=size(A);
[I J]=ndgrid(1:m,1:n);
accumarray(I(:)+J(:)-1,A(:))

% Bruno

[m n]=size(A);
K=hankel(1:m,m+(0:n-1))
accumarray(K(:),A(:))

% Bruno

Thanks to everyone who gave solutions, good to have a range of options to pick from!!