Необходимо из синусоидального сигнала выделить гармонику 50Hz, или наоборот удалить. БИХ фильтры ВЧ, режекторный - составляющая приобретает сдвиг по фазе. Вопрос в том, как получить фильтр который эти 50Hz в результате сдвинет на 5ms (косинус, который мне понадобится для дальнейших вычислений). Т.е. чтобы зная эту задержку, иметь возможность параллельно сдвигать основной сигнал, и далее выполнять действия между ним и отфильтрованным сигналом.

Подскажите какой тип фильтра можно использовать? С помощью чего это все можно хорошо промоделировать, чтобы определить ГВЗ в 50Hz?
Простите за нубство, в ЦОС пока новичок.

может вот - http://www.dsplib.ru/content/allpass/allpass.html
всепропускающий фильтр

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

А полное название книги?

Вы не поверите "Цифровая обработка сигналов"

герцелем просто посчитайте амплитуду и фазу 50Гц, потом вычитайте из основного сигнала и двигайте как хотите.

Скажите, а если сигнал окажется не 50, а скажем может плавно и медленно изменяться в пределах 45-55 Hz ?

Как в таком случае поступают, ведь насколько понимаю, получим растекание спектра и искажение амплитуды?
Похоже придется постоянно измерять частоту, а уже потом герцелем, при чем как-то в realtime его перестраивать, что скорее всего не легко.

Есть ли какие-то методы, которые позволят избежать искажения амплитуды, если частота "отойдет" от предполагаемых 50Hz ?

сделайте окно меньше чем характерное время изменения, считайте в нескольких точках от 45 до 55, ищите максимум и на этой частоте смотрите фазу.

Вариантов может быть несколько.
На каком интервале вы считаете амплитуду (т.е. сколько периодов частоты 50Гц подвергается анализу)? Какую частоту дискретизации используете?

Моделировал Герцеля с динамическим пересчетом спектральных отсчетов на каждом такте дискретизации, использовал Fs=50Khz, N=4000.
Посмотрел sqrt(Re^2+im^2), почему-то возникают "пульсации", это так и есть, или я в чем-то ошибся?

ФАПЧ с астатизмом второго порядка даст вам синусоиду, которая точно совпадает по фазе и частоте со входным сигналом.

...

Склоняюсь к Герцелю, но пока не понял насчет пульсаций, которые на графике при 50Hz.
Попробовал изменить входной сигнал всего на 0.1Hz, получил приличное искажение амплитуды ну и сдвиг по фазе естественно получил.
А получать спектр для гармоники с "небольшим" шагом и выбирать по максимальной амплитуде не получится, потребуется существенно увеличить число точек или уменьшить частоту дискретизации.



Полученная синусоида ведь не совпадет по амплитуде с основной гармоникой?

точное значение частоты через автокорреляцию будет проще найти.

или просто наименьшими квадратами ищите в сигнале A*Sin(2*Pi*f*t+phi) подгоняя амплитуду частоту и фазу.