function [DFT1,DFT2] = tworealDFTs(x1,x2)
    N = length(x1);
    x = x1 + j*x2;
    DFT = fft(x, N);
    % DFT_circfold = conj(circfold(DFT, N));
    DFT_circfold = circfold(conj(DFT), N);
    DFT1 = 0.5*(DFT + DFT_circfold);
    DFT2 = (0.5*(DFT - DFT_circfold))/j;
end

