help-octave
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Real-data DFT in Octave


From: Sergei Steshenko
Subject: Re: Real-data DFT in Octave
Date: Sat, 6 Oct 2012 09:42:18 -0700 (PDT)







>________________________________
> From: Doug Stewart <address@hidden>
>To: Sergei Steshenko <address@hidden> 
>Cc: Torbjörn Rathsman <address@hidden>; "address@hidden" <address@hidden> 
>Sent: Saturday, October 6, 2012 3:46 PM
>Subject: Re: Real-data DFT in Octave
> 
>
>
>
>
>On Sat, Oct 6, 2012 at 9:35 AM, Sergei Steshenko <address@hidden> wrote:
>
>
>>
>>
>>
>>----- Original Message -----
>>> From: Torbjörn Rathsman <address@hidden>
>>> To: address@hidden
>>> Cc:
>>> Sent: Saturday, October 6, 2012 3:18 PM
>>> Subject: Real-data DFT in Octave
>>>
>>> Is it possible to without much effort call the real to complex dft routins 
>>> from
>>> octave. It is easier to manipulate the spectrum that way, since i am 
>>> guarantied
>>> that the mirror part becomes correct. The purpose in this case is to compute
>>> filtered derivative of second order.
>>> _______________________________________________
>>> Help-octave mailing list
>>> address@hidden
>>> https://mailman.cae.wisc.edu/listinfo/help-octave
>>
>>
>>I am not sure what the problem is, but dealing with real input data I use 
>>'conj' and 'fliplr' in the end when I need to produce mirrored spoectrum for 
>>inverse FFT.
>>
>>Regards,
>>  Sergei.
>>
>>
>>_______________________________________________
>>Help-octave mailing list
>>address@hidden
>>https://mailman.cae.wisc.edu/listinfo/help-octave
>>
>fft works for me
>
>
>a=randn(100,1);
>b=fft(a);
>c=ifft(b);
>d=real(c);
>e=a-d;
>e(1:10)
>
>
>this gives close to zero,which is good.
>
>
>
>
>
>-- 
>
>DAS
>
>https://linuxcounter.net/user/206392.html
>
>

I didn't mean 'fft', 'ifft' were not working.

Your example does not include application of filter. Application of filter 
makes sense only on left half of the spectrum since the right half is complex 
conjugate of the left half.

Regards,
  Sergei.



>


reply via email to

[Prev in Thread] Current Thread [Next in Thread]