Phase shift in domain from FFT for multiple signals

I have multiple signals that have the peak at different locations so I want to align them up accordingly and then add them together.

Let say the signal A has dimension of 201 x 1, B has the same dimension as well as C, and so on. Each signal has peak at different locations so I want to shift them so that I can add them. Note that data is noisy and the phase shift is unknown.

I was tried to use FFT function from matlab.

I combined A,B,and C in one single matrix that has dimension off 201 x 3, which I called D.

% Output fast fourier transform for each signal

 X= FFT(D, [], 1),

% magnitude

 X= abs(X)

 and then I add them up.

To be honest, I dont know whether it is wrong or right. I wonder if anyone is familiar with this problem.

Thanks.