High bandwidth real-time oscilloscope

Abstract

A method and apparatus for digitizing a data signal. An input analog data signal, is received and split into a plurality of split signals. At least one of the split signals is mixed with a predetermined periodic function with a predetermined frequency. The split signals are then digitized and combined mathematically to form a single output data stream that is a substantially correct representation of the original input signal.

Description

FIELD OF THE INVENTION

The present invention relates to a high bandwidth real-time digital sampling oscilloscope (DSO) incorporating mixing (or heterodyning) to increase the bandwidth of a typical oscilloscope design with limited bandwidth.

BACKGROUND OF THE INVENTION

A digital sampling oscilloscope (DSO) is the primary tool utilized by engineers to view signals in electronic circuitry. As signals get ever faster, it is very beneficial to have DSOs capable of digitizing, displaying and analyzing these faster signals. The capability of a DSO to digitize fast signals is determined by its bandwidth and sample rate. The sample rate is the number of samples points taken of a waveform in a given amount of time and is inversely proportional to the sample period—the time between samples.

If a sinusoidal frequency sweep is performed from DC up to higher frequencies, the bandwidth is defined as the frequency at which the signal displayed on the DSO screen is approximately 30% smaller than the input sine-wave.

Since one of the uses of the DSO is to design and analyze new electronic devices, high end DSOs must operate at speeds much higher than the present state of the art in electronics. These speeds are generally unachievable through brute-force methods, such as simply providing ever-faster sampling chips, and many methods are employed to overcome this situation. One of the most common methods is a method inherent in the design of the original oscilloscope—that of triggering repeatedly on a periodic event. If an event is frequently, periodically repeating, the waveform at the time of the event can be repeatedly displayed on the screen. Furthermore, data from multiple trigger events average together to provide a good view of the waveform. This technique is the underlying method of a conventional sampling scope. A sampling scope repeatedly triggers on an event and acquires only a few points of the waveform (sometimes only one point of the waveform) on each trigger event. After repeated triggers, the points are reassembled according to the sampling algorithm to form a very high “effective” sample rate version of the waveform. Relatively low sample rates are utilized for each trigger event, and very high bandwidth samples may be generated. Furthermore, the repeated trigger events enable averaging, which can be utilized to increase the signal-to-noise ratio (SNR) and therefore enable further bandwidth increases. However, such a sampling scope presupposes a repetitive input signal so that the representation of the waveform can be generated over many triggers.

A common problem in complex signal analysis is that a signal that is to be analyzed is often not repetitive. In fact, it is very often the case that a non-repetitive event is the cause of some failure in an electronic system. It is the function of the test equipment to help the user identify the cause of the failure. Therefore, a piece of test equipment that requires repetitive signals is of limited usefulness. For example, sometimes the trigger event happens only once, such as in the analysis of bomb blasts. Frequently, however, the trigger event happens repeatedly, but the signal around the trigger event is different. Situations like this require a DSO capable of high bandwidth and sample rate with only a single trigger event. A DSO with these characteristics is called a real-time scope, and acquisitions taken utilizing only a single trigger event are called single-shot acquisitions. The distinction between the sampling scope and the real-time scope is an important one because the tricks that can be utilized to digitize a repetitive waveform are not available to the real-time DSO designer. In general, a real-time DSO is more useful because it does not require the input signal to be repetitive. However, the primary limitation is that the bandwidth of the real-time scope is limited.

In real-time DSO design, the method in common use for overcoming sample rate limitations is the method of interleaving. This method utilizes multiple digitizing elements that sample the same waveform at different points in time such that the waveform resulting from combining the waveforms acquired on these multiple digitizers forms a high sample rate acquisition. Most high-end real-time DSOs have very high sample rates achieved through the use of interleaving and most are capable of “oversampling” an input waveform.

Oversampling is defined as sampling a waveform at a rate whereby virtually no amount of signal content is present at a frequency above one half the sample rate. For example, a DSO with a bandwidth of 6 GHz that does not allow any signal in with a frequency at or above 10 GHz would be sufficiently sampling the waveform at a sample rate of 20 GHz. Any sampling of the waveform above this sample rate would result in an oversampled waveform. Oversampling is not inherently bad, just unnecessary because much more elegant methods can be utilized to produce the highly sampled waveform. The criteria for sufficiency of sample rate, outlined by Nyquist Criterion, states that if a waveform is sampled at a sufficient rate, than the exact analog waveform can be reconstructed. In other words, once the waveform is sufficiently sampled, the waveform can be reconstructed as if physically digitized at any sample rate.

Generally, in real-time DSOs, the interleaving is controlled through a method called channel combination. Combining channels means that the digitizing resources of multiple channels are utilized together to digitize a single waveform. Most often, channel combination is utilized to interleave multiple digitizers for the purpose of increasing the sample rate, but as the acquisition memory is generally connected to individual digitizers, this method is also utilized sometimes to increase the length of the acquisition.

While techniques are generally available for designing high sample rate systems, bandwidth is another issue. Bandwidth is typically dealt with through direct application of very high-speed electronics. In situations where electronics are simply not fast enough, usually the attainment of high bandwidth is achieved by making tradeoffs that the customer simply must accept. For example, higher bandwidth is achievable by removing protection circuitry at the front-end of the scope, thereby making it more susceptible to damage from static discharge or signal overdrive. Also, limitations are placed on the user regarding the range of input signals (to allow for attenuators and active gain components to be eliminated). The tradeoffs foisted on the scope user are often unpalatable, but endured reluctantly by the user with high-bandwidth requirements.

Despite this situation, the fact remains that attempts made to reach high bandwidths are often done at the expense of the overall usability of the scope. In other words, a high bandwidth scope can often not be utilized in a general-purpose manner. Finally, the fact remains that even with every possible trade-off, the bandwidth needs of the real-time scope user are sometimes so high as to be unachievable with the current state of the art.

SUMMARY OF THE INVENTION

