MATLAB: Extract phase of a signal using fft

anglefftphaseunwrap

I have two cerebral signals represented by a 2×10000 matrix and registered for 10 minutes. I have to extract the phase of the signals and calculate the average phase displacement using FFT. Can someone help me with the matlab code? I know I can use the functions "angle" and "unwrap" but I don’t know how to write the code. Thank you in advance!