Re: not equivalent?

[Albert F. Niessner]

Makes perfect sense now. Never knew that transpose did different things.

Thank you very much.

Al

On Thu, 2002-11-07 at 16:02, Paul Kienzle wrote:
> On Thu, Nov 07, 2002 at 03:46:15PM -0500, Albert F. Niessner wrote:
> > 
> > I have two vectors (series and spacing) of 500x1 elements each. If I
> > then do the following (N = length(series);):
> > 
> >   k=0:N-1;
> >   W = e.^(-j*2*pi*(spacing*k)/N);
> >   spectrum = W' * series;
>              ^^
> This is the complex conjugate.
> 
> > 
> > I get some answer. Because of the size of W, it may be worth while to
> > loop over very large vectors so I can also do it this way:
> > 
> >   w = -j*2*pi*spacing/N;
> >   for k=1:N
> >     W = e.^(w*(k-1));
> >     spectrum(k) = sum (series .* W);
>                                    ^
> 
> This is not.
> 
> >   endfor;
> 
> Use W.' to transpose without conjugate.
> 
> Does this help?
> 
> Paul Kienzle
> address@hidden
> 
> 
> 
> -------------------------------------------------------------
> Octave is freely available under the terms of the GNU GPL.
> 
> Octave's home on the web:  http://www.octave.org
> How to fund new projects:  http://www.octave.org/funding.html
> Subscription information:  http://www.octave.org/archive.html
> -------------------------------------------------------------
> 




-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------