function [P, M] = fourier_parts(F)
    P = fftshift(angle(F));
    M = fftshift(abs(F));
end
