Dominic F. Vecchia
Jack C.M. Wang
Statistical Engineering Division, CAML
Paul D. Hale
Optoelectronics Division, EEEL
The design of low cost lightwave communications systems requires accurate measurements of the response of optical to electrical converters in both magnitude and phase. The frequency range of interest is about 1 MHz to 50 GHz or more. To meet this need NIST is investigating methods to calibrate the frequency response of equivalent time sampling devices (both optical and electrical) with impulse or sinusoidal stimuli. Different methods will be used to cross check these calibrations.
In this work, a high-speed sampling oscilloscope automatically can produce
histograms comprising thousands of quasi-random-time samples from
input waveforms swept over many frequencies and power levels.
The model for N random-time samples from a signal generator under
test is given by
where are independently and uniformly distributed on [0, 1/f], and the ej's denote white Gaussian noise with variance . The amplitude a1 is the main parameter of interest, but the amplitude and phase, ak and , are to be estimated if the harmonic term is detected. We have obtained the first twelve moments of the sampling distribution of V for the most likely situations of a second or third harmonic term (k = 2 or 3). For deriving method-of-moments estimates of the parameters, we convert from moments to cumulants, since the latter quantities are simpler expressions than the former. For instance, the first five cumulants for the second harmonic model are: , , , , . These expressions show that, in the event that a2=0, an appropriate estimate of , when tranformed, will provide an estimate of a1, the primary parameter of interest. An estimate of the noise variance can then be obtained from an estimate of , though it need not exist for each sample.
More generally, we use unbiased estimates of the corresponding cumulants to estimate all of the parameters and, by propagation-of-error, their approximate standard errors. Estimates are obtained following an approximate test for the existence of the harmonic term (the procedure accounts for the possibility that, if the harmonic is in phase (), testing for nonzero values of is not sufficient to detect the harmonic). In the worst case ( ), the sixth cumulant is also needed in the derivations. Estimates of higher-order cumulants are required because the standard deviation of kj is a function of and lower order cumulants.
Figure 9: Histograms of 40,000 random time-samples from a 3 GHz waveform. The histogram in the top figure, by its symmetry, would suggest that a second harmonic, if present, is in phase with the fundamental. Asymmetry of the bottom histogram, which was obtained from a signal at a higher power level, is consistent with a second harmonic term. Measurable harmonic content is more likely as the power level is increased.
Date created: 7/20/2001