Товарищи ученые, доценты с кандидатами... Пока мало что понятно, но я стараюсь.
Реь идет о электрораспределительныой сети.
В сети всегда есть основная гармоника и гармоники до 13-й, но основная всегда около 50 Гц, например 45-55 Гц.
Все гармоники в сумме всегда значительно меньше основной.
50 Гц - 90%
100 Гц 0%
150 Гц 3%
200 Гц 0%
250 Гц 2%
...
Хочется выделить основную гармонику (цифровым фильтром) и измерить программно её частоту с точностью 0,01 Гц.

Что я могу:
С помощью ДПФ посчитать амплитуду основной гармоники.
С помощью БПФ посчитать амплитуды всех гармоник.
С помощью цифрового фильтра выделить полосу около 50 Гц и получить отфильтрованную выборку из 64-х точек.

Почему бы Вам не поинтересоваться, как это реализовано в регистраторах показателей качества электроэнергии?
Парма, к примеру.

Надо еще учесть, что в сети всегда присутствует АМ и ЧМ, фликкер и проч. Откройте ГОСТ 13109 и забудьте про цифровой фильтр. Нужен образцовый сигнал.

Pridnya

Выделяйте комплексным полосовым КИХ фильтром интересующий диапазон, АЧХ в полосе пропускания должна быть как можно более равномерная, гармоники(в том числе на отрицательных частотах) должны быть задавлены как можно сильнее, затем вычисляете несколько бинов ДПФ, далее по статье считаете частоту:

http://electronix.ru/forum/index.php?s=&am...t&p=1141831

Советую модельку в симулинке сделать и погонять с различными параметрами.

Удалено, не понял задачи и сморозил чушь.

Сообщение отредактировал Милливольт - Sep 10 2015, 16:02

Обязательно поинтересуюсь, спасибо за ссылку на прибор, который вы считаете образцом.
В сети много чего есть (в неизвестно какой электросети), они же все разные.
Упоминаний о наличии АМ и ЧМ фликкера сетях 6-10 кВ нигде не встречал, даже слово такое впервые слышу,
никогда не нужно было бороться с АМ и ЧМ фликкером в электросетях.
Исходя из используемой вами терминологии, ссылки на ГОСТ (электромагнитная совместимость) и предложением отказаться от цифровой обработки сигналов, предположу, что вы являетесь специалистом в области аналоговых систем связи.
Картинки из ГОСТ 13109 "колебания напряжения", "временное перенапряжение" (в вашей терминологии АМ-фликкер) и "перепады частоты сигнала"
(такой картинки нет, но в вашей терминологии ЧМ-фликкер), это как раз то, что мне нужно и чем я занимаюсь.
Перепады напряжений могут быть в диапазоне от 0 до 2-х кратного номинального, а отклонение частоты не более 10%.
Чтобы все это считать, используется аппаратная часть (в первом моем сообщении указано), микроконтроллер с ядром ARM и модулем FPU осуществляет обработку данных (это о доступной вычислительной мощности, ну, то есть, не компьютер под Windows с MatLab-ом это все обсчитывает).
ГОСТы на ЭМС необходимы только для сертификации готового изделия на ЭМС (как работает изделие в условиях помех в питающей сети ~220 вольт, помех в линиях дискретных входов, линий связи...), анализируются и регистрируются (всякие там аварийные журналы с параметрами аварийных режимов, векторные диаграммы токов и напряжений основной гармоники, ток нулевой, прямой, обратной последовательности основной частоты, напряжение нулевой последовательности основной частоты, углы между всеми этими векторами, среднеквадратичное значение суммы верхних гармоник тока нулевой последовательности...осциллограмма аварийного режима для аналоговых и цифровых входов в формате COMTRADE). Поэтому никак без цифровой фильтрации нельзя, везде фильтрация, свертка...И кстати, качество электроэнергии - это не основная функция, которая меня интересует.


