systems that sample a continuous analog time-domain signal.The invention is particularly applicable to systems that sample a continuous analog time-domain signal and whose analog components have a bandwidth limit below that desired or specified. The invention has been created to address particular problems in the design of digital sampling oscilloscopes (DSOs) that require more bandwidth than that which is easily achievable through traditionally analog techniques. A method and apparatus are provided in the form of a digital filter that is capable of surgically increasing the bandwidth of the system beyond the bandwidth achievable in an analog system. Furthermore, it is demonstrated that this system can perform this bandwidth increase without degradation in the time-domain performance of the system such as pulse or step response. In some cases, the time-domain performance is improved by flattening of the frequency response. Additionally, the system, while boosting the bandwidth, is capable of simultaneously removing noise and therefore producing a digitized signal of higher fidelity than that obtained without the filter in place.

Patent
   RE42809
Priority
Sep 01 2000
Filed
Mar 30 2005
Issued
Oct 04 2011
Expiry
Sep 26 2020
Assg.orig
Entity
Large
0
41
all paid
16. A digital filter for use with a digital sampling oscilloscope, comprising:
a first response portion providing no adjustment to a signal of the sampled system;
a second response portion providing a gain rising substantially log-linearly to said signal of the sampled system; and
a third response portion providing a predetermined attenuation to said signal of said sampled system.
8. A digital filter for increasing the bandwidth of a sampled system, comprising:
a first response portion providing no adjustment to a signal of the sampled system;
a second response portion providing a gain rising substantially log-linearly to said signal of the sampled system; and
a third response portion providing a predetermined attenuation to said signal of said sampled system.
1. A digital filter for increasing the bandwidth of a sampled system, comprising:
a first response portion providing no adjustment to a signal of the sampled system;
a second response portion providing a gain rising substantially log-linearly to said signal of the sampled system;
a third response portion providing a substantially steady gain to said signal of the sampled system;
a fourth response portion providing a gain dropping substantially log-linearly to said signal of the sampled signal; and
a fifth response portion providing predetermined attenuation to said signal of said sampled system.
0. 37. A digital oscilloscope having an upper bandwidth limit, said oscilloscope comprising:
a digital filter applied proximate the upper bandwidth limit, said filter having a frequency response in a frequency range proximate the upper bandwidth limit, the frequency response of the digital filter comprising:
a first response portion providing rising gain in a first frequency range,
a second response portion at a second frequency range higher than said first frequency range, said second response portion providing substantially constant gain, and
a third response portion at a third frequency range higher than said second frequency range, said third response portion providing decreasing gain.
0. 19. A method for increasing the bandwidth of a digital oscilloscope, comprising:
providing the digital oscilloscope having a frequency response with an upper bandwidth limit; and
applying a digital filter to the digital oscilloscope frequency response to obtain a filtered frequency response, the digital filter having a frequency response in a frequency range proximate the upper bandwidth limit, the frequency response of the digital filter comprising:
a first response portion providing rising gain in a first frequency range,
a second response portion at a second frequency range higher than said first frequency range, said second response portion providing substantially constant gain, and
a third response portion at a third frequency range higher than said second frequency range, said third response portion providing decreasing gain;
wherein the filtered frequency response has an increased upper bandwidth limit.
2. The digital filter of claim 1, wherein said bandwidth of said sampled system is increased to approximately 2 GHz.
3. The digital filter of claim 1, wherein said first response portion extends substantially from 0 to 1.8 GHz.
4. The digital filter of claim 1, wherein said second response portion extends from approximately 1.8 GHz to 2.0 GHz.
5. The digital filter of claim 1, wherein said fifth response portion extends substantially above 3 GHz.
6. The digital filter of claim 5, wherein said predetermined attenuation of said fifth response portion is −5 dB.
7. The digital filter of claim 1, wherein said third and fourth response portions are unspecified.
9. The digital filter of claim 8, wherein said digital filter further flattens the frequency response of said sampled system.
10. The digital filter of claim 8, wherein said digital filter further reduces out-of-band noise.
11. The digital filter of claim 8, wherein said digital filter controls a pulse/step response of the sampled system.
12. The digital filter of claim 11, comprising multiple poles and zeros.
13. The digital filter of claim 11, wherein the magnitude in the second response portion is specified as the inverse frequency response of the original system being boosted.
14. The digital filter of claim 8, wherein the digital difference equation is implemented as biquad sections or any other filter topology.
15. The digital filter of claim 8, wherein said digital filter is represented as a Finite Impulse response (FIR) filter, either by sampling the impulse response, or any other similar method.
17. The digital filter of claim 16, wherein said digital filter is implemented in hardware realizations including FPGAs, Gate Arrays and IP Cores.
18. The digital filter of claim 16, wherein a stabilizing zero is removed.
0. 20. The method of claim 19, further comprising generating filter coefficients for the digital filter based at least upon the inverse of the digital oscilloscope frequency response proximate the upper bandwidth limit.
0. 21. The method of claim 19, wherein the digital filter is a finite impulse response filter.
0. 22. The method of claim 19, wherein the digital filter is combined with one or more other digital filters to form a composite digital filter.
0. 23. The method of claim 19, wherein the digital filter further comprises a fourth response portion at a fourth frequency range higher than said third frequency range, said fourth response portion providing substantially constant magnitude response.
0. 24. The method of claim 23, wherein the substantially constant magnitude response is attenuating.
0. 25. The method of claim 19, wherein the digital filter provides substantially no gain or attenuation below the frequency range of the first response portion.
0. 26. The method of claim 19, wherein a magnitude response of the digital filter is the substantial inverse of the digital oscilloscope amplitude response at one or more frequencies proximate the upper bandwidth limit.
0. 27. The method of claim 19, wherein the gain in the first response portion rises substantially log-linearly.
0. 28. The method of claim 19, wherein the second response portion includes a first region in which the gain rises and a second region in which the gain falls.
0. 29. The method of claim 19, wherein the third response portion includes a first region that provides gain and a second region that provides attenuation.
0. 30. The method of claim 19, wherein the step response of the digital oscilloscope after application of the digital filter is substantially consistent with the upper bandwidth limit of the filtered frequency response.
0. 31. The method of claim 19, wherein the digital filter attenuates out-of-band noise.
0. 32. The method of claim 19, wherein the digital filter substantially flattens the frequency response of the digital oscilloscope below the increased upper bandwidth limit.
0. 33. The method of claim 19, further comprising the step of providing one or more samples to the digital filter in addition to samples corresponding to a waveform to be displayed.
0. 34. The method of claim 33, wherein the additional samples allow the digital filter to substantially settle prior to filtration of the samples corresponding to a waveform to be displayed.
0. 35. The method of claim 19, further comprising the step of providing samples to the digital filter in addition to samples representing a waveform to be displayed, wherein the additional samples are excluded from one or more processing steps after application of the digital filter.
0. 36. The method of claim 35, further comprising the step of excluding these additional samples from further processing after filtering.
0. 38. The oscilloscope of claim 37, wherein the magnitude response or phase response of the digital filter is the substantial inverse of an unboosted frequency response of the oscilloscope proximate the upper bandwidth limit.
0. 39. The oscilloscope of claim 37, wherein the digital filter is a finite impulse response filter.
0. 40. The oscilloscope of claim 37, wherein the digital filter is part of a composite filter having additional frequency response portions.
0. 41. The oscilloscope of claim 37, wherein the digital filter further comprises a fourth response portion at a fourth frequency range higher than said third frequency range, said fourth response portion providing substantially constant magnitude response.
0. 42. The oscilloscope of claim 41, wherein the substantially constant magnitude response is attenuating.
0. 43. The oscilloscope of claim 37, wherein the digital filter provides substantially no gain or attenuation below the frequency range of the first response portion.
0. 44. The oscilloscope of claim 37, wherein a magnitude response of the digital filter is the substantial inverse of an unboosted frequency response of the oscilloscope proximate the upper bandwidth limit.
0. 45. The oscilloscope of claim 37, wherein the gain in the first response portion rises substantially log-linearly.
0. 46. The oscilloscope of claim 37, wherein the second response portion includes a first region in which the gain rises and a second region in which the gain falls.
0. 47. The oscilloscope of claim 37, wherein the third response portion includes a first region that provides gain and a second region that provides attenuation.
0. 48. The oscilloscope of claim 37, wherein digital filter does not substantially degrade the step response of the oscilloscope.
0. 49. The oscilloscope of claim 37, wherein the digital filter attenuates out-of-band noise.
0. 50. The oscilloscope of claim 37, wherein the digital filter substantially flattens the frequency response of the oscilloscope below the upper bandwidth limit.
0. 51. The oscilloscope of claim 37, further comprising an acquisition system to provide one or more samples to the digital filter in addition to samples corresponding to a waveform to be displayed.
0. 52. The oscilloscope of claim 51, wherein the additional samples allow the digital filter to substantially settle prior to filtration of the samples corresponding to a waveform to be displayed.
0. 53. The oscilloscope of claim 51, wherein the additional samples are excluded from one or more processing steps after application of the digital filter.

