CN103162901A  Nonlinear calibrating method for multiple temperature points of pressure sensor  Google Patents
Nonlinear calibrating method for multiple temperature points of pressure sensor Download PDFInfo
 Publication number
 CN103162901A CN103162901A CN 201310103444 CN201310103444A CN103162901A CN 103162901 A CN103162901 A CN 103162901A CN 201310103444 CN201310103444 CN 201310103444 CN 201310103444 A CN201310103444 A CN 201310103444A CN 103162901 A CN103162901 A CN 103162901A
 Authority
 CN
 China
 Prior art keywords
 pressure
 value
 adc
 temperature
 interpolation
 Prior art date
Links
 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
 238000005259 measurement Methods 0.000 claims abstract description 11
 238000005070 sampling Methods 0.000 claims abstract description 6
 230000000875 corresponding Effects 0.000 claims description 28
 238000000034 method Methods 0.000 claims description 7
 238000004364 calculation method Methods 0.000 claims description 6
 239000007788 liquid Substances 0.000 claims description 6
 241001269238 Data Species 0.000 claims description 2
 230000003321 amplification Effects 0.000 claims description 2
 238000000151 deposition Methods 0.000 claims description 2
 238000009795 derivation Methods 0.000 claims description 2
 239000011159 matrix material Substances 0.000 claims description 2
 238000003199 nucleic acid amplification method Methods 0.000 claims description 2
 230000000630 rising Effects 0.000 claims description 2
 238000003860 storage Methods 0.000 claims description 2
 238000006467 substitution reaction Methods 0.000 claims description 2
 239000000523 sample Substances 0.000 claims 1
 238000010586 diagram Methods 0.000 description 2
 238000005516 engineering process Methods 0.000 description 2
 238000004519 manufacturing process Methods 0.000 description 2
 XLYOFNOQVPJJNPUHFFFAOYSAN water Substances data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTAwLjUwMScgeT0nMTcwJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+SDwvdGV4dD4KPHRleHQgeD0nMTI2LjExNCcgeT0nMTg2JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjI2cHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+MjwvdGV4dD4KPHRleHQgeD0nMTM4JyB5PScxNzAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID5PPC90ZXh0Pgo8L3N2Zz4K data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTMuMzEyOCcgeT0nNTMuNTkwOScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToyM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6I0U4NDIzNScgPkg8L3RleHQ+Cjx0ZXh0IHg9JzI4LjE1NjknIHk9JzYyLjg2MzYnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MTVweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiNFODQyMzUnID4yPC90ZXh0Pgo8dGV4dCB4PSczNS4wNDU1JyB5PSc1My41OTA5JyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjIzcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojRTg0MjM1JyA+TzwvdGV4dD4KPC9zdmc+Cg== O XLYOFNOQVPJJNPUHFFFAOYSAN 0.000 description 2
 238000001514 detection method Methods 0.000 description 1
