[SciPy-user] Fourier series

[Ivo Maljevic]

I think if you spent some time in deriving formulas it should be possible to
find it:

1) Since the function f(t) is periodic, its spectrum is discrete anyway, so
FFT will do the job, provided that you sample the frequency at multiple of
harmonics.

2) DFT coefficients are complex, but similar to c_n complex coefficients for
the fourier series, there is a relationship:
c_n = a_n -j*b_n, n > 0
c_n = a_n+j*b_n, n < 0

So, for positive m: a_m = 0.5*(c_m + c_{-m}), and similar approach goes for
b_m.

Sorry I cannot give you more info at this time. Hopefully this sketchy
approach will give you some ideas.

Ivo



2008/10/22 Nils Wagner <nwagner at iam.uni-stuttgart.de>

> Hi all,
>
> Is there a function in scipy to compute the Fourier
> coefficients
> a_0, a_1, b_1, a_2, b_2 of a periodic function f(t)=f(t+T)
>
> http://en.wikipedia.org/wiki/Fourier_series
>
> An example would be appreciated.
>
> Nils
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20081022/5ed26fd3/attachment.html>