Интерполяция сигнала, вопрос... Опции

[TigerSHARC]

Хочу минимизировать эфект "растекания" спектра. Сигнал содержит 40 гармоник кратных основной частоте.
Известен интервал изменения основной гармоники (50+ 1Гц)
Частота дискретизации задана жёстко 16000Гц.
По расчётам кол-во точек сигнала равно N = 8192.
Затем по некоторому алгоритму получаем частоту основной гармоники (пусть частота равна fосн. = 50.55Гц).
Далее, по идее, для минимизации "растекания" нужно интерполировать сигнал с учётом полученной частоты основной гармоники.

Расчётные данные показали, что нужно пересчитать на частоту fd = (327.68*fосн) = 16564,224 Гц.
________________________________________________
Вот здесь самое интересное начинается!!!

Чтобы идеализировать интерполяцию просто изначально задаю частоту 16564,224 Гц на том же интервале времени, который равен интервалу времени выборки 8192 отчётов на частоте 16000Гц (т.е. (1/16000)*8192 = 0,512 с.).
Другими словами окно временное оставляю тем же, что и для 16000Гц.

Конечно же в спектре сигнала наблюдается сильное искажение...

Вот такой вопрос:

Я сделал вывод что при интерполяции каким-то образом меняется величина временного окна... иначе никак((
КАК ЭТО ДЕЛАЕТСЯ????

Встал на этом вопросе прирасчётах, очень нужна помощь...

Заранее спасибо всем ответившим!

[thermit]

Хочу минимизировать эфект "растекания" спектра

Такого эффекта на сегодняшний день наука не знает...

Известен интервал изменения основной гармоники (50+ 1Гц)
Частота дискретизации задана жёстко 16000Гц.
По расчётам кол-во точек сигнала равно N = 8192.
Затем по некоторому алгоритму получаем частоту основной гармоники (пусть частота равна fосн. = 50.55Гц).
Далее, по идее, для минимизации "растекания" нужно интерполировать сигнал с учётом полученной частоты основной гармоники.

Расчётные данные показали, что нужно пересчитать на частоту fd = (327.68*fосн) = 16564,224 Гц.

угу. Вы хотите сделать передискретизацию так, чтоб гармоника была кратна бину дпф.

Чтобы идеализировать интерполяцию просто изначально задаю частоту 16564,224 Гц на том же интервале времени, который равен интервалу времени выборки 8192 отчётов на частоте 16000Гц (т.е. (1/16000)*8192 = 0,512 с.).
Другими словами окно временное оставляю тем же, что и для 16000Гц.

Конечно же в спектре сигнала наблюдается сильное искажение...

Вот такой вопрос:

Я сделал вывод что при интерполяции каким-то образом меняется величина временного окна... иначе никак((

Гы. Вы повышаете частоту дискретизации. На интервале 0.512 c получится больше отсчетов.
В Вашем случае их станет 8480,882688. Никакой метафизики тут нету...

КАК ЭТО ДЕЛАЕТСЯ????

Что именно? Как вычислить дпф на 8к точек по 8480,882688-и точкам? Дык вариантов много.
Можно вычислить несколько дпф на 8к точек со смещением по выборке из 8480, а потом их усреднить. Есть еще более сложные способы.
Искажения спектра скорее всего связаны с горбатой интерполяцией.