CN108240836B - 多维折线段测量定标法 - Google Patents
多维折线段测量定标法 Download PDFInfo
- Publication number
- CN108240836B CN108240836B CN201711329155.XA CN201711329155A CN108240836B CN 108240836 B CN108240836 B CN 108240836B CN 201711329155 A CN201711329155 A CN 201711329155A CN 108240836 B CN108240836 B CN 108240836B
- Authority
- CN
- China
- Prior art keywords
- measurement
- broken line
- line segment
- points
- calculating
- 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.)
- Active
Links
Classifications
-
- 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 apparatus or arrangements provided for in groups G01D1/00 - G01D15/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
- G06F17/175—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method of multidimensional data
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种多维折线段测量定标法,包括如下步骤:S1、确定被测量、直接测量参数和影响测量的所有因素;S2、确定每种影响测量因素的测量点数和直接测量参数的测量点数,通过标准的被测量值测得直接测量参数的值,即定标样本数据;直接测量参数构成一维定标样本数据,每增加一种影响测量的因素,就相对多一维定标样本数据;S3、基于定标样本数据进行折线段定标测量;折线段的定标由最高维向一维逐步归一推进。本发明解决了多维折线段的定标问题,突破了折线段定标使用的局限性,特别是当测量点数越细致时,能够大幅提高测量精度。本发明在多个测试仪器的定标中使用,取得了很好的效果。
Description
技术领域
本发明涉及测试测量领域,具体涉及一种多维折线段测量定标法。
背景技术
在测试测量领域中,被测量常受单个或多个因素的影响,并且被测量与影响因素可能存在非线性关系。目前常采用折线段法或查表法定标,由于折线段定标法是把两个测量点之间的非测量点用折线段近似,而查表定标法对非测量点仅用测量点简单代替,所以折线段定标法比查表定标法测量精度更高,使用更广泛。但是当影响测量的因素越多,即定标样本数据维数越多时,查表定标法可以使用;对于折线段定标法,二维定标样本数据如何使用存在难度,特别是当影响被测量的因素越多时,定标样本数据维数会越多,折线段定标法的使用更成问题。
发明内容
为解决上述问题,本发明提供了一种多维折线段测量定标法。
本发明是通过以下技术方案实现的。
多维折线段测量定标法,包括如下步骤:
S1、确定被测量、直接测量参数和影响测量的所有因素;
S2、确定每种影响测量因素的测量点数和直接测量参数的测量点数,通过标准的被测量值测得直接测量参数的值,即定标样本数据;直接测量参数构成一维定标样本数据,每增加一种影响测量的因素,就相对多一维定标样本数据;
S3、基于定标样本数据进行折线段定标测量;折线段的定标由最高维向一维逐步归一推进。
本发明解决了多维折线段的定标问题,突破了折线段定标使用的局限性,特别是当测量点数越细致时,能够大幅提高测量精度。本发明在多个测试仪器的定标中使用,取得了很好的效果。
具体实施方式
下面对本发明的实施例作详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
对比例
(1)测量案例1:被测量为y,直接测量参数为x,y与x成非线性关系,没有其它影响测量的因素。
测量案例1为一维情况。根据测量精度要求,首先确定测量定标样本数据的点数(测量精度要求越高,划分的越细,测量的定标样本数据点数越多,下同),假定为n个测量点,再利用标准的Y值,测得对应的X,X为一维数组,即X[]={X1,X2,X3。。。。。。Xn},定标样本数据见表1。
表1一维定标样本数据
Y | Y<sub>1</sub> | Y<sub>2</sub> | Y<sub>3</sub> | 。。。。。。 | Y <sub>n</sub>-1 | Y<sub>n</sub> |
X | X<sub>1</sub> | X<sub>2</sub> | X<sub>3</sub> | 。。。。。。 | X<sub>n-1</sub> | X<sub>n</sub> |
对于查表定标法,直接根据实测的x值在表中的范围,选择Y赋给y,最终计算出被测量 y。查表法按条件赋值如下:
y=Y1(X1≤x<X2)
y=Y2(X2≤x<X3)
y=Yn-1(Xn-1≤x<Xn)
对于折线段定标法,利用两点法确定n-1个一次函数,y为因变量,x为自变量,测量时,根据实测的x值,选择不同的一次函数,计算出被测量y。折线段方法按条件选择公式如下:
y=(Y2-Y1)/(X2-X1)×(x-X1)+Y1(X1≤x<X2)
y=(Y3–Y2)/(X3–X2)×(x-X2)+Y2(X2≤x<X3)
y=(Yn–Yn-1)/(Xn–Xn-1)×(x-Xn-1)+Yn-1(Xn-1≤x)
(2)测量案例2:被测量为y,直接测量参数为x,y与x成非线性关系,并且a因素影响测量。
测量案例2为两维情况。按照测量精度要求,首先确定因素a的测量点数,假定为n个测量点,即A1、A2....An,然后在每种a因素的情况下,利用标准的Y值,测得m个对应的X, X为两位维数组,即X[][]={XA1Y1,XA1Y2,XA1Y3。。。。。。XA1Ym;XA2Y1,XA2Y2,XA2Y3。。。。。。XA2Ym;。。。。。。;XAnY1,XAnY2,XAnY3。。。。。。XAnYm}。定标样本数据见表2。
对于查表定标法,直接根据实测的a和x的值在表中的范围,选择Y赋值给y,最终计算出被测量y。查表法按条件赋值如下:
y=Y1(A1≤a<A2,XA1Y1≤x<XA1Y2)
y=Y2(A1≤a<A2,XA1Y2≤x<XA1Y3)
y=Y3(A1≤a<A2,XA1Y3≤x<XA1Y4)
y=Y1(An-1≤a<An,XA1Y1≤x<XA1Y2)
y=Y2(An-1≤a<An,XA1Y2≤x<XA1Y3)
y=Y3(An-1≤a<An,XA1Y3≤x<XA1Y4)
本发明实施例提供了一种多维折线段测量定标法,包括如下步骤:以上述的二维举例说明:
第一步,确定被测量、直接测量参数和影响测量的所有因素。
本例被测量为y,直接测量参数为x,y与x成非线性关系,并且a因素影响测量。
第二步,确定每种影响测量因素的测量点数和直接测量参数的测量点数,通过标准的被测量值测得直接测量参数的值,即定标样本数据。直接测量参数构成一维定标样本数据,每增加一种影响测量的因素,就相对多一维定标样本数据,比如影响测量的因素有两种,则最终的定标样本数据为三维。
本例影响测量因素a的测量点数为n,即A1、A2....An,直接测量参数x的测量点数为m,由于测量因素只有一种,所以定标样本数据为二维。利用标准的Y值测得对应的X,二维定标样本数据见表2。
第三步,基于定标样本数据进行折线段定标测量。折线段的定标由最高维向一维逐步归一推进,实质是折线段从最高维向一维转化的过程,并且每维的转换过程包含子步骤。
假定本例a的实际测量值为A',且A1≤A'≤A2,x的实际测量值为X',且(XA1Y1≤X'≤XA1Y2)&&(XA2Y1≤X'≤XA2Y2),计算Y'的步骤如下:
第1步,采用两点[(A1,XA1Y1),(A2,XA2Y1)]确定一条线的折线段法计算xA'Y1。xA'Y1=(XA2Y1-XA1Y1)/(A2-A1)*(A'-A1)+XA1Y1。
第2步,采用两点[(A1,XA1Y2),(A2,XA2Y2)]确定一条线的折线段法计算xA'Y2。xA'Y2=(XA2Y2-XA1Y2)/(A2-A1)*(A'-A1)+XA1Y2。
第3步,采用两点[(xA'Y1,Y1),(xA'Y2,Y2)]计算A'情况下yA'与xA'之间的折线段,yA'=(Y2-Y1)/(xA'Y2-xA'Y1)*(xA'-xA'Y1)+Y1。
第4步,利用第三步的折线段计算Y'。Y'=(Y2-Y1)/(xA'Y2-xA'Y1)*(X'-xA'Y1)+ Y1。
当然X'可以任意,假定a的实际测量值为A',且A1≤A'≤A2,x的实际测量值为X',且 (XA1Y1≤X'≤XA1Y2)&&(XA2Y3≤X'≤XA2Y4),计算Y'的步骤如下:
第1步,采用两点[(A1,XA1Y1),(A2,XA2Y3)]确定一条线的折线段法计算xA'Y1-3。xA'Y1-3=(XA2Y3-XA1Y1)/(A2-A1)*(A'-A1)+XA1Y1。
第2步,采用两点[(A1,XA1Y2),(A2,XA2Y4)]确定一条线的折线段法计算xA'Y2-4。xA'Y2-4=(XA2Y4-XA1Y2)/(A2-A1)*(A'-A1)+XA1Y2。
第3步,采用两点[(xA'Y1-3,Y1),(xA'Y2-4,Y2)]计算A'情况下yA'与xA'之间的折线段,yA'=(Y2-Y1)/(xA'Y2-4-xA'Y1-3)*(xA'-xA'Y1-3)+Y1。
第4步,利用第三步的折线段计算Y'。Y'=(Y2-Y1)/(xA'Y2-4-xA'Y1-3)*(X'-xA'Y1-3) +Y1。
实施例1
被测量为音频电压,记为v,直接测量参数为ADC采样数据,记为d,温度是影响测量的因素,记为t。
第一步,确定被测量、直接测量参数和影响测量的所有因素。本例被测量为音频电压v,直接测量参数为ADC采样数据d,影响测量的因素为温度t。
第二步,确定每种影响测量因素的测量点数和直接测量参数的测量点数,通过标准的被测量值测得直接测量参数的值,即定标样本数据。本例影响测量因素t的测量点数为n,即T1、T2....Tn,直接测量参数d的测量点数为m,由于测量因素只有一种,所以定标样本数据为二维。利用标准的V值测得对应的D,定标样本数据见表3。
第三步,基于定标样本数据进行折线段定标测量。折线段的定标由最高维向一维逐步归一推进,实质是折线段从最高维向一维转化的过程,并且每维的转换过程包含子步骤。
假定本例t的实际测量值为T',且T1≤T'≤T2,d的实际测量值为D',且(DT1V1≤D'≤DT1V2)&&(DT2V1≤D'≤DT2V2),计算V'的步骤如下:
第1步,采用两点[(T1,DT1V1),(T2,DT2V1)]确定一条线的折线段法计算dT'V1。dT'V1=(DT2V1-DT1V1)/(T2-T1)*(T'-T1)+DT1V1。
第2步,采用两点[(T1,DT1V2),(T2,DT2V2)]确定一条线的折线段法计算dT'V2。dT'V2=(DT2V2-DT1V2)/(T2-T1)*(T'-T1)+DT1V2。
第3步,采用两点[(dT'V1,V1),(dT'V2,V2)]计算T'情况下vT'与dT'之间的折线段,vT'=(V2-V1)/(dT'V2-dT'V1)*(dT'-dT'V1)+V1。
第4步,利用第三步的折线段计算V'。V'=(V2-V1)/(dT'V2-dT'V1)*(D'-dT'V1)+ V1。
当然D'可以任意,假定t的实际测量值为T',且T1≤T'≤T2,d的实际测量值为D',且 (DT1V1≤D'≤DT1V2)&&(DT2V3≤D'≤DT2V4),计算V'的步骤如下:
第1步,采用两点[(T1,DT1V1),(T2,DT2V3)]确定一条线的折线段法计算dT'V1-3。dT'V1-3=(DT2V3-DT1V1)/(T2-T1)*(T'-T1)+DT1V1。
第2步,采用两点[(T1,DT1V2),(T2,DT2V4)]确定一条线的折线段法计算dT'V2-4。dT'V2-4=(DT2V4-DT1V2)/(T2-T1)*(T'-T1)+DT1V2。
第3步,采用两点[(dT'V1-3,V1),(dT'V2-4,V2)]计算T'情况下vT'与dT'之间的折线段,vT'=(V2-V1)/(dT'V2-4-dT'V1-3)*(dT'-dT'V1-3)+V1。
第4步,利用第三步的折线段计算V'。V'=(V2-V1)/(dT'V2-4-dT'V1-3)*(D'-dT'V1-3) +V1。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (1)
1.多维折线段测量定标法,其特征在于,包括如下步骤:
步骤S1、确定被测量、直接测量参数和影响测量的所有因素;
被测量为y,直接测量参数为x,y与x成非线性关系,并且a因素影响测量;
步骤S2、确定每种影响测量因素的测量点数和直接测量参数的测量点数,通过标准的被测量值测得直接测量参数的值,即定标样本数据;直接测量参数构成一维定标样本数据,每增加一种影响测量的因素,就相对多一维定标样本数据;
影响测量因素a的测量点数为n,即A1、A2....An,直接测量参数x的测量点数为m,由于测量因素只有一种,所以定标样本数据为二维.利用标准的Y值测得对应的X,二维定标样本数据;
步骤S3、基于定标样本数据进行折线段定标测量;折线段的定标由最高维向一维逐步归一推进;
所述步骤S2按照测量精度要求,首先确定因素a的测量点数,假定为n个测量点,即A1、A2....An,然后在每种a因素的情况下,利用标准的Y值,测得m个对应的X,X为两位维数组,即X[][]={XA1Y1,XA1Y2,XA1Y3......XA1Ym;XA2Y1,XA2Y2,XA2Y3......XA2Ym;......;XAnY1,XAnY2,XAnY3......XAnYm};建立一个或者多个二维定标样本数据表;对于查表定标法,直接根据实测的a和x的值在表中的范围,选择Y赋值给y,最终计算出被测量y.查表法按条件赋值如下:
y=Y1(A1≤a<A2,XA1Y1≤x<XA1Y2);
y=Y2(A1≤a<A2,XA1Y2≤x<XA1Y3);
y=Y3(A1≤a<A2,XA1Y3≤x<XA1Y4);
......
y=Y1(An-1≤a<An,XA1Y1≤x<XA1Y2);
y=Y2(An-1≤a<An,XA1Y2≤x<XA1Y3);
y=Y3(An-1≤a<An,XA1Y3≤x<XA1Y4);
所述步骤S3包括:
步骤S3.1:假定a的实际测量值为A',且A1≤A'≤A2,x的实际测量值为X',且(XA1Y1≤X'≤XA1Y2)&&(XA2Y1≤X'≤XA2Y2),计算Y';
所述步骤S3.1包括:
步骤S3.1.1:采用两点[(A1,XA1Y1),(A2,XA2Y1)]确定一条线的折线段法计算xA'Y1.xA'Y1=(XA2Y1-XA1Y1)/(A2-A1)*(A'-A1)+XA1Y1;
步骤S3.1.2:采用两点[(A1,XA1Y2),(A2,XA2Y2)]确定一条线的折线段法计算xA'Y2.xA'Y2=(XA2Y2-XA1Y2)/(A2-A1)*(A'-A1)+XA1Y2;
步骤S3.1.3:采用两点[(xA'Y1,Y1),(xA'Y2,Y2)]计算A'情况下yA'与xA'之间的折线段,yA'=(Y2-Y1)/(xA'Y2-xA'Y1)*(xA'-xA'Y1)+Y1;
步骤S3.1.4:利用第三步的折线段计算Y'.Y'=(Y2-Y1)/(xA'Y2-xA'Y1)*(X'-xA'Y1)+Y1;
步骤S3.2:当X'取任意值时,假定a的实际测量值为A',且A1≤A'≤A2,x的实际测量值为X',且(XA1Y1≤X'≤XA1Y2)&&(XA2Y3≤X'≤XA2Y4),计算Y';
所述步骤S3.2包括:
步骤S3.2.1,采用两点[(A1,XA1Y1),(A2,XA2Y3)]确定一条线的折线段法计算xA'Y1- 3.xA'Y1-3=(XA2Y3-XA1Y1)/(A2-A1)*(A'-A1)+XA1Y1;
步骤S3.2.2,采用两点[(A1,XA1Y2),(A2,XA2Y4)]确定一条线的折线段法计算xA'Y2- 4.xA'Y2-4=(XA2Y4-XA1Y2)/(A2-A1)*(A'-A1)+XA1Y2;
步骤S3.2.3,采用两点[(xA'Y1-3,Y1),(xA'Y2-4,Y2)]计算A'情况下yA'与xA'之间的折线段,yA'=(Y2-Y1)/(xA'Y2-4-xA'Y1-3)*(xA'-xA'Y1-3)+Y1;
步骤S3.2.4,利用第三步的折线段计算Y'.Y'=(Y2-Y1)/(xA'Y2-4-xA'Y1-3)*(X'-xA'Y1-3)+Y1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711329155.XA CN108240836B (zh) | 2017-12-13 | 2017-12-13 | 多维折线段测量定标法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711329155.XA CN108240836B (zh) | 2017-12-13 | 2017-12-13 | 多维折线段测量定标法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108240836A CN108240836A (zh) | 2018-07-03 |
CN108240836B true CN108240836B (zh) | 2020-11-10 |
Family
ID=62700194
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711329155.XA Active CN108240836B (zh) | 2017-12-13 | 2017-12-13 | 多维折线段测量定标法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108240836B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6358164A (ja) * | 1986-08-28 | 1988-03-12 | Shimadzu Corp | デ−タ処理装置 |
CN103632059A (zh) * | 2013-12-04 | 2014-03-12 | 中国电子科技集团公司第五十研究所 | 测试测量中多维十字定标法 |
CN104501996A (zh) * | 2014-12-16 | 2015-04-08 | 珠海拓普智能电气股份有限公司 | 一种基于内空型温控件的多温区定标方法 |
-
2017
- 2017-12-13 CN CN201711329155.XA patent/CN108240836B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6358164A (ja) * | 1986-08-28 | 1988-03-12 | Shimadzu Corp | デ−タ処理装置 |
CN103632059A (zh) * | 2013-12-04 | 2014-03-12 | 中国电子科技集团公司第五十研究所 | 测试测量中多维十字定标法 |
CN104501996A (zh) * | 2014-12-16 | 2015-04-08 | 珠海拓普智能电气股份有限公司 | 一种基于内空型温控件的多温区定标方法 |
Non-Patent Citations (1)
Title |
---|
复合型两点多段非均匀性校正算法研究;郑军等;《计算机仿真》;20060830(第8期);第106-109页以及157页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108240836A (zh) | 2018-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110286345B (zh) | 一种矢量网络分析仪在片s参数的校准方法、系统及设备 | |
US20070007942A1 (en) | Automatic non-linear phase response calibration and compensation for a power measurement device | |
CN109115257B (zh) | 传感器特性曲线的修正方法、装置、设备及存储介质 | |
Wahab | Interpolation and extrapolation | |
CN108664069B (zh) | 带隙基准电路的校准方法及装置 | |
CN110726475A (zh) | 红外测温标定方法、装置、红外热成像设备及存储装置 | |
CN112946560A (zh) | 电能表校准方法、装置、电能表及电能表系统 | |
JP2017058371A (ja) | 自動車の圧力センサの非線形誤差の減少 | |
CN115616266B (zh) | 一种波形幅度值测量方法、装置、终端设备以及存储介质 | |
CN112165328A (zh) | 一种adc数模转换误差的校正方法、装置及介质 | |
CN113048871A (zh) | 一种基于谐波分量的电容位移检测非线性实时校准方法 | |
CN108240836B (zh) | 多维折线段测量定标法 | |
CN112415458A (zh) | 一种电流传感器线性度测试系统及校准方法 | |
CN103063233B (zh) | 一种采用多传感器降低测量误差的方法 | |
TWI803916B (zh) | 一種晶片測試機測量晶片內阻的方法 | |
CN108107392B (zh) | 多线trl校准方法及终端设备 | |
CN112420066B (zh) | 降噪方法、装置、计算机设备和计算机可读存储介质 | |
CN113804232A (zh) | 传感器校准方法、装置以及终端设备 | |
CN109471151A (zh) | 信号处理方法、装置及用户终端 | |
CN109655749B (zh) | 采集数据矩阵智能跟踪校准方法及装置 | |
CN103076035B (zh) | 一种基于双支持向量机的传感器测量方法 | |
CN116128770A (zh) | 一种图像生成方法、装置、电子设备及存储介质 | |
CN104467701A (zh) | 一种功率放大器的电压校正方法及电子终端 | |
CN106788436B (zh) | 应用于saradc中的pip电容阵列的电压系数校准方法 | |
CN102095668A (zh) | 乳胶漆粘度检测方法以及温-粘修正系数的测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |