CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Main CFD Forum (http://www.cfd-online.com/Forums/main/)
-   -   Derivation by spectral method (http://www.cfd-online.com/Forums/main/8570-derivation-spectral-method.html)

 daniel January 19, 2005 17:50

Derivation by spectral method

I want to evaluate the derivative of a function f=f(x) periodical, by a spectral method. Let's say

(1) f(j)=1/N*SUM(with k=from 1 to N)

( Ck*exp(-i*2pi*(k-1)*(j-1)/N )

where j= index of point N= number of points k is the wave number f is the fundamental frequency N is the number of points Ck are the spectral coefficients

The (1) is the Discrete Fourier decomposition. Now, I calculate the Ck coefficients by the inverse trasformation and I can use the (1) to re-build my function. At the same time I can use the Ck to calculate the derivative of the expression (1), which should be:

(2) df(j)=1/N*SUM(with k=from 1 to N) (-i*2*pi*(k-1)*Ck/(N*DX))*exp(-i*2pi*(k-1)*(j-1)/N )

where DX is the sample. What happens to me is that The re-built function is correct, but not the derivative! Do you have any hint?