CN108240836B - 多维折线段测量定标法 - Google Patents

多维折线段测量定标法 Download PDF

Info

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
Application number
CN201711329155.XA
Other languages
English (en)
Other versions
CN108240836A (zh
Inventor
杨文举
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Microwave Technology CETC 50 Research Institute
Original Assignee
Shanghai Institute of Microwave Technology CETC 50 Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Microwave Technology CETC 50 Research Institute filed Critical Shanghai Institute of Microwave Technology CETC 50 Research Institute
Priority to CN201711329155.XA priority Critical patent/CN108240836B/zh
Publication of CN108240836A publication Critical patent/CN108240836A/zh
Application granted granted Critical
Publication of CN108240836B publication Critical patent/CN108240836B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING 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/00Testing or calibrating apparatus or arrangements provided for in groups G01D1/00 - G01D15/00
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • G06F17/175Function 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。
Figure RE-GDA0001672378180000031
对于查表定标法,直接根据实测的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。
Figure RE-GDA0001672378180000051
第三步,基于定标样本数据进行折线段定标测量。折线段的定标由最高维向一维逐步归一推进,实质是折线段从最高维向一维转化的过程,并且每维的转换过程包含子步骤。
假定本例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。
CN201711329155.XA 2017-12-13 2017-12-13 多维折线段测量定标法 Active CN108240836B (zh)

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)

* Cited by examiner, † Cited by third party
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 珠海拓普智能电气股份有限公司 一种基于内空型温控件的多温区定标方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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