|
|
Ну, во-первых, ничего особенно антинаучного в использовании компаратора и т.д нет. Вопрос только в точности.
Во-вторых, активная мощность
Pa=(1/T)*Integral(u(t)*i(t)*dt)=P*COS(fi).
Имеются ввиду синусоидальные напряжение и ток. Интеграл, соответственно, по периоду. P - полная мощность=Um*Im/2. Um,Im- амплитуды. Т.е. активная мощность находится напрямую, ну может с предварительной фильтрацией напряжения и тока одинаковыми фильтрами.
Это я к тому, что может проще не COS(fi) вычислять, а амплитуды (или действующие значения) напряжения и тока, т.е. фактически полную мощность?
Кому-то уже писал, повторюсь.
Простой способ оценки действ. значения.
Сумма квадратов всех отсчетов синусоиды, деленная на число отсчетов.
Действующее=корню из полученного значения. Число отсчетов лучше, конечно, брать таким, чтобы на этом числе укладывалось целое число периодов.
Еще способ. Позволяет определить как амлитуду, так и фазу.
Множите напряжение и ток на sin и cos частоты сети.
Фильтруете каждый из полученных 4 сигналов ФНЧ. В качестве ФНЧ может оказаться достаточным сумматора, то бишь интегратора.
В результате получаете Re,Im напряжения и Re,Im тока.
Корни из суммы квадратов Re и Im дадут оценку амплитуд (с учетом, естественно, коэф-та передачи фильтра), а arctg(Im/Re) - оценку фаз.
Косинус разности фаз напряжения и тока и будет косинус fi.
Наконец можно все определить, взяв Фурье над напряжением и током.
В спектре легко определяются и амплитуда и фаза нужной частоты.
Правда могут быть маленькие неприятности, связанные с некратностью частот дискретизации и сигнала.
Что такое преобразование Гильберта, а можно конечно и с его помощью, я кому-то отвечал. Вроде бы на предыдущей странице.
Только большого смысла для Вашей задачи в нем не вижу.
По сути оно позволяет определить из исходного сигнала ему сопряженный. Т.е. те же самые Re и Im.
А сведется этот преобразователь Гильберта либо к КИХ фильтру, возможно довольно большого порядка, либо к тому же Фурье.
E-mail: 
info@telesys.ru
