Демодуляция RTL SDR IQ AM

Проект, над которым я сейчас работаю, настраивается на частоты PMR и сохраняет звук в файл WAV, однако, когда я получаю образцы IQ от устройства SDR и демодулирую в AM - я слышу только шум.

Изменить, вот RAW IQ, записанный в файл WAV:https://voca.ro/iSRELps4JBg

Для демодуляции я использовал такой подход:

https://github.com/simonyiszk/minidemod/blob/master/minidemod-wfm-atan.c

Однако звук кажется шумным...

  while(get_samples(buffer, buffer_size) == 1)
    {


        // Demodulate to AM
        uint8_t uiAMBuff[buffer_size];

        double i1=0, i2, q1=0, q2;
        short s; // 2 bytes
        double sum=0;
        double dbAmp=0;
        for (int i = 0; i < buffer_size; i+=2)
        {
            // I / Q

            i2=((unsigned char)buffer[i]-127);
            q2=((unsigned char)buffer[i+1]-127);


            double phi1=get_phase(i1,q1);
            double phi2=get_phase(i2,q2);

            double dphi=phi2-phi1;
            if(dphi<-PI) dphi+=2*PI; //dphi below -180°? Let's interpret it in the other direction!
            if(dphi>PI) dphi-=2*PI;  //dphi above 180°? Let's interpret it in the other direction!

            // Now let's use signed 16 bit output rather than 8 bit unsigned.
            s=((SHRT_MAX-1)/PI)*dphi; //Okay so let's get the actual sample. dphi should be between -PI and +PI.

            // Store AM in the buffer that gets written to a WAV file
            uiAMBuff[i]=s;
            uiAMBuff[i+1]=s>>8;

            i1=i2;
            q1=q2;
            // This is to calculate DB
            double y = s/32768.0;
            sum += y * y; 

        }

Как правильно демодулировать IQ в AM?

Источник

2 ответа

Решение

Я сам хочу опубликовать ответ.

Вот как я демодулирую IQ (RTLSDR) в AM, взяв образцы из SDR.

uint8_t uiAMBuff[размер_буфера];

float sum=0;

for (int i = 0; i < buffer_size; i+=2)
{

    // Get 0 - 255
    int _i=buffer[i];
    int _q=buffer[i+1];

    // This will range : 0 - 32768 (short) 2 bytes
    short amplitude = sqrt((_i*_i)+(_q*_q));
    if(amplitude > 32768) amplitude = 32768;

    // Store in separate buffer (serialize) 2 bytes
    uiAMBuff[i]=amplitude>>0;
    uiAMBuff[i+1]=amplitude>>8;

    // Calculate AM Amplitude which is always 0..1 thats why division by 32768.0
    sum += amplitude / 32768.0;

}
Источник

Ваш код предназначен для FM-демодулятора (частота как производная фазы или дельта фазы в дискретном случае). Для AM-демодуляции вам нужно только получить величину (abs()) компонентов IQ и фильтр нижних частот.

Источник
Другие вопросы по тегам