This invention pertains generally to systems that digitize waveforms; and more specifically systems that convert an analog input signal to a digital output signal whereby the digital signal consists of an array of numbers that represent the amplitude of the analog waveform at known times. This invention also pertains to systems with limited bandwidth where there is a need for higher bandwidth. This invention therefore addresses systems incapable of accurately digitizing very rapidly changing signals.

The most specific application of this invention is to the high-end real-time DSO where extremely high demands are placed on the speed (and bandwidth) of signals digitized in a single-shot acquisition.

It is an object of the invention to demonstrate a method and apparatus whereby the bandwidth of a digitizing system can be increased.

As mentioned previously, channels are often combined within a DSO for the purpose of increasing the sample rate and acquisition memory length. This invention puts forth a method of utilizing channel combination for the purpose of increasing bandwidth.

Still other objects and advantages of the invention will in part be obvious and will in part be apparent from the specification and the drawings.

The invention accordingly comprises the several steps and the relation of one or more of such steps with respect to each of the others, and the apparatus embodying features of construction, combination(s) of elements and arrangement of parts that are adapted to effect such steps, all as exemplified in the following detailed disclosure, and the scope of the invention will be indicated in the claims.

BRIEF DESCRIPTION OF THE DRAWINGS

For a more complete understanding of the invention, reference is made to the following description and accompanying drawings, in which:

FIG. 1 is a block diagram depicting a digitizing system constructed in accordance with the invention;

FIG. 2 is a block diagram showing one possible method of the extension of this technique to 4 channels using a mixing frequency that is at the low side of the frequency band of interest (low side conversion);

FIG. 3 is a block diagram showing another possible method of the extension of this technique to 4 channels using a mixing frequency that is at the high side of the frequency band of interest (high side conversion); and

FIGS. 4–57 are waveform diagrams of signals input to and output from the processing circuits described herein, and are useful in explaining the advantages obtained by the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

FIG. 1 is a block diagram showing a high bandwidth digital oscilloscope architecture according to the present invention. FIG. 1 shows two channels of a DSO combined to digitize waveforms in a manner that effectively doubles the system bandwidth. It should be understood that the bandwidth can be tripled, quadrupled etc. by utilizing three, four or more channels in combination.

An input signal is provided at the input 100. Viewed from the frequency-domain perspective, the input signal might have a frequency content shown as 103. In a standard configuration input signal applied to input 100 directly enters a first channel, CH1 at 104. This signal passes through an analog front end 105, and on to an ADC 106 which digitizes the waveform. The channel has a finite bandwidth, as shown by 107 which results in a digitized waveform of finite bandwidth 108. For the purpose of future explanation, the cutoff frequency at which the bandwidth is limited is designated as F. In a conventional digitizer the CH1 channel (104) and a CH2 channel (109) are combined utilizing preferably an additional output of the CH1 channel (104) front-end (105) connected to an additional input of a CH2 channel ADC (111-connection not shown) for the purpose of doubling the sample-rate and memory length of the acquisition. A less preferable connection combines the channels through a 50 Ohm power splitter at the two channel inputs 104 and 109. Neither of these methods used in the current state of the art can increase the bandwidth. If there is any effect at all, it is to decrease the bandwidth somewhat. This is a side effect and not generally desired. The effect can be minimized with careful design.

The present invention involves the addition of additional analog circuitry 102 between the input 100 and the two channels 104 and 109 and downstream processing of the digital data to account for this additional hardware. In accordance with the present invention, the signal at the input 100 with example frequency content 103 enters a 50 Ohm power splitter 112. The splitter 112 provides the 50 Ohm termination to the input signal and provides at its two outputs, ideally the same signal attenuated. One output of splitter 112 directly connects to CH1 channel (104). The signal proceeds through front-end 105 and is digitized in the normal manner by ADC 106. Since the combination of front-end 105 and ADC 106 is bandwidth limited, as shown in 107, this results in a bandwidth limited acquisition with frequency content shown by 108. CH1 channel (104) is designated as containing the low frequency (LF) portion of the signal. The path through the other output of splitter will be described below and is designated as the high frequency (HF) path.

The addition of the splitter is the only additional component directly in the signal path of CH1 channel (104). In other words, the low frequency signal path with the splitter removed looks identical to the signal path of the DSO not utilizing this invention. Very high quality splitters with very high bandwidth are readily commercially available and as such, do not serve to degrade the signal path, except for possibly decreasing the signal strength.

The HF signal path will now be explained in greater detail. It is well known that frequencies can be shifted through the use of a process called mixing or heterodyning. This process is called “frequency translation”. Mixing is achieved through the time-domain multiplication of a signal with another sinusoidal signal. It is well known that if a sinusoid with a frequency f₀ is mixed with another sinusoid with a frequency f₁, the result is two sinusoids at sum and difference frequencies (i.e. sinusoids at frequencies f₀+f₁ and f₀−f₁) with each sinusoid being half the amplitude of the product of the amplitudes of f₀ and f₁:

$\left[A_0·\cos \left(2·\pi ·f_0\right)\right]·\left[A_1·\cos \left(2·\pi ·f_1\right)\right]\to \frac{A_0·A_1}{2}·\left[\cos \left(2·\pi ·\left(f_0+f_1\right)\right)+\cos \left(2·\pi ·\left(f_0-f_1\right)\right)\right]$ In order to utilize these principles, the signal from the second output 101 of splitter 112 enters a high pass filter 113. High-pass filter 113 is designed to reject to the greatest extent possible all frequencies below frequency F as shown in 114. The result of high pass filtering the input signal with frequency content as shown in 103 is shown in 115. The output of 113 is mixed with a sinusoid at frequency F utilizing mixer 116. The result of mixing the signal with frequency content shown by 115 with the mixing frequency F shown in 117 is the frequency content shown in 118. 118 shows that two images of the content shown in 115 are produced at the sum and difference frequencies, as noted above. In cases where the cutoff of the high-pass filter is inadequate, the mixing frequency could be chosen slightly higher such that a dead-band is utilized to prevent the low frequency edges of the high-pass filter output from folding back into the pass-band. Any increase in mixing frequency, while providing margin, will also serve to degrade the maximum bandwidth achievable.

