Thread Subject: trying to avoid using a for loop

Hi.

I am trying to compute S.
S is a m by m matrix.

To create S, I use the following information.

G is a m by N matrix.
G(:,i) is the ith column of G.
g_bar is a m by 1 vector.

This for-loop can compute S.
S=[];
for i=1:N
     x=G(:,i)-g_bar;
     S=S+(x*x');
end

I am wondering if there is a way to compute S without the for-loop. Thanks.

"Ka Lee" <kaloklee@gmail.com> wrote in message <gdvpsh$26q$1@fred.mathworks.com>...
> Hi.
>
> I am trying to compute S.
> S is a m by m matrix.
>
> To create S, I use the following information.
>
> G is a m by N matrix.
> G(:,i) is the ith column of G.
> g_bar is a m by 1 vector.
>
> This for-loop can compute S.
> S=[];
> for i=1:N
> x=G(:,i)-g_bar;
> S=S+(x*x');
> end

> I am wondering if there is a way to compute S without the for-loop. Thanks.


What is the problem with the for loop?
Does it give a wrong result?
Is it not pretty?

"Ka Lee" <kaloklee@gmail.com> wrote in message <gdvpsh$26q$1@fred.mathworks.com>...
> .......
> for i=1:N
> x=G(:,i)-g_bar;
> S=S+(x*x');
> end
>
> I am wondering if there is a way to compute S without the for-loop.

 G2 = G-repmat(g_bar,1,N);
 S = G2*G2';

Roger Stafford