404 2

404

9. Metody Fouriera

SUBROUTTNE FFT (A, M)

COMPLEX A (1), U, W, T

N=2**M

KV2=N/2

NMl=N—1

J»1

DO 7 1=1, NM1 IF(I.GEJ) GO TO 5 T=A(J)

A(J)=A(I)

A(I)=T

5    K = NV2

6    IF (K.GEJ) GO TO 7 J=J-K

K=K/2 GO TO 6

7    J=J + K

DO 20 L-l, M

LE=2**L

LEL-LE/2

V=(U 0.)

ANG = 3.141592653 58979/LEl W=CMPLX (COS (ANG), SIN (ANG»

DO 20 J=l, LEI DO 10 I=J, N, LE IP=r+LEl T=A(IP)*U A(IP)=A(I)-T 10 A(])=A(I)+T 20 U~U*W RETURN END

3. Napisać podprogram szybkiego przekształcenia Fouriera dla przypadku, gdy jest potęgą 2.

9.4. Przedłużanie okresowe funkcji nieokresowej

Im bardziej regularna jest funkcja, tym szybciej jest zbieżny jej szereg Fouriera, funkcja / i jej pochodne aż do /c-tej włącznie są okresowe i wszędzie ciągłe, a fdnkcj /^<k+1> jest przedziałami ciągła, to

(9.4.1)