The output of mixer 116 connects to CH2 channel (109)—the high frequency channel. The signal passes through a front-end 110 and is digitized by an ADC 111. Since like CH1 channel (104), the combination of front-end 110 and ADC 111 is bandwidth limited, as shown in 119 this results in a bandwidth limited acquisition with frequency content shown by 120.

The acquisition through ADCs 106 and 111 occur simultaneously in parallel, so channels CH1 and CH2 are simultaneously acquired.

To summarize, the signals seen at the inputs to each of the channels are as follows. The LF CH1 channel (104) sees the input waveform directly. The HF CH2 channel (109), sees only the frequency content in the frequency band extending upward from F. Through the use of heterodyning, it sees the frequency content of the input signal at F+ΔF at the frequency location ΔF. Thus, the input signal has been mixed down. Although there were two images (seen in 118), the second image (and part of the first image) were rejected due to the finite channel bandwidth 119. Said differently, the LF CH1 channel acquires the low frequency content of the input signal from 0–F, while the HF CH2 channel acquires the frequency content of the input signal from F−2·F. This signal is mixed down from frequency band F→2*F to the range of 0→F so it “fits” into the bandwidth of the front end. It can be seen that twice the frequency content of the signal has been made to “fit” into the bandwidth of the scope.

Both the LF and HF signals are digitized by the scope in the normal manner. It is assumed that both channels are sufficiently sampling with respect to the channel bandwidth. In other words, each channel is sampling at a sample rate (Fs) such that virtually no frequency content can get through the channels above Fs/2. This might be accomplished through the use of internally interleaving channel digitizers or through the combination of other channels, as mentioned previously. Since a sufficiently sampled channel allows for the complete reconstruction of the signal at any sample rate, each channel is upsampled to a sample rate that at least sufficiently samples the resulting acquisition utilizing upsamplers 121 and 122. In the case where a channel is just sufficiently sampling based on the frequency response of the channel, the acquisition on each channel is upsampled by a factor of 2 when two channels are combined, since the system bandwidth will be effectively doubled in the end. This upsampling is performed utilizing a method such as SinX/X interpolation to interpolate every other sample point. The method and validity of this method of interpolation is well known to those skilled in the art.

The data from the HF CH2 channel is mixed at a mixer 123 digitally (i.e. numerically using a software program and floating point arithmetic) with a sinusoid with the same frequency F (124) as utilized by analog mixer 116 in the analog HF signal path. The result of mixer 123 is two images of the HF signal shown by 125. Each image appears at frequencies from 0–F and from F−2·F. The first, low frequency image is mirrored about F and is unusable. The second image is a replica of the actual high frequency content of the input signal. The output of digital mixer 123 passes through an image reject filter 126, which has a frequency response shown by 27. The result is the frequency content shown in 128.

The result at this point is two digital waveforms, one representing the low frequency portion of the input signal 131 and the other representing the high frequency portion of the input signal 128.

Since both signals passed through an imperfect channel, they are equalized separately to compensate for non-ideal magnitude and phase characteristics of the front-end and digitizing systems. The equalizer for the CH1 channel (129) is shown with its response 130 being an ideal low-pass filter. This results in no change between the frequency content shown in 108 and the equalized content shown in 131 with the understanding that this would not necessarily be the case if there were imperfections in the signal 131. The equalizer for the CH2 channel is shown integrated with the image reject filter 126. Both equalizers also have the difficult job of preparing the signals to accommodate the cross-over from LF to HF.

Finally, the waveform resulting from the LF channel equalizer 129 is added to the waveform resulting from the HF channel equalizer 126 by the summer 132. The result of this addition is shown by adding the LF frequency content in 131 to the HF content in 128, shown graphically in 133. This forms a high-bandwidth, high-sample rate acquisition at the output 135 with the frequency content as shown in 134.

A detailed numerical example stepping through this process is provided below.

To summarize the effect, the input signal with frequency content 103 normally would be digitized by one channel to form an acquisition with frequency content shown in 108. Instead, as a result of this invention, two channels were utilized and the resulting acquisition as twice the bandwidth as demonstrated by the signal frequency content shown in 134. Note that the two other benefits of interleaving—that of doubling the sample rate and memory length are still achieved by this invention.

Multiple channels could be combined in similar arrangements whereby the system bandwidth is increased by a factor equal to the number of combined channels. An example showing how this technique may be used to extend the bandwidth four times using four channels is shown in FIG. 2. FIG. 2 depicts a low side conversion. In this case, each frequency band in 201 (A, B, C & D) is translated down to be digitized by channels 1, 2, 3 and 4 respectively. The filter blocks shown (202) are optional. The goal is to sufficiently isolate the desired frequency band. This can be done, in this example, using a highpass filter, a bandpass filter or no filter at all. In the case of no filter, there will be “images” of the adjacent frequency band that will be digitized by the channel. These can be removed using DSP techniques in the recombination DSP block (203). This recombination technique is shown for a high side conversion (using three channels) below.

Another example of how this technique can be applied is shown in FIG. 3. FIG. 3 depicts a high side conversion. Careful inspection of FIG. 3 will show that the major difference between FIG. 2 and FIG. 3 is the frequency used to “translate” the frequency band of interest into the frequency band of the acquisition channel. FIG. 2 uses a frequency on the low side of the band of interest (F₁ to translate frequency band “B” to the frequency band of the acquisition system) and FIG. 3 uses a frequency on the high side of the band of interest (F₂ to translate frequency band “B” to the frequency band of the acquisition system). Note that the translated frequency band is “reversed” if high side conversion is used. That is, the highest frequency in frequency band “B” (F₂) becomes the lowest frequency in the translated band. This will be corrected in the reconstruction by using F₂ as the up-conversion frequency in the DSP reconstruction of the signal. This recombination technique is shown for a high side conversion (using two channels) below.

