taajuudet FFT-spektrissä

Tyhmyri mutta ahkera

Minulla on paha ongelma spektrianalyysin kanssa.

Minulla on digitaalista reaalista dataa x, joka on saatu näytteenottotaajuudella 48 kHz.

Lasken sille periodogrammispektrin FFT:llä N pisteessä käyttäen käskyä

y=abs(fft(x)).^2.

Spektri on siis N-pisteinen y. Sitten lajittelen spektriarvot Matlabin käskyllä

[value, index]=sort(y).

value kuvaa spektriarvoja ja index niitä vastaavia taajuusarvoja.

Taajuudethan on normalisoitu välille [0,1] tai [-0.5,0.5], missä [0,0.5] vastaa positiivisia taajuuksia.

Nyt tulee kysymys: millä laskuoperaatiolla nämä index-arvot muutetaan alkuperäisiksi Hz-arvoiksi?

1

88

    Vastaukset

    Anonyymi (Kirjaudu / Rekisteröidy)
    5000
    • 4+1

      Kaksipuolinen spektri normalisoidaan välille [-0.5, 0.5]. Huomaa, että jos otat vain välin [0, 0.5], niin joudut kertomaan kahdella jotta signaalin energia säilyy. Reaalisen signaalin tapauksessa toinen puolikas riittää.

      Alkuperäiseen kysymykseen vastaus on, että nuo taajuudet -0.5, 0.5 vastaavat Nyquistin taajuutta. Jos "alkuperäinen Hz-arvo" tarkoittaa esimerkiksi mittalaitteen mittataajuutta, niin kerrot tuon matlabin linspacen suoraan sillä mittalaitteen taajuudella. Laskuoperaation nimi siis on peruskoulusta tuttu kertolasku.

    Ketjusta on poistettu 0 sääntöjenvastaista viestiä.

    Luetuimmat keskustelut

    1. Pupuhuhdasta löytyi lähes sadan kilon miljoonalasti huumeita

      Pupuhuhdasta löytyi lähes sadan kilon miljoonalasti huumeita – neljä Jyväskylän Outlaws MC:n jäsentä vangittu: "Määrät p
      Jyväskylä
      55
      1796
    2. Persut petti kannattajansa, totaalisesti !

      Peraujen fundamentalisteille, vaihtkaa saittia. Muille, näin sen näimme. On helppo luvata kehareille, eikä ne ymmärrä,
      Maailman menoa
      43
      1583
    3. Ei luottoa lakko maahan

      Patria menetti sovitun ksupan.
      Suomen Keskusta
      49
      1539
    4. Nähtäiskö ylihuomenna taas siellä missä viimeksikin?

      Otetaan ruokaöljyä, banaaneita ja tuorekurkkuja sinne messiin. Tehdään taas sitä meidän salakivaa.
      Ikävä
      4
      1485
    5. Sinäkö se olit...

      Vai olitko? Jostain kumman syystä katse venyi.. Ajelin sitten miten sattuu ja sanoin ääneen siinä se nyt meni😅😅... Lis
      Ikävä
      5
      1468
    6. Housuvaippojen käyttö Suomi vs Ulkomaat

      Suomessa housuvaippoja aletaan käyttämään vauvoilla heti, kun ne alkavat ryömiä. Tuntuu, että ulkomailla housuvaippoihin
      Vaipat
      5
      1384
    7. Hyvää yötä ja kauniita unia!

      Täytyy alkaa taas nukkumaan, että jaksaa taas tämän päivän haasteet. Aikainen tipu madon löytää, vai miten se ärsyttävä
      Tunteet
      7
      1298
    8. Lepakot ja lepakkopönttö

      Ajattelin tehdä lepakkopöntön. Tietääkö joku ovatko lepakot talvella lepakkopöntössä ´vai jossain muualla nukkumassa ta
      10
      1260
    9. Revi siitä ja revi siitä

      Enkä revi, ei kiinnosta hevon vittua teidän asiat ja elämä. Revi itte vaan sitä emborullaas istuessas Aamupaskalla
      Varkaus
      3
      1132
    10. Kello on puoliyö - aika lopettaa netin käyttö tältä päivältä

      Kello on 12, on aika laittaa luurit pöydälle ja sallia yörauha kaupungin asukkaille ja työntekijöille. It is past midni
      Hämeenlinna
      3
      1113
    Aihe