[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: not equivalent?
From: |
Paul Kienzle |
Subject: |
Re: not equivalent? |
Date: |
Thu, 7 Nov 2002 16:02:37 -0500 |
User-agent: |
Mutt/1.2.5.1i |
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
-------------------------------------------------------------
- not equivalent?, Albert F. Niessner, 2002/11/07
- Re: not equivalent?,
Paul Kienzle <=