This application claims the benefit of U.S. Provisional Application Ser. No. 60/229,856, filed Sep. 1, 2000, the entire contents of which are incorporated herein by reference.

This invention related generally to the digital manipulation of a continuous time domain sample that is to be sampled in a digital oscilloscope, and more particularly to a digital filter that is capable of increasing the bandwidth of the sampling system beyond the bandwidth range achievable in an analog system.

The present state of the art deals with an attempt to increase bandwidth based upon the assumption that only analog manipulation techniques for modifying a signal to improve the bandwidth characteristics of an apparatus are possible. Other digital techniques are seen as manipulations of the signal that change the output result of the system. This results in a design methodology in which analog design engineers painstakingly design to the best of their ability analog circuitry that has high bandwidth, flat frequency response, good pulse response and is noise-free.

In many cases, these designs are extremely complicated, particularly in the design of a digital oscilloscope. Some reasons for this difficulty are:

What makes things worse is that even upon observing and confirming the existence of problems with the bandwidth, flatness, pulse-response, and noise performance in the system, little can often be done to rectify the situation. This is because circuits designed to fix such problems are often not practically realizable.

Therefore, it would be beneficial to provide an improved Digital Signal Processing (DSP) method and apparatus capable of surgically dealing with lack of bandwidth, while offering some additional control of the pulse-response and flatness, and thereby decreasing the overall noise of the system as well.

