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