https://habr.com/ru/post/460445/
https://habr.com/ru/post/196374/
В постах вроде была литература для "въезда".

Ну да ДПФ врет при отклонении частоты, но в пределах 45-55 Гц не сильно (около процента на герц). Поэтому у некоторых производителей вы можете увидеть что "дополнительная погрешность при изменении частоты в пределах 45-55 не более 5%", что говорит о том, что результаты Фурье фильтрации не корректируются в зависимости от частоты.

Если хотите точнее, тогда надо измерять частоту и далее вносить коррективы (либо обеспечивать постоянство точек на периоде через ресемплинг, либо интегрировать столько сколько нужно, либо колдовать уже над полученными оценками вектора). Так что фильтр Фурье это самый распространенный способ.

Далее можно использовать МНК и посчитать вектор быстрее чем за период, например, но частотные свойства будут хуже. Можно как в ГОСТе по фазорам использовать треугольную взвешивающую функцию. Погрешность от отклонения частоты будет гораздо меньше, но и времени уже 2 периода нужно.

Если нужна только амплитуда без угла, то RMS или вообще средневыпрямленное значение. Какая нибудь супер быстрая ТО может работать чуть ли не по мгновенным значениям (с небольшой фильтрацией). Все зависит от целей алгоритма.

Посмотрите Шнеерсон Цифровая РЗ, Никитин Цифровая РЗ. Основы синтеза измерительной части микпропроцессорных реле.
Или гуглить phasor estimation прям в гугле или на сайте IEEE. Будет вам миллиард статей.

Спасибо Вам за ответ и рекомендации по литературе, подскажите, а про какой ГОСТ идёт речь ?

В МКПА для ПА. Действующие значения токов и напряжение производится через вычисление гармонических составляющих синусоидального сигнала.
В результате алгоритмы устройств ПА будут иметь весь набор составляющих: акт. и реакт. токи фаз и напряжений, а из них уже вычисляются по необходимости всё что угодно: фазы, мощности, сопротивления, частота ....

Нет.
всё проще.
измеряется текущая частота входного сигнала А=Аm*SIN(wt+fi) и подаётся на генератор, который формирует синус a=am*SIN(wt)  и косинус a=am*COS(wt), на которые умножается  входной сигнал и получим:

A1=0,5Am*aSIN(wt+fi)*SIN(wt)  = 0,5Аm*a[COS(fi) - COS(2wt+fi)]
A2=0,5Am*aSIN(wt+fi)*COS(wt) = 0,5Аm*a[SIN(fi) - SIN(2wt+fi)]
видно, что есть постоянные составляющие и переменные двойной частоты. Переменная отфильтровывается (например ФНЧ Чебышева 2-го порядка, частота среза 65 Гц) и остаются
A1=0,5Аm*aCOS(fi)
A2=0,5Аm*aSIN(fi)
первое - действительное значение, второе - мнимое, кв. кор. из суммы вк. получим величину пропорциональную действующему значению (кор. из 2 и знач. а известны), и так обрабатываются все напряжения и токи, всё.
Дальше уже можно  вычислить весь "зоопарк"

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

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

В этом плане определение частоты по переходам через 0 проще... Выходные значения его проще фильтруются... Единственный недостаток этого способа это длч правильного дифференцировария нужно хранение моментов переходов через 0

посмотрите формулы отфильтровывается сигнал двойной частоты

сам проверял на статике от РЕТОМА, когда отлаживал алгоритм АЧР.
напряжение менял в диапазоне 2...100 В, и частоту от 40 до 60, погрешности были в диапазоне 1..2%. например на 5В - 49,96 ... 50,05 Гц. так же следил за значениями т.к. нужно было ещё за одно отлаживать и алгоритм АОСН. Там как бы вообще ни чего не отловил 1...2% как было заявлено.
С измерениями в динамике конечно хуже, но тут вообще нет ни какого способа или метода выполнить более менее точные измерения.
Например, добился, что мой алгоритм  измерения скорости изменения частоты работает не хуже 2...5% при 10...12 Гц/c., а при скорости в диапазоне 1..3 Гц/c вообще 0,3...0,5%.

с временем расчётов ни чего не скажу, т.к. один тик в МКПА - 8 мс, в это время укладывается всё: измерение-перобразование всех аналоговых сигналов - может быть 32 шт., обработка алгоритмов, приём и выдача дискретов - может быть в сумме 128 шт.
как правило МКПА не загружен и на половину, поэтому всё происходит гораздо быстрей 3...5 мс...