Szybka transformacja Fouriera

Opis funkcji FFT

Nagłówek funkcji
bullet void FFT( int p,double *y,double *a,double *b,int *st);
bullet Pobierz funkcję
Zastosowanie
bullet Funkcja FFT oblicza rzeczywiste współczynniki wielomianu trygonometrycznego.

 Dane
bulletp - liczba naturalna taka, że n = 2^p,
bullety - wskazanie na tablicę zawierającą wartości wielomianu w punktach pi*k/n, (k=0,1,2,...,2n-1).
Wyniki
bulleta - wskazanie na tablicę współczynników Aj (j = 0,1,2,...,n),
bulletb - wskazanie na tablicę współczynników Bj (j = 0,1,2,...,n), b[0] = b[n] = 0,
bulletst = 1 jeśli p < 1 lub p > 14, st = 0 w przeciwnym przypadku.

 Przebieg ćwiczenia

Napisać program dokonujący transformacji Fouriera funkcji y(x)
bulletfunkcja y(x) jest funkcją języka C o nagłówku postaci: double y(double x);
bulletpróbki tej funkcji, w punktach pi*k/n, (k=0,1,2,...,2n-1), umieszczone w tablicy, stanowią dane dla funkcji FFT
bulletwynikiem działania programu jest widmo amplitudowe funkcji y(x), czyli wartości sqrt(Aj*Aj+Bj*Bj), j = 0,1,2,...,n
bulletwidmo amplitudowe funkcji należy przedstawić w postaci wykresu prążkowego ChartFX
Przetestować program dla następujących funkcji
bullety(x) = sin(10x)*cos(3x),
bullety(x) = (1+0.5*sin(3x))*sin(40x),
bulletsamodzielnie napisane funkcje reprezentujące sygnały zmodulowane amplitudowo i częstotliwościowo, dla różnych głębokości (indeksów) modulacji.