How to get phase difference between two signals?

[Einstein Mcfly]

Hello all. I have two functions in timef_1(t) and f_2(t) composed of the same set of frequencies such that, say, f_1(t) = \sum a_n \cos(\omega t + d_1(t)) and f_2(t) = \sum b_n cos(\omega t) and I would like to find out the value of the phase difference (I've set the phase in f_2 equal to zero) at each point in time.

Can I get this by comparing the complex parts of the two signal's Fourier transforms? What is the "meaning" of the complex part anyway? Whenever I've used the FT in the past, it's been the absolute value that represents the amplitude of the field at a given frequency, so what are the "meanings" of the real and imaginary parts separately? What is the "meaning" of the phase constructed from the arctan of the ratios of the real and imaginary parts?

Thanks for any help you can give me.

 

[berkeman]

Hello all. I have two functions in timef_1(t) and f_2(t) composed of the same set of frequencies such that, say, f_1(t) = \sum a_n \cos(\omega t + d_1(t)) and f_2(t) = \sum b_n cos(\omega t) and I would like to find out the value of the phase difference (I've set the phase in f_2 equal to zero) at each point in time.

Can I get this by comparing the complex parts of the two signal's Fourier transforms? What is the "meaning" of the complex part anyway? Whenever I've used the FT in the past, it's been the absolute value that represents the amplitude of the field at a given frequency, so what are the "meanings" of the real and imaginary parts separately? What is the "meaning" of the phase constructed from the arctan of the ratios of the real and imaginary parts?

Thanks for any help you can give me.

I could be wrong, but it sure seems like your a_n will have to equal your b_n for all n, in order for "phase difference" to have a stationary meaning.

 

[Einstein Mcfly]

I could be wrong, but it sure seems like your a_n will have to equal your b_n for all n, in order for "phase difference" to have a stationary meaning.

Really? I was thinking that if I calculated the phase for f_1(t) by writing the Fourier transform as \sqrt(Re(F_1(\omega))^2+Im(F_1\ometa)^2)e^{i\theta} where \theta_1=arctan(Im/Re) and then the phase difference at time t would just be
d(t)=\theta_1(t)-\theta_2(t). Wouldn't the problem of the differences in amplitudes be removed by either the fact that we're looking at a ratio for each signal or the fact that the arctan is bounded?

 

[Bob S]

Do your two out-of-phase signals look like the signals in Fig. 4 on page 2?

http://www.eleceng.adelaide.edu.au/personal/wlsoong/documents/pebn005BHcurveandironlossmeasurementsformagneticmaterials.pdf

If both signals, including all harmonics, are known to be cosine-like (like your f₁(t) and f₂(t)), , then just measuring the time delay of the zero-crossings is sufficient to determine the phase difference.

Bob S

 

[Einstein Mcfly]

No, they're a lot more complicated than that. They're actually two independently shaped UV laser fields that form an elliptically polarized pulse (with ellipiticity and phase changing with time) when you form the resultant field. I can plot them both together vs. time and see how the polarization changes, but I want to be able to get the numerical value of the phase at each point in time just from the real amplitudes.

 

[Okefenokee]

Are we talking about a phase difference for two signals at one particular frequency? You could do that easily by using integration with cos(f2pi*t) to single out specific frequency components independently in each signal. Then multiply those frequency components and use trig identities to calculate the phase difference.

The FT of a signal has an entire spectrum of frequency components and the phase of a single component can be different at any point on the spectrum. That means the phase difference between the two signals must be a function of frequency. Is that what you want?

 

[Bob S]

The FT of a signal has an entire spectrum of frequency components and the phase of a single component can be different at any point on the spectrum. That means the phase difference between the two signals must be a function of frequency. Is that what you want?

Although there are many frequency harmonics, they are all exactly cosine-like and have the same phase delay, as stated in the OP. Doesn't this imply that if the phase delay is determined for the fundamental harmonic, it is known for all?

Einstein:
Do you really mean

f₁(t) = ∑a_ncos[(nω₀t) + θ_n] or ∑a_ncos[(nω₀{t-t_n})]

with a different phase shift for every harmonic? Please rewrite your equations for f₁(t) and f₂(t).

Bob S

 

[Einstein Mcfly]

Okay, I think what I was trying isn't going to work. Maybe it would help if I was more general.

What I have is two sets of data that are the real amplitudes of two shaped laser fields at each time point. What I want to do is find the phase difference between the two, which pretty much means (I think) that I have to find the phase of each of them. I only have the real part of E(t)=E^{re}+iE^{im} and so I can't construct the phase of either field from the Fourier Transform. What I HOPE will now work is that I can get it from the Kramers-Kronig relation. That is, that because the real part of the field (which I know) isn't independent from the imaginary part (which I don't) that I can use the Kramers-Kronig relation to calculate it.

Is this the correct use of the Kramers-Kronig relations or am I barking up the wrong tree here?

 

[Bob S]

I have used the Kramers Kroning (Dispersion relations) relations and real part sufficiency to calculate the imaginary part in the frequency domain, but not time domain. It may work. There is a complete page of Kramers-Kronig relations for EE circuits in Bode's book Network Analysis and Feedback Amplifier Design (1945).

Bob S