Задумал я сделать модулятор, в котором ЦАП работает на частоте не целочисленно кратной символьной. В принципе ничего сложно, спектр сигнала сформирован на частоте 16sps, что при скруглении 0.25 соответствует нормированной крайней частоте спектра ~0.16, по идее этого с запасом хватает что бы уложиться в диапазон линейности интерполятора лагранжа 3 го порядка. Но в представленной модельке видны большие искажения спектра сигнала.

Уважаемые гуру ресамплинга подскажите я правильно понимаю что нелинейность АЧХ/ФЧХ интерполятора на высокой частоте работает как смеситель, который рождает гармоники более высоких частот, которые и видно в представленной модели? В таком случае интерполятор 3 го порядка не вариант? нужно либо уходить на интерполятор более высокого порядка (например 5 го как в модели) либо сужением спектра сигнала в 2 раза.
Спасибо.

Да очевидно что фильтр с переменными параметрами это уже нелинейщина.

в таком случае при использовании sinc интерполяторов или FIR интерполяторов этот эффект также будет проявляться? Несмотря на то, что они более широкополосные чем лагранж? Те же Provigent использует в своих чипах такую технику, как они тогда борятся с внеполосными излучениями вызванными нелинейностью?

В принципе любые интерполяторы-дециматоры это нелинейные устройства, борются большими порядками.

понятно, рассчитал модифицированный лагранж 4 го порядка на тот же диапазон mu[-1 : 0] увел нелинейность ниже 80 ти дБ, мне этого хватит. Всего 1500 плиток, зато сколько удовольствия %)

ЗЫ. Для такой интерполяции разрядность mu требуется много больше чем для интерполятора в демодуляторе, в передатчике желательно не менее 12 бит, а в приемнике мне 7ми бит хватало %)

Нифига.
sinc-интерполятор линеен.
(если я правильно понял о чем речь)

Для нелинейности нужно, чтобы в формулы входили произведения или иные нелинейные функции отсчетов преобразуемого сигнала. Просто произвольный момент времени выборки не порождает нелинейность.

Нелинейны сами операции изменения частоты дискретизации, вот и имеем например наложение спектров при децимации, т. е. появление новых спектральных составляющих которых не было в исходном сигнале.

И тем не менее.
F(x+y) = F(x) + F(y)
F(a*x) = a * F(x)



Для сохранения спектра нужна еще стационарность, очевидно.

Да похоже я ошибался.

тогда как объяснить наблюдаемый в модели эффект ?

К сожалению, эффекты модели я не наблюдаю - не запускается.
Но в модели у вас ведь Фарроу? Про Фарроу я ничего не писал.

я думал я один сижу на древнем матлабе 2007а, вот моделька для R14 матлаба. В модели все блоки стандартные, ИМХО проблем быть не должно %)

Насколько я понял литературу Фарроу это реализация интерполятора, я использую интерполятор построенный на многочленах Лагранжа.

Всё верно, только Oldring про интерполятор Лагранжа тоже не писал. Его замечание относилось к идеальному sinc интерполятору. Это полиномиальные интерполяторы нелинейны.
Но Вам же не хочется реализовать sinc-интерполятор))

собрал модельку с управляемым генератором синуса, если поиграться с частотой синуса и кратностью преобразования частоты, видно :
1. чем выше частота синуса, тем больше продуктов нелинейности в выходном спектре
2. чем больше значений дробной задержки используется, тем так же больше продуктов нелинейности (например посмотрите вместо 14/15 кратности 1, 1/2, 1/4)
у меня только одно объяснение, не смотря на то, что при фиксированном mu интерполятор элемент линейный, при разных mu коэффициент передачи изменяется (см рис) поэтому и возникает сия нелинейность на спектре


почему не хочется то? просто под рукой был готовый RTL ный 3 го порядка, он и пошел в дело. Потом уже увидев сей эффект на понарамнике, полез более подробно смотреть. Не подскажите литературу по sinc интерполяторам ?

Да прямо по теореме Котельникова интерполируйте истинным синком

Под аппартными "синк-интерполяторами" обычно понимают такую-же кусочно полиномиальную интерполяцию, только изображающую из себя синк. Такие же как интерполятор Лагранжа сплайны, только более длинные. Там тоже будут продукты? но меньше.
Литературы по общей теории интерполяторов не подскажу, но практический вариант улучшения интерполятора Лагранжа я когда-то давно помещал в форум - там приделывают хвосты к Фарроу и получается несколько лучше
http://electronix.ru/forum/index.php?act=A...st&id=31882