Спасибо за информацию! Чувствуется, что человек понимает о чем пишет. Приятно читать.
Вот симулинком не пользовался никогда. Наверное, придется попробовать. И за симулинк спасибо!


В сетях 6-10 киловольт в основном присутствует основная гармоника (50 Гц), поэтому по ней и работают. Это в сетях 0,4 киловольт (внутридомовые...) гармоник полно...
Вариант синхронного фильтра посмотрю. Тоже не встречал и не пользовался таким фильтром.

Сообщение отредактировал Pridnya - Sep 10 2015, 07:13

Прошу не приписывать мне то, что не было сказано.


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


Понятия "амплитудная модуляция", "частотная модуляция", "доза фликкера" применяются раздельно. Речь о борьбе с ними не шла.


ГОСТ 13109-97 вообще-то называется ЭЛЕКТРИЧЕСКАЯ ЭНЕРГИЯ. НОРМЫ КАЧЕСТВА ЭЛЕКТРИЧЕСКОЙ ЭНЕРГИИ В СИСТЕМАХ ЭЛЕКТРОСНАБЖЕНИЯ ОБЩЕГО НАЗНАЧЕНИЯ, то есть речь никак не об ЭМС. Терминология взята в том числе из него.

Предложение отказаться от ЦОС не вносилось, а было предложение, альтернативное использованию цифрового фильтра. И с областью моих интересов Вы промахнулись.


Это Ваша терминология.

Прошу прощения, если вам показалось, что я в своих рассуждениях что-то вам приписал. Проблема восприятия информации.
Вот и теперь мне кажется, что вы сотрудник С.Пб. фирмы Парма и при первой возможности стараетесь отослать к приборам, к которым вы имеете отношение, выдать их за образец (можно же было и к импортному оборудованию отослать, а у них там фирмы по 100 лет и цифровая обработка сигналов у них зародилась, и сигнальные процессоры...а не у нас). Может у меня что-то с головой, но мне так кажется. Извините. А иначе зачем было отсылать к конкретному изделию конкретной фирмы из вашего города и предлагать мне забыть про цифровую фильтрацию (я так понял, что вы даете мне установку, как Анатолий Кашпировский когда-то "забудьте про...иди за нами, сынок"), а я без неё никак. Не сочтите меня мастером пустословия, все мы разные. Еще раз прошу прощения.
Если будут конкретные предложения с ссылками к алгоритмам, а не к готовым изделиям и их техописаниям, то с удовольствием ознакомлюсь.
Уже не раз сталкивался с тем, что расписан способ один, а используется круче.

http://www.parma.spb.ru/catalog/equipment/...tvo/parma_rk301
Что-то как-то не ярко, похоже очень старая разработка:
угла сдвига фаз между каналами напряжения (только РК03.01ПТ)
Оснащены интерфейсом Centronics для подключения матричного принтера и интерфейсом RS-232 для подключения ПК
Частота f Гц от 45 до 55 Гц ±0,02 Гц (абсолютная, при интервале усреднения 20 секунд)
4.4.2.1 Регистратор обеспечивает измерение текущих значений ПКЭ:
− частоты входного сигнала f от 40 до 70 Гц;
...
Таблица 1 – Нормируемые метрологические характеристики регистраторов
Частота f Гц от 45 до 55
...
4.4.2.1 Время установления рабочего режима – не более 60 с.
4.4.2.2 Потребляемая мощность регистратора не более 15 ВА.

Я считаю фазы всех векторов, направление мощности...
Использую изолированный интерфейс USB, Ethernet, т.е. к любому ноутбуку без всяких переходников могу подключиться.
А мне нужно круче чем ±0,02 Гц, чтобы 50,00 Гц "показометром" показывать (без усреднения за минуту), ну, т.е., чтобы быстро среагировать на изменение частоты.
Измеряем в диапазоне 40-70, а погрешность нормиркется в диапазоне 45-55 Гц.
И куда ж там 15 Вт на энергопотребление уходит, микрокомпьютер что-ли засунули с крякнутым MatLab-ом.

