Демодуляция 1FSK

[dmt0768]

Добрый день!
При демодуляции 1FSK сигнала по следующему алгоритму:

возникла проблема при склеивании arctg c помощью функции unwrap (Python3):

angle = np.arctan2(y_sin, y_cos)
angle = np.unwrap(angle)

в итоге получилось вот что (после unwrap):

при большем приближении можно заметить явную пилу:

И вот что делать с этой пилой я не знаю, понятно во что она превращается при последующем дифференцировании. Можно, конечно, эту пилу успешно отфильтровать, но правильный ли это выход?

[Бахурин Сергей]

1) если FSK, то лучше делать частотный детектор (рисунок 3 этого раздела)
2) Схема не соответствует коду. По схеме арктангенс отношения, по коду арктангенс 2.
3) Графики приведены в попугаях. Выход арктангенса должен быть от -pi до pi. Функция unwrap будет аккумулировать скачки 2pi и раскрывать периодичность. Из ваших графиков непонятно пила это периодичность или полная фаза меандра, который пропустили через FSK модулятор, потому что значения по шкале непонятно в каких единицах (догадываюсь, что их перевели в вольт по какому-то правилу, но unwrap про это ничего не знает и ищет скачки 2pi).

[dmt0768]

Спасибо, да, не разберусь попробую другой алгоритм.

Поторопился с отправлением, забыл переподписать оси, там ось Y в радианах, конечно, всё после раскрытия периодичности. Производная данного графика должна восстановить меандр.

P.S.
Заменил графики.