CFO OFDM Опции

[Ivan55]

Всем доброго дня!
есть вот такая вот инфа:

получают они вот такие красивые картинки:

решил написать код и посмотреть:

CODE

clear all; close all

M=4;
CFO = 1.1;
Nfft = 32;

H_QAM4=modem.qammod('M',M,'PhaseOffset', 0, 'SymbolOrder',...
'binary', 'InputType', 'bit');

In_Data=randi([0 1],Nfft*log2(M),1);
Mod_QAM=modulate(H_QAM4, In_Data);
scatterplot(Mod_QAM); grid on;
y = ifft(Mod_QAM');

nn=0:length(y)-1;
x = y.*exp(1i*2*pi*CFO*nn/Nfft);

figure(2)
plot(angle(x)); hold on; grid on;
plot(angle(y), 'r');
figure(3);
plot(angle(x)-angle(y), 'k'); grid on;

но картинки такими красивыми не получаются, из-за перехода на pi появляются случайным образом провалы в angle(x)-angle(y)
у них же вся фаза расположена в пределах трех и за pi не переваливает

в чем у них шаманство?