вопрос про интерполяцию Опции

[jorikdima]

Здравствуйте.
Возник вопрос, который быть может, надуман, но хочется разобраться.
Есть сигнал - результат НЧ фильтрации исходного сигнала с целью определить его постоянную составляющую. То есть полоса сигнала в общем то составляет доли Герца, но среднее значение (или нулевая гармоника) может варьироваться в довольно большом диапазоне. Частоту дискретизации этого сигнала нужно поднять в 25 раз, для чего требуется провести интерполяцию. В книжках прочитал, что делается это с помощью добавления нужного числа нулей между исходными отсчетами и пропусканием через ФНЧ для удаления размноженных гармоник. В принципе все работает, но появляются артефакты. Дело в том, что во всех книжках, где рассказывается про интерполяцию, нарисован сигнал с околонулевой постоянной составляющей, ну либо нулевая гармоника много меньше других гармоник в сигнале. И тут все прекрасно, добавляются нули, затем они сглаживаются ФНЧ. Но в случае когда, нулевая гармоника является преобладающей добавление в сигнал отсчетов, сильно отличающихся от текущего значения сигнала (а именно нулевые отсчеты) каждый раз вызывают всплеск на выходе фильтра, пропорциональный разнице между нулем и реальным значением сигнала и зависящим от ИХ фильтра. Этот всплеск у меня много больше шума исходного сигнала, являющейся собственно постоянным током и недопустим.
Что делают гуру в таких ситуациях?
У самого есть идея следующая. Не добавлять нулевых отсчетов и не ставить ФНЧ, а вписывать какой-нибудь полином, тот же сплайн, в реальные отсчеты сигнала, а затем, найдя коэффициенты сплайна, рассчитывать 24 значения между исходными выборками. Минус такого решения, что это не риал тайм, нужно дождаться следующего отсчета прежде чем получить 24 точки до него, но мне не страшно это.
Справедливо ли решение или есть какие-то более мудрые действия?
Спасибо.

[des00]

У самого есть идея следующая. Не добавлять нулевых отсчетов и не ставить ФНЧ, а вписывать какой-нибудь полином, тот же сплайн

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

а уйти от около нулевой ПЧ никак ?

[Stanislav]

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

а уйти от около нулевой ПЧ никак ?

При чём здесь нулевая ПЧ, не могли бы просветить?

[des00]

Под зеркалами, хотя это и не правильно, я понимаю то что автор назвал "ФНЧ для удаления размноженных гармоник". Ошибки интерполяции приведут к тому, что они хотя и будут подавлены 30-50дБ, но все равно будут (модель в симулинке, которая демонстрирует этот эффект собирается на раз-два-три)

При чём здесь нулевая ПЧ, не просветите?

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

[Stanislav]

Под зеркалами, хотя это и не правильно, я понимаю то что автор назвал "ФНЧ для удаления размноженных гармоник". Ошибки интерполяции приведут к тому, что они хотя и будут подавлены 30-50дБ, но все равно будут

Это неверно.
Полиномиальная интерполяция способна подавить то, что Вы называете "зеркалами", на любую заданную величину.

PS. Речь идёт о сплайн-интерполяции, конечно. Как того и хочет Автор темы.

У автора проблема с обработкой сигнала на околонулевой ПЧ, то пусть рассмотрит вариант ухода на нормальную, пусть и низкую ПЧ.

Это он Вам сообщил, или сами догадались?

[des00]

Это неверно.
Полиномиальная интерполяция способна подавить то, что Вы называете "зеркалами", на любую заданную величину.

модель в Simulink е в студию, спектр со скруглением RRC = 0.15, интерполяция ну положим на 16ть. Модель должна демонстрировать выходной спектр сигнала. Интерполятор возьмите любой. Нет модели, нет разговора.

Это он Вам сообщил, или сами догадались?

Вам хочется пофлудить, после длительного молчания, или вы сознательно не прочитали слова автора топика "Есть сигнал - результат НЧ фильтрации исходного сигнала с целью определить его постоянную составляющую. То есть полоса сигнала в общем то составляет доли Герца, но среднее значение (или нулевая гармоника) может варьироваться в довольно большом диапазоне."

[Stanislav]

модель в Simulink е в студию, спектр со скруглением RRC = 0.15, интерполяция ну положим на 16ть. Модель должна демонстрировать выходной спектр сигнала. Интерполятор возьмите любой. Нет модели, нет разговора.

Вы про какие-то модели в Симулинке разговор начали - Вам и предъявлять. А я поправлю, если что. Поскольку Симулинком не пользовался сроду.

Вам хочется пофлудить, после длительного молчания,

Что за чушь?
Держите себя в рамках, уважаемый.
Или за время моего отсутствия модераторы электроникса окончательно мутировали в попсу?
Что ни тема - то россыпь перлов. И модераторы в "застрельщиках".

...или вы сознательно не прочитали слова автора топика "Есть сигнал - результат НЧ фильтрации исходного сигнала с целью определить его постоянную составляющую. То есть полоса сигнала в общем то составляет доли Герца, но среднее значение (или нулевая гармоника) может варьироваться в довольно большом диапазоне."

Итак, Вы утверждаете, что, совершив НЧ фильтрацию исходного ПЧ-сигнала, Автор темы получил полосу в доли герца, я правильно Вас понял?

 

Перемежают нулями и выполняют фильтрацию. Наличие артефактов означает неправильно посчитанный фильтр-интерполятор.

Не то, чтобы неправильно, скорей всего... Тут условия не совсем понятны. Видимо, нужно получать результат сразу же по приходу первых отсчётов. А такое действо предполагает некую искусственную процедуру инициализации памяти фильтра - интерполятора.