《一种测量值的粗大误差判别方法.pdf》由会员分享,可在线阅读,更多相关《一种测量值的粗大误差判别方法.pdf(7页珍藏版)》请在专利查询网上搜索。
1、(10)申请公布号 CN 104121877 A (43)申请公布日 2014.10.29 C N 1 0 4 1 2 1 8 7 7 A (21)申请号 201410387859.2 (22)申请日 2014.08.07 G01B 21/00(2006.01) (71)申请人天津航天长征火箭制造有限公司 地址 300462 天津市滨海新区经济技术开发 区西区新业六街19号西区投资服务中 心301号 申请人中国运载火箭技术研究院 (72)发明人李新友 赵瑞峰 乔志峰 李强 陈乃玉 冯叶素 李刚 张志博 杨中宝 熊占兵 马彦辉 徐方强 (74)专利代理机构天津滨海科纬知识产权代理 有限公司 12。
2、211 代理人李震勇 (54) 发明名称 一种测量值的粗大误差判别方法 (57) 摘要 本发明创造提供一种测量值的粗大误差判别 方法,目的是解决在大量测量值中,利用激光跟踪 仪测量薄壁件位姿所测量的数据一般不超过10 个,在剔除粗大误差时,现有的3判断准则已经 不适用的问题。本发明创造提供的一种测量值的 粗大误差判别方法,算法简单、结果可靠,可以查 找小数据量范围内的粗大误差。除此之外,该算 法易于用计算机程序实现,只需要把所需数据输 入程序,即可运行程序得到分析结果,应用领域广 泛。 (51)Int.Cl. 权利要求书1页 说明书4页 附图1页 (19)中华人民共和国国家知识产权局 (12)。
3、发明专利申请 权利要求书1页 说明书4页 附图1页 (10)申请公布号 CN 104121877 A CN 104121877 A 1/1页 2 1.一种测量值的粗大判别方法,其特征在于:该方法包括以下步骤: 步骤一、在刚体上选取n个测量点A i (i1,2n),要求n个测量点不在一条直线上, 通过测量n个点在空间的坐标就可以确定刚体的空间位姿; 步骤二、记刚体上n个测量点A i 在动坐标系o-xyz下的坐标向量为:a i x Ai y Ai z Ai T (i1,2n),点A i 在定坐标系o 1 -x 1 y 1 z 1 中的实际测量值记为A i ,以向量a i x Ai y Ai z A。
4、i T 表示; 步骤三、刚体的实际位姿用转换矩阵R和位置向量T表示,那么点A i 的实测值与通 过坐标转换得到的坐标值之间的误差可以表示为: a i a i -(Ra i +T) (i1,2n) (1) 步骤四、用向量2-范数来评价该误差的大小,可以写成: a i 2 a i -(Ra i +T) 2 (i1,2n) (2) 步骤五、取向量2-范数的平方和: 步骤六、设测量仪器的测量误差范围为-,+,则式(1)的最大值为: a i a i -(Ra i +T) T (4) 步骤七、每个位姿下的目标函数都对应一个最小值,这些最小值是有一个取值范围的, 这些最小值可能达到的最大值为: 步骤八、式(。
5、3)的取值范围是03n 2 ,如果式(3)的值超出此范围,则认 定测量值中至少有一个是粗大误差,为找出该点,依次剔除测量值A j (j1,2 n),然后,利用n-1个测量值计算目标函数,如果只有某次剔除后的n-1个测量值 A 1 ,A 2 .A j-1 ,A j+1 .A n 对应的目标函数值在03(n-1) 2 范围内,而其余 次剔除后n-1个测量值对应的目标函数值超出03(n-1) 2 ,则认定该测量值A j 为 粗大误差,需要剔除; 步骤九、如果n次剔除计算所对应目标函数均超出范围,则认定数据中有1个以上的粗 大误差,此时利用n-2个测量值计算目标函数,判别方法与只有1个粗大误差的方法相。
6、同, 以此类推,直到剔除所有粗大误差。 权 利 要 求 书CN 104121877 A 1/4页 3 一种测量值的粗大误差判别方法 技术领域 0001 本发明创造属于计量测试领域,尤其是涉及一种测量值的粗大误差判别方法。 背景技术 0002 火箭箭体部段、飞机机身和机翼等都属于大型弱刚性薄壁件,在自身重力和外压 力作用下会产生变形,致使两对接面的定位基准很难对正,成为火箭大部件对接装配的瓶 颈问题。如何评估火箭大部段的姿态成为解决上述问题的关键。目前,抑制大型弱刚性薄 壁件变形的方法主要是在薄壁件端面安装刚性端框,对于大型薄壁件的测量主要采用激光 跟踪仪等仪器,通过测量刚性端框表示薄壁件的位置。
7、和姿态。在测量中,会出现人为的测量 误差,即粗大误差。在大量测量值中,比较成熟的剔除粗大误差的方法是3准则。而利用 激光跟踪仪测量薄壁件位姿所测量的数据一般不超过10个,在剔除粗大误差时3准则已 经不适用。 0003 目前,国内航空航天领域所使用的调整机构有手动支架、桁架;国外的航空航天领 域已经开始使用自动化的调姿设备,使用比较多的是POGO柱群。支撑机翼的POGO柱多达 10个,通过联动调解实现机翼的位置和姿态调整,然后实现对接。自动化调姿设备已经在国 外的多种机型中使用,波音737、787等,空客A320等。国内也开始了对自动化对接设备的 研究,国内的浙江大学、中航工业研究所在飞机大部件。
8、装配中引进国外思想,采用了三自由 度POGO柱机构,通过多个POGO柱协调运动进行大部段的调姿装配。 0004 航空航天领域内的大型薄壁件尺寸大,很多是通过螺栓连接的,连接面上少则数 十个连接孔,作则几百,甚至上千个连接孔。连接之前需要对薄壁件进行位置和姿态的调 整,使得连接面对正,才能进行连接。如果有粗大误差存在的话,将导致错误的计算结果。目 前,国内大型薄壁件调整机构在使用过程中,都需要从测量仪器读取测量值,然后根据测量 值计算驱动量,如果有粗大误差存在的话,则会导致不能正确调姿,甚至损伤薄壁件。 发明内容 0005 本发明创造要解决以上大型薄壁件测量时的遇到的问题,提供一种测量值的粗大 。
9、误差判别方法,该方法可以查找测量数据中的误差。 0006 为解决上述技术问题,本发明创造采用的技术方案是:一种测量值的粗大判别方 法,该方法包括以下步骤: 0007 步骤一、在刚体上选取n个测量点A i (i1,2n),要求n个测量点不在一条直线 上,通过测量n个点在空间的坐标就可以确定刚体的空间位姿; 0008 步骤二、记刚体上n个测量点A i 在动坐标系o-xyz下的坐标向量为:a i x Ai y Ai z Ai T (i1,2n),点A i 在定坐标系o 1 -x 1 y 1 z 1 中的实际测量值记为A i ,以向量a i x Ai y Ai z Ai T 表示; 0009 步骤三、。
10、刚体的实际位姿用转换矩阵R和位置向量T表示,那么点A i 的实测值 与通过坐标转换得到的坐标值之间的误差可以表示为: 说 明 书CN 104121877 A 2/4页 4 0010 a i a i -(Ra i +T) (i1,2n) (1) 0011 步骤四、用向量2-范数来评价该误差的大小,可以写成: 0012 a i 2 a i -(Ra i +T) 2 (i1,2n) (2) 0013 步骤五、取向量2-范数的平方和: 0014 0015 步骤六、设测量仪器的测量误差范围为-,+,则式(1)的最大值为: 0016 a i a i -(Ra i +T) T (4) 0017 步骤七、每个。
11、位姿下的目标函数都对应一个最小值,这些最小值是有一个取值范 围的,这些最小值可能达到的最大值为: 0018 0019 步骤八、式(3)的取值范围是03n 2 ,如果式(3)的值超出此范围,则 认定测量值中至少有一个是粗大误差,为找出该点,依次剔除测量值A j (j1,2 n),然后,利用n-1个测量值计算目标函数,如果只有某次剔除后的n-1个测量值 A 1 ,A 2 .A j-1 ,A j+1 .A n 对应的目标函数值在03(n-1) 2 范围内,而其余 次剔除后n-1个测量值对应的目标函数值超出03(n-1) 2 ,则认定该测量值A j 为 粗大误差,需要剔除; 0020 步骤九、如果n次。
12、剔除计算所对应目标函数均超出范围,则认定数据中有1个以上 的粗大误差,此时利用n-2个测量值计算目标函数,判别方法与只有1个粗大误差的方法相 同,以此类推,直到剔除所有粗大误差。 0021 本发明创造具有的优点和积极效果是:一种测量值的粗大误差判别方法,算法简 单、结果可靠,可以查找小数据量范围内的粗大误差。除此之外,该算法易于用计算机程序 实现,只需要把所需数据输入程序,即可运行程序得到分析结果,应用领域广泛。 附图说明 0022 图1是刚性件位置姿态变化示意图。 0023 图中:左侧代表位姿1,位姿1是刚性件在变换之前的位姿,右侧代表位姿2,位姿 2是刚性件在变换之后的位姿;A1、A2、A。
13、3、A4、A5是刚性件上5个测量点。 具体实施方式 0024 下面结合附图对本发明创造的具体实施例做详细说明。 0025 如图1所示,一种粗大误差判别方法的具体步骤为: 0026 步骤一、在刚体上选取n个测量点A i (i1,2n),最佳实施例中取n5。要求 n个测量点不在一条直线上,通过测量n个点在空间的坐标就可以确定刚体的空间位姿。 0027 步骤二、记刚体上n个测量点A i 在动坐标系o-xyz下的坐标向量为:a i x Ai y Ai z Ai T (i1,2n)。点A i 在定坐标系o 1 -x 1 y 1 z 1 中的实际测量值记为A i ,以向量a i x Ai y Ai z A。
14、i T 表示。 说 明 书CN 104121877 A 3/4页 5 0028 步骤三、刚体的实际位姿用转换矩阵R和位置向量T表示。那么点A i 的实测 值与通过坐标转换得到的坐标值之间的误差可以表示为: 0029 a i a i -(Ra i +T) (i1,2n) (1) 0030 步骤四、用向量2-范数来评价该误差的大小,可以写成: 0031 a i 2 a i -(Ra i +T) 2 (i1,2n) (2) 0032 步骤五、取向量2-范数的平方和: 0033 0034 步骤六、设测量仪器的测量误差范围为-,+,则式(1)的最大值为: 0035 a i a i -(Ra i +T) 。
15、T (4) 0036 步骤七、每个位姿下的目标函数都对应一个最小值,这些最小值是有一个取值范 围的。这些最小值可能达到的最大值为: 0037 0038 步骤八、式(3)的取值范围是03n 2 ,如果式(3)的值超出此范围,则 认定测量值中至少有一个是粗大误差。为找出该点,依次剔除测量值A j (j1,2 n),然后,利用n-1个测量值计算目标函数,如果只有某次剔除后的n-1个测量值 A 1 ,A 2 .A j-1 ,A j+1 .A n 对应的目标函数值在03(n-1) 2 范围内,而其余 次剔除后n-1个测量值对应的目标函数值超出03(n-1) 2 ,则认定该测量值A j 为 粗大误差,需要。
16、剔除。 0039 步骤九、如果n次剔除计算所对应目标函数均超出范围,则认定数据中有1个以上 的粗大误差,此时利用n-2个测量值计算目标函数,判别方法与只有1个粗大误差的方法相 同。以此类推,直到剔除所有粗大误差。 0040 一个实施例的过程如下: 0041 测量仪器的测量误差是0.005mm,则3n 2 3.7510 -4 ,根据步骤六、步骤 七计算得到,目标函数的取值范围是03.7510 -4 。 0042 目标函数值超出取值范围,从而可判断:5个测量数据中有含粗大误差的点。 0043 利用4个测量点计算目标函数值。以A 1 ,A 2 ,A 3 ,A 4 这4个点对应的测量 数据为例进行说明。
17、,对应的目标函数值为:因为利用4个测 量点计算实际位姿,即n4,所以此时目标函数值的范围是0310 -4 ,以上目标函数 值超出此范围,所得位姿仍为错误结果。同理,利用测量点的其余组合计算实际位姿,相应 的目标函数值如下: 0044 0045 说 明 书CN 104121877 A 4/4页 6 0046 0047 0048 以上5个目标函数值中,只有A 1 ,A 3 ,A 4 ,A 5 对应的目标函数在允许范围 内,而其余包含测量值A 2 的4个目标函数值均超出允许范围,说明A 2 为含粗大误差的 点。 0049 一种测量值的粗大误差判别方法,算法简单、结果可靠,可以查找小数据量范围内 的粗大误差。除此之外,该算法易于用计算机程序实现,只需要把所需数据输入程序,即可 运行程序得到分析结果,应用领域广泛。 0050 以上对本发明创造的一个实施例进行了详细说明,但所述内容仅为本发明创造的 较佳实施例,不能被认为用于限定本发明创造的实施范围。凡依本发明创造申请范围所作 的均等变化与改进等,均应仍归属于本发明创造的专利涵盖范围之内。 说 明 书CN 104121877 A 1/1页 7 图1 说 明 书 附 图CN 104121877 A 。