У меня есть 129 образцов какого-то сигнала. Я должен выяснить, что это за сигнал (из каких компонентов он состоит) и попытаться восстановить его. поэтому после того, как я сделал БПФ, я хотел найти фазу сигнала, поэтому я пытаюсь:
phs=unwrap(angle(y));
pha_wek=(0:length(y)-1)'/length(y)*129;
plot(pha_wek,phs)
Результат я получил, это не то, что я ожидал. Можете ли вы указать мне, что не так?
Сюжет выглядит совершенно правильно для меня, в зависимости от исходного сигнала. Разрывы фазового угла возникают из-за пересечения нуля амплитудой. Однако скачки фазы меньше, чем число пи (что можно было бы ожидать, если бы входной сигнал был изначально реальным и был подвергнут оконной обработке при применении БПФ). Итак, я предполагаю, что этот сигнал представляет собой некоторую форму сложной модуляции с некоторыми дискретными частотами.
Обычно просмотр фазы говорит вам о сигнале гораздо меньше, чем его амплитуда. Это важно для реконструкции сигнала, но его гораздо сложнее интерпретировать (особенно потому, что он подвержен резким изменениям при простом сдвиге во временной области).
Вы могли бы лучше построить фазу без использования команды unwrap. В этом случае проблема заключается в команде unwrap — ее цель состоит в том, чтобы ограничить скачки по фазе меньше, чем пи, поэтому конечная фаза составляет около -4pi.
Я не уверен, в чем заключается ваша проблема, но неясно, что вы имеете в виду под реконструкцией образцов. Они у вас уже есть, так в чем проблема? (Вас просят реконструировать интерполированные образцы)?
