Rekursywno - adaptacyjna dyskretna transformacja Fouriera

Witaj ,

Musisz poczytać sobie o DSP. Nie wiem jakie masz podstawy matematyczne w temacie analizy Fourierowskiej. Sądząc po pytaniu , niestety mizerne.

Bez podstaw matematyki , dziedziny liczb zespolonych - odradzam branie się za pisanie programu. G...o z tego wyjdzie.

Jest wiele publikacji w sieci na ten temat , ale wszystkim nowicjuszom odpowiadam :

Co książka , to książka ! Głównie jest to matematyka. Jest wiele gotowców w sieci , ale dopóki nie zrozumiesz tematu od podstaw - odpuść sobie.

Pozdrawiam ,

Mariusz Hajduk

Adam Miazga napisał(a):

w dużym uproszczeniu SQRT(-1) to symbol nieoznaczony wprowadzony sztucznie, dzięki niemu j^2 może wynosić -1 i matematycznie wszystko się zgadza Romek

Witam Dokładnie tak. W matematyce jest "i", a w teorii obwodów "j", bo "i" myliło by się z oznaczeniem prądu i to cała tajemnica :). W sposób strawny (dla mnie) teoria obwodów jest przedstawiona w - Athanasios Papoulis "Obwody i układy" WKiŁ 1988. Polecam.

[...]

jesli chcesz tylko moduł widma dla danej harmonicznej to e^(-j*KĄT) = arccos(KĄT) gdzie (KĄT) jest w radianach.

2*PI*f to pulsacja (prędkość kątowa) gdzie leży f w tym wzorze ? ano n/N (indeks aktualnej próbki/ilość próbek) k oznacza aktualną harmoniczną k*f to wielokrotność częstotliwości podstawowej 2*PI/N fachowo nazywa się binem ( najmniejszym odstępem między sąsiednimi harmonicznymi)

oops może wystarczy bo mnie zaraz sclerozis omyli ... :)

jeśli nie potrzebujesz widma fazy, to j ci niepotrzebny

Pozdr el es

Paweł Pracki napisał(a):

Lepsza jest ostatnio dostępna w księgarniach Richard G. Lyons "Wprowadzenie do cyfrowego przetwarzania sygnałów". Znacznie bardziej praktyczna od Papoulisa, nie mówiąc o wykładowcach PW ;)

Adam Miazga napisał(a):

Funkcję exp(ix) możesz rozłożyc na cos(x) + i*sin(x). Dzięki temu możesz uniknąć w programie implementacji liczb urojonych. Podstawy tego masz w w książkach do analizy matematycznej a praktyczne rozwiązania w książkach do metod numerycznych.