观测圆上多个任意点及其圆心坐标拟合.docxVIP

观测圆上多个任意点及其圆心坐标拟合张应芳，张爱国，苏积聪，黄海明，韦霄（厦门理工学院计算机与信息工程学院，福建厦门３６１０２４）摘要：基于最小二乘法的原理，推算圆上观测若干坐标点拟合圆心坐标的公式及其圆心坐标精度评定，并通过实地测量标准圆形构件验证圆心坐标拟合计算公式的正确性。同时，运用ｍａｔｌａｂ编程语言实现圆心坐标拟合计算程序，在大型机械设备的制造与安装等方面具有一定的应用价值。关键词：最小二乘法；圆心坐标拟合；精度评定；间接平差中图分类号：Ｐ２０７文献标志码：Ａ文章编号：１６７１－４６７９（２０１４）０２－００２７－０５ＯｂｓｅｒｖａｔｉｏｎｏｆｍｕｌｔｉｐｌｅａｒｂｉｔｒａｒｙｐｏｉｎｔｓｏｎｔｈｅｃｉｒｃｌｅａｎｄｉｔｓｃｅｎｔｅｒｃｏｏｒｄｉｎａｔｅｓｆｉｔｔｉｎｇＺＨＡＮＧＹｉｎｇ－ｆａｎｇ，ＺＨＡＮＧＡｉ－ｇｕｏ，ＳＵＪｉ－ｃｏｎｇ，ＨＵＡＮＧＨａｉ－ｍｉｎｇ，ＷＥＩＸｉａｏ（ＣｏｌｌｅｇｅｏｆＣｏｍｐｕｔｅｒａｎｄＩｎｆｏｒｍａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇ，ＸｉａｍｅｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｘｉａｍｅｎ３６１０２４Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｂａｓｅｄｏｎｔｈｅｐｒｉｎｃｉｐｌｅｏｆｌｅａｓｔ－ｓｑｕａｒｅｓｍｅｔｈｏｄ，ｉｔｄｅｒｉｖｅｓｔｈｅｆｏｒｍｕｌａｏｆｔｈｅｃｏｏｒｄｉｎａｔｅｏｆｔｈｅｃｉｒｃｌｅｃｅｎｔｅｒａｎｄｅｖａｌｕａｔｅｓｉｔｓａｃｃｕｒａｃｙ．Ｔｈｅｎｔｈｅｃｏｒｒｅｃｔｎｅｓｓｏｆｔｈｅｃｏｏｒｄｉｎａｔｅｏｆｔｈｅｃｉｒｃｌｅｃｅｎｔｅｒｉｓｖｅｒｉｆｉｅｄｂｙｔｈｅｒｅａｌｍｅａｓｕｒｅｍｅｎｔｏｆｃｉｒｃｌｅｃｅｎｔｅｒ．Ｍｅａｎｗｈｉｌｅ，ｕｓｉｎｇｔｈｅｍａｔｌａｂｐｒｏｇｒａｍｍｉｎｇｌａｎｇｕａｇｅ，ｔｈｅｃｏｏｒｄｉｎａｔｅｏｆｔｈｅｃｉｒｃｌｅｃｅｎｔｅｒｆｉｔｔｉｎｇｃａｌｃｕｌａｔｉｏｎｐｒｏｇｒａｍｉｓｉｍｐｌｅｍｅｎｔｅｄ，ｗｈｉｃｈｃａｎｂｅｕｓｅｄｉｎｍａｎｕｆａｃｔｕｒｉｎｇａｎｄｉｎｓｔａｌｌａｔｉｏｎｏｆｌａｒｇｅｍｅｃｈａｎｉｃａｌｅｑｕｉｐｍｅｎｔ．Ｋｅｙｗｏｒｄｓ：ｌｅａｓｔ－ｓｑｕａｒｅｓｍｅｔｈｏｄ；ｔｈｅｃｏｏｒｄｉｎａｔｅｏｆｔｈｅｃｉｒｃｌｅｃｅｎｔｅｒｆｉｔｔｉｎｇ；ａｃｃｕｒａｃｙｅｖａｌｕａｔｉｏｎ；ａｄｊｕｓｔｍｅｎｔｏｆｏｂｓｅｒｖａｔｉｏｎｅｑｕａｔｉｏｎｓ在当前机械企业的大型圆形构件检测问题上，各种圆形大构件的中心精度对于构件能否成功安装以及正常工作具有至关重要的作用。本文提出通过全站仪或ＧＰＳ等测量仪器获得圆形大构件上的若干坐标点，利用最小二乘原理推算其圆心坐标计算公式，并将计算到的圆心坐标与观测到的圆心坐标进行一致性检测，从而检定圆形构件的制作工艺是否合格。同时，通过ｍａｔｌａｂ编程开发圆心三维坐标拟合软件程序，加强原理公式的实际应用，提升其在大型机械设备的制造与安装等方面的应用价值。单个圆的圆心坐标拟合解决思路１．１二维情况下单个圆圆心坐标拟合１．１．１设定圆形大构件的圆心坐标为（ｘ０，ｙ０），半径为Ｒ。圆形构件上均匀分布了ｎ个二维坐标点，且坐标为（ｘｉ，ｙｉ），ｉ＝１，２，…，ｎ，推算公式步骤如下所示。最小二乘拟合圆曲线［１］的基础方程为Ｒ２＝（ｘ－ｘ０）２＋（ｙ－ｙ０）２．展开可得（１）Ｒ２＝ｘ２－２ｘｘ０＋ｘ２２２（）０＋ｙ－２ｙｙ０＋ｙ０．２令ａ＝２ｘ０，ｂ＝２ｙ０，ｃ＝ｘ２＋ｙ２－Ｒ２，即ｘ０＝００ａ，ｙ０＝ｂ，Ｒ＝１槡２＋ｂ２－４ｃ。圆心坐标拟合单个圆的圆心坐标拟合又分二维和三维情况讨论。ａ１２２２则式（２）变换之后有ｘ２＋ｙ２－ａｘ－ｂｙ＋ｃ＝０．（３）将该式推广到ｎ个点后，写成矩阵形式可得收稿日期：２０１３－１０－１７作者简介：张应芳（１９９２－），女，本科生，研究方向：测绘工程．黑龙江工程学院学报第２８卷·２８·２２熿ｘ１熿ｘ１＋ｙ１燄熿１燄ｙ１ｙｎ－１燄熿ａ燄ｘｙ１ｙｎｚＡ熿１１燄熿燄（４）１ｂ＝．（９）Ｂ－＝０．燀ｘｎ－１燅燀ｃ燅燀ｘ２＋ｙ２燅ｎｎ燀ｘｎｚｎ燅燀Ｃ燅燀１燅Ｔ熿ｘ１ｙ１－１燄设定在式（４）两边同乘，由矩阵变熿ｘ１ｚ１燄ｙ１ｙｎ燀ｘｎｙｎ－１燅，Ｍ＝换［２］可得燀ｘｎｚｎ燅２熿∑ｘｉ∑ｘｉｙｉ－∑ｘｉ燄熿ａ燄熿１燄熿Ａ燄２∑ｘｉｙｉｂ＝∑ｙｉ－∑ｙｉ１，ｌ＝Ｘ０＝Ｂ．燀ｃ燅燀－∑ｘｉｎ－∑ｙｉ燅３×３燀Ｃ燅燀１燅∑ｘｉ（ｘ２＋ｙ２ｉｉ）燄熿与二维情况下单个圆圆心坐标拟合相同，经矩阵变∑ｙｉ（ｘ２＋ｙ２ｉｉ）．（５）换可解算出Ａ，Ｂ，Ｃ的值，即可得到圆形大构件表面所在的平面方程。燀－∑（ｘ２＋ｙ２）燅ｉｉＡ则熿燄Ｂ＝熿ａ燄燀Ｃ燅ｂ＝－１２熿∑ｘｉ∑ｘｉｙｉ－∑ｘｉ燄熿∑ｘｉ燄燀ｃ燅－１２∑ｘｉ（ｘ２＋ｙ２ｉｉ）燄熿∑ｘｉ∑