Ну так делайте частоту дискретизации такую чтобы на на 8K выборке помещалось целое число периодов основного тона. Хотите тактовую частоту АЦП меняйте, хотите ресамплинг делайте. Только как-то довольно сомнительно это. ИМХО туже информацию можно извлечь без передискретизации.

p.s. А в дисертациях не обязательно умные вещи пишут...

Укажите какие.
О квалификации Сергиеко А.Б. не Вам судить. Лично с ним знаком.

Да я-то откудова знаю, как по-научному? Рассуждения про растекания спектров ваще из области фантастики. Есть дпф, которое определено для периодической дискретной последовательности периода N. те если выполнить дпф над огрызком синусоиды с периодом, не кратным N, то и получим оценку спектра этого огрызка, а не гармонической последовательности на интервале +-inf. Вот и все растекание...

По теме:
Ознакомился с бумажкой. Потом ознакомился с гостом и приложением. И не очень понял, в чем же проблема? Сделаете Вы передискретизацию, выполните дпф на нужное число точек. Все нужные частоты будут точно совпадать с частотами бинов. По сути вы и измените длину окна. Окно 0.512 c на 16000Гц станет 0.4945...с на 16564,224Гц.

А как объяснить такое.
Читаю:
"... С учётом близкого расположения гармоник на частотной оси ( через каждые 50 Гц) и малых отклонений частоты основной: гармоники (до сотых долей Гц) от номинала требуется высокое разрешение БПФ (вплоть до 0.03 Гц). Таким образом, размер БПФ (длина выборки) превышает 533 тысячи отсчетов прни f=16 кГц. но при выполнении условна отсутствия искажений размер БПФ можно существенно уменьшить - до порвлка 10 тысяч отсчетов при той жe частоте дискретизации)"

Получается не соответствует действительности?

P.S. условие отсутсвия искажений fd*w = N*fосн.
fd - частота дискретизации;
w - некоторое целое число;
N - размер БПФ;
fосн. - частота основной гармоники (49-51Гц)



наконец-то!!! спасибо это и хотел услышать... просто в MATLAB не додумаюсь как это сделать... - изменить длинну окна...
как-то нужно время пересчитать наверное...
мы же изначально на интерполятор пускаем сигнал из отсчётов в окне 0.512 с частотой 16000, преобразую в сигнал с частотой 16564,224 - получается больше отчётов на этом же окне(((
наверное сам принцип "интерполяции" не понимаю...

Сообщение отредактировал TigerSHARC - Dec 1 2009, 15:48

Рабинер Голд Теория и применение цос. Годится не только как базовый учебник, но и как справочник.
Правда, матлабских скриптов там нет. Эт, канешна, серьезный минус...




Нет. О квалификации Сергиенко А.Б. я судить не могу и не хочу. Но как читатель, вполне могу покритиковать его книгу.



Завидую.

Например, меняете частоту АЦП в соответствии и берёте всегда 8192 отсчёта. Тогда "условие отсутствия искажений"
будет выполняться и все гармоники точно лягут на бины ДПФ.

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

Хотя может Фарроу и будет работать за счёт большого оверсамплинга, это вообще-то вопрос. 50гц*40=2 кгц
8 отсчетов период может мало для "прецизионной" интерполяции, но можно пробовать ресамплинг сделать на повышеной частоте, потом прорежать и делать ДПФ (или БПФ)

Да нет все наоборот, нужно получить 8K (== размер FFT) отчетов, из k*8K входных отсчетов (k- коэффициент вроде k=0.99.. 1.01). После ресамплера длина окна постоянна, а фактически, во времени длина окна будет кратна периоду основного тона.

всё верно! Так и надо делать... никакого аппаратного изменения частоты.
мне сам принцип изменения окна временнного при интерполяции интересует.
thermit сказал что окно меняется. Теперь ищу информацию, где доступно написано как именно.

Сообщение отредактировал TigerSHARC - Dec 1 2009, 16:16

Интерполяция в контексте цос - повышение частоты дискретизации. те отсчеты сигнала берутся чаще, и как следствие - их становится больше в единице времени.

Да не нужен Вам никакой ресемплер, наложите Гауссовское окно и посчитайте правильно. Получите очень даже прецизионно. И чему равно Ваше прецизионно? У нас прибор уже работает на этом принципе. Точность по напряжению 0.01% - это полностью, а не только точность вычислений, по частоте и фазе еще существенно выше.

мой алгоритм предполагает реализацию для приборов контроля качества электрической трёхфазной сети в реальном времени.
Интересно:
1) у вас - это у кого? (если не секрет)
2) 0.01% Это с чётом погрешности первичных датчиков?
3) какие параметры системы(чатсота дискретизации, длинна БПФ...)?

Единственное что приходит на ум при применении сглаживающих окон это введение како-то поправки...


Но как вы ответите на такой тезис:
"... несмотря на рекомендации, приведенные в ГОСТ 13109-97 (обязательное приложение Е6), по выбору типа и длины окна, четких правил выбора этих параметров не существует, А поскольку любое окно искажает спектр полезного сигнала» так как произведение функций во временной области соответствует свертке их образов в частотной области, то погрешность вычисления спектра зависит от длины окна и вида спектра анализируемого сигнала. Выбор типа окна должен основываться на предполагаемом спектральном составе сигналов, но спектральный состав сигналов от датчиков тока и напряжения может изменяться во времени, что определяется характером включенной нагрузки. Поэтому единственный правильный вариант - обеспечение условия отсутствия искажений при применении прямоугольного окна. Тем более что это не противоречит ГОСТ 13109-97"

Сообщение отредактировал TigerSHARC - Dec 1 2009, 18:00

Привет SHARK. По Твоей тематике можешь посмотреть работы VU Jiekang (Ziejiang Univ.) . Сейчас работает в Америке в энергетической компании. Может и прибор уже создал. Никакого БПФ! Я использую его метод, но в другой области.

Книга бесусловно хорошая, но есть проблемы с примерами.

Мы знакомы???)

Что за метод?
Где найти его материалы??? поиск ничего не дал(мне по крайней мере)...

Прошу прощения за фамильярность. Навалом ссылок.
Вот журнал: IEEE transactions on power delivery ISSN 0885-8977, 2005, vol. 20, no1, pp. 366-374 [9 page(s) (article)] (22 ref.)
Еще: http://www.springerlink.com/content/r113673647v40251/
Там и почта его есть.
Похоже это та статья. В бесплатном виде пока не нашел.
Мой экземпляр погиб вместе с жестким, на бумаге осталось. Может подскажете где эти transactions почитать подешевле можно?