Вычисление фазы Опции

[Alexey_pashinov]

Здравствуйте! Возникли трудности со следующим вопросом.

Иммется антенна доплеровского локатора K-MC4, и для вычисления направления движущегося объекта мне необходимо вычислить разность фаз сигналов I и Q, идущих с антенны. Рассчет делаю в частотной области, пользуясь формулой tg(fi) = Im/Re, т.е. fi = atan2f(Im, Re); разница фаз I и Q в моем случае должна быть 90 градусов, но в результате получаю никак не анализируемые данные. Вопрос можно ли таким образом, переводя в частотную область, правильно высчитать фазу? Где то слышал что так не делают, к тому же stm в arm_math.h для построения АЧХ предоставляется функция magnitude, и она в моем случае правильно вычисляет бин и используя его я вычисляю частоту, а для вычисления фаз такой функции нет. Совпадение ли или таким образом в частотной области действительно фазу не получить? Можете тогда подсказать какие есть методы и что предпочтительней будет в моем случае? (Весь цос реализуется в stm32)

[Alexey_pashinov]

Смотрите, чтобы вычислить частоту моего синуса, я как раз вычисляю какой частоте соответствует найденный бин (погрешность может быть примерно в 29 Гц, учитывая что Fs = 60 кГц, N = 2048 отсчетов). Таким образом получаю частоту синуса и теперь нужно найти его фазу. Ищется ли фаза этой частоты в том же бине что и частота, если пока не рассматривать, что частота не попадает в бин, т.е. не кратна Fs/N?

[petrov]

Смотрите, чтобы вычислить частоту моего синуса, я как раз вычисляю какой частоте соответствует найденный бин (погрешность может быть примерно в 29 Гц, учитывая что Fs = 60 кГц, N = 2048 отсчетов). Таким образом получаю частоту синуса и теперь нужно найти его фазу. Ищется ли фаза этой частоты в том же бине что и частота, если пока не рассматривать, что частота не попадает в бин, т.е. не кратна Fs/N?

А вы посмотрите внимательно как вычисляется один бин DFT, это всего лишь выход полосового FIR фильтра с прямоугольной огибающей импульсной характеристики.