Алгоритмы измерения ФЧХ Опции

[maxis]

Есть микроконтроллер, к которому по i2s подключен динамик и микрофон (с одним мастер клоком). Необходимо снять АЧХ и ФЧХ системы ЦАП-динамик-среда-микрофон-АЦП на диапазоне 20 - 20000Гц.
С АЧХ всё более менее понятно, генерим белый шум, от принимаемого и передаваемого сигнала считаем Фурье, делим и строим график.
А вот с ФЧХ вопрос. Первое решение которое пришло в голову - генерим сигнал из синусов с дискретом изменения частоты по одному герцу на период, на входящем сигнале считаем время всех переходов через ноль и соответственно определяем таким образом задержку(фазу) каждой частоты.
Может существуют более быстрые и точные алгоритмы снятия ФЧХ?

[khach]

Или квадартурный самплинг на частоте в 4 раза большей чем частота сигнала, тогда математика определения фазы очень простая, или sine fit (метод наименьших квадратов) для двух синусов -опорного и измеренного, при этом частота- общий параметр, а амплитуды, фазы и сдвиг по постоянке- для каждого синуса свои. Метод 4 или 7 параметров (Four-parameter sine wave fitting, dual sine wave fitting)