No niestety nie radzę sobie :/
Funkcja do obliczenia DFT wygląda następująco:
Kod:
//wyliczanie DFT, parametry:
void dft(volatile int N, volatile double samp[], volatile double spect[], volatile double ph[]) {
//zmienne pomocnicze
volatile int n = 0, m = 0 ;
volatile double real = 0.0, imagine = 0.0 ;
for (m = 0 ; m < N ; m++) {
real = 0.0;
imagine = 0.0 ;
//wylicz czesc urojona i rzeczywista
//for(n=0 ; n<N ; n++){
n = 0;
while (n++ < N) {
real += samp[n] * cos(M_TWOPI * (double)(n * m) / (double)(N)) ;
imagine -= samp[n] * sin(M_TWOPI * (double)(n * m) / (double)(N)) ;
}
//usuwanie bledow zaokraglen
if ((int)(imagine * 100000.0) == 0) {
imagine = 0.0;
}
if ((int)(real * 100000.0) == 0) {
real = 0.0;
}
//wyznacz modul wartosci zespolonej
spect[m] = sqrt(real * real + imagine * imagine) ;
ph[m] = atan(imagine / real) / M_PI * 180.0;
}
}
Testuję nie na całym moim buforze 5000 próbek tylko wypełniam bufor z przykładu 128 pierwszymi próbkami:
Kod:
for (index_res = 0 ; index_res < 128 ; index_res++) {
samples[index_res] = measre_audio[index_res];
}
i uruchamiam ją jak w przykładzie:
Kod:
dft(128, samples, spectrum, phase) ;
Funkcja zwraca tablicę trzech wartości:
Częstotliwość, wartość spektrum (? nie wiem czy dobrze rozumiem), i phase (jeszcze nie wiem co to
)
W każdym razie wszystko kompiluje się i uruchamia, ale wyniki są kompletnie bezsensu i tak dla przykładu, kiedy podałem na ADC atmegi sygnał z karty dźwiękowej 500Hz uzyskuję taki wynik:
dla 1000Hz
dla 100Hz
Myślałem, że coś robię nie tak przy samplowaniu wiec wypełniłem tablicę próbek na sztywno wyliczeniami:
(nie wiem ile to warte, było w przykładzie)
Kod:
//wyznacz sygnal do analizy (fs=100Hz)
index = 0 ;
for(time=0.0 ; time<0.1 ; time += 0.01){
samples[index++] = 5.0 + 10.0 * sin(M_TWOPI * 20.0 * time + M_PI/4) + 15.0 * sin(M_TWOPI * 30.0 * time) ;
}
i dostałem coś takiego:
Ktoś ma jakiś pomysł dlaczego?
Może analizuje za mało próbek?
I największa zagadka z tego co wyczytałem to wartość wyliczeń zależy (i pewnie słusznie) z częstotliwości próbkowania, ale skąd u licha w/w funkcja ma wiedzieć że moje sample są próbkowane 1000Hz?
Kanał - ATNEL tech-forum
Zaloguj
Zarejestruj
Informacje o cookies
FAQ
Szukaj
