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

[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Гц.

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

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

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

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

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

[TigerSHARC]

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

... ну и как это звучит по научному?

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

да

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

совершенно верно

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

так я же пишу, что здесь не интерполяция, пока, а просто частоту дискретизации сменил на том же интервале...

Вопрос в том что при интерполяции мы имеем дело с с некоторым количеством отчётов на некотором жестко заданом интервале (в частности 0.512с). А я думаю что при интерполяции длина окна должна меняться.
вобщем нужно менять длинну окна... как-то...

что бы не быть голословным приведу вырезку из диссертации...
на этой странице говорится об особенностях применения ДПФ для анализа 40 гармоник промышленной частоты (50Гц)...

Эскизы прикрепленных изображений