Как шумит радио (цифровая обработка и слуховой приём)

[alex_m]

фаза
Я не увидел в каком месте вы передаете последнее значение фазы из предыдущего вызова процедуры в последующий.

Вот:

        private double last_phase;            // тут хранится значение фазы между вызовами

        private const double TWO_PI = 2D * Math.PI;


        public void Demodulate(Complex* iq, float* audio, int length)
        {
            var gain = FMGain;
            if (_mode == FmMode.Narrow)
                gain *= NarrowAFGain;
            for (var i = 0; i < length; i++)
            {
                // Angle estimate
                var phase = Math.Atan2(iq[i].Imag, iq[i].Real);

                var dphase = phase - last_phase;
                if (dphase < -Math.PI) dphase += TWO_PI;
                if (dphase > Math.PI) dphase -= TWO_PI;
                dphase /= Math.PI;
                last_phase = phase;


                // Scale
                audio[i] = (float)(dphase * gain);
            }
        }

[EU1SW]

просто я явно передавал его в процедуру, равно как и возврат из процедуры, так вернее и красифше КМК

float fmdemod(complexf* input, float *output, int input_size, float last_phase)
{
    бла бла
    {
бла бла бла
    }
    return last_phase;
}

ну и деемфазис после демодуляции ЧМ крайне желателен, для компенсации треугольного спектра шумов

[Genadi Zawidowski]

Тут с удивлён ем обнаружил, что полосу пропускания до демодулятора nfm ставил как удвоенную верхнюю нч.. 9 кГц правильно?

[EU1SW]

С моей точки зрения все верно, делаю так же, для индексов модуляции до 1 (узкополосный ЧМ) можно принимать полосу сигнала до демодулятора как удвоенную верхнюю частоту модуляции, ну если индекс будет немного более 1, то можно взять немного с запасом, 9 кГц КМК вполне
Другое дело что иногда и 15, и 25 кГц не помешали бы, но я остановился на максимальной полосе для АМ/ЧМ в 10 кГц, столкнусь с реальной проблемой нехватки полосы до детектора - перекомпилирую, тут не долго )))

Добавлено через 9 минут(ы):

если конкретно энергетику просчитывать, то при М=1, полосе НЧ 3 кГц, результирующий спектр ЧМ будет иметь полосу 12 кГц по уровню -34 дБ,
полосы второго порядка имеют уровень по амплитуде в 0,11, т.е. даже с 6 кГц фильтром мы потеряем только ~2,5% от мощности сигнала, а через 9 кГц фильтр и подавно

[Genadi Zawidowski]

Вот так оно получается... ПЧ=9 кГц, НЧ=300..3400

[alex_m]

хм... сегодня послушал эфир с новым вариантом FM демодулятора - а ведь действительно лучше работает
На сигналах средней мощности пропал дребезг, звук стал мягким

[EU1SW]

Ну и отлично! Я завтра/послезавтра доберусь до исходников, запощу сюда код для постобработки деемфазиса. Будет еще лучше )

[Genadi Zawidowski]

Я уже забыл... А с чего пришлось вам делать демодулятор? Вроде код откуда-то взят... Там что-то не так было?

[EU1SW]

Насколько я понял это с сдр шарп код... Ну похож )

[alex_m]

да, с сдр шарп экспериментирую.
Думаю свой сдр приемник сделать, собираю информацию и алгоритмы, разбираюсь как оно работает