US7532992B2  Measuring apparatuses and methods of using them  Google Patents
Measuring apparatuses and methods of using them Download PDFInfo
 Publication number
 US7532992B2 US7532992B2 US11337011 US33701106A US7532992B2 US 7532992 B2 US7532992 B2 US 7532992B2 US 11337011 US11337011 US 11337011 US 33701106 A US33701106 A US 33701106A US 7532992 B2 US7532992 B2 US 7532992B2
 Authority
 US
 Grant status
 Grant
 Patent type
 Prior art keywords
 polynomial
 terms
 correction
 stage
 data
 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.)
 Expired  Fee Related, expires
Links
Images
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
 G01F23/00—Indicating or measuring liquid level, or level of fluent solid material, e.g. indicating in terms of volume, indicating by means of an alarm
 G01F23/14—Indicating or measuring liquid level, or level of fluent solid material, e.g. indicating in terms of volume, indicating by means of an alarm by measurement of pressure
 G01F23/18—Indicating, recording or alarm devices actuated electrically

 G—PHYSICS
 G01—MEASURING; TESTING
 G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
 G01D18/00—Testing or calibrating of apparatus or arrangements provided for in groups G01D1/00  G01D15/00
 G01D18/008—Testing or calibrating of apparatus or arrangements provided for in groups G01D1/00  G01D15/00 with calibration coefficients stored in memory

 G—PHYSICS
 G01—MEASURING; TESTING
 G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
 G01D3/00—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups
 G01D3/02—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups with provision for altering or correcting the law of variation
 G01D3/022—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups with provision for altering or correcting the law of variation having an ideal characteristic, map or correction data stored in a digital memory

 G—PHYSICS
 G01—MEASURING; TESTING
 G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
 G01D3/00—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups
 G01D3/028—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure
 G01D3/036—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure on measuring arrangements themselves
 G01D3/0365—Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure on measuring arrangements themselves the undesired influence being measured using a separate sensor, which produces an influence related signal

 G—PHYSICS
 G01—MEASURING; TESTING
 G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
 G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through the meter in a continuous flow
 G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through the meter in a continuous flow by using mechanical effects
 G01F1/34—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through the meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure
 G01F1/36—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through the meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure the pressure or differential pressure being created by the use of flow constriction
 G01F1/363—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through the meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure the pressure or differential pressure being created by the use of flow constriction with electrical or electromechanical indication

 G—PHYSICS
 G01—MEASURING; TESTING
 G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
 G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through the meter in a continuous flow
 G01F1/05—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through the meter in a continuous flow by using mechanical effects
 G01F1/34—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through the meter in a continuous flow by using mechanical effects by measuring pressure or differential pressure
 G01F1/50—Correcting or compensating means