In either case shown in FIGS. 2 and 3, it is required that the phase of the translation frequency is known in order to reconstruct the original signal. This can be accomplished, by example, by summing a pilot tone into the signal channel, or locking the mixer phase to the sample clock.

EXAMPLES

The first example shows how a step can be digitized using two 5 GHz bands and low side downconversion.

rt := .045                       Risetime of edge specified (ns).
f                      bw                                        :=                                          0.344                      rt fbw = 7.644 Bandwidth of criticallydamped second order system.
ω0 := 1.5542 · π · fbw           Calculate the center frequencyfor the system.
ω                        ⁢                                                                                                  ⁢                        0                                                                    2                        ·                        π                                                              =                    11.879 Center frequency (GHz).
TD := 5                          Time delay for step edge (ns).
H                      ⁡                                              (                        s                        )                                                              =                                                                                            ω                          ⁢                                                                                                          ⁢                                                      0                            2                                                                                                                                (                                                                                          s                                2                                                            +                                                                                                                                    ω                                    ⁢                                                                                                                                                  ⁢                                    0                                                                    Q                                                                ·                                s                                                            +                                                              ω                                ⁢                                                                                                                                  ⁢                                                                  0                                  2                                                                                                                      )                                                    ·                          s                                                                    ·                                              e                                                                              -                            s                                                    ·                          TD Laplace transform of the stepspecified.

The inverse Laplace transform provides the time-domain step waveform.f(t):=if[t<TD,0,[−1−ω0·(t−TD)]·e^{[−ω0·(t−TD)]}+1]To simulate the behavior of the analog components, it is modeled digitally with an extremely high sample rate.

FShi := 1000                     Sample rate for simulating analog system (GHz).
KH := 10000                      kh := 0 . . . KH − 1
th                      kh                                        :=                                          kh                                              FS                        hi Time of each point (ns).

Utilize a raised cosine window to minimize effects of the FFT.

wh                      kh                                        :=                                                                  1                        2                                            -                                                                        1                          2                                                ·                                                  cos                          ⁡                                                      (                                                          2                              ·                              π                              ·                                                              kh                                                                  KH                                  -                                  1                                                                                                                      )
xhkh := f(thkh) · whkh           Calculate the windowed step.
t                      10                                        :=                                                                  .53181160838961202015                                                  ω                          ⁢                                                                                                          ⁢                          0                                                                    +                      TD t10 = 5.007
t                      90                                        :=                                                                  3.8897201698674290579                                                  ω                          ⁢                                                                                                          ⁢                          0                                                                    +                      TD t90 = 5.052
t90 − t10 = 0.045                Verify that risetime is correct.

FIG. 4 shows a picture of the simulated 45 ps step.

	Xh := CFFT(xh)	Calculate the FFT.
	$\mathrm{NH}:=\frac{\mathrm{KH}}{2}$	nh := 0 . . . NH	${\mathrm{fh}}_{\mathrm{nh}}:=\frac{\mathrm{nh}}{\mathrm{NH}}·\frac{{\mathrm{FS}}_{\mathrm{hi}}}{2}$

FIG. 5 shows the frequency content of the simulated step.

As we know, the scope does not have the bandwidth to digitize this signal. Therefore, we apply the method of this invention. First, we will utilize a system bandwidth of 5 GHz. then, we develop bandpass filters that select 5 GHz bands of the signal. Note that because the system is bandlimited, it is not actually necessary to utilize bandpass filters—only high pass filters need be utilized, but bandpass filters are used to simplify the discussion. Furthermore, the first band does not even need a filter—the scopes limited bandwidth will do this for us. (inside the scope, a digital low pass filter would be utilized to provide the hard bandwidth limiting).

• BW:=5 System bandwidth utilized for each (GHz).Make low pass and bandpass filters for each band.nn:=1 . . . NH−1Mfl_nh:=if(fh_nh≦BW,1,0)Mfh_nh:=if(BW<fh_nh≦2·BW,1,0)Mfl_NH+nn:=Mfl_NH−nn Mfh_NH+nn:=Mfh_NH−nn

FIG. 6 shows the response of the low frequency band filter.

FIG. 7 shows the response of the high frequency band filter.

Apply these filters to the input waveform.Xfl:={right arrow over ((Xh·Mfl))}Xfh:={right arrow over ((Xh·Mfh))}Xf:={right arrow over ((Xh·M))}

FIG. 8 shows the frequency domain result of applying the low frequency band filter to the simulated step.

FIG. 9 shows the frequency domain result of applying the high frequency band filter to the simulated step.

Calculate the inverse FFT of these filtered bands.xfl:=ICFFT(Xfl)xfh:=ICFFT(Xfh)xf:=ICFFT(Xf)

FIG. 10 shows the input waveform.

FIG. 11 shows the input waveform bandlimited to 10 GHz.

FIG. 12 shows the time domain result of applying the low frequency band filter to the simulated step.

FIG. 13 shows the time domain result of applying the high frequency band filter to the simulated step.

It is useful to add these two signals together and compare them to the input waveform. FIG. 14 shows this. You will note the sum is not identical to the input because the system has limited the bandwidth at 10 GHz. The 10 GHz bandwidth limited signal is the best that we will be able to provide.

It is also useful to compare the low frequency and actual input waveforms directly. FIG. 15 shows this. The point of FIG. 15 is to demonstrate the problem that this invention is designed to solve. The limited bandwidth slows the edge of the step. This simulates the analog waveform that gets sampled by a digitizer with a front-end bandwidth of 5 GHz. Our goal is to digitize the actual waveform with a much higher bandwidth.

First, the high frequency band is applied to the mixer.F_mixer0:=1·BW Φ_mixer0:=rnd(2·π)The frequency of the high frequency mixer is at the cutoff frequency of the first band. Apply the mixer.xfhm_kh:=xfh_kh·2·cos(2·π·F_mixer0th_kh+Φ_mixer0)Look at the frequency content.Xfhm:=CFFT(xfhm)Low pass filter the mixer outputs.Xfhml:={right arrow over ((XfhmMfl))}

