есть ли методы дизеринга (dithering) для линейной интерполяции при повышении тактовой частоты (upsampling), хочется что-то аналогичное дизерингу при генерации синуса Опции

[Krys]

Здравствуйте. Подскажите, пожалуйста, есть ли методы дизеринга (dithering) для линейной интерполяции при повышении тактовой частоты (upsampling, digital up-converter), хочется что-то аналогичное дизерингу при генерации синуса в DDS, когда периодичность ошибки фазы вызывает группировку шума на каких-то частотах спектра, а дизеринг этот шум размазывает, улучшая SFDR. Мне надо что-то такое же.
Поясню на картинках.

На этой картинке синим спектр исходного сигнала (ЛЧМ), до передискретизации. Красным - спектр после повышения тактовой и линейной интерполяции. Зелёным - спектр после повышения тактовой и фильтрации обычным КИХ-фильтром.


Эта картинка из временного домена, после повышения тактовой (увеличено вблизи конца временной выборки, когда частота ЛЧМ близка к максимальной). Синим - результат линейной интерполяции, красным - результат интерполяции через КИХ-фильтр (зелёный спектр на предыдущей картинке), а зелёным - разница между синей и красной кривыми.

По первой картинке (красная кривая) видно, что присутствуют зеркальные спектры от ЛЧМ. Хочется их маленько ослабить, размазав, т.к. есть запас по частоте в половину диапазона (полоса ЛЧМ - четверть частоты дискретизации). Т.е. аналогично размазывают дизерингом составляющие от DDS. Вот я подумал, нельзя ли что-то подобное сделать для линейной интерполяции? Куда тут только шум влепить... в узловые точки интерполяции?

[Krys]

Почему вы думаете, что метод перераспределения внеполосной энергии, который подходит для случая 1 подойдет для случая 2?

"у меня нет ресурсов на то, чтобы реализовать фильтр с умножителем, сплайн или даже CIC звенья, но, чтобы стало лучше, я могу куда-нибудь добавить шум"

Отвечу на оба сообщения: я не хотел сказать, что мне подойдёт именно один в один метод, который используется в DDS. Я лишь хочу верить, что возможно существуют простые в реализации способы разравнивания спектра и для линейной интерполяции, просто я о них не знаю. Поэтому и спросил.
Offtopic: Дополнительно 2 Fat Robot: почему вы позволяете себе поясничать и выставлять себя выше других? Нельзя просто нормально разговаривать?

да, кстати, а как на лету считать коэффициент наклона для интерполяции между точками x(i), x(i+1) : k = (x(i+1)-x(i))/K, это же деление) не очень вяжется с крайне ограниченным ресурсом.

Это не стыкуется с:
Нужно делить на 10.. "(грубо)". ;)

Для простоты положим, что апсемплинг на величину степени двойки. Более того, за 10 тактов посчитать один коэффициент наклона всё же можно, даже при минимуме ресурсов и без степени двойки. Тем более за 100 тактов.

Еще подумайте сразу над тем, как в своей полосе скомпенсировать неравномерность в 5-10 дБ без умножителей или объяснить заказчику, что это "ничего страшного".

Не понял, откуда неравномерность? У меня по картинке это порядка 1 дБ (я приблизил в матлабе).