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

[maxis]

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

[andyp]

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

Можно использовать широкополосный сигнал с хорошими автокорреляционными свойствами для получения оценки импульсной характеристики, скажем, оптимальной в среднеквадратическом смысле. АЧХ и ФЧХ можно будет получить обратным преобразованием Фурье от ИХ.

[maxis]

Можно использовать широкополосный сигнал с хорошими автокорреляционными свойствами для получения оценки импульсной характеристики, скажем, оптимальной в среднеквадратическом смысле. АЧХ и ФЧХ можно будет получить обратным преобразованием Фурье от ИХ.

Ну импульсную характеристику можно вероятно получить как реакцию на дельта-функцию Дирака. При этом частота дискретизации у нас 48кГц, так что её спектр будет равномерный до 22кГц, тут всё нормально. Но если после мк будет стоять фильтр с бесконечной импульсной характеристикой то по идее импульсная характеристика так же будет бесконечной. Что мы получим в этом случае? Просто меньшую точность?

Не напомните как именно получить АЧХ и ФЧХ из ИХ?

[petrov]

Ну импульсную характеристику можно вероятно получить как реакцию на дельта-функцию Дирака. При этом частота дискретизации у нас 48кГц, так что её спектр будет равномерный до 22кГц, тут всё нормально.

С дельта-функцией возникает проблема формирования импульса достаточной энергии. Лучше использовать CAZAC последовательность, которая позволяет точно измерить конечную ИХ.

Но если после мк будет стоять фильтр с бесконечной импульсной характеристикой то по идее импульсная характеристика так же будет бесконечной. Что мы получим в этом случае? Просто меньшую точность?

Вот и проверьте, обрежьте ИХ БИХ фильтра, постройте ЧХ полученного КИХ фильтра, обрежьте дальше повторите и сравните.


Ещё вариант решения - по шуму адаптивным фильтром идентифицировать.