Note again that the typical manner of low pass filtering the mixer outputs would be to use the scope front-end. This filtering is being shown here as actual low pass filters applied. FIG. 16 shows the effect at the high frequency band mixer and image reject filter output in the frequency domain.

Take the inverse FFT to generate the analog mixer output signals—the analog signals input to the channel digitizers.xfhml:=ICFFT(Xfhml)

FIGS. 17 and 18 shows the low frequency and high frequency band signal due to the processing of the simulated step input signal.

It is interesting to see what the sum of these two waveforms is—their sum does not produce anything useful. This is shown in FIG. 19. This is shown to exemplify that further processing of the digitized waveforms is required to achieve the objective.

At this point, the waveforms are digitized. The waveforms must be sampled at a rate sufficient to satisfy Nyquist Criterion. For this example, this means that they must be sampled at at least 2 times BW, or 10 GS/s. After the waveforms have been digitized, they are immediately upsampled using SinX/x interpolation. This is possible because all digitized waveforms are bandlimited. It is useful to upsample the waveforms to a sample rate capable of meeting Nyquist Criterion for the system bandwidth—I have chosen 40 GS/s. The upsampling is trivial and for the purpose of this example, I simply use a 40 GS/s digitizer with the understanding that the exact same waveform would result from sampling the waveform at 10 GS/s and upsampling by a factor of 4.

FS := 40                         Upsampled digitizer sample rate.
D                    :=                                                                  FS                        hi                                            FS D = 25 Upsampling factor for analog waveform model.
K                    :=                                          KH                      D k := 0 . . . K − 1

Sample the waveforms.

$\begin{array}{c}{t}_k:=\frac{k}{\mathrm{FS}}\\ x_{l_k}:={\mathrm{xfl}}_{k·D}\\ x_{h_k}:={\mathrm{xfhml}}_{k·D}\\ x_k:={\mathrm{xh}}_{k·D}\\ w_k:={\mathrm{wh}}_{k·D}\end{array}$

Generally, at this point, we would apply the sharp cutoff filter. If a sharp cutoff analog filter was not used, we'd have to satisfy Nyquist Criterion such that any extra frequency content would not fold back into the 5 GHz band. I've already applied a sharp cutoff filter to the analog signal, so this is not necessary.

Also, at this point, some magnitude and phase compensation would probably be necessary to account for non-deal channel frequency response characteristics. This example shows the signal digitized with ideal digitizers with ideal frequency response characteristics. Next, the high and very high frequency waveforms are mixed up to there appropriate frequency location and digitally bandpass filtered.

Note that these digital mixers know the phase of the analog mixers—some mechanism must be provided for determining this—either through a pilot tone or locking of the mixer phase to the sample clock.

Apply digital mixers.x_hmk:=x_hk·(2·cos(2·π·F_mixer0t_k+Φ_mixer0))Bandpass filter the mixer outputs.

$\begin{array}{c}N:=\frac{K}{2}\\ n:=0\dots N\\ f_n:=\frac{n}{N}·\frac{\mathrm{FS}}{2}\\ \mathrm{Xhm}:=\mathrm{CFFT}\left(x_{\mathrm{hm}}\right)\\ \mathrm{Xlm}:=\mathrm{CFFT}\left(x_l\right)\\ {\mathrm{Xfhm}}_n:=\mathrm{if}\left(f_n<1·\mathrm{BW},0,{\mathrm{Xhm}}_n\right)\\ {\mathrm{Xfhm}}_n:=\mathrm{if}\left(f_n>2·\mathrm{BW},0,{\mathrm{Xfhm}}_n\right)\\ \mathrm{nn}:=1\dots N-1\\ {\mathrm{Xfhm}}_{N+\mathrm{nn}}:=\stackrel{\_}{{\mathrm{Xfhm}}_{N-\mathrm{nn}}}\\ X_h:=\mathrm{CFFT}\left(x_h\right)\\ X_l:=\mathrm{CFFT}\left(x_l\right)\end{array}$

FIG. 20 shows the low frequency band waveform frequency content.

FIG. 21 shows the high frequency band waveform frequency content before and after digital mixing and filtering.

FIG. 22 shows the results of summing the output waveforms. We have acquired the waveform with a 10 GHz bandwidth utilizing two 5 GHz bandwidth channels.

Now let's see how the time domain waveforms compare.xfhm:=Re(ICFFT(Xfhm))

FIG. 23 shows a comparison of the analog low frequency portion of the input waveform to the digitized and processed low frequency waveform result. FIG. 23 shows that these are identical.

FIG. 24 shows a comparison of the analog high frequency portion of the input waveform to the mixed, digitized and digitally remixed processed high frequency waveform result. FIG. 24 shows that these are identical.

FIGS. 25 and 26 show how successive addition of the processed band outputs causes the resulting waveform to more and more closely approximate the input waveform.

FIG. 26 shows that the result is identical to the input waveform band limited to 10 GHz and shows that the 10 GHz bandwidth limited step is fully recreated.

The second example shows how a step can be digitized using three 5 GHz bands and high side downconversion.

rt = .035                        Risetime of edge specified (ns).
f                      bw                                        =                                          0.344                      rt fbw = 9.829
ω0 = 1.5542 · π · fbw            Calculate the center frequency forthe system.
ω                        ⁢                                                                                                  ⁢                        0                                                                    2                        ·                        π                                                              =                    15.274 Center frequency (GHz).
TD = 5                           Time delay for step edge (ns).
H                      ⁡                                              (                        s                        )                                                              =                                                                                            ω                          ⁢                                                                                                          ⁢                                                      0                            2                                                                                                                                (                                                                                          s                                2                                                            +                                                                                                                                    ω                                    ⁢                                                                                                                                                  ⁢                                    0                                                                    Q                                                                ·                                s                                                            +                                                              ω                                ⁢                                                                                                                                  ⁢                                                                  0                                  2                                                                                                                      )                                                    ·                          s                                                                    ·                                              e                                                                              -                            s                                                    ·                          TD Laplace transform of the stepspecified.

