Abstract:
A test and measurement system includes a signal creation tool to generate a complex-valued stimulus signal, at least one waveform generator to receive the stimulus signal and produce at least one pair of baseband signals, a test instrument to capture the at least one pair of baseband signals and produce captured baseband signals, a pre-compensation coefficients estimation block to receive the captured baseband signals, characterize the captured baseband signals and to generate pre-compensation coefficients, and a pre-compensation block to apply the pre-compensation coefficients to the complex-valued stimulus signal when there is a device under test. 
     A method of characterizing a test system includes generating a multi-tone stimulus signal, producing at least one pair of baseband signals form the multi-tone stimulus signal, capturing the at least one pair of baseband signals with a test instrument, characterizing the at least one pair of baseband signals to generate pre-compensation coefficients, and applying the pre-compensation coefficients to signals applied to a device under test.

Description:
CROSS REFERENCE TO RELATED APPLICATIONS 
       [0001]    This application claims benefit of, and priority to, U.S. Provisional Application No. 62/257,612, “JOINT ESTIMATION OF COEFFICIENTS FOR SKEW, GAIN IMBALANCE AND CHANNEL RESPONSE FOR SIGNAL SOURCES,” filed Nov. 19, 2015, which is incorporated by reference herein in its entirety. 
     
    
     TECHNICAL FIELD 
       [0002]    This disclosure relates to baseband testing on devices under test, more particularly to pre-distortion coefficients used to filter the baseband signal. 
       BACKGROUND 
       [0003]    With the advancements in Digital Sampling and Semiconductor technology, Arbitrary Waveform Generators (AWGs) with high sampling rate Digital to Analog Converters (50GS/s and 25GS/s) have become available on the market. These high bandwidth AWGs can create wide bandwidth IQ (baseband) signals required for coherent optical signal communication receiver testing, where I represents the in-phase signal and Q represents the quadrature signal. 
         [0004]    AWGs that generate test signals for wide band applications such as optical communication are ideally required to have a flat magnitude and linear phase frequency response in the band of interest. In addition to flatness, skew and gain between the I and Q channels must be matched when testing wideband IQ signals. The quality of a digitally modulated signal is measured using Error Vector Magnitude (EVM). Any distortion in the signals due to skew or gain mismatches between the I and Q channels will increase the signal&#39;s EVM, indicating that the digitally modulated signal has poor quality. 
         [0005]    The test signal must undergo de-embedding or pre-compensating for any magnitude distortions, phase distortions, or skew and gain imbalance introduced by the test setup before being used to qualify a Device Under Test (DUT). Typically, pre-distortion coefficients are determined before testing a DUT, and may be used to de-embed or pre-compensate the test signal during testing. 
         [0006]    One method used to obtain the pre-distortion coefficients is a two-step process, wherein the frequency response of each of I and Q channels are measured independently. The frequency response is then inverted and two real-valued pre-distortion coefficients are obtained. These pre-distortion coefficients can be used by correction filters which are applied on the I and Q signals separately. 
         [0007]    Another method obtains complex pre-distortion coefficients by directly measuring the complex frequency response of the I and Q signals. 
         [0008]    The main disadvantage of the above methods lies in that they obtain frequency responses without considering the skew and magnitude imbalance between the channels. The skew and magnitude have a significant role in degradation of EVM of the signals. For example, measuring the complex frequency response without correcting for skew imbalances between the channels may result in incorrect phase measurements. 
         [0009]    To obtain the skew and gain imbalance, another prior art method may be used. For example, a known signal such as a step signal or a sinusoidal signal may be generated simultaneously from both the I and Q channels. The skew and gain imbalances between the channels are measured manually and corrected by either mathematically shifting and scaling one of the signals, or changing the clock phase between the AWGs and the amplitude control of the channels. 
         [0010]    One drawback of this prior art method results from the user having to manually obtain the skew and gain correction values. This often becomes a time consuming method of trial-and-error, and prone to user measurement errors. 
         [0011]    Automated methods exist to synchronize and adjust the skew between AWG channels, but they do not account for distortions caused by cables running to the Device Under Test (DUT). Currently, no methods are known for automatically adjusting gain imbalances between the channels. 
       SUMMARY 
       [0012]    An embodiment consists of a test and measurement system that has a signal creation tool to generate a complex-valued stimulus signal, at least one waveform generator to receive the stimulus signal and produce at least one pair of baseband signals, a test instrument to capture the at least one pair of baseband signals and produce captured baseband signals, a pre-compensation coefficients estimation block to receive the captured baseband signals, characterize the captured baseband signals and to generate pre-compensation coefficients, and a pre-compensation block to apply the pre-compensation coefficients to the complex-valued stimulus signal when there is a device under test. 
         [0013]    Another embodiment consists of a method of characterizing a test system that includes generating a multi-tone stimulus signal, producing at least one pair of baseband signals form the multi-tone stimulus signal, capturing the at least one pair of baseband signals with a test instrument, characterizing the at least one pair of baseband signals to generate pre-compensation coefficients, and applying the pre-compensation coefficients to signals applied to a device under test. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0014]      FIG. 1  shows a prior art baseband test setup. 
           [0015]      FIG. 2  shows an embodiment of a system to jointly calculate skew and gain imbalances and obtaining a frequency response of each channel. 
           [0016]      FIG. 3  shows a flowchart of an embodiment of jointly calculating skew and gain imbalances and obtaining a frequency response for baseband signals. 
           [0017]      FIG. 4  shows an embodiment of a pre-compensation coefficients estimation block. 
           [0018]      FIG. 5  shows an embodiment of a skew estimation block. 
       
    
    
     DETAILED DESCRIPTION OF THE EMBODIMENTS 
       [0019]      FIG. 1  shows an exemplary prior art system  10  for performing baseband testing on a device under test (“DUT”)  18  using two AWG channels. The AWGs receive a generated signal from the signal creation tool  12  and produce two signals one for the I channel or signal, and one for the Q channel or signal. Although  FIG. 1  shows separate AWGs  14  and  16  for each channel, a single AWG with two channels may also be used. AWGs that generate test signals for wide band applications such as optical communication are ideally required to have a flat magnitude and linear phase frequency response in the band of interest. In addition to flatness, skew and gain between the I and Q channels must be matched when testing wideband IQ signals. The quality of a digitally modulated signal is measured using Error Vector Magnitude (EVM). Any distortion in the signals due to skew or gain mismatches between the I and Q channels will increase the signal&#39;s EVM. 
         [0020]    The current embodiments disclose a method to jointly obtain the frequency response of each channel, and calculate the skew and the gain imbalances between the channels by using a single complex multi-tone signal as a stimulus. Though the rest of this document describes using one or more AWGs as the signal source, the disclosed method can be extended to any device capable of generating a baseband signal.  FIG. 2  shows an exemplary system  20  for jointly calculating skew and gain imbalances and obtaining a frequency response of each channel, and  FIG. 3  shows a flowchart of an embodiment of a method to determine these measures to perform pre-compensation. 
         [0021]    In  FIG. 2 , the system  20  has a signal creation tool  22 , and two AWGs, one 14 for the I channel or signal, and one 16 for the Q channel or signal, a test instrument  28 , such as a wideband oscilloscope, a pre-compensation coefficients estimation block  30  and a pre-compensation block  32 . 
         [0022]      FIG. 3  shows an embodiment of a process to jointly calculate skew and gain imbalances and obtaining a frequency response of each channel using a system such as that in  FIG. 2 . The signal creation tool  22  of  FIG. 2 , which may consist of one or more devices capable of generating baseband signals creates a complex-valued stimulus signal at  40 . A wideband oscilloscope or other test instrument such as  28  in  FIG. 2  then captures the I and Q signals generated by the AWGs  24  and  26  at  42 . If necessary, the test instrument adjusts the sampling rate for the captured signals at  44 . The pre-compensation coefficients estimation block  30  then estimates the skew imbalance between the channels at  46  and the gain imbalance between the channels at  48 . It also obtains the channel frequency response at  50 . These measures are then used to generate coefficients for the pre-compensation filter or block  32  of  FIG. 2 . The pre-compensation block  32  then operates on the incoming signal during operation of the system when it performs testing. When the DUT then undergoes testing, the signal used by the AWGs to produce the I and Q baseband signals will have been pre-compensated for the testing set up, allowing more accurate testing. 
         [0023]    Although  FIG. 2  depicts two separate AWGs, a single AWG with two separate I and Q channels could also be used. In some embodiments, more than two AWGs could also be used, for example, such as when multiple I and Q channels exist. Each of the signal creation tool  22 , the pre-compensation block  32 , and the pre-compensation coefficients estimation block  30  may be implemented by one or more of the AWGs  24  and  26 , by the scope  28 , or by another instrument or a device such as a general-purpose computer. 
         [0024]    As discussed at  40  in  FIG. 3 , a signal creation took generates a complex-valued multi-tone signal as a stimulus signal. The phases of the tones within the stimulus signal are chosen to have a good Peak to Average Power Ratio (PAPR). This is to make sure the full dynamic range of the AWG&#39;s Digital to Analog Converter (DAC) is used. In one embodiment, the signal creation tool uses either Newman phases or a combination of pseudo-random binary sequences (PRBS) used in standard wireless communication for channel estimation. An example stimulus signal may be described by the following equations: 
         [0000]    
       
         
           
             
               s 
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   m 
                   = 
                   1 
                 
                 Nt 
               
                
               
                 exp 
                  
                 
                   ( 
                   
                     
                       j 
                        
                       
                         ( 
                         
                           
                             2 
                             * 
                             π 
                             * 
                             n 
                             * 
                             
                               f 
                               m 
                             
                           
                           Fs 
                         
                         ) 
                       
                     
                     + 
                     
                       φ 
                       m 
                     
                   
                   ) 
                 
               
             
           
         
       
       
         
           
             
               
                 s 
                  
                 
                   ( 
                   n 
                   ) 
                 
               
               = 
               
                 
                   ∑ 
                   
                     m 
                     = 
                     1 
                   
                   Nt 
                 
                  
                 
                   ( 
                   
                     
                       cos 
                        
                       
                         ( 
                         
                           
                             cos 
                              
                             
                               ( 
                               
                                 2 
                                 * 
                                 π 
                                 * 
                                 n 
                                 * 
                                 
                                   f 
                                   m 
                                 
                               
                               ) 
                             
                           
                           Fs 
                         
                         ) 
                       
                     
                     + 
                     
                       j 
                        
                       
                           
                       
                        
                       
                         sin 
                          
                         
                           ( 
                           
                             
                               2 
                               * 
                               π 
                               * 
                               n 
                               * 
                               
                                 f 
                                 m 
                               
                             
                             Fs 
                           
                           ) 
                         
                       
                     
                   
                   ) 
                 
               
             
             ; 
           
         
       
     
         [0000]    where Nt is the number of tones, and Fs is the AWG&#39;s sampling rate. If two or more baseband signal generating devices are used, Fs represents the common sampling rate for all generating devices. 
         [0025]    In this case, s(n) is a complex-valued waveform, wherein the real part of s(t) is generated from an I channel, and the imaginary part of s(t) is generated from a Q channel. 
         [0026]    The ‘I’ &amp; ‘Q’ outputs of the AWG/AWGs are then captured on two or more channels of an instrument such as a wideband oscilloscope, where both the channels are triggered at the same time. To ensure an accurate measurement, the instrument must be well calibrated and any skew or gain imbalances between the instrument&#39;s channels must have been corrected prior to characterization. 
         [0027]      FIG. 4  shows an embodiment of characterizing the baseband signal. Before performing an FFT (Fast Fourier Transform) on the captured signals, it is important to make sure the frequency bins of the generated tones fall on the integer bins of the FFT output. If the duration of the signal required for frequency analysis (1*n/frequency resolution of the tones, where n is an integer) is not an integer number of samples, the bins will fall on a non-integer index. This is known as the picket fence effect. 
         [0028]    One common method to remedy the picket fence effect is to use time domain windowing such as a Hanning or flat top window, on the captured waveform  62 . This works well with the magnitude response, but distorts the phases. Using adjacent bins to calculate the power and correct for phase distortion does not help because of the intermodulation distortion (IMD) components in the frequency response. To solve this issue, the embodiments here resample the captured waveforms at  62  to obtain integer samples and an integer index for frequency bins of interest. This is what is meant by adjusting the sampling rate as necessary. 
         [0029]    For example, if the duration of the generated waveform from the AWG/AWGs is represented by D, and the sampling rate of the test instrument or scope is represented by F scope , then the number of samples required for analysis (Na) would be: 
         [0000]        Na=D*F   scope . 
         [0030]    When Na is an integer, the sampling rate used for analysis (F analyzer ) is equal to the sampling rate of the test instrument: 
         [0000]        F   analyzer   =F   scope . 
         [0031]    But, if Na is not an integer, then a new sampling rate can be calculated as, F analyzer =Round (Na)/D. The captured waveform can then be resampled based on F analyzer , such as by a factor of F analyzer /F scope . One should note that the re-sampler might add its own artifacts. But if the captured sampling rate (F analyzer ) is much larger the input bandwidth of the signal, the artifacts are minimized Next, the skew and the gain imbalances are estimated and corrected, before the channel response is estimated. 
         [0032]    Skew between the channels can be obtained at  64  by using the resampled real (I) and imaginary (Q) signals that were captured from the scope. The captured signals are independently cross correlated with the real and imaginary components of a reference signal. By observing the peak energy of the cross-correlated data, the lag for each channel can be determined with respect to the reference signal&#39;s I or Q component. The skew between the channels can be estimated by determining the difference in the lag between the I channel and the Q channel. 
         [0033]    To obtain a higher resolution of skew estimation, the captured signals may first be resampled again, using a higher sampling rate referred to as F skewanalysisrate . For example, to obtain a picosecond of resolution, the captured signal theoretically should be sampled at a rate of 1,000 GS/s. The maximum skew that can be estimated is also limited by the frequency resolution between the individual tones in the multi-tone signal, since the cross correlation peaks would repeat every cycle for the duration of the waveform. Increasing the frequency resolution, such as by lowering the minimum frequency, will increase the maximum skew that can be estimated. The reference signal for cross correlation is created using same phase values as the stimulus signal, but with the same sampling rate as the resampled captured signal, F skewanalysisrate . 
         [0034]    If x(n) represents the complex signal represented by the captured I and Q signals, which were captured at a sample rate of F analyzer  and x′(m) is the resampled captured signal, with a new sampling rate of F skewanalysisrate . x′ r  (m) and x′ i  (m), be the real and imaginary part of the resampled captured signal. Let r(m) be the reference signal created with sampling rate of F skewanalysisrate  and let r r  (m) and r i  (m), be the real and imaginary part of r(m). 
         [0035]    I channel index: 
         [0000]        C   xr,rr ( k )=Σ x′   r ( m )* r   n ( m−k )
 
         [0000]        C   peakreal =max(abs( C   xr,rr ( k )) 
         [0000]      Index real =Index of  C   peakreal  in  C   xr,rr ( k ) 
         [0036]    Q channel index: 
         [0000]        C   xi,ri ( k )=Σ x′   i ( m )* r   i ( m−k )
 
         [0000]        C   peakimag =max(abs( C   xi,ri ( k )) 
         [0000]      Index imag =Index of  C   peakimag  in  C   xi,ri ( k ). 
         [0000]    The estimated skew between the I and Q channels (Δt)=(Index imag −Index real ) F skewanalysisrate . Based on the estimated channel skew, the skew in the index (ΔIndex) can be estimated as (Δt)*F analyzer . Once the skew has been estimated, it can be used to adjust the imaginary component of x(n), to obtain skew-adjusted x″(n) at  66 : 
         [0000]      If  x ( n )= x   r ( n )+ jx   i ( n ), 
         [0000]        x ″( n )= x   r ( n )+ jx   i ( n −((ΔIndex))
 
         [0000]    By using this approach, fractional delay filters can be applied to correct for the skew. This results in a better resolution than the signal sampled at F analyzer . 
         [0037]    Gain imbalance refers to a difference in gain between the I and Q paths of the signals. If the ratio of gain between the I channel and the Q channel is not equal to 1, this will result in distorted constellation plot. The gain imbalance is estimated at 68 using the skew-corrected multi-tone signal. The gain imbalance is estimated by comparing the gain ratio between the I and Q components in the reference signal to the gain ratio of the I and Q signals in the skew corrected signal. 
         [0000]    
       
         
           
             
               IQ 
                
               
                   
               
                
               Gain 
                
               
                   
               
                
               ratio 
                
               
                   
               
                
               in 
                
               
                   
               
                
               the 
                
               
                   
               
                
               reference 
                
               
                   
               
                
               signal 
             
             = 
             
               
                 E 
                  
                 
                   ( 
                   
                     real 
                      
                     
                       ( 
                       
                         normalized 
                          
                         
                             
                         
                          
                         reference 
                          
                         
                             
                         
                          
                         signal 
                       
                       ) 
                     
                   
                   ) 
                 
               
               
                 E 
                  
                 
                   ( 
                   
                     imaginary 
                      
                     
                       ( 
                       
                         normalized 
                          
                         
                             
                         
                          
                         reference 
                          
                         
                           
                               
                           
                            
                           
                               
                           
                         
                          
                         signal 
                       
                       ) 
                     
                   
                   ) 
                 
               
             
           
         
       
       
         
           
             
               IQ 
                
               
                   
               
                
               Gain 
                
               
                   
               
                
               ratio 
                
               
                   
               
                
               in 
                
               
                   
               
                
               the 
                
               
                   
               
                
               skew 
                
               
                   
               
                
               corrected 
                
               
                   
               
                
               signal 
             
             = 
             
               
                 E 
                  
                 
                   ( 
                   
                     real 
                      
                     
                       ( 
                       
                         normalized 
                          
                         
                             
                         
                          
                         skew 
                          
                         
                             
                         
                          
                         corrected 
                          
                         
                             
                         
                          
                         signal 
                       
                       ) 
                     
                   
                   ) 
                 
               
               
                 E 
                  
                 
                   ( 
                   
                     imaginary 
                      
                     
                       ( 
                       
                         normalized 
                          
                         
                             
                         
                          
                         skew 
                          
                         
                             
                         
                          
                         corrected 
                          
                         
                             
                         
                          
                         signal 
                       
                       ) 
                     
                   
                   ) 
                 
               
             
           
         
       
     
         [0000]    where E is the expected value. In this context, the expected value represents the average value over the length of the samples. 
         [0000]    
       
         
           
             
               Gain 
                
               
                   
               
                
               imbalance 
                
               
                   
               
                
               
                 ( 
                 
                   G 
                   imb 
                 
                 ) 
               
             
             = 
             
               
                 { 
                 
                   
                     { 
                     
                       
                         ( 
                         
                           IQ 
                            
                           
                               
                           
                            
                           Gain 
                            
                           
                               
                           
                            
                           ratio 
                            
                           
                               
                           
                            
                           in 
                            
                           
                               
                           
                            
                           skew 
                            
                           
                               
                           
                            
                           corrected 
                            
                           
                               
                           
                            
                           signal 
                         
                         ) 
                       
                       
                         ( 
                         
                           IQ 
                            
                           
                               
                           
                            
                           Gain 
                            
                           
                               
                           
                            
                           ratio 
                            
                           
                               
                           
                            
                           in 
                            
                           
                               
                           
                            
                           reference 
                            
                           
                               
                           
                            
                           signal 
                         
                         ) 
                       
                     
                     } 
                   
                   - 
                   1 
                 
                 } 
               
               . 
             
           
         
       
     
         [0038]    Gain imbalance can now be corrected at  70  by applying G imb  to the real and imaginary portions of x″(n) to obtain the gain and skew-corrected signal x′″(n), as shown in the following equation: 
         [0000]    
       
         
           
             
               
                 x 
                 ′′′ 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   
                     x 
                     r 
                     ″ 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
                 * 
                 
                   { 
                   
                     1 
                     + 
                     
                       { 
                       
                         
                           G 
                           imb 
                         
                         
                           ( 
                           
                             2 
                             + 
                             
                               G 
                               imb 
                             
                           
                           ) 
                         
                       
                       } 
                     
                   
                   } 
                 
               
               + 
               
                 j 
                 * 
                 
                   
                     x 
                     i 
                     ″ 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
                 * 
                 
                   
                     ( 
                     
                       1 
                       - 
                       
                         { 
                         
                           
                             G 
                             imb 
                           
                           
                             ( 
                             
                               2 
                               + 
                               
                                 G 
                                 imb 
                               
                             
                             ) 
                           
                         
                         } 
                       
                     
                     } 
                   
                   . 
                 
               
             
           
         
       
     
         [0039]    Channel frequency response is now estimated at 72 using x′″(n). A time domain window is applied on x′″(n) to smooth its frequency response. Frequency domain analysis is done by performing a Fast-Fourier Transform (FFT) on the time-windowed version of x′″(n). 
         [0040]    The relative magnitude between all of the tones within x′″(n) is obtained. Group delay distortion is obtained by comparing the phases of the tones within x′″(n) with the phases in the stimulus signal and removing the constant group delay or linear component in the phase values. Pre-compensation coefficients are obtained by taking the inverse response. 
         [0000]    
       
         
           
             
               X 
                
               
                 ( 
                 w 
                 ) 
               
             
             = 
             
               F 
                
               
                 ( 
                 
                   
                     x 
                     ′′′ 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
                 ) 
               
             
           
         
       
       
         
           
             
               H 
                
               
                 ( 
                 w 
                 ) 
               
             
             = 
             
               
                 1 
                 
                   X 
                    
                   
                     ( 
                     w 
                     ) 
                   
                 
               
               . 
             
           
         
       
     
         [0000]    The impulse response (h(m)) is obtained from H(w) using a frequency sampling method. 
         [0041]    A user defined baseband (IQ) waveform, y(n) can be corrected for skew by using the estimated skew Δt, obtained using the method described above at  74  in  FIG. 4 , and shown in more detail in  FIG. 5 . Skew adjustment can be performed by circularly shifting the imaginary (Q) component of y(n) (represented as y i (n)) by Δt. Fractional delay filters can be used to shift samples by less than a whole sample. As a result, skew-corrected signal y′(n) can be obtained, as represented in the following equation: 
         [0000]        y ′( n )= y   r ( n )+ jy   i ( n −((−Δ t ))
 
         [0042]    Gain imbalance G imb  is corrected in both the real and imaginary parts of the signal y′(n) by using following equation: 
         [0000]    
       
         
           
             
               
                 y 
                 ″ 
               
                
               
                 ( 
                 n 
                 ) 
               
             
             = 
             
               
                 
                   
                     y 
                     r 
                     ′ 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
                 * 
                 
                   { 
                   
                     1 
                     + 
                     
                       { 
                       
                         
                           - 
                           
                             G 
                             imb 
                           
                         
                         
                           ( 
                           
                             2 
                             + 
                             
                               ( 
                               
                                 - 
                                 
                                   G 
                                   imb 
                                 
                               
                               ) 
                             
                           
                           ) 
                         
                       
                       } 
                     
                   
                   } 
                 
               
               + 
               
                 j 
                 * 
                 
                   
                     y 
                     i 
                     ′ 
                   
                    
                   
                     ( 
                     n 
                     ) 
                   
                 
                 * 
                 
                   
                     ( 
                     
                       1 
                       - 
                       
                         { 
                         
                           
                             ( 
                             
                               - 
                               
                                 G 
                                 imb 
                               
                             
                             ) 
                           
                           
                             ( 
                             
                               2 
                               + 
                               
                                 ( 
                                 
                                   - 
                                   
                                     G 
                                     imb 
                                   
                                 
                                 ) 
                               
                             
                             ) 
                           
                         
                         } 
                       
                     
                     } 
                   
                   . 
                 
               
             
           
         
       
     
         [0043]    The channel frequency response is corrected by applying the inverse complex valued impulse response filter h(m) on y″(n), as shown in the following equation: 
         [0000]        y ′″( n )= y ″( n ) convolve with  h ( m ).
 
         [0044]    In this manner, the system used to test a device can be characterized and the signals generated can be pre-compensated using a more accurate and automatic process than previously available. 
         [0045]    Although specific embodiments of the invention have been illustrated and described for purposes if illustration, it will be understood that various modifications may be made without departing from the spirit and scope of the invention.