[Scilab-users] discrete Fourier transform

Stefan Du Rietz sdr at durietz.se
Thu Mar 21 18:46:14 CET 2013


On 2013-03-21 17:13, Serge Steer wrote:
--------------------
> Le 21/03/2013 16:41, haasejos a écrit :
>> hello,
>> for signalanalysis I would like to use discrete Fourier transform
>> (dft). To
>> see, how it works, I use the simple example below. Why is *XfA =
>> abs(Xf)*2/n* respectively why is  XfA = abs(Xf) wrong?
> Why do you say that  XfA = abs(Xf) is wrong
> Note however. It is much more efficient using fft instead of dft.
>
> Serge Steer

Maybe haasejos meant and should have used XfA == abs(Xf)

Stefan

>>
>> clear; clc; xdel;
>>
>> function y = f(x);
>>      y = sin(x);
>> endfunction;
>> n = 200;
>> x=linspace(0,2*%pi,n);
>> y=f(x);
>>
>> mat = [x',y'];
>> //disp(mat);
>> //plot2d(x , y);
>>
>> //xtitle('DATA','n''','y''');
>>
>> Xf=dft(y,-1);
>> XfA = abs(Xf)*2/n;
>> plot2d3([1:n/2]',XfA(1:n/2));
>>
>>
>>
>>
>> --
>> View this message in context:
>> http://mailinglists.scilab.org/discrete-Fourier-transform-tp4026318.html
>>
>> Sent from the Scilab users - Mailing Lists Archives mailing list
>> archive at Nabble.com.
>> _______________________________________________
>> users mailing list
>> users at lists.scilab.org
>> http://lists.scilab.org/mailman/listinfo/users
>>
>
> _______________________________________________
> users mailing list
> users at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users
>





More information about the users mailing list