The inverse Laplace transform provides the time-domain step waveform.f(t)=if[t<TD,0,[−1−ω0·(t−TD)]·e^{[−ω0·(t−TD)]}+1]To simulate the behavior of the analog components, it is modeled digitally with an extremely high sample rate.

FShi = 1000                      Sample rate for simulating analog system (GHz).
KH = 10000                       kh = 0 . . . KH − 1
th                      kh                                        =                                          kh                                              FS                        hi Time of each point (ns).

Utilize a raised cosine window to minimize effects of the FFT.

wh                      kh                                        =                                                                  1                        2                                            -                                                                        1                          2                                                ·                                                  cos                          ⁡                                                      (                                                          2                              ·                              π                              ·                                                              kh                                                                  KH                                  -                                  1                                                                                                                      )
xhkh = f(thkh) · whkh            Calculate the windowed step.
t                      10                                        =                                                                  .53181160838961202015                                                  ω                          ⁢                                                                                                          ⁢                          0                                                                    +                      TD t10 = 5.006
t                      90                                        =                                                                  3.8897201698674290579                                                  ω                          ⁢                                                                                                          ⁢                          0                                                                    +                      TD t90 = 5.041
t90 − t10 = 0.035                Verify that risetime is correct.

FIG. 27 shows a picture of the simulated 35 ps step.

	Xh = CFFT(xh)	Calculate the FFT.
	$\mathrm{NH}=\frac{\mathrm{KH}}{2}$	nh = 0 . . . NH	${\mathrm{fh}}_{\mathrm{nh}}=\frac{\mathrm{nh}}{\mathrm{NH}}·\frac{{\mathrm{FS}}_{\mathrm{hi}}}{2}$

FIG. 28 shows the frequency content of the simulated step.

As we know, the scope does not have the bandwidth to digitize this signal. Therefore, we apply the method of this invention. First, we will utilize a system bandwidth of 5 GHz. then, we develop bandpass filters that select 5 GHz bands of the signal. Note that because the system is bandlimited, it is not actually necessary to utilize bandpass filters—only high pass filters need be utilized, but bandpass filters are used to simplify the discussion. Furthermore, the first band does not even need a filter—the scopes limited bandwidth will do this for us. (inside the scope, a digital low pass filter would be utilized to provide the hard bandwidth limiting).

• BW:=5 System bandwidth utilized for each (GHz).Make low pass and bandpass filters for each band.nn:=1 . . . NH−1Mfl_nh=if(fh_nh≦BW,1,0)Mfl_NH+nn=Mfl_NH−nn

FIG. 29 shows the response of the low frequency band filter.Mfh_nh=if(BW<fh_nh≦2·BW,1,0)Mfh_NH+nn=Mfh_NH−nn

FIG. 30 shows the response of the high frequency band filter.Mfhh_NH+nn=Mfhh_NH−nn Mfhh_nh=if(2·BW<fh_nh≦3·BW,1,0)

FIG. 31 shows the response of the high frequency band filter. Apply these filters to the input waveform.Xfl={right arrow over ((Xh·Mfl))}Xfh={right arrow over ((Xh·Mfh))}Xfhh={right arrow over ((Xh·Mfhh))}Xf={right arrow over ((Xh·M))}

FIG. 32 shows the frequency domain result of applying the low frequency band filter to the simulated step.

FIG. 33 shows the frequency domain result of applying the high frequency band filter to the simulated step

FIG. 34 shows the frequency domain result of applying the very high frequency band filter to the simulated step.

Calculate the inverse FFT of these filtered bands.

xfl=ICFFT(Xfl)xfh=ICFFT(Xfh)xfhh=ICFFT(Xfhh)xf=ICFFT(Xf)

FIG. 35 shows the band limited input waveform.

FIG. 36 shows the time domain result of applying the low frequency band filter to the simulated step.

FIG. 37 shows the time domain result of applying the high frequency band filter to the simulated step.

FIG. 38 shows the time domain result of applying the very high frequency band filter to the simulated step.

It is useful to add these three signals together and compare them to the input waveform. FIG. 39 shows this. You will note the sum is not identical to the input because the system has limited the bandwidth at 15 GHz. The 15 GHz bandwidth limited signal is the best that we will be able to provide.

It is also useful to compare the low frequency and actual input waveforms directly. FIG. 40 shows this. The point of FIG. 40 is to demonstrate the problem that this invention is designed to solve. The limited bandwidth slows the edge of the step. This simulates the analog waveform that gets sampled by a digitizer with a front-end bandwidth of 5 GHz. Our goal is to digitize the actual waveform with a much higher bandwidth. First, the high frequency and very high frequency bands are applied to the mixers.F_mixer0=BWΦ_mixer0=rnd(2·π)

The frequency of the high frequency mixer is at the twice the cutoff frequency of the first band. The frequency of the very high frequency mixer is three times the cutoff frequency of the first band.F_mixer1=2·BW Φ_mixer1=rnd(2·π)Apply the mixers.xfhm_kh=xfh_kh·2·cos(2·π·F_mixer0th_kh+Φ_mixer0)xfhhm_kh=xfhh_kh·2·cos(2·π·F_mixer1th_kh+Φ_mixer1)Look at the frequency content.Xfhm=CFFT(xfhm)Xfhhm=CFFT(xfhhm)Low pass filter the mixer outputs.Xfhml={right arrow over ((XfhmMfl))}Xfhhml={right arrow over ((XfhhmMfl))}

Note again that the typical manner of low pass filtering the mixer outputs would be to use the scope front-end. This filtering is being shown here as actual low pass filters applied.

FIG. 41 shows the effect at the high frequency band mixer and image reject filter output in the frequency domain.

FIG. 42 shows the effect at the very high frequency band mixer and image reject filter output in the frequency domain.

Take the inverse FFT to generate the analog mixer output signals—the analog signals input to the channel digitizers.xfhml=ICFFT(Xfhml)xfhhml=ICFFT(Xfhhml)

FIGS. 43, 44 and 45 show the low frequency, high frequency, and very high frequency band digitizer input signal due to the processing of the simulated step input signal.