It is therefore an object of the invention to provide an improved Digital Signal Processing (DSP) method and apparatus capable of surgically dealing with lack of bandwidth, while offering some additional control of the pulse-response and flatness.

Another object of the invention is to provide an improved Digital Signal Processing (DSP) method and apparatus capable of surgically dealing with lack of bandwidth, while offering some additional control of the pulse-response and flatness, and in which the overall noise of the system can be decreased.

A still further object of the invention is to provide an improved Digital Signal Processing (DSP) method and apparatus capable of surgically dealing with lack of bandwidth, while offering some additional control of the pulse-response and flatness by increasing the bandwidth in a very controlled manner.

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.

Generally speaking, in accordance with the invention, consider a system whose frequency response is shown in FIG. 2. As is apparent in FIG. 2, the response violates the limits on a specified 2 GHz frequency response by dropping below ωz0=2·π·fz0ωp0=2·π·fp0










Considering the filter in Equation 1, it is clear that the magnitude response of such a system is described as follows:

M ( f ) = n = 0 N - 1 20 · log [ ( j · 2 · π · f ) 2 + 2 · π · fz n Qz n · ( j · 2 · π · f ) + ( 2 · π · fz n ) 2 ( j · 2 · π · f ) 2 + 2 · π · fp n Qp n · ( j · 2 · π · f ) + ( 2 · π · fp n ) 2 · ( 2 · π · fp n ) 2 ( 2 · π · fz n ) 2 ] Equation 3
The solution of the variables in Equation 2 is performed by finding these set of variables in which the magnitude response given by Equation 3 best matches the filter design criteria in the least-squares sense. Unfortunately, Equation 3 represents a non-linear function of these variables therefore requiring methods of non-linear fitting. The method used for this invention is the Levenberg-Marquardt algorithm. In order to use this algorithm, several items must be provided. First, the partial derivatives of must be provided with respect to each of the variables being solved for:

fp n M ( f ) = 20 ln ( 10 ) · [ 2 · ( Qp n ) 2 · f 4 - 2 · ( Qp n ) 2 · ( fp n ) 2 · f 2 + ( fp n ) 2 · f 2 ] [ fp n · ( Qp n ) 2 · f 4 - 2 · ( Qp n ) 2 · ( fp n ) 3 · f 2 + ( Qp n ) 2 · ( fp n ) 5 + ( fp n ) 3 · f 2 ] Equation 4 Qp n M ( f ) = 20 ln ( 10 ) · ( fp n ) 2 · f 2 [ ( Qp n ) 3 · f 4 - 2 · ( Qp n ) 3 · ( fp n ) 2 · f 2 + ( Qp n ) 3 · ( fp n ) 4 + Qp n · ( fp n ) 2 · f 2 ] Equation 5 fz n M ( f ) = - 20 ln ( 10 ) · [ f 2 · ( fz n ) 2 - 2 · f 2 · ( Qz n ) 2 · ( fz n ) 2 + 2 · f 4 · ( Qz n ) 2 ] [ fz n · ( Qz n ) 2 · f 4 + ( fz n ) 3 · f 2 - 2 · f 2 · ( Qz n ) 2 · ( fz n ) 3 + ( Qz n ) 2 · ( fz n ) 5 ] Equation 6 Qz n M ( f ) = - 20 ln ( 10 ) · ( fz n ) 2 · f 2 [ f 4 · ( Qz n ) 3 + Qz n · f 2 · ( fz n ) 2 - 2 · f 2 · ( Qz n ) 3 · ( fz n ) 2 + ( Qz n ) 3 · ( fz n ) 4 ] Equation 7
Next, the vectors containing the frequencies and the responses desired must be created:

