[Aqvis]

Dlaczego przy korzystaniu z algorytmu fft, gdy na poczatku okresla sie zakres dziedziny czasu, odrzuca sie ostatnia mozliwa probke?

#studbaza #matematyka #nauka

[groman43]

@Aqvis: Przyznam, że pierwszy raz o czymś takim słyszę. Nie widzę też powodu, dla którego, wartość jednej próbki mogła coś zmienić.

Jedyne sensowne wytłumaczenie, jakie przychodzi mi do głowy, to że po prostu musiałeś zadbać, żeby liczba próbek była potęga dwójki.

[Aqvis]

@groman43 https://www.mathworks.com/help/matlab/ref/angle.html
Chodzi o dziedzine czasu w zakladce fft phase

[Aqvis]

@groman43 i faktycznie jest to konieczne tylko do stworzenia widma fazowego

[sezzart]

@Aqvis: a to nie jest tak ze 100Hz przez 1 sekunde to znaczy ze potrzebujesz 100 time-stampow? jak policzyc od 0 do 1 wlacznie to masz 101 probek wiec odrzucasz ostatnia?

[groman43]

@Aqvis: Widmo fazowe też się rozmazuje (przynajmniej w Octavie, nie mam teraz MATLABa pod ręką) jeśli liczba próbek nie jest dokładną wielokrotnością częstotliwości próbkowania.
Jestem w pracy i mam ważniejsze rzeczy do roboty, ale wieczorem postaram się to policzyć na kartce.

[Aqvis]

@sezzart tak tez o tym myslalem, ale nadal nie sprawdza sie, gdy czestotliwosc nie jest liczba calkowita

[Aqvis]

@groman43 chyba wlasnie o to chodzi

[sezzart]

@Aqvis: to nie wiem :)

[Aqvis]

@groman43 w kazdym razie tak jak wspominalem wyzej probowalem z czestotliwoscia sygnalu, ktora nie jest liczba calkowita i wtedy widmo( przy zachowaniu reszty parametrow) tez sie rozjezdza

[groman43]

@Aqvis: Chyba (słowo klucz) w końcu skojarzyłem o co chodzi. Od razu jak zobaczyłem Twoje pytanie, to ta odpowiedź przyszła mi do głowy, ale nie pamiętałem szczegółów.

Najprawdopodobniej to jest typowy przykład spectral leakage. DFT (a więc również i FFT) ze względu na swoją dyskretną naturę, może reprezentować tylko skończony zbiór częstotliwości. Bardziej konkretnie, częstotliwości, które możemy obserwować „bez problemu” muszą być wielokrotnością fs / L (L reprezentuję długość wektora zwróconego

[Aqvis]

@groman43 dziekuje bardzo, chyba juz rozumiem