Abstract
Description
This invention relates to measuring apparatuses and methods of using them such as for example apparatuses for measuring physical parameters, such as pressure, flow rates of liquids, chemical organic vapor concentrations and temperature.
Commonly, measuring instruments are affected by noise, nonlinearity and interference. Noise is related to random error—differences between the measured physical parameter value and the actual physical parameter value that cannot be corrected by additional information. Nonlinearity is a systematic error that arises from assuming a linear response of the instrument. Many transducers have an approximate mathematical linear relationship between its electrical signal voltage output and the measured physical parameter value. This linear relationship can greatly simplify the data processing of the instrument. There are trade offs between ease of data processing and precision in the design of such instruments. If accuracy and precision are more important than ease of data processing for a given application, then these nonlinear errors must be considered and treated. Additional information concerning these nonlinearities can be used to correct the measured physical parameter value. Interference is another systematic error that arises from the influences of other physical parameters on the measured signal. Information on these other physical parameters can be used to correct the original measured physical parameter value. This patent is mainly concerned with treatment of these last two error contributions: nonlinearity and interference.
For example, in one type of measuring apparatus, the depth of water is measured by sensing the pressure at the bottom of the water with a pressure sensor. The pressure sensors in these instruments are transducers that convert pressure to voltage. This pressure sensor has a voltage signal output which can be converted to a pressure signal value which can be further translated into a depth measurement by converting the pressure into units of depth. In some instruments, the depth measurement is used to determine the velocity of flow as in U.S. Pat. No. 5,275,042 or average velocity as in U.S. Pat. No. 5,371,686. Due to natural electrical “white noise” which is not predictable from additional information, the measured pressure parameter value has a noise contribution to the measurement error. These sensors are further assumed to have a linear mathematical relationship between the actual pressure and signal voltage output. Since this relationship is only approximate, the measured pressure parameter value has a nonlinear contribution to the measurement error. These sensors are also affected by temperature, another physical parameter. The measured pressure parameter value has an interference contribution to the measurement error.
In another type of sensor, chemical organic vapor concentrations are sensed by a tin oxide bead that varies its resistance to current as a function of organic vapor level changes due to competition between the organic vapor and oxygen in the air. The organic vapors reduce the tin oxide to metallic tin; whereas, the oxygen oxidizes the tin back to tin oxide. Tin and tin oxide have different electrical resistances. In this type of transducer, the measurement is affected by temperature and humidity. Temperature and water vapor in the air influence the resistance of the tin oxide bead. This effect is used in several other types of instruments such as for example in an analyzer of water for organic impurities as described in U.S. Pat. No. 6,123,904. The above two examples are provided for illustration since there are many different sensors in many different types of apparatuses that are affected by noise, nonlinearity and interference. The accuracy and precision of these sensors is reduced because the output signal voltage is also affected by other physical conditions such as temperature or humidity.
It is known to improve the precision of measurements by using higherorder, multivariant polynomial calibration curves to correct the measurements for nonlinearity and interference. It is also known to obtain the optimum coefficients of terms in the polynomial calibration curve by any of several methods including the least squares regression method. Commonly, the calibration curve is applied to measurements through a microcontroller.
In the prior art use of polynomial calibration curves, the general form of the polynomial such as the number of terms and the degree of the terms must be selected before the coefficients can be determined. Although the general form of the polynomial greatly influences the precision obtained from the use of the calibration polynomial, no completely satisfactory automatic approach for some calibration needs is known.
Accordingly, it is an object of the invention to provide a novel measuring system and method of using it.
It is a further object of the invention to provide a novel technique for calibrating measuring systems.
It is a still further object of the invention to provide a novel calibration system.
It is a still further object of the invention to provide a systematic technique for providing a calibration polynomial.
In accordance with the above and further objects of the invention, a measuring instrument for measuring at least a first value includes at least one sensor or transducer for generating a first signal representing the first value or parameter and at least one other sensor for generating a second signal representing a value or parameter that may interfere with the measurement of the first value, create error or otherwise cause the first value to be inaccurate or imprecise. There may be one or a plurality of such sensors that sense values to be measured and one or a plurality of sensors or transducers to measure values that may interfere, create error or otherwise cause the value or values being measured to be inaccurate or imprecise. The sensor or sensors for sensing the value or values being measured will hereinafter be referred to as the primary sensor or transducer or primary sensors or transducers and the value or values being measured will hereinafter be referred to as the primary value or values. The sensor or sensors measuring the value or values that may interfere, create error or otherwise cause the value being measured to be inaccurate hereinafter is referred to as a secondary sensor or sensors or a secondary transducer or transducers. A signal provided by the primary sensor or transducer is referred to in this specification as a primary measurement signal. The signal resulting from correction of a primary measurement signal in accordance with the methods described herein is referred to hereinafter as a calibrated signal.
In this specification, the word “value” and the word “parameter” each mean any physical characteristic that may be sensed such as for example the depth of water or pressure or temperature or intensity of energy. The word “value” herein includes not only values that are to be used or understood by a person but also values of interference that may interfere with values that are to be used. For example, the value to be measured may be pressure as measured by a pressure sensor and an interference related value may be temperature which alters the signal provided by the pressure sensor.
The measuring instrument may include a microcontroller for correcting the primary measurement signal to provide the calibrated signal. Hereinafter, the values that may interfere with this signal are referred to as interference related values and the signals generated by a secondary sensor or transducer from the interference related values are referred to as interference related signals. The microcontroller or other apparatus for correcting the signals includes a correction polynomial or correction curve hereinafter referred to as a criteriaoptimized correction polynomial or curve or as a criteriaoptimized calibration curve or polynomial. In this specification, the words “criteriaoptimized correction polynomial” or “criteriaoptimized correction curve” or “criteriaoptimized calibration curve” or “criteriaoptimized calibration polynomial” means a relationship formed using a special procedure described hereinafter The criteriaoptimized correction polynomial or criteriaoptimized correction curve includes the interference related signals and the uncorrected value of interest as independent variables and the corrected value to be measured as indicated by the calibrated signal as a dependent variable.
The microcontroller is connected to receive the primary measurement signal and one or more interference related signals. This enables the microcontroller to correct the primary measurement signal for the interference to provide the calibrated signal representing the corrected value of interest. In one embodiment, the measuring instrument is an apparatus for determining the volumetric rate of flow of a liquid within a flow bed. For this purpose, it includes apparatus for determining the average rate of flow of the liquid in the flow bed and the primary sensor is a pressure sensor positioned at the bottom of a flow path for the liquid. Thus the primary measurement signal is related to the depth of the liquid. An interference related sensor is a temperature measuring sensor in this embodiment.
In this embodiment, the criteriaoptimized correction polynomial has the signals from the pressure sensor and temperature sensor as independent variables and depth as a dependent variable. The volumetric flow rate is calculated by the microcontroller by multiplying the depth by the average flow velocity of the water. In the preferred embodiment, the average flow velocity is determined by transmitting an ultrasonic signal into the liquid in the flow bed; receiving reflected ultrasonic signals; and utilizing the signals to calculate an approximate average velocity. The approximate average velocity is calculated by performing a Fourier transform on the digital signals and averaging certain of the coefficients of the Fourier transform as described in U.S. Pat. No. 5,371,686; the disclosure of which is incorporated herein by reference.
In one embodiment, the procedure for forming the criteriaoptimized correction polynomial or criteriaoptimized correction curve includes the steps of forming a preparatory calibration polynomial with a plurality of independent variables. One of the independent variables referred to hereinafter as a primary independent variable is the primary measurement signal. The other independent variables are the interference related signals. The dependent variable refers to the sought after measured value such as for example the primary signal corrected for interference and nonlinearity from a pressure sensor. This signal is referred to as the calibrated signal in this specification although the calibrated signal may be further modified such as by being converted to a digital signal from an analog signal or by being operated upon by other components of the circuit to make it suitable for other operations significant to the operation of the instrument. If this pressure sensor is used to indicate the depth of a liquid in which the pressure sensor is immersed, then the dependent variable is the depth of the liquid since the pressure multiplied by a constant reflecting the water density in the units being utilized provides a signal indicating the depth of the liquid.
The preparatory calibration polynomial includes as its terms the independent variables raised to a degree that is selected by criteria relating to the accuracy and precision desired or needed in the dependent variable for the application of the measuring instrument. The accuracy and precision may be selected because it is important to the usefulness of the instrument or the needs of the final procedure using the instrument. For example, a surgical instrument may require a certain procedure or not be usable for safety reasons. In this case, the safety of the patient is the criteria determining the precision of the instrument. On the other hand, the precision may be selected in view of the equipment that is to be used. The higher the degree of the independent variables and the number of cross terms, the larger the memory of the microcontroller or other apparatus or techniques used. Accordingly, if the size of the memory is limited then the degree of the independent variables will be limited even though that will have an adverse effect on the accuracy of the final measurement. In the preferred embodiment, the criteria can be the size of the memory or speed of a microcontroller.
Once the highest degree of the independent variables has been selected, a preparatory or working polynomial is formed including a plurality of terms. Each term has only one of the independent variables in it but the preparatory correction polynomial includes lowerdegree variables. In the preferred embodiment, the preparatory correction polynomial prepared at this point of development includes terms having at least some of the lower degree variables. In the preferred embodiment, one or all of the independent variables of all degrees from the highest selected degree to the first degree are included. For example, if the fourth degree were to be selected for temperature then a variable, C_{4}T^{4}, would have a coefficient, C_{4}, and all of the lower degree of temperature such as T^{3}, T^{2 }and T have corresponding coefficients C_{3}, C_{2 }and C_{1}.
At this point in the development of the correction polynomial, no terms have multiple independent variables as part of them. Terms that include the product of two or more independent variables are hereinafter referred to as cross terms. In this specification, the degree of the variable refers to the power to which the variable is raised so that for example, the degree of T^{4 }is 4.
After the dependent variable, the independent variables and their degree have been chosen, some cross terms are added to the preparatory correction polynomial and the coefficients are evaluated using any known method such as the well known method referred to as the least squares regression method. The coefficients are determined using data obtained from a test arrangement that includes primary and secondary sensors. The cross terms that are added at this point preferably include only lower degree independent variables. Any arbitrary number of cross terms may be selected including none but it is desirable to select no more than half of the possible cross terms at this point. After the general form of the variables has been obtained including the dependent, independent variables and the cross terms with the selected degree of the variables, the coefficients are evaluated.
The coefficients are evaluated using data obtained from a test arrangement. An instrument designed to use the correction polynomial to obtain calibrated signals may serve as the test arrangement or a special test arrangement may be used. In either case, a primary and at least one secondary transducer are used to sense the primary value and the interference related value or values as the primary value is varied in a measurable manner to permit the calibrated signal to be determined. Each set of data comprising the primary value, interference related value and corrected value is obtained and collectively they are compiled into a data base used in testing the terms of the preparatory correction polynomials. For example, a pressure transducer may be the primary transducer and a thermistor may be a secondary transducer and the depth of the primary transducer may be the dependent variable. Data may be gathered by locating the transducers at a known depth in a body of water and the coefficients of a correction polynomial for depth may be found using any suitable method such as the least squares regression method. The preparatory correction polynomial at this point in its development is referred to in this specification as a firststage preparatory correction polynomial.
After the firststage preparatory correction polynomial has been formed, secondstage and thirdstage preparatory correction polynomials are formed leading to a final step of arriving at the correction polynomial that is to be used to calculate the calibrated signal. The secondstage preparatory correction polynomials result from the elimination of terms that provide the least improvement to the precision of the measurement and the thirdstage correction polynomials result from the addition of cross terms that provide the greatest improvement on the precision of the measurement.
In the preferred embodiment, the second and thirdstage preparatory correction polynomials are alternately formed until secondstage and thirdstage preparatory correction polynomials are substantially the same. At this point, a term that was eliminated from a secondstage preparatory correction polynomial is added to form a thirdstage preparatory polynomial or a term that was added to form a thirdstage preparatory correction polynomial is removed to form a secondstage correction polynomial. At this point in the development of the correction polynomial, the polynomial is said to be selfconsistent. In the preferred embodiment, this selfconsistent correction polynomial is used in the measuring instrument to form the calibrated signal from the values sensed by the instrument.
To form a secondstage preparatory correction polynomial, each term of the firststage preparatory polynomial is tested and the term that improves the precision of the measurement the least using the data from the data base, a plurality of data sets is eliminated from the secondstage correction polynomial. A series of secondstage correction polynomials are formed by eliminating terms that improve the precision of the measurement the least. In the preferred embodiment, terms are eliminated until the selected criteria are met. For example, terms may be eliminated until the polynomial is usable in the size of the memory chosen for the microcontroller or until the thirdstage preparatory correction polynomial formed by eliminating terms that improve the precision of the instrument the least and adding terms that improve the precision the most becomes selfconsistent. In the preferred embodiment, the series of secondstage correction polynomials are not consecutive but alternates with the formation of thirdstage correction polynomials.
To determine which term improves the precision of the measurement the least when forming a secondstage preparatory correction polynomial, a series of trial correction polynomials are formed. A trial polynomial is formed by eliminating a term from the firststage correction polynomial, determining the coefficients with the term removed to obtain a trial correction polynomial and determining that the trial correction polynomial improved the precision of the measurement less than other trial correction polynomials being considered. These trial polynomials are compared to determine which trial polynomial provides the most precise and accurate determination of the dependent variable. That trial correction polynomial becomes a secondstage preparatory correction polynomial which may result in a final secondstage preparatory correction polynomial except for one or more possible thirdstage correction polynomials formed by adding terms if the criteria is met or if the correction polynomial becomes self consistent.
In the preferred embodiment, each of the secondstage preparatory correction polynomials is selected from a comparison of every trial secondstage polynomial formed from elimination of every cross term from one other first or thirdstage preparatory polynomial. However, it is possible to obtain good results by eliminating terms that include only one independent variable in addition to cross terms and to form and consider trial correction polynomials formed by eliminating terms to other secondstage correction polynomials and/or to form fewer trial correction polynomials for comparison by not forming trial polynomials by elimination of each cross term but only select some of the cross terms to be eliminated to form a trial polynomial.
In the preferred embodiment, secondstage trial correction polynomials are compared by calculating the dependent variable several times for each trial correction polynomial using a plurality of different data sets. Deviation between the calculated value and the actual value is determined for each trial polynomial at each of the plurality of data sets and a comparison made using these results. In the preferred embodiment, several measurements are made at each preset value for the dependent variable. For example, in the case in which pressure and temperature are the independent variables and the depth under water is the dependent variable, several measurements of pressures and temperatures are made with the test apparatus at known depths in generating the data base. Each trial correction polynomial is used to calculate depth from the values of pressure and temperature in the data base and the calculated values of depth are compared to the measured values from the data base at a plurality of data sets from the data base. A data set in this specification is a set of measurements at different depths, temperatures and pressures. Several sets are obtained for each correction polynomial and each trial correction polynomial is tested with several sets.
In the preferred embodiment, the maximum deviation from the plurality of data sets for each trial correction polynomial is taken as the standard to be used in determining which term will be discarded although any other value directly related to the disagreement error and reflecting the effect of the term on improving the precision or the measurement such as a value within ten percent of the maximum deviation could be used. The average of the maximum deviations for each trial correction polynomial is calculated and this is used to determine the trial correction polynomial that causes the lowest average maximum deviation. The secondstage trial correction polynomial that causes the lowest average maximum deviation is taken as a secondstage correction polynomial. Although the average maximum deviation is used as the comparison standard in the preferred embodiment, any other indicator of the central tendency of the deviations such as for example the mean deviation could be used for this selection provided it results in a correction polynomial that provides a more precise and accurate measurement.
To increase the precision of measurement, new cross terms are added to the secondstage preparatory correction polynomials one by one at a time to form trial thirdstage preparatory correction polynomials. Each trial thirdstage preparatory correction polynomial includes all of the terms of a secondstage preparatory polynomial with one added cross term. The trial thirdstage correction polynomials are each evaluated and the cross term that improves the precision of the measurement the most are left in the final thirdstage preparatory polynomial. Each time a cross term is added, the coefficients of all of the terms are optimized and the disagreement errors are found. When all of the disagreement errors have been found, a central tendency value such as the average maximum deviation is determined for the disagreement errors. The cross term that provides the lowest average deviation in a trial thirdstage correction polynomial is added permanently. The polynomials obtained with the above procedure may be used to calibrate a measuring instrument. Thus, measuring instruments may be economically made and nevertheless provide high precision.
From the above summary it can be understood that the instrument of this invention and the method of using it has several advantages, such as for example: (1) it can be more easily calibrated; (2) some aspects of it are more easily automated; and (3) it can provide more precise operation.
The above noted and other features of the invention will be better understood from the following detailed description when considered in connection with the accompanying drawings, in which:
In
The primary transducer 18 is intended to generate a signal representing a value of interest in some applications. The secondary transducer 16 is an interference related value. It is used to correct the reading from the primary transducer 18 to correct for interference related influences. The value that is sensed by the primary transducer 18 is the primary value. The electrical signal generated by the primary transducer 18 is referred to as the primary measurement signal and after it is corrected for nonlinearity and interference by the correction polynomial, it is referred to as the calibrated signal. Because changes in environment affect the primary measured signal, these changes are interference values. The primary values and interference values may be considered points of stimulation. The interference related signals from the secondary transducer 16 are referred to as secondary because they are not used directly for the purpose of obtaining measurement signals but are only used for removing error from the reading of the primary transducer 18, which is used for these purposes of obtaining measurement signals.
For example in the preferred embodiment, the system 10 generates a calibration curve for depth of a liquid to be used in a measuring instrument that measures depth and average velocity through a known flow bed and determines volumetric flow rate. In this instrument and in the system 10, the primary transducer 18 is a pressure sensor that responds to pressure by generating an electrical voltage. The primary transducer 18 is positioned at the bottom of the flow bed in the instrument and is positioned at different known depths of liquid 22 in a container 20 in the system 10. The secondary transducer 16 is a thermistor that generates voltage in response to temperature of the liquid 22. Temperature is interference in the instrument since it changes the voltage response to pressure of the primary transducer 18.
In the preferred embodiment, the system 10 generates a plurality of sets of data correlating depth of liquid, temperature and voltage output of the primary transducer 18. The depth measurement is a primary value that results in a primary measurement signal. Each set of data includes a plurality of readings of voltage output from the primary transducer 18, which are primary measurement signals and temperature of the liquid 22 which is an interference related value at different depths of the primary transducer 18. For example, the liquid 22 can be brought to a particular temperature by a temperature control device 36 and the depth of the primary transducer 18 in the liquid 22 can be changed by removing liquid from the container 20 or the liquid above the primary transducer 18 can be changed by changing the depth of the primary transducer 18. The voltage output from the primary transducer 18 can be read and recorded in the calibration microcontroller 14 at each depth of the primary transducer 18 and then the liquid 22 brought to a different temperature and the depth changed again to take readings to generate and record a set of depth, temperature voltage data at different temperatures and depths. This process can be repeated to generate and record several sets of data.
The data sets are obtained by obtaining a series of output values from one of a first of a primary transducer 18 and a secondary transducer 16 while the one of the first of a primary transducer 18 and secondary transducer 16 is at a first stimuli point and the other of the primary transducer and secondary transducer is subject to said series of other stimulating points. A series of output values, each at a different stimuli value of a series of stimuli values are obtained from the primary and secondary transducers. The primary independent variable is obtained from the primary transducer 18. In this specification, the words “stimulating point” mean a value of the independent variable. For example, in forming a data set in which water level is the dependant variable, the different levels at which the pressure sensor 18 is located are each stimulating points and the different temperatures for measurement are stimulating points.
The data sets are used to eliminate the least significant terms of the preparatory calibration polynomial by eliminating one term at a time of the preparatory polynomial. The coefficients of the remaining terms are determined and the maximum deviation determined until the maximum deviation has been determined for several of the terms. The terms that increase the maximum deviation the most are permanently eliminated. After the least significant terms have been eliminated from the calibration polynomial, cross terms are added one at a time, the coefficients for the polynomial each time a cross term is added are determined and the maximum deviation is determined. The cross terms having the most effect on the maximum deviation are permanently added.
The primary value is the value sensed by the primary transducer 18. In this example, the primary value is depth of the liquid. The electrical signal generated by the primary transducer 18 is referred to as the primary measurement signal and after it is corrected for nonlinearity and interference by the correction polynomial it is referred to as the calibrated signal. Because changes in temperature affect the pressure readings, temperature changes are interference values. The secondary transducer 16 measures the temperature and the signals from the thermistor 16 are interference related signals. Signals from the secondary transducer 16 are referred to as secondary because they are not used directly for the purpose of obtaining depth of the water nor average velocity of the water in the flow stream but for removing interference from the reading of the primary transducer 18, which is used for these purposes.
The calibration microcontroller 14 obtains readings of output voltage from the pressure sensor 18 at a plurality of data points which, in the preferred embodiment, are levels of water or other pressure inducing means. It also obtains readings of temperatures at which the readings of the level of water are obtained. It also receives or generates a general form of polynomial to serve as the framework of a model. While a calibration microcontroller is referred to in
In the preferred embodiment, the polynomial consists of terms having temperature as an independent variable, terms having pressure as an independent variable and cross terms that are the product of temperature and pressure. The degree of the variables (powers or exponents) is selected in accordance with the precision needed and the microcontroller memory available. Coefficients of the starting polynomial are determined from the data sets using any appropriate technique. In the preferred embodiment, the well known least squares regression method is used as described in “Data Reduction and Error Analysis For The Physical Sciences”, Chapters Eight and Nine, Pgs. 134186. McGrawHill Book Company, New York, 1969. After these coefficients are found, the least significant cross terms are eliminated and the most significant cross terms that are within the design criteria are added in a manner described hereinafter. This polynomial is used to calibrate the measuring instrument.
In
With this arrangement, the inputs to the control and computation system 48 result in a more precise signal that may be read out by the information input/output system 49. In the preferred embodiment, the criteriaoptimized instrument 30 is a depth measuring system in which the first sensor 18 is a pressure sensor. This pressure sensor at the bottom of a flow stream is subject to pressure from the head of liquid in the flow stream and converts it to an electrical signal that is transmitted to the control and computation system 48. In the preferred embodiment, only one interference sensor is used. This sensor is the second sensor 16, which detects temperature. The temperature is used to correct for variations and the signal output from the first sensor 18 through the use of the criteriaoptimized polynomial recorded in the microprocessor within the control and computation system 48.
In
The step 66 of forming third stage preparatory calibration polynomials is performed only if a step 118 of determining whether a term was eliminated that had been previously added. If a term was eliminated that had been previously added, the polynomial is selfconsistent as shown at step 78 and the formation of the calibration polynomial is complete. Similarly, the program 24 returns to the step 56 of eliminating one term at a time from the first or third stage polynomial to produce a set of trial second stage polynomials only if a step 120 of determining that a term was not eliminated that had been previously added. If a term was eliminated that had been previously added, the polynomial is self consistent as shown at step 78 and the calibration curve has been determined.
Since the firststage polynomial formed in step 26 may have many terms, a set of many trial secondstage preparatory polynomials is formed in step 56—one for every term eliminated. For example, if the firststage polynomial includes three terms, then the set of trial secondstage preparatory polynomials includes three polynomials. As a further example, if the firststage polynomial were: y=C_{1}X^{2}+C_{2}X^{2}+C_{3}XT, then the complete set of trial secondstage preparatory polynomials would be: y=C_{2}X^{2}+C_{3}XT, y=C_{1}X+C_{3}XT and y=C_{1}X+C_{2}X^{2}.
To select the most significant cross terms to be added to the preparatory secondstage polynomial, the process includes the step 66 of adding one new cross term from a list generated in step 84 (
More specifically, the precision of measurements in a measuring instrument is increased while staying within a selected criteria by adding new cross terms to form thirdstage preparatory correction polynomials. The cross terms that are added are selected from a set of all possible cross terms that can be formed from the independent variables in the firststage preparatory correction polynomials except cross terms that were initially selected to be part of the firststage preparatory correction polynomial. Terms that were selected to be part of the original firststage preparatory correction polynomial remain in the second and third preparatory correction polynomials unless they are eliminated in forming a secondstage preparatory polynomial as described above.
Each of the cross terms that are considered for addition to the firststage preparatory polynomial is added to a secondstage preparatory correction polynomial as shown in step 66 of
In the preferred embodiment, the cross terms that are added to form trial correction polynomials are: (1) cross terms not in the firststage preparatory correction polynomial; and (2) cross terms that include an independent variable of lower or equal degree to the highest degree of the independent variable in the firststage preparatory correction polynomial. For example, if the highest degree on one of the independent variables is T^{4}, cross terms are added to the extent necessary to provide a cross term that includes one or more of the independent variables T^{4}, T^{3}, T^{2 }and T but not T^{5 }in the preferred embodiment although it is possible under some circumstances to include an independent variable of a higher degree. The selection of lower degrees of the independent variables to form cross terms is done for each independent variable in the preferred embodiment. In the preferred embodiment, only one new cross term is included in the trial polynomial although more than one could be included.
Each time a cross term is added, the coefficients of all the terms are evaluated and the disagreement errors are found as described in connection with the development of secondstage preparatory correction polynomials. When all of the disagreement errors have been found, a central tendency value such as the average maximum deviation is determined for the disagreement errors. The cross term that provides the lowest average deviation in a trial thirdstage correction polynomial is added permanently. Thus, cross terms having independent variables in them lower than the highest order independent variable in a term having only one independent variable and which improves the precision of the measurement are permanently added to the preparatory correction polynomial. Cross terms are added until the preparatory correction polynomials are self consistent or all of the cross terms with independent variables having degrees equal to or lower than the highest degree of an independent variable in the firststage preparatory correction polynomial have been considered.
While the preferred embodiment only includes cross terms having independent variables with degrees lower than the independent variables in the firststage preparatory correction polynomial, this is not absolutely necessary. If the design criteria is better met by adding cross terms with higher degree variables or terms with only one independent variable, these terms may be evaluated by the procedure herein, selected after evaluation and added if they meet the selected criteria. Thus, if the criteria are a limit on the amount of memory required by a correction polynomial and a larger capacity memory is provided, new terms for addition to the correction polynomial may be evaluated by this procedure.
In the preferred embodiment, the process 24 continues the steps 56, 60, 118, or 66, 67 and 120 until a term is eliminated that was previously added, indicating that the polynomial is self consistent. While the language optimized calibration curve is used in this specification, the calibration curve need not be perfectly or completely developed into the most effective form to be considered an optimized calibration polynomial to be within the scope of this invention. The steps 60 and 66 need not be carried out to any fixed point since they make incremental improvements at each step. Moreover, the elimination of the least significant terms need not be alternated with the step of adding significant terms nor need the same number of terms be eliminated as added and only one of the processes of eliminating terms of lower significance of adding terms of more significance may be used. However, in the preferred embodiment, the steps of eliminating a term of low significance and adding a term of high significance are alternated and the repetitive process is continued until a term is eliminated that was added during the step in which terms are added or a term is added that had been earlier eliminated. At this point, the calibration curve is said to be self consistent.
In
The selection of the maximum degree of the independent variables and the number of terms of the polynomial generally involves a compromise between the complexity of the equipment, the size of the memory to be used and the precision of the measurements to be obtained. The precision of measurement increases with an increase in the number of terms, the degree of the variables and the number of cross terms. On the other hand, in embodiments implemented by a microcontroller, an increase in the number of terms, the degree of the variables, and the number of cross terms require an increase in the complexity and the size of memory. In the preferred embodiment, a microcontroller is used although the processes described herein in this specification can be performed by hand without the use of a microcontroller.
Once the highest degree of the independent variables has been selected, a firststage preparatory or working polynomial is formed including a plurality of terms. Some of these terms include powers of independent variables and others are cross terms including a product of two or more independent variables as shown at step 84. Terms that include the product of two or more independent variables are herein referred to as cross terms. The coefficients for the terms are selected using any known curve fitting method such as the well known least squares regression method. In the preferred embodiment, one or all of the independent variables of all degrees from the highest selected degree to the first degree are included. For example, if the fourth degree were to be selected for temperature, then a variable, C_{4}T^{4 }would have a coefficient, C_{4 }and all of the lower powers of temperature such as T^{3}, T^{2 }and T^{1 }have corresponding coefficients C_{3}, C_{2 }and C_{1}.
In
After completing the analysis of the data base, the average value of all of the maximum deviations for all data sets within the data base is determined in step 64. The value of the average maximum deviation is related to the rank of significance for that particular trial polynomial. The steps 58, 61, 62 and 64 are repeated as shown at step 63 for the next trial polynomial in order to determine its rank of significance. This process is continued until all trial polynomials have been ranked in significance. After all of the trial polynomials have been ranked, the trial polynomial with the lowest value for the average maximum deviation is selected as the best polynomial as shown at step 65. This trial polynomial has the least deviation or disagreement error with the data base. The process 60 applies the same to both secondstage and thirdstage trial preparatory calibration polynomials.
To determine which term improves the precision of the measurement the least when forming a secondstage preparatory correction polynomial, a series of trial correction polynomials are formed. A trial polynomial is formed by eliminating a term from the firststage correction polynomial, determining the optimum coefficients with the term removed to obtain a secondstage trial correction polynomial and eliminating secondstage trial correction polynomials that improved the precision of the measurement less than other secondstage trial correction polynomials being considered. These trial polynomials are compared to determine which trial polynomial provides the most precise determination of the dependent variable. The secondstage trial correction polynomial that provides the most precise measurement becomes a secondstage preparatory correction polynomial, which may result in a final secondstage preparatory correction polynomial except for one or more possible thirdstage correction polynomials formed by adding terms if the criteria is met or the correction polynomial becomes self consistent.
In the preferred embodiment, each of the secondstage preparatory correction polynomials is selected from a comparison of every trial secondstage polynomial formed from elimination of every cross term from one other first or thirdstage preparatory polynomial. However, it is possible to obtain good results by eliminating terms that include only one independent variable in addition to cross terms and to form and consider trial correction polynomials formed by eliminating terms from other secondstage correction polynomials and/or to form fewer trial correction polynomials for comparison by not forming trial polynomials by elimination of each cross term but only select some of the cross terms to be eliminated to form a trial polynomial.
In the preferred embodiment, secondstage trial correction polynomials are compared by calculating the dependent variable several times for each trial correction polynomial using a plurality of different data sets. Deviation between the calculated value and the actual value is determined for each trial polynomial at each of the plurality of data sets and a comparison made using these results. In the preferred embodiment, several measurements are made at each preset value for the dependent variable. For example, in the case in which pressure and temperature are the independent variables and the depth under water is the dependent variable, several measurements of pressures and temperatures are made with the test apparatus at known depths in generating the data base. Each trial correction polynomial is used to calculate depth from the values of pressure and temperature in the data base and the calculated values of depth compared to the measured values from the data base at a plurality of data sets from the data base. A data set in this specification is a set of measurements at different depths, temperatures and pressures. Several sets are obtained for each correction polynomial and each trial correction polynomial is tested with several sets.
In the preferred embodiment, the maximum deviation from the plurality of data sets for each trial correction polynomial is taken as the standard to be used in determining which term will be discarded although any other value directly related to the disagreement error and reflecting the effect of the term on improving the precision or the measurement such as a value within ten percent of the maximum deviation could be used. The average of the maximum deviations for each trial correction polynomial is calculated and this is used to determine the trial correction polynomial that causes the lowest average maximum deviation. The secondstage trial correction polynomial that causes the lowest average maximum deviation is taken as a secondstage correction polynomial. Although the average maximum deviation is used as the comparison standard in the preferred embodiment, any other indicator of the central tendency of the deviations such as for example the mean deviation could be used for this selection provided it results in a correction polynomial that provides a more precise measurement.
In
In
Steps 92, 94 and 96 in repetition generate only one data set. The data base is considered a collection or plurality of data sets. In other words, a data set is a subset of the data base. The data set only includes calibration data for particular transducers and a particular day. The data base includes calibration data for all transducers and all days. Calibration data collected on a particular transducer on the first day would constitute a complete data set. Calibration data collected on the same transducer on another day would constitute a separate but complete data set. Calibration data collected on another transducer would constitute yet another separate but complete data set. All three data sets would be contained in the data base. Elsewhere, the data base is also referred to as the plurality of data sets.
With the above procedure, an instrument is calibrated by forming a calibration polynomial to provide a calibrated signal indicating a measured value corrected for interference and nonlinearity. To prepare the calibration polynomial, a firststage preparatory correction polynomial including the primary independent variable, at least one interference related independent variable and some or no cross terms is first prepared. A plurality of data sets, that is a data base, of the dependent variable, the primary independent variable and the at least one interference related independent variable is obtained. The data sets are used to eliminate the least significant terms of the firststage preparatory correction polynomial and to add the most significant cross terms as described above.
In making an instrument, a design criteria is selected such as the precision needed or the size of memory that can be used to calibrate the instrument and a criteriaoptimized calibration curve is prepared that will satisfy this criteria. The primary transducers 18 and at least one secondary transducer 16 are selected to sense the values being measured and the interference related signals or factors that might reduce the precision of the measuring instrument. The microcontroller 14 is programmed to correct the value measured by the primary transducer 18. The instrument may use this corrected value to make further calculations and/or may provide a display of the corrected value.
In
The average velocity sensing system 114 includes an ultrasonic Doppler transmitter and receiver under the control of an automatic range and threshold setting system. The velocity meter transmits sound through a representative section of a flow stream or through the entire cross section of the flow stream and receives a complex signal back which is digitized and analyzed using a fast Fourier transform analyzer. With this arrangement, receive and transmit transducers 34 and 32 are positioned at an angle to the horizontal to radiate a beam to and receive reflections from a representative portion of the flow stream.
The resolution of the measurement depends on the number of ranges of frequencies selected for each term of the Fourier transform analyzer across the full range of frequency shifts caused by the range of possible velocities in the flow stream. The expected velocity range is determined in the preferred embodiment and 256 bands of frequencies are selected for positive and negative terms of the Fourier transform analyzer.
The words, “representative portion”, in this specification means a portion of the total flow stream which has a volume that includes within it smaller portions of fluid streams at each velocity flowing in the total flow stream with the fluid streams for each of the velocities of the smaller portions of flow streams having a reflective portion that is in the same proportion to the size of the reflective portion of the total flow stream having the same velocity as any other reflective smaller portion with a different velocity in the representative portion. That proportionality can be achieved in part by reflecting signals from a volume of the liquid rather than from an imaginary plane cutting the flow stream.
In this definition of a representative portion, each unit area flowing at a particular velocity in the representative portion has a ratio to the area of liquid of the total flow stream flowing at that rate which is the same ratio as every other cross sectional area flowing at that flow rate. Thus, this representative portion truly reflects the average flow rate of the entire flow stream.
In practice, some inaccuracy always occurs because of the failure to properly sample either the entire cross sectional area of the flow stream or a portion that is precisely a representative portion. Because it is easier to utilize a representative portion than the total cross sectional area, the preferred embodiment utilizes a representative portion and preferably arrives at this representative portion by selecting an angle at which the ultrasonic sound is reflected and selecting an angle at which it is received so that proportional amounts of the fluid flowing at each velocity reflect signals to the receiving transducer or transducers 34. One such error is caused by a lack of symmetry in the reflected signal with respect to an axis perpendicular to the transducer when the reflected signal is considered as a cone. It can be easily corrected, however. The lack of symmetry can be corrected by using a fixed factor such as two percent or other value to account for the discrepancy.
The average velocity sensing system 114 further includes an input circuit 40, a reflection processing circuit 42, a timecontrol and computation system 48A, the information input output system 49A, a timing circuit 44 and a transmitting signal generator 46. The input circuit 40 is electrically connected to the receiving transducer or transducer array 34 through a conductor 65 to receive signals therefrom, amplify them with automatic gain control and transmit the signals to the reflection processing circuit 42 through a conductor 74.
The reflection processing circuit 42 is electrically connected to: (1) the timing circuit 44 through conductors 50, 52, 54 and 58 which control the scanning of amplitudes and setting of a threshold value; and (2) the time control and computation system 48A through conductors 70 and 72 through which it transmits data for use by the time control and computation system 48A and through the conductor 68 from the time control and computation system 48A from which it which receives signals which control the time of transmission of data to the time control and computation system 48A.
The time control and computation system 48A is electrically connected to the timing circuit 44 through conductors 63 and 76 to control the synchronization of the entire flow meter 30A and to establish sampling rates and frequency cutoff points to the input circuit 40 to adjust the amplitude level of an automatic gain control circuit to obtain an adequate signal. The timing circuit 44 is electrically connected to the transmit signal generator 46 through conductor 61 through which it transmits signals to control the time at which the transmit signal generator 46 transmits signals to the transmitting transducer or transducer array 32 through a conductor 67. These signals control the sampling time and the repetition rate of the transmitted ultrasonic signals for the purpose of scanning across a range of sample times and rates for increased precision.
A pressure sensor 18A is electrically connected to the time control and computation system 48A to transmit depth information thereto and the information input output system 49A, which includes a computer keyboard and other input devices, supplies information to the time control and computation system 48A, such as a crosssectional area of the flow stream. With these values, the time control and computation system 48A is able to calculate the area of flow in the flow stream and the average velocity, and from that, calculate the rate of flow of liquid in the flow stream in a manner known in the art.
In general, an attempt is made to receive Doppler shift information from the entire crosssection of the stream. The signals are intended to represent all of the actual velocities and the crosssectional area of each of the velocities. These velocities and crosssectional areas may be represented in a curve, with the velocities being represented along the abscissa and the amount of area of the crosssection having each velocity or small range of velocities as the ordinates when viewed graphically. These values are measured with the Doppler frequency shift representing the velocity and the amplitude of the received ultrasonic signal as the area having that velocity. The average velocity system is explained in greater detail in the aforementioned U.S. Pat. No. 5,777,892, the disclosure of which is incorporated herein by reference.
While a somewhat detailed explanation of criteria optimized polymers has been provided, the words “criteriaoptimized polynomial” is not intended to be limited to all of the features of this explanation. The words “criteriaoptimized polynomial” refers to any procedure which includes the steps of forming a preparatory calibration polynomial including the primary independent variable and at least one interference related independent variable, obtaining a plurality of data sets of the dependent variable, the primary independent variable and at least one interference related independent variable and using the data sets to eliminate the least significant terms of the preparatory calibration polynomial or to add the most significant cross terms to the preparatory calibration polynomial until a criteria arbitrarily selected for a general purpose has been met or to add terms that provide a correction polynomial that improves the precision of the measurement.
In this specification, the words “criteriaoptimized” when applied to an instrument means an instrument using a criteriaoptimized polynomial to improve its precision and when applied to a calibration curve or polynomial means that the calibration curve or polynomial has been formed by first forming a preparatory or preparatory calibration polynomial including the primary independent variable and at least one interference related independent variable, obtaining a plurality of data sets of the dependent variable, the primary independent variable and the at least one interference related independent variable either before or after forming the polynomial and then using the data sets to eliminate the least significant terms of the preparatory calibration polynomial or to add the more significant terms.
From the above description, it can be understood that the instrument of this invention and the method of using it has several advantages, such as for example: (1) it can be more easily calibrated; (2) some aspects of it are more easily automated; and (3) it can provide more precise operation.
While a preferred embodiment of the invention has been described with some particularity, many modifications and variations in the system are possible without deviating from the invention. Therefore, it is to be understood that, within the scope of the appended claims, the invention may be practiced other than as specifically described.
Claims (8)
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

US11337011 US7532992B2 (en)  20060120  20060120  Measuring apparatuses and methods of using them 
Applications Claiming Priority (8)
Application Number  Priority Date  Filing Date  Title 

US11337011 US7532992B2 (en)  20060120  20060120  Measuring apparatuses and methods of using them 
KR20087020215A KR20090005293A (en)  20060120  20070117  Measuring apparatus and methods of using them 
JP2008551330A JP2009524055A (en)  20060120  20070117  Measuring device and how to use it 
CA 2637188 CA2637188A1 (en)  20060120  20070117  Measuring apparatuses and methods of using them 
PCT/US2007/001148 WO2007087199A3 (en)  20060120  20070117  Measuring apparatuses and methods of using them 
CN 200780006126 CN101473196A (en)  20060120  20070117  Measuring apparatus and methods of using them 
EP20070762617 EP1974289A4 (en)  20060120  20070117  Measuring apparatuses and methods of using them 
US12435155 US20090216477A1 (en)  20060120  20090504  Measuring apparatuses and methods of using them 
Related Child Applications (1)
Application Number  Title  Priority Date  Filing Date 

US12435155 Division US20090216477A1 (en)  20060120  20090504  Measuring apparatuses and methods of using them 
Publications (2)
Publication Number  Publication Date 

US20070174013A1 true US20070174013A1 (en)  20070726 
US7532992B2 true US7532992B2 (en)  20090512 
Family
ID=38286574
Family Applications (2)
Application Number  Title  Priority Date  Filing Date 

US11337011 Expired  Fee Related US7532992B2 (en)  20060120  20060120  Measuring apparatuses and methods of using them 
US12435155 Abandoned US20090216477A1 (en)  20060120  20090504  Measuring apparatuses and methods of using them 
Family Applications After (1)
Application Number  Title  Priority Date  Filing Date 

US12435155 Abandoned US20090216477A1 (en)  20060120  20090504  Measuring apparatuses and methods of using them 
Country Status (7)
Country  Link 

US (2)  US7532992B2 (en) 
EP (1)  EP1974289A4 (en) 
JP (1)  JP2009524055A (en) 
KR (1)  KR20090005293A (en) 
CN (1)  CN101473196A (en) 
CA (1)  CA2637188A1 (en) 
WO (1)  WO2007087199A3 (en) 
Cited By (3)
Publication number  Priority date  Publication date  Assignee  Title 

US20060253253A1 (en) *  20050509  20061109  Reynolds Christopher I  Environmental characteristic determination 
US20090216477A1 (en) *  20060120  20090827  Teledyne Isco, Inc.  Measuring apparatuses and methods of using them 
US20150107367A1 (en) *  20131021  20150423  Continental Automotive Systems, Inc.  Dual range high precision pressure sensor 
Families Citing this family (5)
Publication number  Priority date  Publication date  Assignee  Title 

US8152366B2 (en) *  20080123  20120410  University Of Delaware  Estimation of subsurface thermal structure using sea surface height and sea surface temperature 
EP2652456A1 (en) *  20101215  20131023  Jaguar Land Rover Limited  Wading vehicle depth measurement apparatus 
US9007377B2 (en) *  20110527  20150414  Molecular Devices, Llc  System and method for displaying parameter independence in a data analysis system 
FR3034863A1 (en) *  20150410  20161014  Cie Ind D'applications Thermiques  Methods for determining a volume of liquid at a predetermined temperature, of amounts of thermal energy calculation and determination of a hot water consumption profile 
US20170082471A1 (en) *  20150923  20170323  Texas Instruments Incorporated  Ultrasonic flow meter autotuning for reciprocal operation of the meter 
Citations (30)
Publication number  Priority date  Publication date  Assignee  Title 

US4481596A (en) *  19811102  19841106  Kaye Instruments Inc.  Method of and apparatus for automatically compensating for variations in output response characteristics of sensors and the like 
US5275042A (en)  19901004  19940104  Isco, Inc.  Variable gate flow analyzing method and apparatus 
US5371686A (en)  19920330  19941206  Isco, Inc.  Ultrasonic fluid velocity measuring system 
US5731995A (en)  19960412  19980324  Mc Donnell Douglas Corp  Method for determining pressure 
US5777892A (en)  19920330  19980707  Isco, Inc.  Doppler shift velocity measuring system with correction factors 
US5783836A (en)  19941228  19980721  Hoechst Aktiengesellschaft  Optical sensor apparatus for detecting vapor of organic solvent 
US5940780A (en)  19950929  19990817  Advanced Thermal Solutions, Inc.  Universal transceiver 
US6123904A (en)  19950210  20000926  Isco, Inc.  Water impurity analyzer 
US6308553B1 (en)  19990604  20011030  Honeywell International Inc  Selfnormalizing flow sensor and method for the same 
US6311136B1 (en)  19971126  20011030  Invensys Systems, Inc.  Digital flowmeter 
US6314380B1 (en)  19990603  20011106  Robert Bosch Corporation Corp Of Delaware  Ultrasound transducer temperature compensation methods, apparatus and programs 
US6352001B1 (en)  19990830  20020305  General Electric Company  Noniterative method for obtaining mass flow rate 
US6369712B2 (en)  19990517  20020409  The Goodyear Tire & Rubber Company  Response adjustable temperature sensor for transponder 
US6401045B1 (en)  19960604  20020604  Warren Rogers Associates, Inc.  Method and apparatus for monitoring operational performance of fluid storage systems 
US6449571B1 (en)  19990709  20020910  Mykrolis Corporation  System and method for sensor response linearization 
US6505519B2 (en)  20000323  20030114  Invensys Systems, Inc.  Correcting for twophase flow in a digital flowmeter 
US6681625B1 (en)  20000119  20040127  Lockheed Martin Corporation  Constanttemperaturedifference bidirectional flow sensor 
US6681787B2 (en)  19990709  20040127  Mykrolis Corporation  System and method of operation of a digital mass flow controller 
US6691061B1 (en)  19960604  20040210  Warren Rogers Associates, Inc.  Method and apparatus for monitoring operational performance of fluid storage systems 
US6701274B1 (en)  19990827  20040302  Rosemount Inc.  Prediction of error magnitude in a pressure transmitter 
US6769307B1 (en)  19971121  20040803  Perceptron, Inc.  Method and system for processing measurement signals to obtain a value for a physical parameter 
US6775632B1 (en)  19991214  20040810  The Goodyear Tire & Rubber Company  Calibration of a transponders for a tire pressure monitoring system 
US6777940B2 (en)  20021108  20040817  Ultima Labs, Inc.  Apparatus and method for resistivity well logging 
US6796173B1 (en)  19981009  20040928  Fti Flow Technology, Inc.  Fuel flowmeter 
US6850857B2 (en)  20010713  20050201  Honeywell International Inc.  Data fusion of stationary array sensor and scanning sensor measurements 
US6871551B2 (en)  20000628  20050329  Endress + Hauser Wetzer Gmbh + Co. Kg  Apparatus for generating and conducting a fluid flow, and method of monitoring said apparatus 
US6879926B2 (en)  20010629  20050412  National Instruments Corporation  Measurement system software architecture for easily creating highperformance measurement applications 
US6877387B1 (en)  19980710  20050412  Faure Herman  Method for measuring displacement of a fluid in a conduit 
US6892572B2 (en)  19940509  20050517  Automotive Technologies International, Inc.  Method and apparatus for measuring the quantity of a liquid in a vehicle container 
US20070288213A1 (en) *  20040722  20071213  Rainer Schantl  Method For Analyzing The Behavior Of Complex Systems, Especially Internal Combustion Engines 
Family Cites Families (6)
Publication number  Priority date  Publication date  Assignee  Title 

JP2579143B2 (en) *  19840202  19970205  ハネウエル・インコーポレーテッド  The method of digital correction of process variable sensors and process variable transmitter therefor 
US4754641A (en) *  19870210  19880705  Schlumberger Technology Corporation  Method and apparatus for measurement of fluid flow in a drilling rig return line 
US5691914A (en) *  19960724  19971125  American Sigma, Inc.  Fluid flow measurement correcting system, and methods of constructing and utilizing same 
JP2000180230A (en) *  19981217  20000630  Tokimec Inc  Flowmeter and method for measuring flow rate 
US7401497B2 (en) *  20051031  20080722  Honeywell International Inc.  Microdischarge detector method and apparatus 
US7532992B2 (en) *  20060120  20090512  Teledyne Isco, Inc.  Measuring apparatuses and methods of using them 
Patent Citations (35)
Publication number  Priority date  Publication date  Assignee  Title 

US4481596A (en) *  19811102  19841106  Kaye Instruments Inc.  Method of and apparatus for automatically compensating for variations in output response characteristics of sensors and the like 
US5275042A (en)  19901004  19940104  Isco, Inc.  Variable gate flow analyzing method and apparatus 
US5777892A (en)  19920330  19980707  Isco, Inc.  Doppler shift velocity measuring system with correction factors 
US5371686A (en)  19920330  19941206  Isco, Inc.  Ultrasonic fluid velocity measuring system 
US6892572B2 (en)  19940509  20050517  Automotive Technologies International, Inc.  Method and apparatus for measuring the quantity of a liquid in a vehicle container 
US5783836A (en)  19941228  19980721  Hoechst Aktiengesellschaft  Optical sensor apparatus for detecting vapor of organic solvent 
US6123904A (en)  19950210  20000926  Isco, Inc.  Water impurity analyzer 
US5940780A (en)  19950929  19990817  Advanced Thermal Solutions, Inc.  Universal transceiver 
US5731995A (en)  19960412  19980324  Mc Donnell Douglas Corp  Method for determining pressure 
US6691061B1 (en)  19960604  20040210  Warren Rogers Associates, Inc.  Method and apparatus for monitoring operational performance of fluid storage systems 
US6401045B1 (en)  19960604  20020604  Warren Rogers Associates, Inc.  Method and apparatus for monitoring operational performance of fluid storage systems 
US6769307B1 (en)  19971121  20040803  Perceptron, Inc.  Method and system for processing measurement signals to obtain a value for a physical parameter 
US6311136B1 (en)  19971126  20011030  Invensys Systems, Inc.  Digital flowmeter 
US6754594B2 (en)  19971126  20040622  Invensys Systems, Inc.  Digital flowmeter 
US6877387B1 (en)  19980710  20050412  Faure Herman  Method for measuring displacement of a fluid in a conduit 
US6796173B1 (en)  19981009  20040928  Fti Flow Technology, Inc.  Fuel flowmeter 
US6369712B2 (en)  19990517  20020409  The Goodyear Tire & Rubber Company  Response adjustable temperature sensor for transponder 
US6470286B1 (en)  19990603  20021022  Robert Bosch Corporation  Ultrasound transducer temperature compensation methods, apparatus and programs 
US6314380B1 (en)  19990603  20011106  Robert Bosch Corporation Corp Of Delaware  Ultrasound transducer temperature compensation methods, apparatus and programs 
US6782332B2 (en)  19990603  20040824  Robert Bosch Corporation  Ultrasound transducer temperature compensation methods, apparatus and programs 
US6308553B1 (en)  19990604  20011030  Honeywell International Inc  Selfnormalizing flow sensor and method for the same 
US6715339B2 (en)  19990604  20040406  Honeywell International Inc.  Selfnormalizing flow sensor and method for the same 
US6449571B1 (en)  19990709  20020910  Mykrolis Corporation  System and method for sensor response linearization 
US6681787B2 (en)  19990709  20040127  Mykrolis Corporation  System and method of operation of a digital mass flow controller 
US6701274B1 (en)  19990827  20040302  Rosemount Inc.  Prediction of error magnitude in a pressure transmitter 
US6352001B1 (en)  19990830  20020305  General Electric Company  Noniterative method for obtaining mass flow rate 
US6775632B1 (en)  19991214  20040810  The Goodyear Tire & Rubber Company  Calibration of a transponders for a tire pressure monitoring system 
US6681625B1 (en)  20000119  20040127  Lockheed Martin Corporation  Constanttemperaturedifference bidirectional flow sensor 
US6758102B2 (en)  20000323  20040706  Invensys Systems, Inc.  Correcting for twophase flow in a digital flowmeter 
US6505519B2 (en)  20000323  20030114  Invensys Systems, Inc.  Correcting for twophase flow in a digital flowmeter 
US6871551B2 (en)  20000628  20050329  Endress + Hauser Wetzer Gmbh + Co. Kg  Apparatus for generating and conducting a fluid flow, and method of monitoring said apparatus 
US6879926B2 (en)  20010629  20050412  National Instruments Corporation  Measurement system software architecture for easily creating highperformance measurement applications 
US6850857B2 (en)  20010713  20050201  Honeywell International Inc.  Data fusion of stationary array sensor and scanning sensor measurements 
US6777940B2 (en)  20021108  20040817  Ultima Labs, Inc.  Apparatus and method for resistivity well logging 
US20070288213A1 (en) *  20040722  20071213  Rainer Schantl  Method For Analyzing The Behavior Of Complex Systems, Especially Internal Combustion Engines 
NonPatent Citations (2)
Title 

Data Reduction and Error Analysis for the Physical Sciences, Philip R. Bevington, 1969, Chapters 8 and 9, pp. 134186. 
Eric W. Weisstein's "Least Squares FittingPolynomial"; Math Worldhttp://mathworld.wolfram.com/leastsquaresfittingpolynomial.html, Sep. 22, 2008. 
Cited By (4)
Publication number  Priority date  Publication date  Assignee  Title 

US20060253253A1 (en) *  20050509  20061109  Reynolds Christopher I  Environmental characteristic determination 
US7650238B2 (en) *  20050509  20100119  Northrop Grumman Corporation  Environmental characteristic determination 
US20090216477A1 (en) *  20060120  20090827  Teledyne Isco, Inc.  Measuring apparatuses and methods of using them 
US20150107367A1 (en) *  20131021  20150423  Continental Automotive Systems, Inc.  Dual range high precision pressure sensor 
Also Published As
Publication number  Publication date  Type 

EP1974289A4 (en)  20120926  application 
CN101473196A (en)  20090701  application 
CA2637188A1 (en)  20070802  application 
EP1974289A2 (en)  20081001  application 
WO2007087199A2 (en)  20070802  application 
US20090216477A1 (en)  20090827  application 
WO2007087199A3 (en)  20081204  application 
US20070174013A1 (en)  20070726  application 
KR20090005293A (en)  20090113  application 
JP2009524055A (en)  20090625  application 
Similar Documents
Publication  Publication Date  Title 

US5877997A (en)  Pulse echo distance measurement  
US5511425A (en)  Flaw detector incorporating DGS  
Stoffelen  Toward the true near‐surface wind speed: Error modeling and calibration using triple collocation  
US5827963A (en)  System and method for determining a density of a fluid  
US5544531A (en)  Flowmeter having active temperature compensation  
US6234019B1 (en)  System and method for determining a density of a fluid  
US5792951A (en)  Measurement system with self calibrating probe  
US20040258127A1 (en)  Ultrasound transducer fault measurement method and system  
US5072416A (en)  Method and apparatus for calibrating a flowmeter using a master meter and a prover  
US6701265B2 (en)  Calibration for vector network analyzer  
EP0807804A2 (en)  Method for calibrating a differential pressure fluid flow measuring system  
US5716272A (en)  Moisture/yield monitor grain simulator  
US6290652B1 (en)  Correction of blood pressure measurements in invasive liquidfilled systems  
Ludwig et al.  Calibration of hydrophones based on reciprocity and time delay spectrometry  
US20080034893A1 (en)  Method and Apparatus for Determining Flow Pressure Using Density Information  
WO1982003121A1 (en)  Bridge circuit compensation for environmental effects  
US7984637B2 (en)  System and method for field calibration of flow meters  
US6643610B1 (en)  Process transmitter with orthogonalpolynomial fitting  
US5596508A (en)  High resolution measurement of a thickness using ultrasound  
US6085577A (en)  Surface tension measurement in a pressurized environment  
JPH10111152A (en)  Flowmeter and massflowrate control apparatus  
US6257049B1 (en)  Ambient humidity measurement using microwaves  
US20030074953A1 (en)  Device for determining the change in the density of a medium  
US6345683B1 (en)  System and method for an improved device for measuring water levels  
US6490535B1 (en)  Method and apparatus for calibrating an instrument 
Legal Events
Date  Code  Title  Description 

AS  Assignment 
Owner name: TELEDYNE ISCO, INC., NEBRASKA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:LIESCHESKI, PHILLIP B.;REEL/FRAME:017514/0380 Effective date: 20060120 

REMI  Maintenance fee reminder mailed  
LAPS  Lapse for failure to pay maintenance fees  
FP  Expired due to failure to pay maintenance fee 
Effective date: 20130512 