f spec := [ 0 F sbs 4 F sbs 2 F hs F sbs F sbs + 1 4 · ( F pbs - F sbs ) F sbs + 1 2 · ( F pbs - F sbs ) F sbs + 3 4 · ( F pbs - F sbs ) F pbs F pbe F sbe 4 ] Equation 8
Note the introduction of a frequency Fhs. It is useful to control this frequency in the specification, as can be seen later. For now, assume it is three quarters of Fsbs.

The specified frequencies given are calculated specifically to provide enough points in the flat region to ensure flatness, enough points in the boost ramp region to provide a controlled boost. In short, these are the most important points in the design and therefore there are more points specified in this region.

Next, the response vector is calculated. The response vector contains the desired response at each of these frequency points. If controlling of flatness in the system is desired, frequencies and responses may be contrived which are the negative of the actual, unboosted response. I have chosen to use the specification as shown in FIG. 3 which can be described programmatically as:

[ [ H des ( f ) := return 0 if f < F sbs [ return ( log ( f F sbs ) log ( F pbs sbs ) ) · B ] if F pbs > f F sbs return B if F pbe f F pbs [ return B - ( log ( f F pbe ) log ( F pbe F sbe ) ) · - ( B - A ) ] if F sbe > f F pbe return A if f F sbe - 100 ] ] H des ( f ) := return 0 if f < F sbs [ return ( log ( f F sbs ) log ( F pbs F sbs ) ) · B ] if F pbs > f F sbs return B if F pbe f F pbs [ return B - ( log ( f F pbe ) log ( F pbe F sbe ) ) · - ( B - A ) ] if F sbe > f F pbe return A if f F sbe - 100 Equation 9
The response vector is generated as:
Mspec=Mdes(fspec)  Equation 10
Finally, the Levenberg-Marquardt algorithm requires a vector of guesses at the values of the variables being solved for. Since Table 2 outlined the approximate pole and zero locations, it seems reasonable to use these approximations in the fitting algorithm. It also seems reasonable to select Q values that are somewhat high because of the sharp changes in the filter criteria. Therefore, the guess vector becomes:

g = ( F pbs 5 F sbs 5 F pbe 5 F sbe 5 ) Equation 11

Note the order of the values in this vector. This vector supplies the initial guesses at the variables, and is altered on each iteration of the Levenberg-Marquardt algorithm until convergence is achieved. Therefore, subsequently, the values of this guess vector are assumed to correspond to the following variables being solved for:

g 0 = fp 0 g 4 = fp 1 g 1 = Qp 0 g 5 = Qp 1 g 2 = fz 0 g 6 = fz 1 g 3 = Qz 0 g 7 = Qz 1 Equation 12
Upon determining these guesses, all that is necessary is to run multiple iterations of the Levenberg-Marquardt algorithm. Typically, iteration is halted once the mean-squared error has become small enough. An iteration of the Levenberg-Marquardt algorithm is shown here:

g i = for k ε0…K - 1 for each point specified R k M ( f spec k , g ) - M spec k Generate a residual W k , k 1 Generate the weights matrix J k , 0 fp 0 M ( f spec k , g ) Fill in a row of the Jacobian J k , 1 Qp 0 M ( f spec k , g ) matrix with the partial derivatives J k , 2 fz 0 M ( f spec k , g ) with respect to each variable J k , 3 Qz 0 J k , 4 fp 1 M ( f spec k , g ) J k , 5 Qp 1 M ( f spec k , g ) J k , 6 fz 1 M ( f spec k , g ) J k , 7 Qz 1 M ( f spec k , g ) H J T · W · J Generate the approximate Hessian for v ε0…V - 1 matrix , and a diagonal matrix whose D v , v H v , v elements are those on diagonal . ms ɛ i k = 0 SP - 1 ( R k ) 2 Calculate the mean - squ ared error Δ P ( H + λ · D ) - 1 · Calculate the delta to th e var iables J T · W · R b eing solved for and apply the delta g g - Δ P ms ɛ f k = 0 SP - 1 ( M ( f spec k , g ) - H spec k ) 2 Calculate the new mean - squ ared error λ λ · 10 if ms ɛ f > ms ɛ i If th e new mean squared error improve d , λ λ · 10 if ms ɛ f > ms ɛ i favor N ewton - G auss con vergence , otherwise favor steepest decent convergence g Equation 13
For the filter specified in FIG. 3, the result of successive iterations of this algorithm result in the following filter values:

