logo elektroda
logo elektroda
X
logo elektroda
REKLAMA
REKLAMA
Adblock/uBlockOrigin/AdGuard mogą powodować znikanie niektórych postów z powodu nowej reguły.

Transformata FFT funkcji autokorelacji - jak wyznaczyć okres

majery 16 Lip 2009 18:57 4006 5
REKLAMA
  • #1 6786584
    majery
    Poziom 2  
    Witam, mam następujący problem.
    Piszę w Matlabie skrypt wyznaczający odstęp częstotliwościowy Δf pewnego sygnału OFDM. Wyznaczyłem widmową gęstość mocy (PSD), następnie funkcję autokorelacji widmowej gęstości mocy Ryy. Maksimum Ryy odpowiada Δf. Kolejne maksima pojawiają się z okresem odpowiadającym Δf, czyli na przykład:
    odczytana z wykresu odległość pomiędzy kolejnymi maksimami w próbkach wynosi ok. 61, cz. próbkowania fs = 8000 Hz, nfft = 4096, zatem 61*fs/nfft = 119.1406, co w przybliżeniu daje Δf (120 Hz). Tylko, że ja chciałbym ten okres wyznaczyć za pomocą FFT, tzn. [peakValue, indeks] = max(abs(FFT(Ryy))) - nie potrafię jakoś zinterpretować (przeskalować) wyniku, tzn. pozycji na której występuje maksimum modułu, aby otrzymać wielkość w przybliżeniu równą 120 Hz.
    Za podpowiedź będę bardzo wdzięczny :).
  • REKLAMA
  • #2 6790866
    Tantalos
    Poziom 18  
    A jaki sygnał kolega chce analizować i jakie otrzymuje wyniki?
  • REKLAMA
  • #3 6791683
    majery
    Poziom 2  
    OFDM :) - szerokość widma ok. 2800 Hz, liczba podnośnych 20 + 1 pilot (na 3300 Hz), Δf = 120 Hz, szybkość manipulacji w podkanale ok. 75 Bd, modulacja DQPSK (wariant B - pseudo-ośmiowartościowa).
    Maksimum transformaty fft funkcji autokorelacji obliczam tak:

    [peakValue, indeks] = max(abs(fft(Ryy, nfft)))

    peakValue = 100.7748, indeks = 4.

    Nie wiem jak teraz wyznaczyć Δf. Myślałem, że tak: indeks*(fs/nfft), ale wynik otrzymuję: 7.8125 :(
  • REKLAMA
  • #4 6802085
    Tantalos
    Poziom 18  
    Z tego co znalazłem to widmo sygnału OFDM jest płaskie, dlatego trudno jest rozróżnić podnośne. Jeżeli jest możliwość podania stałej wartości na modulator DQPSK, tak żeby faza nośnej się nie zmieniała w sposób przypadkowy, to widmo nośnych będzie wąskie i wtedy można będzie rozróżnić minima i maksima, a sygnał będzie bardziej skorelowany.
  • REKLAMA
  • #5 6802495
    majery
    Poziom 2  
    Tantalos dziękuję za odpowiedź. Ja jednak uparcie poszukuję odpowiedzi na pytanie jak wyznaczyć okres funkcji autokorelacji widmowej gęstości mocy. Wydawało mi się, że najprościej wyznaczyć transformatę FFT funkcji autokorelacji, tylko że potem mam problem z interpretacją wyników (przeskalowaniem osi X). Metodę takiego wyznaczania parametru Δf znalazłem w literaturze, ale nie potrafię tego praktycznie zrobić. Jak pisałem wyżej, "ręczna" interpretacja wykresu funkcji autokorelacji daje poprawny wynik - teraz chciałbym "automatycznie" wyznaczyć okres Ryy. A oto moje wyniki graficznie:

    Transformata FFT funkcji autokorelacji - jak wyznaczyć okres
  • #6 6806476
    Tantalos
    Poziom 18  
    Sygnał po autokorelacji ma również składowe o niższych częstotliwościach co pokazuje FFT. To dlatego, że podnośne zachodzą na siebie. Może zamiast autokorelacji użyj funkcji autokowariancji.
REKLAMA