MATLAB Central Newsreader - how to vectorise formula

how to vectorise formula

Richard — Sun, 08 Aug 2010 22:52:05 +0000

Hi, having difficulty working out how to vectorise this formula

for j=1:x
sum over h=1:n of U(h)*G(h,j) * sum over t=1:m of (X(h,t)*P(t))/G(h,t)

Any help much appreciated.
thanks.

Re: how to vectorise formula

Roger Stafford — Sun, 08 Aug 2010 23:51:05 +0000

"Richard " <REMOVETHISrcaldwellie@yahoo.com> wrote in message <i3ncel$s9b$1@fred.mathworks.com>...
> Hi, having difficulty working out how to vectorise this formula
>
> for j=1:x
> sum over h=1:n of U(h)*G(h,j) * sum over t=1:m of (X(h,t)*P(t))/G(h,t)
>
> Any help much appreciated.
> thanks.
- - - - - - - - - - - -
Your expression is somewhat difficult to interpret but does this do what you want?

((U(:).*((X./G)*P(:))).')*G

Roger Stafford

Re: how to vectorise formula

Roger Stafford — Mon, 09 Aug 2010 00:48:04 +0000

"Roger Stafford" <ellieandrogerxyzzy@mindspring.com.invalid> wrote in message <i3nft9$sc0$1@fred.mathworks.com>...
> Your expression is somewhat difficult to interpret but does this do what you want?
>
> ((U(:).*((X./G)*P(:))).')*G
>
> Roger Stafford
- - - - - - - - - - -
I just realized that since the same G matrix appears twice in the above expression that implies that x is equal to m. If that is not so, you would have to modify the formula, perhaps something like this:

((U(:).*((X./G(:,1:m))*P(:))).')*G(:,1:t)

Roger Stafford

Re: how to vectorise formula

Richard — Mon, 09 Aug 2010 09:19:07 +0000

"Roger Stafford" <ellieandrogerxyzzy@mindspring.com.invalid> wrote in message <i3nj84$n2m$1@fred.mathworks.com>...
> "Roger Stafford" <ellieandrogerxyzzy@mindspring.com.invalid> wrote in message <i3nft9$sc0$1@fred.mathworks.com>...
> > Your expression is somewhat difficult to interpret but does this do what you want?
> >
> > ((U(:).*((X./G)*P(:))).')*G
> >
> > Roger Stafford
> - - - - - - - - - - -
> I just realized that since the same G matrix appears twice in the above expression that implies that x is equal to m. If that is not so, you would have to modify the formula, perhaps something like this:
>
> ((U(:).*((X./G(:,1:m))*P(:))).')*G(:,1:t)
>
> Roger Stafford
Hi Roger, thanks. Although the G matrix appears twice, the first time its a specific column vector from G and the second time its the complete matrix. I 'think' I've vectorised it as
(U.*G(:,j))'*(P'*(X(:,:)./G(:,:))')'
the amount of columns in G = t