It is interesting to see what the sum of these three waveforms are—there sums to not produce anything good. This is shown in FIG. 46. This is shown to exemplify that further processing of the digitized waveforms is required to achieve the objective.

At this point, the waveforms are digitized. The waveforms must be sampled at a rate sufficient to satisfy Nyquist Criterion. For this example, this means that they must be sampled at at least 2 times BW, or 10 GS/s. After the waveforms have been digitized, they are immediately upsampled using SinX/x interpolation. This is possible because all digitized waveforms are bandlimited. It is useful to upsample the waveforms to a sample rate capable of meeting Nyquist Criterion for the system bandwidth—I have chosen 40 GS/s. The upsampling is trivial and for the purpose of this example, I simply use a 40 GS/s digitizer with the understanding that the exact same waveform would result from sampling the waveform at 10 GS/s and upsampling by a factor of 4.

FS = 40                          Upsampled digitizer sample rate.
D                    =                                                                  FS                        hi                                            FS D = 25 Upsampling factor for analog waveform model.
K                    =                                          KH                      D k = 0 . . . K − 1

Sample the waveforms.

$\begin{array}{c}{t}_k=\frac{k}{\mathrm{FS}}\\ x_{l_k}={\mathrm{xfl}}_{k·D}\\ x_{h_k}={\mathrm{xfhml}}_{k·D}\\ x_k={\mathrm{xh}}_{k·D}\\ w_k={\mathrm{wh}}_{k·D}\\ x_{{\mathrm{hh}}_k}={\mathrm{xfhhml}}_{k·D}\end{array}$

Generally, at this point, we would apply the sharp cutoff filter. If a sharp cutoff analog filter was not used, we'd have to satisfy Nyquist Criterion such that any extra frequency content would not fold back into the 5 GHz band. I've already applied a sharp cutoff filter to the analog signal, so this is not necessary.

Also, at this point, some magnitude and phase compensation would probably be necessary to account for non-deal channel frequency response characteristics. This example shows the signal digitized with ideal digitizers with ideal frequency response characteristics. Next, the high and very high frequency waveforms are mixed up to there appropriate frequency location and digitally bandpass filtered.

Note that these digital mixers know the phase of the analog mixers—some mechanism must be provided for determining this—either through a pilot tone or locking of the mixer phase to the sample clock.

Apply digital mixers.x_hmk=x_hk·(2·cos(2·π·F_mixer0t_k+Φ_mixer0))x_hhmk=x_hhk·(2·cos(2·π·F_mixer1t_k+Φ_mixer1))Bandpass filter the mixer outputs.

$\begin{array}{c}N=\frac{K}{2}\\ n=0\dots N\\ f_n=\frac{n}{N}·\frac{\mathrm{FS}}{2}\\ \mathrm{Xhm}=\mathrm{CFFT}\left(x_{\mathrm{hm}}\right)\\ \mathrm{Xhhm}=\mathrm{CFFT}\left(x_{\mathrm{hhm}}\right)\\ \mathrm{Xlm}=\mathrm{CFFT}\left(x_l\right)\\ {\mathrm{Xfhm}}_n=\mathrm{if}\left(f_n>\mathrm{BW},{\mathrm{Xhm}}_n,0\right)\\ {\mathrm{Xfhhm}}_n=\mathrm{if}\left(f_n>2·\mathrm{BW},{\mathrm{Xhhm}}_n,0\right)\\ \mathrm{nn}=1\dots N-1\\ {\mathrm{Xfhm}}_{N+\mathrm{nn}}=\stackrel{\_}{{\mathrm{Xfhm}}_{N-\mathrm{nn}}}\\ {\mathrm{Xfhhm}}_{N+\mathrm{nn}}=\stackrel{\_}{{\mathrm{Xfhhm}}_{N-\mathrm{nn}}}\\ X_h=\mathrm{CFFT}\left(x_h\right)\\ X_l=\mathrm{CFFT}\left(x_l\right)\\ X_{\mathrm{hh}}=\mathrm{CFFT}\left(x_{\mathrm{hh}}\right)\end{array}$

FIG. 47 shows the low frequency band waveform frequency content.

FIG. 48 shows the high frequency band waveform frequency content before and after digital mixing and filtering.

FIG. 49 shows the very high frequency band waveform frequency content before and after digital mixing and filtering.

FIG. 50 shows the results of summing the output waveforms. We have acquired the waveform with a 15 GHz bandwidth utilizing three 5 GHz bandwidth channels.

Now let's see how the time domain waveforms compare.xfhm=Re(ICFFT(Xfhm))xfhhm=Re(ICFFT(Xfhhm))

FIG. 51 shows a comparison of the analog low frequency portion of the input waveform to the digitized and processed low frequency waveform result. FIG. 51 shows that these are identical.

FIG. 52 shows a comparison of the analog high frequency portion of the input waveform to the mixed, digitized and digitally remixed processed high frequency waveform result. FIG. 52 shows that these are identical.

FIG. 53 shows a comparison of the analog very high frequency portion of the input waveform to the mixed, digitized and digitally remixed and processed very high frequency waveform result. FIG. 53 shows that these are identical.

FIGS. 54, 55 and 56 show how successive addition of the processed band outputs causes the resulting waveform to more and more closely approximate the input waveform.

FIG. 57 shows that the result is identical to the input waveform band limited to 15 GHz and shows that the 15 GHz bandwidth limited step is fully recreated.

It should be obvious to one skilled in the art that there are many combinations of translation frequencies and filter choices that will accomplish the objective of this invention. Each has different tradeoffs and implementation considerations depending on the specific application.

It will be understood that the above description of the present invention is susceptible to various modifications, changes and adaptations, and the same are intended to be comprehended within the meaning and range of equivalents of the appended claims. The most obvious modification, for example, is the use of more than two channels.

It will thus be seen that the objects set forth above, among those made apparent from the preceding description, are efficiently attained and, since certain changes may be made in carrying out the above method (process) without departing from the spirit and scope of the invention, it is intended that all matter contained in the above description shall be interpreted as illustrative and not in a limiting sense.

