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?
taajuudet FFT-spektrissä
Tyhmyri mutta ahkera
1
90
Vastaukset
- 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
- 972763
Jotain puuttuu
Kun en sinua näe. Et ehkä arvaisi, mutta olen arka kuin alaston koivu lehtiä vailla, talven jäljiltä, kun ajattelen sinu1042350- 912110
- 161874
Hei A, osaatko
sanoa, miksi olet ihan yhtäkkiä ilmestynyt kaveriehdotuksiini Facebookissa? Mitähän kaikkea Facebook tietää mitä minä en441731- 381725
- 781696
Persuilla ja Saksi-Riikalla meni sitten pornon levittämiseksi koko touhu.
Onko kenellekään yllätys?1031676Synnittömänä syntyminen
Helluntailaisperäisillä lahkoilla on Raamatunvastainen harhausko että ihminen syntyy synnittömänä.1291487Mitä tämä tarkoittaa,
että näkyy vain viimevuotisia? Kirjoitin muutama tunti sitten viestin, onko se häipynyt avaruuteen?411294