фрагмент кода частотной модуляции (FM) - PullRequest
Новая
       37

фрагмент кода частотной модуляции (FM)

0 голосов
/ 30 декабря 2011

Я написал следующий код для частотной модуляции аудиосигнала.Сам звук имеет длительность 1 секунда, частота дискретизации составляет 8000 Гц.Я хочу применить FM к этому аудиосигналу, используя синусоидальный сигнал с частотой 50 Гц (выраженный в виде доли от частоты дискретизации).Модулирующий сигнал имеет индекс модуляции 0,25, чтобы создать только одну пару боковых полос.

for (i = 0; i < 7999; i++) {
    phi_delta = 8000 - 8000 * (1 + 0.25 * sin(2* pi * mf * i));
    f_phi_accum += phi_delta; //this can have a negative value
    /*keep only the integer part that'll be used as an index into the input array*/
    i_phi_accum = f_phi_accum;   
    /*keep only the fractional part that'll be used to interpolate between samples*/
    r_phi_accum = f_phi_accum - i_phi_accum; 
    //If I'm getting negative values should I convert them to positive
    //r_phi_accum = fabs(f_phi_accum - i_phi_accum);   
    i_phi_accum = abs(i_phi_accum);
    /*since i_phi_accum often exceeds 7999 I have to add this if statement so as to      prevent out of bounds errors   */  
    if (i_phi_accum < 7999)
        output[i] = ((input[i_phi_accum] + input[i_phi_accum + 1])/2) * r_phi_accum;            
}

1 Ответ

1 голос
/ 30 декабря 2011

Ваш расчет phi_delta отключен с коэффициентом 8000 и смещением - это должно быть 1 +/- небольшое значение, т.е. 

phi_delta = 1.0 + 0.25 * sin(2.0 * pi * mf * i));

, в результате чего phi_delta будет иметь диапазон от 0,75 до 1,25.

...