Calibration Equations

The power spectrum at the input of the correlator when the telescope is pointed at the source of interest, tex2html_wrap_inline531, and at the reference target, tex2html_wrap_inline533, can be written as
eqnarray8
where

tabular23

tex2html_wrap_inline551, tex2html_wrap_inline553, and tex2html_wrap_inline555 are assumed not to vary significantly with frequency over the bandpass of the spectrometer. Taking ratios and doing appropriate algebra,

displaymath36
Since the gain is the unknown that one aims to remove by this calibration, it is assumed not vary with power level. One then obtains the difference between the signal spectrum and the reference spectrum:
 equation49
The ratio of the mean power levels can be obtained from a broadband detector or power meter or similar device. In general the spectrum will have poor baselines because the system gain is invariably non-linear to some small degree and the non-linearity is a function of frequency.

The system temperatures for the reference spectrum and the signal spectrum, respectively, are
eqnarray62
Equation 1 gives the signal, both line and continuum, by which the source differs from the reference:
eqnarray73

Since spectral lines are typically weak, tex2html_wrap_inline557 and tex2html_wrap_inline559, so that
displaymath88
and therefore
displaymath94
This allows us to rewrite equation 1 as
 equation101
and suggests the definition of a normalized calibrated spectrum as
equation109
so that
 equation114

Noise Analysis

The expected r.m.s. noise in tex2html_wrap_inline577 is calculated from equation 2:
eqnarray135
Substituting in the radiometer equation for the normalized spectrum
displaymath154
we get
 equation146
which, assuming that the sig and ref integration times are the same, leads to the result
 equation159
Generally, the sig/ref duty cycle is 50%, and the total integration sig+ref integration time, tex2html_wrap_inline579 is used in the radiometer equation, in which case
 equation171
A tidy way to avoid having to worry about the constant in the radiometer equation, and also to handle cases where the two spectra don't have the same integration times, is to define the integration time for the calibrated spectrum so that the constant in the radiometer equation remains unity:
equation179
so that
 equation187
In the case of a 50% duty cycle leads to tex2html_wrap_inline585.

Combining Spectra

We have two spectra tex2html_wrap_inline591 and tex2html_wrap_inline593 which we want to combine optimally so that the combined spectrum tex2html_wrap_inline595 has the minimum noise. To do this we combine them with weights tex2html_wrap_inline597 and tex2html_wrap_inline599 so that
displaymath209
The proof below shows that this requires that
 equation201
To apply this result, we note from equation 9 above that
 equation209
Since any constant divides out, we can choose tex2html_wrap_inline601 and thus
equation217

Using equation 4 to replace tex2html_wrap_inline600 with tex2html_wrap_inline602 and equation 16 (see below) to replace the denominator, we obtain the equivalent normalized spectrum:
equation231


Proof: The noise in the combined spectrum is
 eqnarray242
We want to minimize this noise by an optimum choice of tex2html_wrap_inline605:
eqnarray256
which proves equation 10.

Noise Propagation

With equation 10, equation 15 reduces to:
eqnarray308
Using equation 9, and cancelling out tex2html_wrap_inline609, we can write this as
 equation342
where tex2html_wrap_inline614 and tex2html_wrap_inline613 are defined as the effective system temperature and integration time of the combined spectrum. If tex2html_wrap_inline615, then proper scaling of the combined spectrum requires that tex2html_wrap_inline620, and thus that
equation357
It then follows from equation 16 that
equation361



Tom Kuiper
Thu Nov 4 20:42:14 UTC 1999