Abstract
The invention relates to the relevant technical field of pressure characteristic value acquisition of pressure sensors, in particular to a method and a system for acquiring a pressure characteristic value of a pressure sensor with a nonlinear temperature compensation algorithm. The method includes the steps: sampling the pressure characteristic value of the pressure sensor to obtain a pressure characteristic value D; calculating a temperature drift correction coefficient delta t related with the outside temperature according to the outside temperature; and calculating a corrected value D+delta Dt of the pressure characteristic value. The sensor is calibrated by an upper computer, the upper computer makes a data sheet according to a series of data measured at different temperatures and under different pressures in calibration as the upper computer can execute largescale computation, a calculating parameter for calculating the pressure value can be rapidly found in an EEProm (electrically erasable programmable read only memory) sheet according to signals measured by a temperature sensor and the pressure sensor when a lower computer is used for measurement, and an accurate pressure value is obtained by calculating. As the temperature drift correction coefficient is added, acquisition accuracy of the pressure characteristic value is remarkably and particularly improved by the aid of a onestep compensation method in a hightemperature area and a lowtemperature area.
Description
Technical field
The present invention relates to pressure transducer pressure characteristic value and gather correlative technology field, particularly a kind of pressure transducer pressure characteristic value that contains the nonlinear temperature compensation algorithm is converted into the method and system of force value.
Background technology
The manufacturer of pressure transducer provides the pressure transducer of strip resistance network building out mode mostly now, and this sensor has stable performance, circuit design characteristics simple, easy to use.Because the limitation of manufacture craft can't satisfy high precision and high conforming requirement, but precision generally can only satisfy 1% left and right to this sensor due to resistor network.
Although some sensor software compensations, because the imperfection of its computing method, precision can only accomplish 5 ‰, and the described algorithm of patent of the present invention can be accurate to 0.5 ‰ to measuring accuracy, totally can improve an order of magnitude.
In the fields such as early warning of military affairs, aviation, earthquake and disaster, the pressure transducer and the liquid level gauge that need, very strict to the requirement of precision, such as detecting the instrument of aircraft flight height on aircraft; Ensure the Boiler Steam control system that naval vessel energy highspeed smooth is travelled on warship and load water yield control system; System of the prepolice observation WATER LEVEL CHANGES of the disasteies such as reply earthquake etc.Use common pressure transducer in these applications, will because of the low precision of sensor, measure out of true inaccurate with control, the result that finally leads to errors.In the higher field of accuracy requirement, many temperature spots gamma correction method of the present invention is more important.
Summary of the invention
The invention provides a kind of pressure transducer pressure characteristic value acquisition method and system that contains the nonlinear temperature compensation algorithm, to solve prior art in the pressure transducer detection system, because the making of algorithm imperfection or compensating circuit is coarse, and the inaccurate technical matters of pressure characteristic value sampling precision that causes.
The technical scheme that adopts is as follows:
Due to the inconsistency of the pressure transducer of actual production, be issued to highprecision pressure at wide temperature range and measure, need to demarcate each pressure transducer.The scaling method that uses is that pressure transducer is placed in temperature control box, and the pressure criteria that use can be traced to the source produces nominal pressure, provides different temperatures by incubator, provides the different pressures value by pressure controller for pressure transducer.Measure the A/D collection value P of the output voltage V p to obtain reflecting actual reception pressure by the A/D converter
_{adc}Obtain simultaneously reflecting the A/D collection value T of the output voltage V t of environment temperature
_{adc}, in conjunction with the standard setting pressure P of pressure source
_{set}Just can reflect the nonlinear and temperature characterisitic of pressure transducer comprehensively.In order to illustrate the process of demarcation, with 5 Temperature Settings, 5 pressure settings are example.
First pressure sensor temperature is set as minimum temperature spot T1, to the pressure transducer pressurization, from the pressure P of minimum
_{Set (1)}Be added to successively maximum pressure P
_{Set (5)}, record simultaneously the A/D collection value P of the Vp of each spot pressure
_{adc}(1,1), P
_{adc}(2,1) ..., P
_{adc}(5,1); Also to record T corresponding to this temperature
_{adc}(1).
Again the pressure transducer environment temperature is set as time low temperature spot T2, to the pressure transducer pressurization, from the pressure P of minimum
_{Set (1)}Be added to successively maximum pressure P
_{Set (5)}, record simultaneously the A/D collection value P of the Vp of each spot pressure
_{adc}(1,2), P
_{adc}(2,2) ..., P
_{adc}(5,2); Record equally T corresponding to this temperature
_{adc}(2).
Adjust successively temperature from T1 until T5 repeats to demarcate and records P
_{adc}(m, n) and T
_{adc}(n), wherein m represent pressure from low to high the order of order (m=1,2 ..., 5).N represent temperature from low to high the order of order (n=1,2 ..., 5).
Just obtained like this three arrays in our hand:
1 " 5 of pressure source standard setting pressure: Pset (m).
2 " 5 of different temperatures collection values:
(n).
3 " under different temperatures, the collection value of different pressures is 5*5=25 altogether:
(m, n).
These nominal datas are kept to prepare against in the host computer hard disk when being used for making the form of downloading to slave computer EEProm when host computer is demarcated slave computer use.
[005] the present invention due in pressure transducer management system, has adopted the correction of temperature drift algorithm, makes pressure characteristic value acquisition precision higher.Owing to having added the correction of temperature drift coefficient, take especially the method for single compensation at high and low warm area, significantly improved the precision that gathers the pressure characteristic value.Simultaneously, due to gaging pressure cell pressure eigenwert more accurately, therefore can better the detected pressures sensor.
Description of drawings
[006] Fig. 1 is the pressure transducer schematic diagram
Fig. 2 is the modular structure figure of the embodiment of the present invention;
Fig. 3 is that host computer of the present invention asks EEProm tables of data and slave computer usage data table to carry out the process flow diagram of calculation of pressure by interpolation;
Fig. 4 is that interpolation of the present invention obtains revised tonogram;
Embodiment
[007] below in conjunction with the drawings and specific embodiments, the present invention is described in further details, as shown in Figure 2.
[008] the computer control pressure controller is delivered to gas on solenoid valve with certain pressure, controlled by solenoid valve controller again and open that way solenoid valve, gas arrives on normal pressure meter and tested pressure transducer, then controls temperature by incubator by solenoid valve, orifice valve.
[009] at certain temperature, add testing pressure to pressure transducer, measure the A/D collection value P of current sensor output voltage Vp by the A/D converter
_{ad}And the A/D collection value T of output voltage V t
_{ad}, the nominal data that the combined standard pressure gauge obtains just can obtain pressure measuring value P by correction, and the computation process of correction was divided into for two steps.
The first step is calculated, and is to the method for nominal data by cubic spline interpolation, from the A/D collection value T of current output voltage Vt
_{ad}Interpolation obtains the A/D collection value P ' of the output voltage V p under all pressure calibration points
_{ad}Process successively the data under the different pressures set point, the explanation as an example of certain spot pressure order force order m' example, m'=1 here, 2,3,4,5.As shown in Figure 3, transverse axis is Tadc (m', the n) data under current pressure point order force order m', and one has 5 temperature spot data, and the longitudinal axis is Padc (m', n) data corresponding to each temperature spot.Obtain P ' (m ') data under current pressure point order m ' according to the Tad numerical interpolation of current measurement.Need to calculate successively corresponding P ' (m ') interpolated data to all 5 pressure setpoint.
Because the variation of the A/D image data of uniform pressure, different temperatures is little, but larger on the correction impact, therefore adopt cubic spline functions method interpolation in the sampling point interval, so both can avoid the unsettled shortcoming of high order interpolation precision, and can guarantee again the slickness of interpolating function curve integral body.Cubic spline functions S (Tadc (m')) the expression formula in [Tadc (m', k), Tadc (m', k+1)] (k=0,1,2,3,4) here is:
Use threemoment method can simplify the calculated amount of cubic spline functions, under natural boundary conditions, a0=a4=0 is arranged, other unknown quantity a1, a2, a3 satisfy following Diagonal Equation group:
Wherein:
${\mathrm{\μ}}_{k}=\frac{{h}_{k}}{{h}_{k1}+{h}_{k}}$
This system of equations matrix of coefficients diagonal dominance is considered the ubiquity of algorithm, can use chasing method to find the solution.Determine between the location of interpolation T ad, with interpolation T ad substitution (1) formula, just can be in the hope of interpolation result P'(m').Interpolation outside sampling point is used linear interpolation.Cubic spline interpolation can obtain higher interpolation precision, but also brings more calculated amount simultaneously.Because calculated amount is very large, the singlechip microcomputer that sensor is used is difficult to realization, and these calculate all and must complete by host computer when sensor being carried out the linearization demarcation.
Second step calculates, and is that the A/D collection value P' to the output voltage V p under all pressure calibration points obtained in the previous step carries out interpolation, obtains revised pressure measuring value P from the A/D collection value Pad interpolation of current output voltage Vp.As shown in Fig. 4, transverse axis is the P'(m' under all spot pressure order m') data, one has 5 groups of data, and the longitudinal axis is each selfcorresponding Pset (m') data of all spot pressure order m'.Carry out linear interpolation according to the Pad numerical value of current measurement and just can obtain the revised pressure P of current measurement data.This step calculates only needs to use linear interpolation can guarantee measuring accuracy.
After we obtain the AD of current temperature value and pressure transducer, by [010], [011], [012], [013] described algorithm, can obtain accurate force value.But the calculating that this process need is a large amount of due to the speed of singlechip microcomputer and the limitation of hardware resource, is completed this calculated amount very difficult.So also need this algorithm is done further processing.
When sensor being carried out linearization calibration, we are by test, the AD collection value when obtaining five temperature that under five spot pressures, each pressure is corresponding: P
_{adc}(m, n)., take first force value as example, from the minimum temperature value, how many AD collection values that computer software utilizes this algorithm can calculate 1 ℃ of correspondence of every rising should be.Calculate maximum temperature always.So we just obtain one when actual pressure is constant, the array of the pressure AD value that every degree temperature is corresponding: P
_{adc}(1, n ').Do same processing to remaining four force value, just obtain one to the new twodimensional array after the amplification of test gained twodimensional array: P
_{adc}(m, n '), n ' is that the value of going forward one by one is 1 ℃ from low temperature being begun the order to maximum temperature here.
Open up 5 data blocks in slave computer EEprom storage space, DD value corresponding to different temperatures when depositing respectively five spot pressures used during for measurement.
Computing machine is passed to sensor to the array that [015] is made by serial ports, and sensor stores these data into EEProm corresponding position.
Slave computer is as follows to the use step of temperature compensation algorithm:
The first step, when beginning to measure, first use the temperature sensor measurement temperature, obtain trying to achieve by tabling lookup after temperature this temperature value be maximum pressure value, minimum pressure values with maximum pressure value and downforce value between after the quartern value: the corresponding P0 of DD_Min(), the corresponding P1 of DD_M1(), the corresponding P2 of DD_M2(), the corresponding P3 of DD_M3(), the corresponding P4 of DD_Max().
Second step, measure DD value current and pressure core body force value, by seeing relatively the DD value is at which section of four sections pressure stages.
The 3rd step, to illustrate for example at minimum one section: i.e. DD_Min＜DD＜DD_M1, curve a: P=a*DD+b who asks at this section.
Asking equational a, b value: P0 and DD_Min, P1 and DD_M1 to bring formula into, to obtain the linear equation in two unknowns formula as follows:
Ｐ０＝ａ＊ＤＤ＿Ｍｉｎ＋ｂ
Ｐ１＝ａ＊ＤＤ＿Ｍ１＋ｂ
P0, P1 are by artificial fixed in temperature compensation, and DD_Min, DD_M1 table lookup and obtain, and solving equations obtains a, b value.
The 4th step, this time, a and b became numerical value known, that the DD value gathers for AD.At this time just can calculate force value according to P=a*DD+b.
The 5th step, because P=ρ h, (ρ is the density of liquid) after obtaining pressure value P, just can go out liquid level according to this derivation of equation: h=P/ ρ.
If measure the DD value lower than DD_Min or higher than DD_Max, test specification has surpassed range, and software can carry out outrange and report to the police.If the DD value can be according to abovementioned computing method within other scope, corresponding force value and the level value asked.
Above said method is just to have gathered AD value corresponding to five temperature spots under five spot pressures, explain orally host computer and make calibration with the algorithm of twodimensional array, in the use of the array of slave computer with used a curve to explain orally to the computing method of force value.Do like this product of producing and accomplished 5/10000ths precision, reach domestic and international advanced level.If increase host computer needed number pressure and temperature number when making the calibration array, adopt quafric curve or cubic curve during slave computer calculating pressure value, just can obtain higher precision.
Claims (7)
1. many temperature spots of pressure transducer of the present invention gamma correction method is characterized in that: the pressure criteria that use can be traced to the source produces nominal pressure, provides different temperatures by incubator, provides the different pressures value by pressure controller for pressure transducer.
2. many temperature spots of pressure transducer of the present invention gamma correction method is characterized in that: measure by the A/D converter and can obtain reflecting that sensor probe bears the A/D collection value P of the output voltage V p of pressure
_{adc}Obtain simultaneously reflecting the A/D collection value T of the output voltage V t of environment temperature
_{adc}, in conjunction with the standard setting pressure P of pressure source
_{set}Just can reflect the nonlinear and temperature characterisitic of pressure transducer, in order to illustrate the process of demarcation, with 5 Temperature Settings, 5 pressure settings are example comprehensively;
First pressure sensor temperature is set as minimum temperature spot T1, to the pressure transducer pressurization, from the pressure P of minimum
_{Set (1)}Be added to successively maximum pressure P
_{Set (5)}, record simultaneously the A/D collection value P of the Vp of each spot pressure
_{adc}(1,1), P
_{adc}(2,1) ..., P
_{adc}(5,1); Also to record T corresponding to this temperature
_{adc}(1);
Again the pressure transducer environment temperature is set as time low temperature spot T2, to the pressure transducer pressurization, from the pressure P of minimum
_{Set (1)}Be added to successively maximum pressure P
_{Set (5)}, record simultaneously the A/D collection value P of the Vp of each spot pressure
_{adc}(1,2), P
_{adc}(2,2) ..., P
_{adc}(5,2); Record equally T corresponding to this temperature
_{adc}(2);
Adjust successively temperature from T1 until T5 repeats to demarcate and records P
_{adc}(m, n) and T
_{adc}(n), wherein m represent pressure from low to high the order of order (m=1,2 ..., 5); N represent temperature from low to high the order of order (n=1,2 ..., 5);
Just obtained like this three arrays in our hand:
1 " 5 of pressure source standard setting pressure: Pset (m);
2 " 5 of different temperatures collection values:
(n);
3 " under different temperatures, the collection value of different pressures is 5*5=25 altogether:
(m, n);
These nominal datas are kept to prepare against in the host computer hard disk when being used for making the form of downloading to slave computer EEProm when host computer is demarcated slave computer use.
3. based on the described algorithm of feature 2, from the A/D collection value T of current output voltage Vt
_{ad}The A/D collection value P ' that interpolation obtains the output voltage V p under all pressure calibration points processes the data under the different pressures set point successively, the explanation as an example of certain spot pressure order force order m' example, m'=1 here, 2,3,4,5; As shown in Figure 4, transverse axis is Tadc (m', the n) data under current pressure point order force order m', and one has 5 temperature spot data, and the longitudinal axis is Padc (m', n) data corresponding to each temperature spot; Obtain P ' (m ') data under current pressure point order m ' according to the Tad numerical interpolation of current measurement; Need to calculate successively corresponding P ' (m ') interpolated data to all 5 pressure setpoint;
Because the variation of the A/D image data of uniform pressure, different temperatures is little, but larger on the correction impact, therefore adopt cubic spline functions method interpolation in the sampling point interval, so both can avoid the unsettled shortcoming of high order interpolation precision, and can guarantee again the slickness of interpolating function curve integral body; Cubic spline functions S (Tadc (m')) the expression formula in [Tadc (m', k), Tadc (m', k+1)] (k=0,1,2,3,4) here is:
Use threemoment method can simplify the calculated amount of cubic spline functions, under natural boundary conditions, a0=a4=0 is arranged, other unknown quantity a1, a2, a3 satisfy following Diagonal Equation group:
Wherein:
${\mathrm{\μ}}_{k}=\frac{{h}_{k}}{{h}_{k1}+{h}_{k}}$
This system of equations matrix of coefficients diagonal dominance, consider the ubiquity of algorithm, can use chasing method to find the solution, determine between the location of interpolation T ad, with interpolation T ad substitution (1) formula, just can be in the hope of interpolation result P'(m'), interpolation outside sampling point is used linear interpolation, cubic spline interpolation can obtain higher interpolation precision, but also brings more calculated amount simultaneously, can use optimized algorithm or operating speed singlechip microcomputer faster.
4. described based on feature 2, the A/D collection value P' of the output voltage V p under all pressure calibration points obtained in the previous step is carried out interpolation, obtain revised pressure measuring value P from the A/D collection value Pad interpolation of current output voltage Vp; As shown in Fig. 4, transverse axis is the P'(m' under all spot pressure order m') data, one has 5 groups of data, the longitudinal axis is each selfcorresponding Pset (m') data of all spot pressure order m', carries out linear interpolation according to the Pad numerical value of current measurement and just can obtain the revised pressure P of current measurement data; This step calculates only needs to use linear interpolation can guarantee measuring accuracy; Open up 5 data blocks in slave computer EEprom storage space, DD value corresponding to different temperatures when depositing respectively five spot pressures used during for measurement.
5. described based on feature 2, feature 3 and feature 4, we are by test, the AD collection value when obtaining five temperature that under five spot pressures, each pressure is corresponding: P
_{adc}(m, n); , take first force value as example, from the minimum temperature value, how many AD collection values that computer software utilizes this algorithm can calculate 1 ℃ of correspondence of every rising should be; Calculate maximum temperature always, so we just obtain one when actual pressure is constant, the array of the pressure AD value that every degree temperature is corresponding: P
_{adc}(1, n '), do same processing to remaining four force value, just obtain one to the new twodimensional array after the amplification of test gained twodimensional array: P
_{adc}(m, n '), n ' is that the value of going forward one by one is 1 ℃ from low temperature being begun the order to maximum temperature here; Computing machine is by serial ports, array P
_{adc}(m, n ') passes to sensor, and sensor stores these data into EEProm corresponding position.
6. many temperature spots of pressure transducer of the present invention gamma correction method is characterized in that:
Slave computer is as follows to the use step of temperature compensation algorithm:
The first step, when beginning to measure, first use the temperature sensor measurement temperature, obtain asking by tabling lookup after temperature this temperature value be maximum pressure value, minimum pressure values with maximum pressure value and downforce value between value after the quartern: the corresponding P0 of DD_Min(), the corresponding P1 of DD_M1(), the corresponding P2 of DD_M2(), the corresponding P3 of DD_M3(), the corresponding P4 of DD_Max();
Second step, measure DD value current and pressure core body force value, by seeing relatively the DD value is at which section of four sections pressure stages;
The 3rd step, to illustrate for example at minimum one section: i.e. DD_Min＜DD＜DD_M1, curve a: P=a*DD+b who asks at this section;
Asking equational a, b value: P0 and DD_Min, P1 and DD_M1 to bring formula into, to obtain the linear equation in two unknowns formula as follows:
Ｐ０＝ａ＊ＤＤ＿Ｍｉｎ＋ｂ
Ｐ１＝ａ＊ＤＤ＿Ｍ１＋ｂ
P0, P1 are by artificial fixed in temperature compensation, and DD_Min, DD_M1 table lookup and obtain, and solving equations obtains a, b value;
The 4th step, this time, a and b became numerical value known, that the DD value gathers for AD; At this time just can calculate force value according to P=a*DD+b;
The 5th step, because P=ρ h, (ρ is the density of liquid) after obtaining pressure value P, just can go out liquid level according to this derivation of equation: h=P/ ρ.
7. described based on feature 6, slave computer uses once curve to come the pressure of calculation pressure sensor, is going out liquid depth with calculation of pressure, if require more highprecision product, can calculate with quafric curve or cubic curve.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN 201310103444 CN103162901A (en)  20130328  20130328  Nonlinear calibrating method for multiple temperature points of pressure sensor 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN 201310103444 CN103162901A (en)  20130328  20130328  Nonlinear calibrating method for multiple temperature points of pressure sensor 
Publications (1)
Publication Number  Publication Date 

CN103162901A true CN103162901A (en)  20130619 
Family
ID=48586172
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN 201310103444 CN103162901A (en)  20130328  20130328  Nonlinear calibrating method for multiple temperature points of pressure sensor 
Country Status (1)
Country  Link 

CN (1)  CN103162901A (en) 
Cited By (21)
Publication number  Priority date  Publication date  Assignee  Title 

CN103837300A (en) *  20140319  20140604  成都千嘉科技有限公司  Pressure sensor calibration method with temperature compensation function 
CN103968998A (en) *  20140412  20140806  沈阳仪表科学研究院有限公司  Method for conducting pressure calibration on diffused silicon sensors through upper computer 
CN104359611A (en) *  20141027  20150218  温州大学  Method for discriminating nonlinearity of pressure sensors under influence of internal stress 
CN104614125A (en) *  20141205  20150513  贵州航天凯山石油仪器有限公司  Pressure gage zero drift solving method 
CN104697711A (en) *  20150215  20150610  广东风华芯电科技股份有限公司  Pressure calibration method and device for tire pressure sensor 
CN104776957A (en) *  20150410  20150715  北京振兴计量测试研究所  Pressure sensor calibration method and pressure sensor calibration device 
CN104964790A (en) *  20150612  20151007  广东电网有限责任公司电力科学研究院  Method for correcting dynamic pressure in combustion chamber measured by pressure guiding tube 
CN105092145A (en) *  20150527  20151125  重庆川仪自动化股份有限公司  Temperature and pressure compensation linear correction method on intelligent transmitter 
CN105258846A (en) *  20151008  20160120  歌尔声学股份有限公司  Method and device for calibrating pressure sensor 
US20160033991A1 (en) *  20140729  20160204  Sensirion Ag  Sensor chip 
CN105527056A (en) *  20140928  20160427  沈阳中科奥维科技股份有限公司  Temperature referencebased pressure compensation calibration method 
CN104359611B (en) *  20141027  20170104  温州大学  A kind of internal stress influence lower pressure sensor nonlinear discriminant method 
CN106610331A (en) *  20151027  20170503  沈阳中科奥维科技股份有限公司  Control method for pressure meter calibration for petroleum transmission pipelines in petroleum industry 
CN106644193A (en) *  20170127  20170510  武汉立易方科技有限公司  Pressure intensity value determination method and system 
CN106840516A (en) *  20170302  20170613  成都信息工程大学  A kind of pressure gauge temperature drift scaling method based on fitting of a polynomial 
CN108362311A (en) *  20180118  20180803  上海兰宝传感科技股份有限公司  A kind of distance measuring sensor segmented compensated distance method 
CN109323796A (en) *  20181027  20190212  哈尔滨工业大学  A kind of full warm area range pressure sensor scaling method 
CN109323792A (en) *  20181115  20190212  航宇救生装备有限公司  A kind of intelligent temperaturecompensation method based on ejector seat program controller pressure sensor 
CN109425461A (en) *  20170905  20190305  上海融德机电工程设备有限公司  Intelligent pressure transmitter and its temperaturecompensation method 
CN111458070A (en) *  20200417  20200728  齐亚斯(上海)物联网科技有限公司  Hightemperature melt pressure transmitter system 
CN111681392A (en) *  20200513  20200918  河南省南阳水文水资源勘测局  River hydrological measurement system 

2013
 20130328 CN CN 201310103444 patent/CN103162901A/en not_active Application Discontinuation
Cited By (27)
Publication number  Priority date  Publication date  Assignee  Title 

CN103837300B (en) *  20140319  20160316  成都千嘉科技有限公司  With the pressure sensor calibration method of temperature compensation 
CN103837300A (en) *  20140319  20140604  成都千嘉科技有限公司  Pressure sensor calibration method with temperature compensation function 
CN103968998A (en) *  20140412  20140806  沈阳仪表科学研究院有限公司  Method for conducting pressure calibration on diffused silicon sensors through upper computer 
US10718639B2 (en)  20140729  20200721  Sensirion Ag  Sensor chip 
US20160033991A1 (en) *  20140729  20160204  Sensirion Ag  Sensor chip 
CN105318893A (en) *  20140729  20160210  盛思锐股份公司  Sensor chip 
CN105527056A (en) *  20140928  20160427  沈阳中科奥维科技股份有限公司  Temperature referencebased pressure compensation calibration method 
CN104359611A (en) *  20141027  20150218  温州大学  Method for discriminating nonlinearity of pressure sensors under influence of internal stress 
CN104359611B (en) *  20141027  20170104  温州大学  A kind of internal stress influence lower pressure sensor nonlinear discriminant method 
CN104614125A (en) *  20141205  20150513  贵州航天凯山石油仪器有限公司  Pressure gage zero drift solving method 
CN104697711A (en) *  20150215  20150610  广东风华芯电科技股份有限公司  Pressure calibration method and device for tire pressure sensor 
CN104776957A (en) *  20150410  20150715  北京振兴计量测试研究所  Pressure sensor calibration method and pressure sensor calibration device 
CN105092145B (en) *  20150527  20171226  重庆川仪自动化股份有限公司  Applied to the temperature and pressure compensation linear correction method on intelligent transducer 
CN105092145A (en) *  20150527  20151125  重庆川仪自动化股份有限公司  Temperature and pressure compensation linear correction method on intelligent transmitter 
CN104964790A (en) *  20150612  20151007  广东电网有限责任公司电力科学研究院  Method for correcting dynamic pressure in combustion chamber measured by pressure guiding tube 
CN104964790B (en) *  20150612  20171212  广东电网有限责任公司电力科学研究院  The modification method of dynamic pressure in combustion chamber is measured using pressure guiding pipe 
CN105258846A (en) *  20151008  20160120  歌尔声学股份有限公司  Method and device for calibrating pressure sensor 
CN106610331A (en) *  20151027  20170503  沈阳中科奥维科技股份有限公司  Control method for pressure meter calibration for petroleum transmission pipelines in petroleum industry 
CN106644193B (en) *  20170127  20200414  武汉立易方科技有限公司  Method and system for measuring pressure value 
CN106644193A (en) *  20170127  20170510  武汉立易方科技有限公司  Pressure intensity value determination method and system 
CN106840516A (en) *  20170302  20170613  成都信息工程大学  A kind of pressure gauge temperature drift scaling method based on fitting of a polynomial 
CN109425461A (en) *  20170905  20190305  上海融德机电工程设备有限公司  Intelligent pressure transmitter and its temperaturecompensation method 
CN108362311A (en) *  20180118  20180803  上海兰宝传感科技股份有限公司  A kind of distance measuring sensor segmented compensated distance method 
CN109323796A (en) *  20181027  20190212  哈尔滨工业大学  A kind of full warm area range pressure sensor scaling method 
CN109323792A (en) *  20181115  20190212  航宇救生装备有限公司  A kind of intelligent temperaturecompensation method based on ejector seat program controller pressure sensor 
CN111458070A (en) *  20200417  20200728  齐亚斯(上海)物联网科技有限公司  Hightemperature melt pressure transmitter system 
CN111681392A (en) *  20200513  20200918  河南省南阳水文水资源勘测局  River hydrological measurement system 
Similar Documents
Publication  Publication Date  Title 

Morris et al.  Measurement and instrumentation: theory and application  
CN100533087C (en)  Error compensation method of highprecision pressure sensor  
CN103398807B (en)  A kind of novel threedimensional stress measurement sensor  
US4437164A (en)  Ridge circuit compensation for environmental effects  
RU2453931C1 (en)  Improved temperature compensation of multivariate pressure gage  
US10422631B2 (en)  Optical fiber and method and apparatus for accurate fiber optic sensing under multiple stimuli  
CN1856700B (en)  Calibration of a process pressure sensor  
CN100554902C (en)  Be applicable to the temperature simple calibrating method of reflectiontype polarizationpreserving fiber temperature sensor  
CN103837300B (en)  With the pressure sensor calibration method of temperature compensation  
CN107728094A (en)  Current calibration coefficient measurement apparatus and method, current sensing means and method  
da Silva Hack et al.  Measurement uncertainty: Literature review and research trends  
CN101936791B (en)  Digital pressure gauge  
CN101793493B (en)  Precision improvement calibrating method for current vortex sensor  
US20130218502A1 (en)  Temperature compensated pressure transducer  
CN102597707B (en)  With the angle detector of compound selfcorrecting function  
CN104953283B (en)  A kind of phased array antenna quick calibration method and calibration system  
CN101943615A (en)  Temperature measuring device and method based on Raman light reflection  
CN102305612A (en)  Displacement/deflection measuring system and method  
CN103364751B (en)  A kind of vector network analyzer Electronic Calibration part and calibration steps  
CN103001583A (en)  Temperature compensation method and crystal oscillator  
CN103575414B (en)  A kind of electric thermocouple temperature measurement Error Compensation method  
JPWO2013100156A1 (en)  Output value correction method for physical quantity sensor device, output value correction method for physical quantity sensor, physical quantity sensor device, and output value correction device for physical quantity sensor  
CN203719811U (en)  Weak voltage signal detection device of microvoltage sensor  
GB2085597A (en)  Method and apparatus for detemining the mass flow of a fluid  
CN106706153B (en)  A kind of processing method of low temperature temperaturesensitive resistancetemperature characteristic data 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
WD01  Invention patent application deemed withdrawn after publication 
Application publication date: 20130619 

C02  Deemed withdrawal of patent application after publication (patent law 2001) 