Ну 3 раза имелось в виду вычисление действующего, а так конечно на период порядка 200 измерений, на 100 мс я тоже делал по такому принципу, но там пришлось синхронизироваться с сетью.

5-6 вычислений RMS
при использовании встроенного АЦП

Кстати, а зачем корень вычислять? Проще уставку возвести в квадрат и с этим значением и сравнивать насчитанное.

Дык и не пользуйтесь плавающей арифметикой. Возьмите fast integer square root
или вот такой, страшно математисский

Касательно "fast integer square root" - может кто-то доступно объяснить подоплёку отого чЮда ?
Всмысле вроде метод нъютона - это где производные, а откуда там взялись эти махинации с битами ?? .....

На первый взгляд - это просто метод последовательных приближений. Как будто в АЦП, только с квадратом внутри.

Если махинации с битами - то это уже не Ньютона :-)
Вот, нашёл в гуглогрупсах старое
(там цитирования надо поразворачивать)
О, тут не вся тема, зато чистым текстом http://files.kpv.nnov.ru/sqrt16b.txt

А для компенсации постоянки отдельно вычислять среднее арифметическое и вычитать его из каждого отсчёта не нужно (заодно массив отсчётов пропадает, экономия ОЗУ).
Делается это так:
Накапиливается в двух накопителях - сумма квадратов и сумма отсчётов.
В конце цикла измерения вычисляем средний квадрат без компенсации постоянки и среднее значение (постоянку).
Из среднего квадрата вычитаем квадрат среднего.
Всё, получили среднее квадрата отклонения от постоянки.

Т.е. (тут нормально формулы не написать - нет ни тегов sup/sub, ни надчерка)
sqr() - это возведение в квадрат

average( sqr( x - average(x) ) ) = average( sqr(x) ) - sqr( average( x ) )

При необходимости из этого извлечь корень.


О, нашёл в RU.EMBEDDED ещё более старую тему, но там есть С-шный 8-битниковый вариант, писался под 8051.

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

Batman
Подскажите где взять инфу про "пояс Роговского"???

Это совершенно правильно. Но только в том случае, если нас устраивает получение значения только 1 раз за период (20 мс при 50 Гц). Чтобы получать значение вдвое чаще - придётся завести ещё 2 накопителя и вычислять всё это со сдвигом в 10 мс. И т.д., по 2 накопителя на каждое дополнительное получение значения за период. Если нужно вычислять значения с максимально возможной частотой (с частотой осчётов АЦП), то эффективнее использовать FIFO отсчётов, FIFO квадратов отсчётов, накопитель отсчётов и накопитель квадратов отсчётов.
При получении очередного отсчёта АЦП:
1. Прибавлять его к накопителю остчётов, квадрат его прибавлять к накопителю квадратов отсчётов.
2. Значения из хвостов соответствующих FIFO вычитать из соответствующих накопителей.
3. Заносить новые значения в головы соответствующих FIFO (при этом значения из хвостов FIFO затрутся).
Ну а дальше как у вас - из накопителя квадратов отсчётов вычитаем квадрат накопителя отсчётов.

Поскольку обычно постоянная составляющая подмешивается осознано, то куда проще ее завести на отдельный вход АЦП и мерять отдельно, а потом вычитать ее из сигнала.

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


Где, где... в ГуглЕ Вот первая же ссылка http://www.phys.nsu.ru/electricity/Labwork4-3.html

Среднее значение синуса равно нулю.

Вы в курсе, что петля Роговского дифференцирует сигнал ?
Я эти проблемы уже решал, поэтому в курсе.


Не иначе как блок защиты к автоматческому выключателю делаете.


Без теории тут никак не обойтись.


А если подумать, то корень не понадобится вообще.


Точно.


Все правильно.
Еще одна мысль. Можно RMS считать не за 20 мс, а на за 10 мс, как все нормальные люди делают.

Я тако понял, что нужно еще выводить текущее значение тока. Так что одним сравнением не обойтись.



Я бы предложил вести накопитель для среднего арифметического с обновлением через несколько полупериодов. И вычитать его из каждого измерения и возводить в квадрат. Насколько я понимаю, то тут постоянная составляющая будет достаточно "постоянна". Массивы тут по поределению не нужны.

RMS можно получить выполнив Дискретное Преобразование Фурье для первой гармоники входного сигнала и умножив ее на 0,7. Делая Фурье вы автоматом избавляетесь от постоянной составляющей, а получая амплитуду первой гармоники (50Гц) вы автоматом избавляетесь от влияния несинусоидальности входного сигнала и можете смело умножать на 0,7.

Кстати, есть еще и интегральное решение True RMS-to-DC Converter от Analog Devices, например AD636, AD737. На выходе дает RMS входного сигнала в виде постоянного напряжения.