Неадекватное отображение фазового спектра

%% Параметры
Tm=0.00025;% Длина сигнала (с)
Fd=1024*1024;% Частота дискретизации (Гц)
Ak=0.5;% Постоянная составляющая (Попугаев)
A1=5;% Амплитуда первой синусоиды (Попугаев)
A2=0.7;% Амплитуда второй синусоиды (Попугаев)
F1=50000;% Частота первой синусоиды (Гц)
F2=100;% Частота второй синусоиды (Гц)
Phi1=0;% Начальная фаза первой синусоиды (Градусов)
Phi2=37;% Начальная фаза второй синусоиды (Градусов)
An=1*A1;% Дисперсия шума (Попугаев)
FftL=2048*1024;% Количество линий Фурье спектра
Epsilon=18000;% Область максимума
 
%% Генерация рабочих массивов
T=0:1/Fd:Tm;% Массив отсчетов времени
 
Noise=An*randn(1,length(T));% Массив случайного шума длиной равной массиву времени
 
 
Signal=A1*sind((F1*360).*T+Phi1);
 
%% Спектральное представление сигнала
FftS=abs(fft(Signal,FftL));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./FftL;% Нормировка спектра по амплитуде
FftS(1)=FftS(1)/2;% Нормировка постоянной составляющей в спектре
%FftSh=abs(fft(Signal+Noise,FftL));% Амплитуды преобразования Фурье смеси сигнал+шум
FftSh=angle(fft(Signal,FftL));% Фазовый спектр
FftSh(1)=FftSh(1)/2;% Нормировка постоянной составляющей в спектре
[C,i] = max(FftS);%Максимум и его аргумент в массиве частот спектра
 
 
%% Построение графиков
subplot(2,1,1);% Выбор области окна для построения
plot(T,Signal);% Построение сигнала
title('Сигнал');% Подпись графика
xlabel('Время (с)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика
subplot(2,1,2);% Выбор области окна для построения
plot(T,Signal+Noise);% Построение смеси сигнал+шум
title('Сигнал+шум');% Подпись графика
xlabel('Время (с)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика
 
%%
F=0:Fd/FftL:Fd/2-1/FftL;% Массив частот вычисляемого спектра Фурье
figure% Создаем новое окно
 
subplot(3,1,1);% Выбор области окна для построения
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
title('Спектр сигнала');% Подпись графика
xlabel('Частота (Гц)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика
 
 
subplot(3,1,2);% Выбор области окна для построения
plot(F(max(1,i-Epsilon):min(length(F),i+Epsilon)),...
FftSh(max(1,i-Epsilon):min(length(F),i+Epsilon)),'-r');% Построение спектра Фурье сигнала
title('Фазовый спектр');% Подпись графика
xlabel('Частота (Гц)');% Подпись оси х графика
ylabel('Фаза');% Подпись оси у графика
 
subplot(3,1,3);% Выбор области окна для построения
plot(   F(max(1,i-Epsilon):min(length(F),i+Epsilon)),...
        FftS(max(1,i-Epsilon):min(length(F),i+Epsilon)),'-r',...
        F(i),FftS(i),'ko');% Посмотрим нагляднее у максимума
title('Спектр сигнала в окрестности максимума');
xlabel('Частота(Гц)');% Подпись оси х графика
ylabel('Ампл.');% Подпись оси у графика