USH1458H - Signal amplitude distribution analyzer - Google Patents
Signal amplitude distribution analyzer Download PDFInfo
- Publication number
- USH1458H USH1458H US08/083,962 US8396293A USH1458H US H1458 H USH1458 H US H1458H US 8396293 A US8396293 A US 8396293A US H1458 H USH1458 H US H1458H
- Authority
- US
- United States
- Prior art keywords
- voltage
- signals
- signal
- amplitude
- window
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 claims description 13
- 230000001186 cumulative effect Effects 0.000 abstract description 20
- 238000005259 measurement Methods 0.000 abstract description 20
- 238000004458 analytical method Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 18
- 230000003252 repetitive effect Effects 0.000 description 8
- 238000012360 testing method Methods 0.000 description 5
- 238000005315 distribution function Methods 0.000 description 4
- 239000003990 capacitor Substances 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/26—Measuring noise figure; Measuring signal-to-noise ratio
Definitions
- This invention relates generally to the field of signal analyzers and, more particularly, to a signal analyzer designed to measure the amplitude probability distribution of electrical signals.
- the invention relates to a signal analyzer designed to measure the amplitude probability density function of noise signals to allow characterization of the signal's statistics in terms of skewness, clipping, kurtosis and other non-Gaussian deviations.
- the Gaussianicity of noise signals was estimated by measuring the signal's cumulative amplitude distribution. This measurement was done through repetitive sampling of signal amplitude in combination with computer analysis, or by a simple comparator circuit such as that described on pages 272-274 of the Dec. 11, 1986 EDN magazine article titled, "Test whether a noise source is Gaussian.”
- the cumulative amplitude distribution function gives the probability that a variable will assume a value equal to or greater than a particular value over a range of values.
- the relative Gaussianicity of noise signals is determined by the degree of conformance of measured data to a theoretical cumulative distribution curve.
- a limitation of the cumulative amplitude distribution method is that measured data do not show the actual distribution of amplitudes. Because of this, departures of data from the theoretical curve do not characterize precisely how a noise signal is different from true Gaussian noise. As will be visually depicted in this disclosure, this method indicates non-Gaussianicity only by a departure of the signal's plotted distribution from the theoretical cumulative amplitude distribution curve. Abnormal characteristics of the noise signal such as skewness, kurtosis, clipping and other non-Gaussian deviations cannot be readily identified.
- the cumulative amplitude distribution method can also be relatively insensitive and give results that are difficult to relate to in Gaussian terms, thereby making it difficult to judge the overall degree of non-Gaussianicity of noise signals.
- Another probability distribution is the amplitude probability density distribution. This distribution gives the probability that a variable will assume a value near any particular value in its range of values.
- the amplitude probability density distribution is determined through an invention described in U.S. Pat. No. 3,581,200.
- This invention produces a probability density function profile through the use of a wave generator, spectrum analyzer, sweep generator and x-y recorder.
- the '200 invention converts signal amplitude distributions into a frequency domain profile for visualization by way of the spectrum analyzer.
- the invention is of relative high complexity, and nature of the design appears to make system calibration difficult.
- Such an invention should be able to readily permit the perception of noise signal skewness, kurtosis, clipping, as well as other non-Gaussian deviations.
- the signal amplitude distribution analyzer of the invention is designed to measure the amplitude probability density function of electrical signals.
- This invention has been particularly devised to measure the amplitude probability distribution of noise signals. Such measurements may be used to determine the "Gaussianicity" of noise signals, that is, a measurement of how closely the amplitude density distribution of noise signals corresponds to theoretical values derived from the Gaussian probability density distribution function (i.e. a normal distribution curve).
- This theoretical density function represents the relative percentage of time that a noise signal is at a given amplitude.
- the invention gives an approximation of this function by measuring the amount of time a noise signal is between a window of two adjustable voltage levels. This is accomplished by producing an output voltage proportional to the amount of time a noise signal amplitude falls within the window of values defined by the two adjustable voltage levels. Both repetitive and non-repetitive signals may be measured.
- the center point of the window is then plotted versus the invention's output voltage, giving the amplitude distribution density. This is then compared to the theoretical density function, plotted on the same graph, to determine signal "Gaussianicity.”
- the invention is calibrated to give an output that is Gaussian when analyzing a known Gaussian input. A continuous resolution measurement of the signal amplitude probability density function is possible, permitting accurate analysis of noise statistics in terms of skewness, clipping, etc. The cumulative amplitude distribution of noise signals can also be measured.
- It is an object of this invention is to provide an improved signal analyzer.
- Another object of this invention is to provide an improved signal analyzer that is simple in operation.
- Yet another object of this invention is to provide an improved signal analyzer that is capable of continuous resolution.
- a further object of this invention is to provide an improved signal analyzer that may be easily calibrated so as to provide a quantitative measurement.
- Still a further object of this invention is to provide an improved signal analyzer that measures amplitude probability density functions/probability density distributions.
- FIG. 1 is a schematic of a representative embodiment of the signal amplitude distribution analyzer of the invention.
- FIG. 2 depicts waveforms useful in describing the calibration procedure of the invention.
- FIG. 3 there are plotted outputs of the invention as compared to a theoretically derived signal amplitude probability density distribution and a cumulative amplitude distribution.
- FIG. 1 a schematic diagram of signal amplitude distribution analyzer 10 is shown.
- all resistors are in ohms and capacitors are in microfarads ( ⁇ F).
- All 3 integrated circuit (IC) power connections are decoupled by 100 ohm, and 0.1 ⁇ F filter networks.
- resistor R1 provides a 50-ohm termination for signal source 12 being tested. The source is fed as an input signal, V in , that is then connected to an inverting input of integrated circuit 1 (ICl), pin 4, and the non-inverting input of IC2 (pin 3) through isolating resistors R2 and R3, respectively.
- Integrated circuits IC1 and IC2 are high speed comparators such as LM361s.
- the reference voltage for IC1 and IC2, V x is set by resistors R4, R5 and R6, and is applied to pins 3 and 4 of IC1 and IC2, respectively.
- Capacitor C1 provides filtering so that V x remains at the voltage set by R5 in spite of any feedback of comparator switching transients.
- the other inputs to IC1 and IC2 are variable voltages that compensate for IC offset voltages and currents, and also set comparator thresholds.
- the voltage (V H ) at pin 4 of ICl is primarily determined by a voltage divider including R9, R10, R2, R1 and the signal source resistance.
- the voltage (V L ) at pin 3 of IC2 is developed similarly by R7, R8, R3, R1 and the signal source resistance.
- the voltage at pin 4 of IC1 (V H ) is normally set to some small negative voltage and that at pin 3 of IC2 (V L ) to a small positive voltage creating a preselected "window" of voltage levels.
- V x set for 0 volts DC (VDC)
- IC1 output will be high (positive) and IC2 will be low (nominally zero volts).
- IC1 and IC2 outputs are applied to pins 1 and 2 of an AND gate IC3, such as a 74HC08, the output at pin 3 of which will thus be low.
- the AND gate output (pin 3 of IC3) is thus high only when V in is between the range between the small negative voltage and the small positive voltage. At all other times the AND gate output is low.
- the "window" of V in voltages that produces a high output is determined by the bias voltages set at pin 4 of IC1 by R9, and at pin 3 of IC2 by R7. It should be noted that when V in decreases, the same result occurs. That is, when V in goes below the window upper level (i.e. to some small positive voltage), pin 4 of IC1 will become more negative than pin 3, giving a high output at IC1, and at IC3. But when V in drops below the window lower level (i.e. beyond a small negative voltage), pin 3 of IC2 will go negative, producing a low IC2 and IC3 output.
- the output of AND gate IC3 is thus positive only when V in transitions through the window of voltages set by R7 and R9.
- DMM 14 effectively integrates the filtered output pulses from pin 3 of IC3. As explained above, IC3 will produce pulse lengths equal to the amount of time that V in is within the preset window voltage range. DMM 14, therefore, reads a voltage that is proportional to the amount of time V in is within this window. Put another way, the DMM indicates the relative amount of time that an input signal is at a given input level range.
- the window will simply be raised or lowered.
- the window voltage range i.e. the difference between the window lower and upper voltage levels
- DMM 14 will indicate the relative amount of time that the signal spends at these levels.
- the window lower voltage, V L , and higher voltage, V H are set symmetrically about V x .
- the digital multimeter therefore indicates the relative amount of time that V in equals V x , where V x is approximated by a window of voltages between V L and V H .
- V out ⁇ 10 is labeled as "V out ⁇ 10" since in this representative embodiment, the invention is normally calibrated by R12 to give a DMM reading of 0.286 VDC with the 140 mv rms noise voltage in, with V x set to 0 VDC. This makes measurements easily related to the theoretical curve 16 (i.e. DMM reading times 10), while not requiring an additional amplifier with a gain of 10.
- the invention is easily calibrated by applying a sinewave at a frequency of about 4 kHz to the V in input, and setting V x equal to 0 VDC.
- the input signal, V in , and the output pulses at pin 3 of IC3, V W are simultaneously viewed on an instrument such as a dual trace oscilloscope.
- Potentiometer R7 is then set so that the V W pulses rise when V in is at a -V L level
- Resistor R12 is then adjusted, as noted above, for 0.286 VDC output with 140 mv rms of noise in.
- measurements can be done as follows. First, the noise source to be tested, signal source 12, is connected to analyzer 10's V in input. The output of the noise source is adjusted for an amplitude of 140 mv rms, as read on a true rms voltmeter. The value of V x is then varied over the range, for example, -0.5 to +0.5 VDC, and V out ⁇ 10 is recorded for each value of V x tested. By closely spacing the values of V x continuous measurement resolution is possible. These measurements are plotted on a graph such as FIG. 3 so that measured curves can be directly compared to theoretical curve 16 for P(v). FIG. 3 shows examples of such test data.
- the measured values for a good noise source fall very closely to theoretical curve 16.
- the curve for a noise source known to be faulty shows a marked deviation from Gaussian theory. This curve (18), indicates that the noise signal contains inordinately high negative peaks, positive peak clipping, and a positively skewed, leptokurtic distribution. All of these non-Gaussian characteristics could be cross-correlated with an oscilloscope display of the signal. It should be noted that none of these characteristics would be readily visible from a cumulative amplitude distribution of the faulty noise source.
- the signal amplitude distribution analyzer of the invention can also be used to measure cumulative amplitude distribution. This is done by setting switch S1 to the CUM position, thereby disconnecting the output from IC1 and connecting the pin 1 AND input of IC3 to a high level of +5 VDC. This in effect makes IC3 just a buffer stage whose output is the same as the input of its pin 2. The circuit is then recalibrated by first terminating the V in input in 50 ohms (no signal), setting V x to some negative value, and setting R12 for a 1 VDC DMM reading.
- the 1 VDC DMM reading is set to indicate a 100% cumulative distribution value.
- V x is varied over a selected range such as 0 to +0.5 VDC and the cumulative distribution is read on the DMM for a number of V x values.
- the recorded values of V x and V out are then co-plotted with the theoretical Gaussian cumulative distribution curve 20 as shown on the left of FIG. 3.
- theoretical curve 20 is defined by P(cum v), which is the integral of the distribution function from V x to infinity.
- the P(cum v) theoretical curve 20 can also be computed using available tables and the formula: ##EQU2## where "erfc" is the complementary error function.
- FIG. 3 in addition to showing distribution density function profiles, also shows cumulative distribution data obtained from the same noise sources tested for the distribution density function. Note that the good noise source (plotted as "X"s) produces test data points that lie closely to theoretical curve 20. The faulty noise source data (shown as dots), on the other hand, departs from theoretical curve 20. However, unlike the data from distribution density function testing, little indication is given as to the degree or character of this departure from Gaussianicity.
- FIG. 3 provides a direct comparison of results from cumulative and distribution density function testing. It shows how much more revealing it is to measure the distribution density function itself, rather than to measure cumulative distribution, as is typically done.
- FIG. 3 also shows the measured amplitude distribution of a sinewave (22) for positive values of V x . This demonstrates that the signal amplitude distribution analyzer of the invention can be used to make amplitude distribution measurements of any waveform type whether they be repetitive or non-repetitive.
- the advantages of the signal amplitude distribution analyzer of the invention over older methods of measuring signal amplitude distribution are many.
- the invention gives the capability of measuring probability density function itself rather than cumulative distribution functions, providing a more sensitive measure as to the degree of conformance of a noise source with theory.
- the invention is therefore capable of characterizing deviations from theory rather than only indicating a degree of nonconformance.
- the invention is easily calibrated from a known input source, precise measurement of a signal source's deviations is possible. These sources may be of a repetitive or non-repetitive type, with the invention providing a continuous, overlapping, resolution to enhance measurement accuracy.
- the invention provides not only a mechanism to sense the presence of sought-after signals, but also provides the relative amount of time that an input signal is within a selected signal range. This latter relationship is exactly what is given by theoretical probability density functions.
- the signal amplitude distribution analyzer of the invention uses a very simple implementation scheme, requiring only three integrated circuit components. 0f course the invention could be made more accurate by using more complex circuitry. It also could be improved by using an analog-to-digital converter, computer interface circuits and computer analysis. The invention could also be expanded to provide automatic threshold or window stepping, so that results could be printed automatically on an X/Y recorder.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
A signal amplitude distribution analyzer measures the amplitude probability density function of electrical signals, and in particular, the amplitude probability distribution of noise signals. Such measurements may be used to determine the "Gaussianicity" of noise signals, that is, a measurement of how closely the amplitude distribution of noise signals corresponds to theoretical values derived from the Gaussian probability distribution density function. This theoretical density function represents the relative percentage of time that a noise signal is at a given amplitude. The invention gives an approximation of this function by measuring the amount of time a noise signal is between a window of two adjustable voltage levels. This is accomplished by producing an output voltage proportional to the amount of time a noise signal amplitude falls within the window of values defined by the two adjustable voltage levels. The center point of this window is then plotted versus the invention's output voltage, giving the amplitude distribution density. This is then compared to the theoretical density function, plotted on the same graph, to determine signal "Gaussianicity." The invention is calibrated to give an output that is Gaussian when analyzing a known Gaussian input. A continuous resolution measurement of the signal amplitude probability density function is possible, permitting accurate analysis of noise statistics in terms of skewness, clipping, etc. The cumulative amplitude distribution of noise signals can also be measured by the analyzer of the invention.
Description
The invention described herein may be manufactured and used by or for the Government of the United States of America for governmental purposes without the payment of any royalties thereon or therefor.
1. Field of the Invention
This invention relates generally to the field of signal analyzers and, more particularly, to a signal analyzer designed to measure the amplitude probability distribution of electrical signals. In greater particularity, the invention relates to a signal analyzer designed to measure the amplitude probability density function of noise signals to allow characterization of the signal's statistics in terms of skewness, clipping, kurtosis and other non-Gaussian deviations.
2. Description of the Related Art
In the past, the Gaussianicity of noise signals was estimated by measuring the signal's cumulative amplitude distribution. This measurement was done through repetitive sampling of signal amplitude in combination with computer analysis, or by a simple comparator circuit such as that described on pages 272-274 of the Dec. 11, 1986 EDN magazine article titled, "Test whether a noise source is Gaussian."
The cumulative amplitude distribution function gives the probability that a variable will assume a value equal to or greater than a particular value over a range of values. By the cumulative amplitude distribution method, the relative Gaussianicity of noise signals is determined by the degree of conformance of measured data to a theoretical cumulative distribution curve.
A limitation of the cumulative amplitude distribution method is that measured data do not show the actual distribution of amplitudes. Because of this, departures of data from the theoretical curve do not characterize precisely how a noise signal is different from true Gaussian noise. As will be visually depicted in this disclosure, this method indicates non-Gaussianicity only by a departure of the signal's plotted distribution from the theoretical cumulative amplitude distribution curve. Abnormal characteristics of the noise signal such as skewness, kurtosis, clipping and other non-Gaussian deviations cannot be readily identified.
The cumulative amplitude distribution method can also be relatively insensitive and give results that are difficult to relate to in Gaussian terms, thereby making it difficult to judge the overall degree of non-Gaussianicity of noise signals.
Another probability distribution is the amplitude probability density distribution. This distribution gives the probability that a variable will assume a value near any particular value in its range of values.
Techniques exist for calculating amplitude probability density distributions. One such technique is described in U.S. Pat. No. 3,626,168 issued to Keith H. Norsworthy. This patent describes an invention capable of a multitude of signal measurements. For measurements of amplitude probability density distribution, this invention outputs an indicating signal when an input signal falls between two known signal levels. The indicating signal is then apparently sent to an averaging bin to provide a digitally constructed density distribution. The '168 patent describes an invention that is highly complex and because of the use of a finite number of averaging bins, it cannot provide continuous resolution capability.
In a second scheme, the amplitude probability density distribution is determined through an invention described in U.S. Pat. No. 3,581,200. This invention produces a probability density function profile through the use of a wave generator, spectrum analyzer, sweep generator and x-y recorder. The '200 invention converts signal amplitude distributions into a frequency domain profile for visualization by way of the spectrum analyzer. The invention is of relative high complexity, and nature of the design appears to make system calibration difficult.
In a related but different area, the invention of U.S. Pat. No. 4,625,283 issued to James R. Hurley describes an invention designed to analyze repetitive signals, e.g. sinewave signals. This invention measures the elapsed times it takes for an alternating current signal to cross predetermined reference values and compares these with known values to determine characteristics of an alternating current signal being analyzed. Signal characteristics such as frequency, size of direct current offset and waveform amplitude apparently can be determined with this invention.
A need thus exists for a simple, calibrated device that permits continuous resolution of signal amplitude probability distribution. Such an invention should be able to readily permit the perception of noise signal skewness, kurtosis, clipping, as well as other non-Gaussian deviations.
The signal amplitude distribution analyzer of the invention is designed to measure the amplitude probability density function of electrical signals. This invention has been particularly devised to measure the amplitude probability distribution of noise signals. Such measurements may be used to determine the "Gaussianicity" of noise signals, that is, a measurement of how closely the amplitude density distribution of noise signals corresponds to theoretical values derived from the Gaussian probability density distribution function (i.e. a normal distribution curve). This theoretical density function represents the relative percentage of time that a noise signal is at a given amplitude. The invention gives an approximation of this function by measuring the amount of time a noise signal is between a window of two adjustable voltage levels. This is accomplished by producing an output voltage proportional to the amount of time a noise signal amplitude falls within the window of values defined by the two adjustable voltage levels. Both repetitive and non-repetitive signals may be measured.
The center point of the window is then plotted versus the invention's output voltage, giving the amplitude distribution density. This is then compared to the theoretical density function, plotted on the same graph, to determine signal "Gaussianicity." The invention is calibrated to give an output that is Gaussian when analyzing a known Gaussian input. A continuous resolution measurement of the signal amplitude probability density function is possible, permitting accurate analysis of noise statistics in terms of skewness, clipping, etc. The cumulative amplitude distribution of noise signals can also be measured.
It is an object of this invention is to provide an improved signal analyzer.
Another object of this invention is to provide an improved signal analyzer that is simple in operation.
Yet another object of this invention is to provide an improved signal analyzer that is capable of continuous resolution.
A further object of this invention is to provide an improved signal analyzer that may be easily calibrated so as to provide a quantitative measurement.
Still a further object of this invention is to provide an improved signal analyzer that measures amplitude probability density functions/probability density distributions.
Other objects, advantages and new features of the invention will become apparent from the following detailed description of the invention when considered in conjunction with the accompanied drawings.
FIG. 1 is a schematic of a representative embodiment of the signal amplitude distribution analyzer of the invention.
FIG. 2 depicts waveforms useful in describing the calibration procedure of the invention.
In FIG. 3 there are plotted outputs of the invention as compared to a theoretically derived signal amplitude probability density distribution and a cumulative amplitude distribution.
Referring now to FIG. 1 a schematic diagram of signal amplitude distribution analyzer 10 is shown. In the representative embodiment of the invention shown, all resistors are in ohms and capacitors are in microfarads (μF). All 3 integrated circuit (IC) power connections are decoupled by 100 ohm, and 0.1 μF filter networks. In FIG. 1 resistor R1 provides a 50-ohm termination for signal source 12 being tested. The source is fed as an input signal, Vin, that is then connected to an inverting input of integrated circuit 1 (ICl), pin 4, and the non-inverting input of IC2 (pin 3) through isolating resistors R2 and R3, respectively. Integrated circuits IC1 and IC2 are high speed comparators such as LM361s. The reference voltage for IC1 and IC2, Vx, about which a noise signal is analyzed, is set by resistors R4, R5 and R6, and is applied to pins 3 and 4 of IC1 and IC2, respectively. Capacitor C1 provides filtering so that Vx remains at the voltage set by R5 in spite of any feedback of comparator switching transients. The other inputs to IC1 and IC2 are variable voltages that compensate for IC offset voltages and currents, and also set comparator thresholds. The voltage (VH) at pin 4 of ICl is primarily determined by a voltage divider including R9, R10, R2, R1 and the signal source resistance. The voltage (VL) at pin 3 of IC2 is developed similarly by R7, R8, R3, R1 and the signal source resistance.
Referring now to FIGS. 1 and 2, the voltage at pin 4 of IC1 (VH) is normally set to some small negative voltage and that at pin 3 of IC2 (VL) to a small positive voltage creating a preselected "window" of voltage levels. With Vx set for 0 volts DC (VDC), then, when Vin is very negative, IC1 output will be high (positive) and IC2 will be low (nominally zero volts). IC1 and IC2 outputs are applied to pins 1 and 2 of an AND gate IC3, such as a 74HC08, the output at pin 3 of which will thus be low. When Vin rises to a negative voltage slightly smaller than the positive bias at pin 3 of IC2, pin 3 becomes positive with respect to pin 4 and IC2 goes high. Both inputs to IC3 are then high, producing a high output. As Vin continues to rise, however, it will become positive enough to overcome the negative bias at pin 4 of IC1. When this occurs, pin 4 will become positive with respect to pin 3, and IC1 output will go low. This in turn will cause AND gate output (pin 3 of IC3) to go low.
The AND gate output (pin 3 of IC3) is thus high only when Vin is between the range between the small negative voltage and the small positive voltage. At all other times the AND gate output is low. The "window" of Vin voltages that produces a high output is determined by the bias voltages set at pin 4 of IC1 by R9, and at pin 3 of IC2 by R7. It should be noted that when Vin decreases, the same result occurs. That is, when Vin goes below the window upper level (i.e. to some small positive voltage), pin 4 of IC1 will become more negative than pin 3, giving a high output at IC1, and at IC3. But when Vin drops below the window lower level (i.e. beyond a small negative voltage), pin 3 of IC2 will go negative, producing a low IC2 and IC3 output. The output of AND gate IC3 is thus positive only when Vin transitions through the window of voltages set by R7 and R9.
The output of IC3 is voltage divided by Rll and potentiometer R12, and filtered by capacitor C2. Output voltage, Vout is then applied to a digital multimeter (DMM) 14. DMM 14 effectively integrates the filtered output pulses from pin 3 of IC3. As explained above, IC3 will produce pulse lengths equal to the amount of time that Vin is within the preset window voltage range. DMM 14, therefore, reads a voltage that is proportional to the amount of time Vin is within this window. Put another way, the DMM indicates the relative amount of time that an input signal is at a given input level range.
By setting R5 to make Vx some positive or negative value, the window will simply be raised or lowered. The window voltage range (i.e. the difference between the window lower and upper voltage levels) remains constant. Therefore, R5 can be set for different values of Vx, and DMM 14 will indicate the relative amount of time that the signal spends at these levels. In practice, the window lower voltage, VL, and higher voltage, VH, are set symmetrically about Vx. The digital multimeter therefore indicates the relative amount of time that Vin equals Vx, where Vx is approximated by a window of voltages between VL and VH.
This relationship is exactly what is given by theoretical probability density functions. That is, for Gaussian noise, the amplitude probability density function is: ##EQU1## where v is the instantaneous voltage amplitude and s is the root mean square (rms) voltage level. This relationship is plotted on FIG. 3 for a value of s=140 mv rms, see plot 16. The 140 mv rms value of s is a relatively arbitrary level, chosen in this case to be small enough for noise sources that were actually tested, but high enough to be relatively insensitive to circuit fluctuations. The right hand ordinate of FIG. 3 is labeled as "Vout ×10" since in this representative embodiment, the invention is normally calibrated by R12 to give a DMM reading of 0.286 VDC with the 140 mv rms noise voltage in, with Vx set to 0 VDC. This makes measurements easily related to the theoretical curve 16 (i.e. DMM reading times 10), while not requiring an additional amplifier with a gain of 10.
In practice the invention is easily calibrated by applying a sinewave at a frequency of about 4 kHz to the Vin input, and setting Vx equal to 0 VDC. The input signal, Vin, and the output pulses at pin 3 of IC3, VW, are simultaneously viewed on an instrument such as a dual trace oscilloscope. Potentiometer R7 is then set so that the VW pulses rise when Vin is at a -VL level, and R9 is set so that the trailing edge of VW occurs when Vin is at a VH level. This defines the sampling window as VH -VL, with Vx =0 VDC as the average value. Resistor R12 is then adjusted, as noted above, for 0.286 VDC output with 140 mv rms of noise in. The area under the P(v) curve of FIG. 3 between any two values of Vx represents the relative probability (or amount of time) that a signal voltage is within the windowed voltage range. Therefore, for greatest accuracy, the measurement window range should be as small as possible. However, very small window ranges can produce very low output voltages as well as high sensitivities to circuit imperfections. For practical purposes, therefore, with Vx =0 V, VL can be set for -25 mv and VH can be set for +25 mv, giving a 50 mv window range. This value has given adequate accuracy for a number of commercial noise sources tested.
Once the invention has been calibrated as described above, measurements can be done as follows. First, the noise source to be tested, signal source 12, is connected to analyzer 10's Vin input. The output of the noise source is adjusted for an amplitude of 140 mv rms, as read on a true rms voltmeter. The value of Vx is then varied over the range, for example, -0.5 to +0.5 VDC, and Vout ×10 is recorded for each value of Vx tested. By closely spacing the values of Vx continuous measurement resolution is possible. These measurements are plotted on a graph such as FIG. 3 so that measured curves can be directly compared to theoretical curve 16 for P(v). FIG. 3 shows examples of such test data. As can be seen, the measured values for a good noise source (shown as "X's) fall very closely to theoretical curve 16. However, the curve for a noise source known to be faulty (plotted with dots) shows a marked deviation from Gaussian theory. This curve (18), indicates that the noise signal contains inordinately high negative peaks, positive peak clipping, and a positively skewed, leptokurtic distribution. All of these non-Gaussian characteristics could be cross-correlated with an oscilloscope display of the signal. It should be noted that none of these characteristics would be readily visible from a cumulative amplitude distribution of the faulty noise source.
As earlier discussed, the signal amplitude distribution analyzer of the invention can also be used to measure cumulative amplitude distribution. This is done by setting switch S1 to the CUM position, thereby disconnecting the output from IC1 and connecting the pin 1 AND input of IC3 to a high level of +5 VDC. This in effect makes IC3 just a buffer stage whose output is the same as the input of its pin 2. The circuit is then recalibrated by first terminating the Vin input in 50 ohms (no signal), setting Vx to some negative value, and setting R12 for a 1 VDC DMM reading. Since IC2 and IC3 are both continuously high under these conditions, simulating the case where Vin is always greater than 0 V, the 1 VDC DMM reading is set to indicate a 100% cumulative distribution value. Next, with a noise source having an amplitude of 140 mv rms connected to the Vin input, and with Vx set for 0 VDC, R7 is set for a 0.5 VDC reading on the DMM. In this case, R7 acts only as an offset nulling adjustment for IC2, and the DMM reading indicates that the signal is above 0 VDC 50% of the time, as it should with Vx =0 VDC.
Once calibrated, cumulative distribution measurement is made in the same way as for other distribution measurements. The value of Vx is varied over a selected range such as 0 to +0.5 VDC and the cumulative distribution is read on the DMM for a number of Vx values. The recorded values of Vx and Vout (the left hand ordinate of FIG. 3) are then co-plotted with the theoretical Gaussian cumulative distribution curve 20 as shown on the left of FIG. 3. In this case theoretical curve 20 is defined by P(cum v), which is the integral of the distribution function from Vx to infinity. The P(cum v) theoretical curve 20 can also be computed using available tables and the formula: ##EQU2## where "erfc" is the complementary error function.
FIG. 3, in addition to showing distribution density function profiles, also shows cumulative distribution data obtained from the same noise sources tested for the distribution density function. Note that the good noise source (plotted as "X"s) produces test data points that lie closely to theoretical curve 20. The faulty noise source data (shown as dots), on the other hand, departs from theoretical curve 20. However, unlike the data from distribution density function testing, little indication is given as to the degree or character of this departure from Gaussianicity.
FIG. 3 provides a direct comparison of results from cumulative and distribution density function testing. It shows how much more revealing it is to measure the distribution density function itself, rather than to measure cumulative distribution, as is typically done. FIG. 3 also shows the measured amplitude distribution of a sinewave (22) for positive values of Vx. This demonstrates that the signal amplitude distribution analyzer of the invention can be used to make amplitude distribution measurements of any waveform type whether they be repetitive or non-repetitive.
The advantages of the signal amplitude distribution analyzer of the invention over older methods of measuring signal amplitude distribution are many. The invention gives the capability of measuring probability density function itself rather than cumulative distribution functions, providing a more sensitive measure as to the degree of conformance of a noise source with theory. The invention is therefore capable of characterizing deviations from theory rather than only indicating a degree of nonconformance.
As the invention is easily calibrated from a known input source, precise measurement of a signal source's deviations is possible. These sources may be of a repetitive or non-repetitive type, with the invention providing a continuous, overlapping, resolution to enhance measurement accuracy. The invention provides not only a mechanism to sense the presence of sought-after signals, but also provides the relative amount of time that an input signal is within a selected signal range. This latter relationship is exactly what is given by theoretical probability density functions.
It should be noted that the signal amplitude distribution analyzer of the invention uses a very simple implementation scheme, requiring only three integrated circuit components. 0f course the invention could be made more accurate by using more complex circuitry. It also could be improved by using an analog-to-digital converter, computer interface circuits and computer analysis. The invention could also be expanded to provide automatic threshold or window stepping, so that results could be printed automatically on an X/Y recorder.
Obviously, many modifications and variations of the invention are possible in light of the above teachings. It is therefore to be understood that within the scope of the appended claims the invention may be practiced otherwise than as has been specifically described.
Claims (10)
1. An apparatus for analyzing signals comprising:
a first comparator for comparing voltage of said signals to a preselected lower voltage and for providing an output signal when the voltage of said signals is above said preselected lower voltage;
a second comparator for comparing voltage of said signals to a preselected higher voltage and for providing an output signal when the voltage of said signals is below said preselected higher voltage; and
an AND circuit coupled to said first and second comparators for providing a pulse upon simultaneously receiving said output signals from said comparators, said pulse having a pulse length corresponding to the elapsed time said signals are within a voltage window bounded by said lower and higher voltages.
2. An apparatus according to claim 1 further including:
an integrator coupled to said AND circuit for integrating said pulse provided by said AND circuit to provide a voltage proportional to said elapsed time said signals are within said voltage window.
3. An apparatus according to claim 2 in which said integrator is a voltmeter.
4. An apparatus according to claim 1 in which said apparatus is calibrated by analyzing a known input signal.
5. An apparatus according to claim 4 in which said known input signal is a Gaussian signal.
6. A method for approximating the amplitude probability density distribution of signals comprising the steps of:
comparing voltage of said signals to a preselected lower voltage in a first comparator;
providing an output signal from said first comparator when the voltage of said signals is above said preselected lower voltage;
comparing voltage of said signals to a preselected higher voltage in a second comparator;
providing an output signal from said second comparator when the voltage of said signals is below said preselected higher voltage;
providing a pulse from an AND circuit operably coupled to said first and second comparators upon said AND circuit simultaneously receiving said output signals from said comparators, said pulse having a pulse length corresponding to the elapsed time said signals are within a voltage window bounded by said lower and higher voltages; and
integrating said pulse provided by said AND circuit in an integrator to provide an output voltage proportional to said elapsed time said signals are within said voltage window.
7. A method according to claim 6 further including a step of
analyzing a known input signal as a calibration of said comparators and said integrator.
8. A method according to claim 7 in which said known input signal is a Gaussian signal.
9. A method according to claim 8 further including steps of:
analyzing unknown input signals by shifting said voltage window over a desired voltage range, each voltage window having a middle voltage point; and
plotting said output voltage versus said middle voltage point for each voltage window examined.
10. A method according to claim 9 further including the step of:
comparing said plot to a theroretically derived probability density distribution to determine the magnitude of deviation of said plot from said theoretically derived probability density function.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US08/083,962 USH1458H (en) | 1993-06-23 | 1993-06-23 | Signal amplitude distribution analyzer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US08/083,962 USH1458H (en) | 1993-06-23 | 1993-06-23 | Signal amplitude distribution analyzer |
Publications (1)
Publication Number | Publication Date |
---|---|
USH1458H true USH1458H (en) | 1995-07-04 |
Family
ID=22181785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US08/083,962 Abandoned USH1458H (en) | 1993-06-23 | 1993-06-23 | Signal amplitude distribution analyzer |
Country Status (1)
Country | Link |
---|---|
US (1) | USH1458H (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5682336A (en) * | 1995-02-10 | 1997-10-28 | Harris Corporation | Simulation of noise behavior of non-linear circuit |
US6509728B1 (en) * | 1998-05-28 | 2003-01-21 | Anritsu Corporation | Spectrum analyzer having function of displaying amplitude probability distribution effectively |
US20030171883A1 (en) * | 2002-03-05 | 2003-09-11 | I-Shih Tseng | Method and system for measuring characteristics of liquid crystal display driver chips |
US8664978B2 (en) * | 2011-08-30 | 2014-03-04 | Taiwan Semiconductor Manufacturing Company, Ltd. | Methods and apparatus for time to current conversion |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2905888A (en) * | 1957-04-08 | 1959-09-22 | Collins Radio Co | Amplitude distribution analyzer |
US3497684A (en) * | 1967-06-30 | 1970-02-24 | Bell Telephone Labor Inc | Direct reading statistical distributor estimator |
US3581200A (en) * | 1969-02-18 | 1971-05-25 | Ralph W Mifflin | Apparatus for deriving the probability density function of a signal with respect to amplitude |
US3626168A (en) * | 1969-07-07 | 1971-12-07 | Boeing Co | Measuring correlation, amplitude probability and probability density distributions, and signal response averaging |
US3665505A (en) * | 1969-04-24 | 1972-05-23 | Ver Flugtechnische Werke | Apparatus and method for measuring and analyzing dynamic processes |
US3742352A (en) * | 1970-12-02 | 1973-06-26 | H Euer | Apparatus and method for measuring and analyzing dynamic processes |
US4129827A (en) * | 1976-05-20 | 1978-12-12 | Southwick R | Amplitude probability detector |
US4342965A (en) * | 1980-09-22 | 1982-08-03 | Baldwin Johnny E | Waveform instability detector |
US4528591A (en) * | 1982-08-16 | 1985-07-09 | Tektronix, Inc. | Raster scan frame digitizer method and apparatus |
US4625283A (en) * | 1982-05-07 | 1986-11-25 | Cooper Industries, Inc. | Method and apparatus for digitally measuring alternating current |
US4985844A (en) * | 1989-05-08 | 1991-01-15 | Tektronix, Inc. | Statistical waveform profiler employing counter/timer |
US5003248A (en) * | 1986-11-10 | 1991-03-26 | Hewlett-Packard Company | Probability density histogram display |
-
1993
- 1993-06-23 US US08/083,962 patent/USH1458H/en not_active Abandoned
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2905888A (en) * | 1957-04-08 | 1959-09-22 | Collins Radio Co | Amplitude distribution analyzer |
US3497684A (en) * | 1967-06-30 | 1970-02-24 | Bell Telephone Labor Inc | Direct reading statistical distributor estimator |
US3581200A (en) * | 1969-02-18 | 1971-05-25 | Ralph W Mifflin | Apparatus for deriving the probability density function of a signal with respect to amplitude |
US3665505A (en) * | 1969-04-24 | 1972-05-23 | Ver Flugtechnische Werke | Apparatus and method for measuring and analyzing dynamic processes |
US3626168A (en) * | 1969-07-07 | 1971-12-07 | Boeing Co | Measuring correlation, amplitude probability and probability density distributions, and signal response averaging |
US3742352A (en) * | 1970-12-02 | 1973-06-26 | H Euer | Apparatus and method for measuring and analyzing dynamic processes |
US4129827A (en) * | 1976-05-20 | 1978-12-12 | Southwick R | Amplitude probability detector |
US4342965A (en) * | 1980-09-22 | 1982-08-03 | Baldwin Johnny E | Waveform instability detector |
US4625283A (en) * | 1982-05-07 | 1986-11-25 | Cooper Industries, Inc. | Method and apparatus for digitally measuring alternating current |
US4528591A (en) * | 1982-08-16 | 1985-07-09 | Tektronix, Inc. | Raster scan frame digitizer method and apparatus |
US5003248A (en) * | 1986-11-10 | 1991-03-26 | Hewlett-Packard Company | Probability density histogram display |
US4985844A (en) * | 1989-05-08 | 1991-01-15 | Tektronix, Inc. | Statistical waveform profiler employing counter/timer |
Non-Patent Citations (2)
Title |
---|
S. R. Michaels, "Test Whether a Noise Source is Gaussian", EDN, 11 Dec. 1, pp. 272-274. |
S. R. Michaels, Test Whether a Noise Source is Gaussian , EDN, 11 Dec. 1986, pp. 272 274. * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5682336A (en) * | 1995-02-10 | 1997-10-28 | Harris Corporation | Simulation of noise behavior of non-linear circuit |
US6509728B1 (en) * | 1998-05-28 | 2003-01-21 | Anritsu Corporation | Spectrum analyzer having function of displaying amplitude probability distribution effectively |
US20030171883A1 (en) * | 2002-03-05 | 2003-09-11 | I-Shih Tseng | Method and system for measuring characteristics of liquid crystal display driver chips |
US6925415B2 (en) * | 2002-03-05 | 2005-08-02 | Chroma Ate Inc. | Method and system for measuring characteristics of liquid crystal display driver chips |
US8664978B2 (en) * | 2011-08-30 | 2014-03-04 | Taiwan Semiconductor Manufacturing Company, Ltd. | Methods and apparatus for time to current conversion |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0242700A2 (en) | AC level calibration method and apparatus | |
US4608657A (en) | Method and apparatus for testing probe calibration | |
US5844412A (en) | Board test apparatus and method for fast capacitance measurement | |
USH1458H (en) | Signal amplitude distribution analyzer | |
Deyst et al. | A fast-pulse oscilloscope calibration system | |
US4517511A (en) | Current probe signal processing circuit employing sample and hold technique to locate circuit faults | |
CA2536083C (en) | Method and apparatus for measuring impedance of electrical component under high interference conditions | |
US6348816B1 (en) | Tracking percent overload signal as indicator of output signal magnitude | |
US20020136337A1 (en) | Method and apparatus for high-resolution jitter measurement | |
KR940002724B1 (en) | Ac evaluation equipment and the mehtod for an ic tester | |
US5194818A (en) | Risetime and falltime test system and method | |
US4719408A (en) | Apparatus for indicating proper compensation of an adjustable frequency compensation network | |
US3475683A (en) | Method and apparatus for measuring signal to noise ratio | |
Mihai | Statistical tools for analysis of the performance of a high resistance measurement bridge | |
Souders et al. | Accurate RF voltage measurements using a sampling voltage tracker | |
EP0077725A1 (en) | Current probe signal processing circuit | |
US3324390A (en) | Apparatus for determining the voltage ratio of a voltage divider that is energized by pulsed signals | |
Procházka et al. | Validation of partial discharge calibrators for small charge rates | |
Hudlicka et al. | Practical aspects of a pulse generator calibration | |
US3710249A (en) | Slideback peak circuits with constant tone indications | |
Kuenzi et al. | Fundamentals of oscilloscope measurements in Automated Test Equipment (ATE) | |
JPH0466878A (en) | Apparatus and method for measuring electrostatic capacitance, resistance and inductance | |
Schoenwetter | Settling time measurements | |
Wright | Calibration of IEC 1000-3-2 harmonic analysers. | |
Claudi | Calibration of digital impulse measuring systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: UNITED STATES OF AMERICA, THE, AS REPRESENTED BY T Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SLACK, ROBERT A.;REEL/FRAME:006620/0331 Effective date: 19930622 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |