Как объясняется увеличение амплитуды фазового спектра на положительной частоте, при изменении начальной фазы

function [] = ...
    matlab2()
 
        a = input('Частота первой синусоиды: ');
        
        maxFrequency = max([a]); % Максимальная частота
        Fs = maxFrequency*10 + 1000; % Частоты дискретизации сигнала
        Tm = 100; % Длительность сигнала (в ед. времени)
        T = 1/min([a]); % Наибольший период
        time = 0 : 1/Fs : Tm; % Массив времени
        y1 = redef_signal(a, time);
 
            NFFT = Fs*Tm; % Количество спектральных линий
            frequency = -Fs/2 : Fs/NFFT : Fs/2-Fs/NFFT; % Массив частот
            Y1 = imag(fftshift(fft(y1, NFFT))) * -1 + fftshift(fft(y1, NFFT)); % Сдвиг результата БПФ
            Y1 = 2*Y1./NFFT;
 
            figure_handler = figure(1);
            set(figure_handler, 'Position', get(0,'ScreenSize'), 'Color', 'w'); 
            clf(figure_handler)
 
            subplot(1, 2, 1), stem(frequency, Y1,'Marker','none'), grid; % График амплитудного спектра
            xlim([-maxFrequency*2 maxFrequency*2]); % Ограничить вывод по оси частот
            ylim([-3 3]); % Ограничить вывод по оси амплитуды
            xlabel('Частота, Гц'); % Название оси X
            ylabel('Амплитуда'); % Название оси Y
            title('Спектр сигнала'); % Название графика          
            
            time = 0 : 1 / (Fs * 10) : Tm; % Массив времени
            y1 = redef_signal(a, time);
            
            subplot(1, 2, 2), plot(time, y1, 'r'), grid; % Вывести график сигнала
            xlim([0 T*5]); % Ограничить вывод по оси времени 5 периодами сигнала
            ylim([-4.5 4.5]); % Ограничить вывод по оси амплитуды
            xlabel('Время'); % Название оси X
            ylabel('Амплитуда'); % Название оси Y
            title('Сигнал'); % Название графика
 
end
 
function [y1] = redef_signal(a, time) %Функция переопределения исходного чистого сигнала
    y1 = sin(a*2*pi*time); % Генерация сигнала
end

function [y1] = redef_signal(a, time) %Функция переопределения исходного чистого сигнала
    y1 = sin(a*2*pi*time+pi/4); % Генерация сигнала
end

Y1 = abs(fftshift(fft(y1, NFFT)))