It is also to be understood that the following claims are intended to cover all of the generic and specific features of the invention herein described and all statements of the scope of the invention which, as a matter of language, might be said to fall therebetween.

Claims

1. A method for digitizing a data signal, comprising the steps of: receiving an input analog data signal spanning a frequency range;splitting the received input analog data signal into a plurality of split signals, each spanning a frequency range smaller than the original frequency range;mixing at least one of said split signals with a predetermined periodic function with a predetermined frequency, thereby translating the frequency range of the at least one split signal;digitizing each of said split signals;converting the frequency range of the at least one split signal back to its original frequency range; andcombining said digitized split signals mathematically to form a single output data stream that is a substantially correct representation of the received input signal over the original frequency range.

2. The method of claim 1, wherein the input analog data signal is split by a 50 Ohm splitter.

3. The method of claim 1, further comprising the step of band limiting each of the split signals to a corresponding one of a plurality of predetermined band ranges before mixing.

4. The method of claim 3, wherein a frequency of said predetermined periodic function is at a low side of said band of said at least one split signal that is mixed.

5. The method of claim 3, wherein a frequency of said predetermined periodic function is at a high side of said band of said at least one split signal that is mixed.

6. The method of claim 1, further comprising the step of passing the mixed split signal through an image reject filter, having a second predetermined frequency.

7. The method of claim 6, wherein said image reject filter comprises an intrinsic bandwidth of a digitized channel used for digitizing said mixed split signal.

8. The method of claim 1, wherein said mixing is accomplished by a sampling action of a digitizer used for digitizing said at least one split signal.

9. The method of claim 1, wherein said predetermined periodic function is a low-distortion sinusoid.

10. A method for digitizing a data signal, comprising the steps of: receiving an input analog data signal;splitting the received input analog data signal into channels having a low frequency split signal and a high frequency split signal respectively;mixing at least one of said low frequency and said high frequency split signals with a predetermined periodic function with a predetermined frequency;digitizing said split signals;band limiting each of the split signals to a predetermined band range;upsampling each digitized split signal to a predetermined sample rate;mixing said at least one of said low frequency and said high frequency split signals with a predetermined periodic function with said predetermined frequency, resulting in two images of said at least one split signal;passing said at least one of said high frequency and said low frequency split signals through an image reject filter, having a second predetermined frequency, to remove an unwanted one of the two images; andcombining said digitized split signals mathematically to form a single output data stream that is a substantially correct representation of the received input signal.

11. The method of claim 10, further comprising the step of passing at least one of said low frequency and high frequency split signals through a band limiting filter.

12. The method of claim 10, further comprising the step of equalizing the low frequency and the high frequency split signals separately to compensate for non-ideal magnitude and phase characteristics of front-end and digitizing systems for each channel.

13. The method of claim 10, wherein the input analog data signal is split by a 50 Ohm splitter.

14. The method of claim 10, wherein said predetermined periodic functions are low-distortion sinusoids.

15. An apparatus for digitizing a data signal, comprising: an input for receiving an input analog data signal spanning a frequency range;a splitter for splitting the received input analog data signal into a plurality of split signals, each spanning a frequency range smaller than the original frequency range;a mixer for mixing at least one of said split signals with a predetermined periodic function with a predetermined frequency, thereby translating the frequency range of the at least one split signal;a digitizer for digitizing each of said split signals;a converter for converting the frequency range of the at least one split signal back to its original frequency range; anda combining unit for combining said digitized split signals mathematically to form a single output data stream that is a substantially correct representation of the original input signal over the original frequency range.

16. The apparatus of claim 15, wherein said splitter is a 50 Ohm splitter.

17. The apparatus of claim 15, further comprising a band limiter for band limiting each of the split signals to a corresponding one of a plurality of predetermined band ranges before mixing.

18. The apparatus of claim 17, wherein a frequency of said predetermined periodic function is at a low side of said band of said at least one split signal that is mixed.

19. The apparatus of claim 17, wherein a frequency of said periodic predetermined function is at a high side of said band of said at least one split signal that is mixed.

20. The apparatus of claim 15, further comprising an image reject filter having a second predetermined frequency for receiving and passing the mixed split signal.

21. The apparatus of claim 20, wherein said image reject filter comprises an intrinsic bandwidth of a digitized channel used for digitizing said mixed split signal.

22. The apparatus of claim 15, wherein said mixing is accomplished by a sampling action of a digitizer used for digitizing said at least one split signal.

23. The apparatus of claim 15, wherein said predetermined periodic function is a low-distortion sinusoid.

24. An apparatus for digitizing a data signal, comprising: an input for receiving an input analog data signal;a splitter for splitting the received input analog data signal into channels having a low frequency split signal and a high frequency split signal respectively;a first mixer for mixing at least one of said high frequency and said low frequency split signals with a predetermined periodic function with a predetermined frequency;a digitizer for digitizing said split signals;a band limiter for band limiting each of the split signals to a predetermined band range;an upsampler for upsampling each digitized split signal to a predetermined sample rate;a second mixer for mixing said at least one of said high frequency and said low frequency split signals with a predetermined periodic function with said predetermined frequency, resulting in two images of said at least one split signal;an image reject filter, having a second predetermined frequency for receiving and passing said at least one of said high frequency and said low frequency split signals to remove an unwanted one of the two images; anda combining unit for combining said digitized split signals mathematically to form a single output data stream that is a substantially correct representation of the received input signal.

25. The apparatus of claim 24, further comprising a band limiting filter for receiving and passing at least one of said high frequency and low frequency split signals.

26. The apparatus of claim 24, further comprising an equalizer for equalizing the high frequency and the low frequency split signals separately to compensate for non-ideal magnitude and phase characteristics of front-end and digitizing systems for each channel.

27. The apparatus of claim 24, wherein said splitter is a 50 Ohm splitter.

28. The apparatus of claim 24, wherein said predetermined periodic functions are low-distortion sinusoids.