Не годится.

Сообщение отредактировал Pridnya - Sep 10 2015, 09:23

Парма - это, как раз, образец безграмотной разработки.

Стандартный способ точного определения частоты основной гармоники сети такой (примерно).
1. Выборка с частотой порядка 25 кГц точек (24 бита) за, примерно, 300 мс (это параметры реального прибора класса 0.05 с частотным диапазоном порядка 8 кГц).
2. На выборку накладывается окно Гаусса с аккуратно подобранными параметрами. Делается FFT (само FFT - float, обработка - double). При этом провал бин между соседними гармониками оказывается порядка 80 дБ.
3. По МНК (парабола или, лучше, гаусс, по 5 - 7 соседним бинам с весами) определяется частота максимума. Если она сильно смещена относительно центра выбранных бин - операция повторяется с коррекцией центра.

При этом точность определения частоты основной гармоники не хуже 0.005 Гц даже при 10% гармониках.

Спасибо!!! Буду пробовать сам метод, не смотря на то, что у меня АЦП 16 бит и частота выборки 3200 Гц (могу поднять до 6400).

Сообщение отредактировал Pridnya - Sep 10 2015, 11:35

Кажется, кажется.. Еще у меня мультиметр Fluke, а осциллоскоп простенький LeCroy. Давайте определяйте по приборам, какой из этих фирм я сотрудник?


Вот и выросло поколение воинствующих инженеров, не отличающих точность от разрешающей способности, мгновенную частоту от частоты основной гармоники

Тепрь-то можно написать хоть про Снеговика с ведром на голове, на содеянное это не повлияет. ))

На данном этапе можно любую терминологию использовать, лишь бы Src и Dst понимали друг друга, а для других - хоть белый шум. Кстати, единого языка с единой терминологией для нучного мира так и не создали.

Сообщение отредактировал Pridnya - Sep 10 2015, 19:58

Точность и разрешающая способность - это совершенно разные вещи, указанная цифра дает именно значение погрешности, т.е. среднеквадратичное отклонение. Разрешающая способность на 2 порядка выше.

Понятия "мгновенная частота" не существует.

При периодическом сигнале частота основной гармоники совпадает с измеренной за любые несколько периодов. Иначе сигнал периодическим ( в строгом смысле этого слова) не является и говорить о частоте основной гармоники бессмысленно.

Дарю. Пользуйтесь.

Спасибо, конечно, но это не метрологическое, а чисто теоретическое понятие, которое применяется только на чистом синусе (ну, в крайнем случае, при узкополосном синале, не содержащем гармоник). Непосредственному измерению не подлежит, поскольку определяется через другой, опять же, непосредственно не измеряемый параметр - скорость изменения фазы, которую опять же невозможно определить по отсчетам АЦП, поскольку, даже при моночастотном сигнале, его амплитуда на момент измерения неизвестна. Для понимания этого полезно почитать литературу, указанную в списке литературы под номером [3] в вашей же ссылке (Финк Л. М. Сигналы, помехи, ошибки…).

Мгновенную частоту можно оценить только задним числом по достаточно большой выборке, но, при этом, само это понятие теряет смысл. И что вы будете делать если есть несколько гармоник? Как вы собираетесь определить понятие фазы такого сигнала?

Именно для того, чтобы избежать этого гемора, и используется метрологически корректное понятие частоты основной гармоники периодического сигнала. Но и оно не является полностью корректным, поскольку есть случайный шум, помехи (частично случайные), время измерения ограничено, сигнал существует не бесконечно долгое время... Но это уже погрешности второго порядка.

f(t) = d/dt(arg(y(t))), простейший частотный дискриминатор, известный со времён старого доброго ЧМ, нет? )

А ещё оказывается в современном ЦОС есть целый класс алгоритмов, т.н. DIFM, например
 GE_47.pdf ( 358.18 килобайт ) Кол-во скачиваний: 386

Вполне себе практические вещи делают.