[ [ [ [ f p = ( 1 , 935 2 , 515 ) Q p = ( 4 , 767 3 21 ) f z = ( 1 85 3 , 049 ) Q z = ( 5 , 557 3 59 ) ] ] f p = ( 1.935 2.515 ) Q p = ( 4.767 3.21 ) f z = ( 1.85 3.049 ) Q z = ( 5.557 3.59 )
The poles and zero locations of this filter are the roots of the four equations of the following form:

s 2 + 2 · π · f Q · s + ( 2 · π · f ) 2 Equation 14
This form is shown in FIG. 6. Note that the magnitudes are shown as frequency, not ω which is sometimes customary.

The response of this filter can be compared to the filter specifications and is shown in FIG. 7. It can be seen that the analog prototype filter matches the specifications quite well.

Some observations regarding FIG. 7 are discussed below. As implied previously, this filter is a valid analog filter design and could be implemented using actual electronic components such as resistors, capacitors and inductors. For the purpose of this invention, this filter will be implemented as a digital filter with the design so far providing an analog prototype filter.

The method chosen here for conversion to a digital filter is the bilinear transformation. In order to make this transformation, the pole and zero locations must be pre-warped to account for the nonlinear frequency mapping enforced by the bilinear transform:

α = - 2 · f s · ( - α f s - 1 ) ( - α f s + 1 ) Equation 15
Where fs is the sampling rate of the system (in GS/s, in this case) and α is the pole or zero being pre-warped. Note that at this point in the filter design, the sampling rate of the system must be known. In the implementation of this invention in a digital oscilloscope where the sample rate is variable, the design steps starting with the application of Equation 15 are performed dynamically within the oscilloscope itself as the sample rate is changed.

After the pre-warping, the new f and Q values are calculated. These are calculated as follows:

f = α 2 · π Q = 1 2 · cos ( π - arg ( α ) ) Equation 16
Note that only one of the complex conjugate pairs of the two sets of poles and zeros need be considered in Equation 16.

In order to convert this prototype into a digital filter, the transfer function described in Equation 1 must be factored in s and placed in the following form:

H ( s ) = n = 0 2 - N a n · s n n = 0 2 - N b n · s n Equation 17
Once this is done, the filter coefficients are calculated as:

ωz 0 = 2 · π · fz 0 ωz 1 = 2 · π · fz 1 D num = ( ωz 0 ) 2 · ( ωz 1 ) 2 a 4 = 1 D num a 3 = ( ωz 0 Qz 0 + ωz 1 Qz 1 ) · 1 D num a 2 = [ ( ωz 0 ) 2 + ωz 0 · ωz 1 Qz 0 · Qz 1 + ( ωz 1 ) 2 ) · 1 D num a 1 = [ ( ωz 0 ) 2 · ωz 1 Qz 1 + ( ωz 1 ) 2 · ωz 0 Qz 0 ] · 1 D num a 0 = 1 ωp 0 = 2 · π · fp 0 ωp 1 = 2 · π · fp 1 D den = ( ωp 0 ) 2 · ( ωp 1 ) 2 b 4 = 1 D den b 3 = ( ωp 0 Qp 0 + ωp 1 ωp 1 ) · 1 D den b 2 = [ ( ωp 0 ) 2 + ωp 0 · ωp 1 Qp 0 · Qp 1 + ( ωp 1 ) 2 ) · 1 D den b 1 = [ ( ωp 0 ) 2 · ωp 1 Qp 1 + ( ωp 1 ) 2 · ωp 0 Qp 0 ] · 1 D den b 0 = 1 Equation 18
Finally, the analog filter coefficients are converted to digital filter coefficients using the bilinear coefficient formulae:

BF ( i , n , N ) := 2 i · k = max ( n - N + i , 0 ) min ( i , n ) i ! · N - i ! k ! · ( i - k ) ! · ( n - k ) ! · ( N - i - n + k ) ! · ( - 1 ) k A n = i = 0 2 - N a i · f s i · BF ( i , n , N ) B n = i = 0 2 - N b i · f s i · BF ( i , n , N ) A := A B 0 B := B B 0 Equation 19
For a sample rate of 50 GS/s, the filter coefficients for the design specified are calculated as:

A = ( 0.746 - 2.737 3.892 - 2.539 0.643 ) B = ( 1 - 3 .705 5 .286 - 3 .437 0 .861 )
The final Z domain transfer function is of the form:

H ( z ) = n = 0 2 - N A n · z - n n = 0 2 - N B n · z - n Equation 20

The affect of this filter on the overall magnitude response of the system is shown in FIG. 9. The digital filter performance is shown in FIG. 8. The best form of the difference equation which implements the filter is found to be:

Equation 21 [ [ y k = A 0 · x k + n = 1 2 · N ( A n · x k - 1 - B n · y k - 1 ) ] ] y k = A 0 · x k + n = 1 2 · N ( A n · x k - n - B n · y k - n )
Where xk is the data sampled by the digitizing system and yk is the data at the output of the boost filter. The filter implementation is Infinite Impulse Response (IIR).

1.1 Noise Boost in the Pass-band

The application of this filter will boost, along with the signal, any noise contained in the boost ramp region and boost region. Therefore, before application of this filter, the unboosted system noise profile must be analyzed to determine the applicability of this filter. It is important to note that while boosting noise in these regions, the filter also attenuates noise beyond the pass-band of the system. This may or may not result in an overall noise performance improvement, depending on the noise sources.

In order to check this, a 260.6 mV rms sine wave at 2 GHz is applied using an RF signal generator and a sine-fit is applied. The sinefit tells the Signal to Noise Ratio (SNR) and also the Effective Number of Bits (ENOB). FIG. 12 and FIG. 13 show a comparison with and without the boost filter in place. FIG. 12, without the boost filter, shows an rms value of 131 mV, an attenuation of approximately 6 dB, and ENOB and SNR of 4.67 and 23, respectively. FIG. 13, with the boost filter, shows an rms value of 226.2 mV, an attenuation of 1.23 dB, and an ENOB and SNR of 5.35 and 32.1, respectively. This means that for this application, the boost filter resulted in an increase in the effective bits by over 0.6 while making the bandwidth specification. In other words, the bandwidth was boosted while simultaneously improving the noise performance of the unit.

1.2 Nyquist Limitation and Stabilizing Zero Placement

As with all digital filters, there is the limitation of sample rate on the system. Because of Nyquist's criteria, a digital filter can perform like an analog system up to ½ the sample rate of the system, after which the filter response repeats over and over. In other words, frequencies above ½ the sample rate appear as aliases at frequencies under ½ the sample rate. Even worse, the conversion of the analog prototype filter has big problems using the bilinear transformation if any pole or zero locations are above the Nyquist rate.

Fc Fe in the filter specification must appear at or below the Nyquist rate. The determination of the location of the stabilizing zero is performed mostly through the specification of the attenuation at Fsbe. The fact that Fsbe in this design has attenuation A specified constrains the stabilizing zero to appear at a frequency between and Fsbe and Fe. While this causes the objectionable decrease in attenuation in the attenuation region, it does help in the realization of the digital filter. The attenuation of Fsbe may be decreased, which will move the stabilizing zero higher in frequency, but the design must keep this zero below the Nyquist rate of the system, otherwise the filter design will fail.

1.3 Filter Startup

Any system employing memory will take some time to stabilize after the signal appears for the first time. This is not generally a problem, and most designs handle this by waiting some time for signals to stabilize. In the case of a digital oscilloscope, this is accomplished through pre-trigger hold-off. When the acquisition system is armed and acquiring, the trigger is held off until enough time has passed for everything to stabilize. In the case of a digital oscilloscope utilizing the present invention, there is an additional problem. There will be some stabilization time associated with the filter, and the system does not get to see the waveform until the point in time that the waveform has been acquired and is being read-out of the acquisition system memory. This means that additional samples must be acquired at the beginning of the waveform so that when the filter is applied, the system will have stabilized prior to the point at which the waveform comes on-screen. The points prior to the left edge of the oscilloscope screen where the filter is stabilizing are simply discarded. Since the filter-taps are loaded with zeros (or more commonly, the first point in the input waveform), the first point entering the filter looks like a step. Therefore, the startup time can be estimated by examining the impulse response of the filter.

The impulse response for the design example provided is shown in FIG. 14. This response was calculated using Equation 21. It is useful to note that this sampled impulse response could be used as the coefficients of a Finite Impulse Response (FIR) filter design.

In order to determine the precise startup time required, FIG. 15 is provided which shows the magnitude of the amplitude with respect to time plotted log-log. From this plot you can see that a filter startup time allowance of 3 ns allows the system to stabilize to within 0.1% of its final value.

1.4 The Importance of the Stabilizing Zero in the Digital Filter Design.

The final zero in the system, the stabilizing zero, has the effect of leveling off the attenuation of the system. Some might consider this objectionable, preferring the attenuation to continue, and thus gaining the maximum noise attenuation. This is possible, however some problems arise in the removal of this zero.

First, Equation 15, which provides the pre-warping equation, works only with an equal number of poles and zeros in the system. This is a huge benefit, because Equation 15 provides a digital filter whose frequency domain performance matches almost identically the analog filter performance, even when the frequencies of interest are up near the Nyquist rate of the system. Therefore, removal of the stabilizing zero will cause difficulties in matching the digital filter to the analog prototype filter.

Second, keep in mind that digital filters repeat above the Nyquist rate—they first fold about the Nyquist frequency up to the sample rate of the system, after which the filter image repeats over and over again. This means that a boost at a frequency at 2 GHz in a system sampling at 8 GS/s, for example, will also boost frequencies at 6 GHz. Typically, only noise is present at this frequency, but the boosting of this noise in conjunction with the obliteration of all other noise might cause artifacts which are objectionable.

1.5 The Placement of Fhs.

Equation 8 makes reference to a frequency in the design specification called Fhs. Note in FIG. 7, that the placement of this frequency has an effect on the performance of the filter between Fhs and Fsbs. The performance of the filter in this range of frequencies is unconstrained, and this lack of constraint is utilized in the design by allowing a dip in the gain of the system during this region. In other words, Fhs must be placed at a frequency such that in between this frequency, and the start of the boost ramp region, this slight attenuation can be tolerated. As shown in FIG. 9, a slight bite has been taken out of the un-boosted response around 1.7 GHz. Over-constraining the system by pushing Fhs up close to effectively enforcing a sharper edge, will generate a system that has higher Q values which is usually undesirable.

1.6 The Compromise Between Noise Performance and Pulse-response.

In this design in accordance with the invention, there are particular frequency response specifications that have significant implications with regard to the time-domain performance and to the noise reduction. Out to the boost ramp region, any smoothening of the system roll-off will improve the pulse response and reduce overshoot while simultaneously reducing noise. Beyond this point, a trade-off must be made regarding noise and pulse-response. Since the purpose of this invention is the boosting of the bandwidth of the system, the implication is that the system frequency response is rolling off around the boost point, and the roll-off rate is generally increasing. In situations like this, the pulse response will only be worsened unless the filter prior to the boost point has already applied some attenuation. In any case, the roll-off rate of the system can be controlled to some extent.

It has already been noted that the stabilizing zero does not allow the system to reach it's full potential with regard to noise performance. This zero does, however, provide the benefit of controlling the roll-off in the boost drop region. This is a good thing with regard to pulse response. It can be shown that with proper selection of the frequency and response for the frequencies Fpbe and Fsbe, and the attenuation at Fe, the trade-off between system roll-off (and thus pulse-response performance) and noise attenuation outside of the pass-band can be made. The placement of these frequencies and their responses must be made with care, however, in order to maintain the ability to fit the analog filter to the design criteria and to provide an analog filter which is realizable as a digital filter with the given sample rate constraints.

In FIG. 9, it can be seen that the roll-off was not completely controlled, but the lack of discussion of this point simplified the explanation of the invention. This lack of control of the roll-off resulted in less than optimal noise attenuation performance and slight worsening of the pulse response as can be seen by comparing FIG. 10 and FIG. 11. Furthermore, the present invention is being applied to a system that requires 3 dB margin at the bandwidth specified in order to meet the bandwidth specification over the temperature range of the oscilloscope.

1.7 The Inverse Response Specification

The filter design specification as shown in FIG. 3 is very easy to understand conceptually and provides a rather generic filter spec. It is possible to tailor the specification to a particular unit's frequency response by substituting the negative of the frequency response for the unit for the response provided in Equation 9. This would enable virtually complete flattening/control of the frequency response. In order to do this, however, there needs to be a complete understanding of the response with regard to the filter order required to flatten the response. In general, the more bumps in the frequency response, the more poles and zeros (and thus a higher order) is required for the boost filter. The ability of the boost filter to flatten the response by specifying the inverse system response is implicit to the present invention.

In the implementation of any high order filter, it is sometimes useful to separate the filter into sections, effectively cascading sections of lower order. Usually the sections are second order biquad sections. This was not deemed necessary for this design in accordance with the invention. The method of separating this filter into multiple, lower order, cascaded sections is well known by those practiced in the art of digital signal processing.

It will thus be seen that the objects set forth above, among those made apparent from the preceding description, are efficiently attained and, because certain changes may be made in carrying out the above method and in the construction(s) set forth without departing from the spirit and scope of the invention, it is intended that all matter contained in the above description and shown in the accompanying drawings 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.

Pupalaikis, Peter J.

Patent Priority Assignee Title
Patent Priority Assignee Title
4788653, Dec 23 1986 General Electric Company Digital filter for power system stabilizer
4875166, Oct 09 1987 Input/Output, Inc.; INPUT OUTPUT, INC Bandwidth enhancing seismic acquisition system and method
5208596, Apr 10 1992 PLYMOUTH DEWITT, INC DAC distortion compensation
5226059, Sep 07 1990 Nortel Networks Limited DSP line equalizer
5239578, May 15 1990 Plantronics, Inc. Noise cancelling apparatus for a telephone handset
5280353, Aug 15 1989 CINTEL INTERNATIONAL LIMITED, A U K CORP Method and apparatus for signal processing by interpolation and filtering with simultaneous frequency response compensation and offset generation
5283483, Jan 27 1993 Fairchild Semiconductor Corporation Slimmer circuit technique
5341177, Sep 04 1991 SAMSUNG ELECTRONICS CO , LTD , A CORP OF THE REPUBLIC OF KOREA System to cancel ghosts generated by multipath transmission of television signals
5388062, May 06 1993 Thomson Consumer Electronics, Inc Reconfigurable programmable digital filter architecture useful in communication receiver
5487023, Feb 14 1994 Tektronix, Inc. Repeatable finite and infinite impulse response integrated circuit structure
5574639, Oct 12 1994 Los Alamos National Security, LLC System and method for constructing filters for detecting signals whose frequency content varies with time
5577117, Jun 09 1994 Nortel Networks Limited Methods and apparatus for estimating and adjusting the frequency response of telecommunications channels
5698984, Jan 30 1996 Fluke Corporation Adaptive digital filter for improved measurement accuracy in an electronic instrument
5754437, Sep 10 1996 Tektronix, Inc.; Tektronix, Inc Phase measurement apparatus and method
5812009, Apr 03 1995 Fujitsu Limited Boost type equalizing circuit
5854756, Feb 01 1996 Racal-Datacom Limited Digital filters
5978742, Apr 04 1997 Tektronix, Inc. Method and apparatus for digital sampling of electrical waveforms
6115418, Feb 09 1998 National Semiconductor Corporation Simplified equalizer for twisted pair channel
6175849, Feb 10 1998 WSOU Investments, LLC System for digital filtering in a fixed number of clock cycles
6184748, Sep 30 1998 AVAGO TECHNOLOGIES INTERNATIONAL SALES PTE LIMITED Magnitude and group delay shaping circuit in continuous-time read channel filters
6219392, Jan 20 1999 Matsushita Electric Industrial., Ltd. Filter device and method of acquiring filter coefficients
6279021, Jan 30 1998 Sanyo Electric Co. Ltd. Digital filters
6289063, Sep 02 1998 RPX CLEARINGHOUSE LLC QAM receiver with improved immunity to crosstalk noise
6304134, Mar 29 2000 Texas Instruments Incorporated High-frequency boost technique
6437932, May 16 1996 AVAGO TECHNOLOGIES GENERAL IP SINGAPORE PTE LTD Decision based time-varying equalizers
6539318, Nov 17 2000 TELEDYNE LECROY, INC Streaming architecture for waveform processing
6542914, Sep 01 2000 TELEDYNE LECROY, INC Method and apparatus for increasing bandwidth in sampled systems
6581080, Apr 16 1999 SONNOX LIMITED Digital filters
6631175, Apr 03 1998 TELECOM HOLDING PARENT LLC Spectrally constrained impulse shortening filter for a discrete multi-tone receiver
20020009135,
20020075267,
20020109496,
20020163959,
20020171408,
20030002574,
20030007583,
20030055855,
20030118094,
20050080831,
JP10282153,
WO38384,
//////
Executed onAssignorAssigneeConveyanceFrameReelDoc
Mar 30 2005LeCroy Corporation(assignment on the face of the patent)
Mar 30 2007LeCroy CorporationMANUFACTURERS AND TRADERS TRUST COMPANYSECURITY AGREEMENT0193310239 pdf
Jul 29 2010LeCroy CorporationMANUFACTURERS AND TRADERS TRUST COMPANYSECURITY AGREEMENT0248920689 pdf
Aug 03 2012LeCroy CorporationTELEDYNE LECROY, INC MERGER SEE DOCUMENT FOR DETAILS 0291620724 pdf
Aug 22 2012RBS CITIZENS, N A TELEDYNE LECROY, INC RELEASE BY SECURED PARTY SEE DOCUMENT FOR DETAILS 0291550478 pdf
Oct 09 2012MANUFACTURERS AND TRADERS TRUST COMPANY, AS AGENTLeCroy CorporationRELEASE BY SECURED PARTY SEE DOCUMENT FOR DETAILS 0291280280 pdf
Date Maintenance Fee Events
Feb 12 2013STOL: Pat Hldr no Longer Claims Small Ent Stat
Oct 01 2014M1553: Payment of Maintenance Fee, 12th Year, Large Entity.


Date Maintenance Schedule
Oct 04 20144 years fee payment window open
Apr 04 20156 months grace period start (w surcharge)
Oct 04 2015patent expiry (for year 4)
Oct 04 20172 years to revive unintentionally abandoned end. (for year 4)
Oct 04 20188 years fee payment window open
Apr 04 20196 months grace period start (w surcharge)
Oct 04 2019patent expiry (for year 8)
Oct 04 20212 years to revive unintentionally abandoned end. (for year 8)
Oct 04 202212 years fee payment window open
Apr 04 20236 months grace period start (w surcharge)
Oct 04 2023patent expiry (for year 12)
Oct 04 20252 years to revive unintentionally abandoned end. (for year 12)