[Octave-bug-tracker] [bug #39566] Suspect interpft is wrong regarding treatment of Nyquist frequency

[Alain Cochard]

Follow-up Comment #3, bug #39566 (project octave):

Mike Miller writes:
 > Update of bug #39566 (project octave):
 > 
 >                   Status:               Confirmed => In Progress           

 >              Assigned to:                    None => mtmiller              

 > 
 >     _______________________________________________________
 > 
 > Follow-up Comment #2:
 > 
 > I've got a patch for both problems you described, please test the
 > attached if you are able to and let me know if this does what you
 > expect.

Yes, this does what I expect.


 > Also if you have access to Matlab, which it seems you do, it would
 > be helpful to post some examples of interpft on complex input,
 > since the Nyquist bin will not be purely real and the FFT is not
 > conjugate symmetric. For example, what do you get with
 > 
 > 
 > > fft (interpft ((1+j)*[1:4]', 6))

With matlab, I get:

  15.0000 +15.0000i
  -6.0000 - 0.0000i
  -1.5000 - 1.5000i
        0          
  -1.5000 - 1.5000i
  -0.0000 - 6.0000i


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?39566>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/