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

[jorikdima]

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

[Stanislav]

 

Этот всплеск у меня много больше шума исходного сигнала, являющейся собственно постоянным током и недопустим.
Что делают гуру в таких ситуациях?

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

...Справедливо ли решение или есть какие-то более мудрые действия?

Линейная фильтрация - лишь частный приближённый способ интерполяции сигнала, основанный на предположении об ограниченности его спектра.

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

ЗЫ. Это на пальцах, простите...
ЗЗЫ. О том, что у любого ФНЧ существует задержка, не забыли?
Кроме того, Вы не привели тип использованного фильтра-интерполятора. Я понял так, что это БИХ или минимально-фазовый, но необходимо уточнить.