Polarization Calibration: Difference between revisions

From EOVSA Wiki
Jump to navigation Jump to search
Line 11: Line 11:
</center>
</center>
In terms of autocorrelation powers, we have the 4 polarization products XX*, YY*, XY* and YX*, where the * denotes complex conjugation.  The quantities RR* and LL* are then
In terms of autocorrelation powers, we have the 4 polarization products XX*, YY*, XY* and YX*, where the * denotes complex conjugation.  The quantities RR* and LL* are then
 
<center>
:<math>
:<math>
\begin{align}  
\begin{align}  
Line 18: Line 18:
\end{align}
\end{align}
</math>
</math>
 
</center>
One problem is that there is generally a non-zero delay in Y with respect to X.  This creates phase slopes in XY* and YX* from which we can determine the delay very accurately.  As a check,
One problem is that there is generally a non-zero delay in Y with respect to X.  This creates phase slopes in XY* and YX* from which we can determine the delay very accurately.  As a check,
 
<center>
:<math>
:<math>
\begin{align}  
\begin{align}  
Line 27: Line 27:
\end{align}
\end{align}
</math>
</math>
 
</center>
For completeness:
For completeness:
 
<center>
:<math>\begin{align}  
:<math>\begin{align}  
Stokes \, Q &= XX^* - YY^* \\
Stokes \, Q &= XX^* - YY^* \\
Line 37: Line 37:
\end{align}
\end{align}
</math>
</math>
 
</center>
When I plot the quantities I, V, R and L as measured (Figure 1) for geosynchronous satellite Ciel-2, the results look reasonable, except that there are parts of the band where R and L are mis-assigned, and others where they do not separate well.
When I plot the quantities I, V, R and L as measured (Figure 1) for geosynchronous satellite Ciel-2, the results look reasonable, except that there are parts of the band where R and L are mis-assigned, and others where they do not separate well.


Line 43: Line 43:


The problem is that residual phase slope of Y with respect to X, caused by a difference in delay between the two channels.  This can be seen in the upper panel of Figure 2, which shows the uncorrected phases of XY* and YX*.  To correct the phases, the RCP phase was fit by a linear least-squares routine, and then the phases were offset by &pi;/2 for both XY* and YX* according to:
The problem is that residual phase slope of Y with respect to X, caused by a difference in delay between the two channels.  This can be seen in the upper panel of Figure 2, which shows the uncorrected phases of XY* and YX*.  To correct the phases, the RCP phase was fit by a linear least-squares routine, and then the phases were offset by &pi;/2 for both XY* and YX* according to:
 
<center>
:<math>\begin{align}
:<math>\begin{align}
XY^{*'} &= XY^*e^{-i(\phi(v) + \frac{\pi}{2})} \\
XY^{*'} &= XY^*e^{-i(\phi(v) + \frac{\pi}{2})} \\
Line 49: Line 49:
\end{align}
\end{align}
</math>
</math>
 
</center>
where &phi;(v) is the phase fit by the linear function. The corrected phases are shown in the lower panel of Figure 2.
where &phi;(v) is the phase fit by the linear function. The corrected phases are shown in the lower panel of Figure 2.


Line 55: Line 55:


When the corrected (primed) quantities are used in
When the corrected (primed) quantities are used in
 
<center>
:<math>\begin{align}
:<math>\begin{align}
RR^* &= XX^* - iXY^{*'} + iYX^{*'} + YY^* \\
RR^* &= XX^* - iXY^{*'} + iYX^{*'} + YY^* \\
Line 61: Line 61:
\end{align}
\end{align}
</math>
</math>
 
</center>
the resulting channel separation is shown in Figure 3, which shows pretty much perfect channel separation, with the correct polarity.
the resulting channel separation is shown in Figure 3, which shows pretty much perfect channel separation, with the correct polarity.



Revision as of 18:38, 27 September 2016

Linear to Circular Conversion

At EOVSA’s linear feeds, in the electric field the linear polarization, X and Y, relates to RCP and LCP (R and L) as:

In terms of autocorrelation powers, we have the 4 polarization products XX*, YY*, XY* and YX*, where the * denotes complex conjugation. The quantities RR* and LL* are then

One problem is that there is generally a non-zero delay in Y with respect to X. This creates phase slopes in XY* and YX* from which we can determine the delay very accurately. As a check,

For completeness:

When I plot the quantities I, V, R and L as measured (Figure 1) for geosynchronous satellite Ciel-2, the results look reasonable, except that there are parts of the band where R and L are mis-assigned, and others where they do not separate well.

Figure 1: Amplitudes of Stokes I, Stokes V, RCP and LCP for Ciel-2 satellite. The red vertical lines show the center of RCP channels, while the green vertical lines show the center of LCP channels. There is pretty good separation of RCP and LCP at lower and higher frequencies in the range, but around 12.4 GHz the polarization has switched and is mis-assigned.

The problem is that residual phase slope of Y with respect to X, caused by a difference in delay between the two channels. This can be seen in the upper panel of Figure 2, which shows the uncorrected phases of XY* and YX*. To correct the phases, the RCP phase was fit by a linear least-squares routine, and then the phases were offset by π/2 for both XY* and YX* according to:

where φ(v) is the phase fit by the linear function. The corrected phases are shown in the lower panel of Figure 2.

Figure 2: Upper panel) Raw phases on XY* (blue) and YX* (green). Lower panel) Phases corrected by removing a delay slope of order 3 ns, and shifting the phase such that the XY* phases on RCP channels is -π/2. This adjustment automatically makes XY* phases on LCP channels become +π/2.

When the corrected (primed) quantities are used in

the resulting channel separation is shown in Figure 3, which shows pretty much perfect channel separation, with the correct polarity.

Figure 3: Amplitudes of Stokes I, Stokes V, RCP and LCP for Ciel-2 satellite, after correcting for delay. The red vertical lines show the center of RCP channels, while the green vertical lines show the center of LCP channels. There is now PERFECT correspondence between RCP channels and the plot, and likewise for LCP channels.

This analysis can be performed on all antennas in a completely automated way, to produce the results of Fig. 3. The curves for each antenna are very similar, except that they show differences in amplitude (gain) as a function of frequency. When the curves are scaled to have the same Stokes I levels, the result is shown in Figure 4. Here the 13 antennas are shown in different colors. It is clear that matching Stokes I also causes an excellent match in the other polarizations. Remaining mismatch is likely due to not subtracting a no-signal background prior to scaling the Stokes I curves.

Figure 4: Same as Fig. 3, but for all 13 atnennas, with the Stokes I curve for each antenna scaled so that they agree. The beacon frequency in marked in the lower panel.

Also marked in Figure 4 is the location of the beacon signal for Ciel-2, which is at 12.209 GHz. The beacon signal appears in precisely the correct bin (frequency subchannel), indicating that we have precise tuning (certainly not unexpected, but nice to see).

Just to show that the approach is general, the results for a different communications satellite, Nimiq-5, is shown in Figure 5.

Figure 4: Same as Fig. 5, but for the Nimiq-5 commun-ications satellite.

Polarization Mixing Correction

Due to relative feed rotation between az-al mounted antennas and equatorial mounted antennas