|
Отправлено
Vallav 18 января 2002 г. 13:40 В ответ на: почитал ветку - интересно отправлено yes 18 января 2002 г. 11:09 |
 |
|
|
"берете коэффициенты Фурье (БПФ) и по ним считаете сумму интегралов в смещенной точке?"
Изложу проблемму еще раз ( алглротмический подход? )
Есть исходная выборка Xn, n от 0 до N-1. Но редкая, хочется значений
в промежуточных точках. Делаем так:
Берем БПФ. Получаем Ck, k от 0 до N-1. Умножаем на 4, разрезаем в
районе N/2 и вставляем в разрез N*3 нулей. Берем обратное БПФ.
Получаем N*4 точек. Точки 0, 4, 8, ... повторяют исходный массив,
остальные - промежуточные значения. Все хорошо, если Ck в районе
k=N/2 близки к нулю, тогда это идеальная интерполяция - восстановление
сигнала. Но если CN/2 и окресность большие?
1. Куда девать CN/2? Оставить слева, сдвинуть вправо, распилить над-
вое, без разницы?
2. Своими манипуляциями мы внесли скачек в частотной области, что
неизбежно приведет к переколебаниям в временной ( эффект Гиббса ).
Но ведь у нас туча нулей! Не записать ли вместо них что нибудь
более путное. Вариант, который обсуждался - записать вместо нулей
периодически продолженный исходный спектр и помножить затем на под-
ходящий фильтр. Но какой?
"кажется чтобы победить Гиббса можно сдвигать окно (на котором считается БПФ) не на N, а на N/2"
Если не ошибаюсь, это из другой оперы - вычисление скользящего спек-
тра на большом массиве.
E-mail: 
info@telesys.ru