Weaver (математика)

[GenaSPB]

Да понятно. Я про fpga.

[GenaSPB]

Кстати, маскиратор с полностью цифровой обработкой тоже не должен свистеть.

[Игорь 2]

Кстати, маскиратор с полностью цифровой обработкой тоже не должен свистеть.

У меня он был полностью аналоговый. Сто лет назад дело было...

[ra0ahc]

ЛАБА
Смеситель

freq = 400;//сигнал герц
freqOpora = 2000;// несушка

Заполнение буферов выглядит так:

        float32_t amp = 600;//амплитуда, макс 2048 для 12бит

        for (int i = 0; i < FRAME_SIZE; i++) {
            float32_t opora = arm_sin_f32(phaseOpora);
            float32_t signal = arm_sin_f32(phase);


            pOutDac[dacTxBankNawNumber][i ] = ((uint32_t) (amp + signal*opora * amp));

            phase += delta;
            if (phase >= pp) phase -= pp;

            phaseOpora += deltaOpora;
            if (phaseOpora >= pp) phaseOpora -= pp;

            opora = arm_sin_f32(phaseOpora);
            signal = arm_sin_f32(phase);

            pOutDac[dacTxBankNawNumber][i ] |= (((uint32_t) (amp + signal*opora *amp))) << 16;
            phase += delta;
            if (phase >= pp) phase -= pp;

            phaseOpora += deltaOpora;
            if (phaseOpora >= pp) phaseOpora -= pp;

        }

[ra0ahc]

а вот тоже самое, но для сигнала 645Гц

[ra0ahc]

После запуска FIR фильтра на 3кГц + шум
Получился на выходе цап такой сигнал

[ra0ahc]

Интересный феномен я нашел. Не знаю что это за глюк в математике. Но вот смотрите (вопрос к математикам)

амплитуду умножаю на 10. (все числа float)
 float32_t amp = 199*10;
.
.
дальше перемножаю сигналы (они тоже float) и делю на теже 10.  цифра от -1 до +1

arrDac[i ] = (signal* opora)/10;  //цифра от -1 до +1

дальше формирую амплитуду для цап от 0...4000
умножаю на амплитуду и смещаю на эту же амплитуду, убираю отрицательный значения.

 pOutDac[dacTxBankNawNumber][y] = (uint32_t) (amp + firOutDac[i ] * amp);

И получается вот такая красивая картинка....
(все дальше 3кГц это перегруз звуковой платы, в программе работает фир фильтр со срезом 3кГц)

[ra0ahc]

А теперь убираю умножения 10 и деление на 10
и получается вот такая картинка

[ra0ahc]

ачх FIR

[ra0ahc]

Нашел причину 
Программер фигов. А всего то надо было выспаться.
Смещение амплитуды надо делать константой.

pOutDac[dacTxBankNawNumber][y ] = (uint32_t) (2048 + firOutDac[i-1] * amp);