Есть статьи Куликова А.Л., Фальшиной В.А. по поводу подавления апериодики. Проанализированы существующие алгоритмы, предложены свои. Большим плюсом, конечно же, является список литературы, где есть ссылки в том числе и на алгоритмы китайских исследователей. Есть возможность давить конкретную апериодику (с конкретной постоянной времени) или же в неком диапазоне.
Если говорить в целом за хороший новый ЦОС, то, наверное, без адаптивной фильтрации не обойтись. Алгоритм сам настроится на сигнал и выделит оттуда и экспоненты и основную гармонику и субгармоники, если они есть. И в зависимости от изменения сигнала будут меняться коэффициенты фильтра.
Об этом написано в "Методы обработки цифровых сигналов энергосистем" за авторством Антонов, Лазарева, Пуляев (+список литературы!!). В принципе об адаптивной фильтрации (структурном анализе) эти авторы (и не только эти) пишут с конца 80-х годов. Например, статья 1988 года "Адаптивная цифровая обработка входных величин релейной защиты" Лямец, Антонов, Арсентьев. А сейчас в сборниках статей ЭКРА, доступных на сайте, опять же эта тема не забыта. Подробно также есть в книге "Цифровой спектральный анализ и его приложения" Марпл-мл.
Идея, по мне, хороша: моделировал такие фильтры невысокого порядка на сигнале sin+exp. Пока экспонента существенна, алгоритм тратит адаптивную часть на ее подавление, потом, по мере затухания экспоненты, алгоритм превращается в обычный Фурье.
но!
Во-первых, нужно довольно много считать, так как, по сути, необходимо решать переопределенную СЛАУ практически неполного ранга
Во-вторых, в момент смены режима ЭС, пока в фильтр (даже адаптивный) не зайдет сигнал нового режима, не нужно ожидать хороших результатов (в отличие от Фурье, который эти переходные области обрабатывает нормально. Например, траектория Z из режима нагрузки в режим КЗ за 20 переходных мс не будет заскакивать в зону срабатывания). Так что нужно фиксировать момент коммутации и начинать работать уже с сигналом нового режима.
В-третьих, а стоит ли игра свеч?
А так, при асинхронном ходе, алгоритм выделяет из сигнала две частоты, например, 50 и 51 Гц. Можно обе показать на дисплее блока)