最小二乘法正规方程组c语言,编程思路构造正规方程组求解.ppt

编程思路构造正规方程组求解

球的最小二乘拟合 学 院：精密仪器与光电子工程学院 组 员：苏智琨(1012202035) 武志超(1012202041) 霍树春(1012202020) 应用背景 三坐标测量机是目前比较通用的接触式精密测量仪器，其测量精度能够达到微米级别。在进行测量之前一般都要进行探针的校准，以确定探针的实际位置和补偿测端球径与探针挠曲变形误差。实际校准过程一般是测量一个已知的实物标准，如标准球。通过测量点拟合出一个球，求得其直径，与标准球的标称直径进行比较，进而得到探针的补偿值，用于实际的测量。 非线性方程的线性化 设Pi(xi，yi，zi)(i=1，2，3，…，n)为标准球面上n个测量点，设理想球面方程为 其中，(x0，y0，z0)为球心坐标，R为球半径。 非线性方程的线性化 引入如下变换 圆的方程成为如下线性方程 线性化： 这是一个关于(x0，y0，z0)和R的二次非线性方程。当(x0，y0，z0)足够小的时候，可以认为三者的平方项是高阶微量，可以忽略不计。 编程思路 构造正规方程组求解 根据极值原理，欲使F取得最小值，必须满足 构造回归方程求解 编程思路 以x0、y0、z0、C为变量，构造回归方程 result=regress(Y，X) 求解后的处理 当( x0、y0、z0 )不足够小时，需要对测量点进行平移，重新求解，直到结果足够小。 实验结果 用一台三坐标测量机对一个标准球和一个未知球进行测量，三坐标测量软件会给出这两个球的直径，该结果当做真值。用上述两种方法在Matlab中求解，求得的结果与真值进行比较，考查两种方法的可靠性。 标称值 三坐标测量结果 正规方程结果 回归方程结果 15.8752mm 15.8748mm 15.8748mm 15.8748mm 三坐标测量结果 正规方程结果 回归方程结果 39.7759mm 39.776mm 39.776mm 其中h为某一点的高度，I(h)为该点的干涉信号的光强，I0为背景光强，γ是一个经验系数，lc与λ0为与光源相关的常数。 从式(2)可以看出，连续小波变换就是使用不同频率的子小波沿着时间方向扫描信号，从而计算出信号和小波之间的相干系数。当信号和小波的频率和空间特性高度匹配时，会得到一个最大的相干系数。由式(3)可以看到，干涉信号是一个受高斯函数调制的平稳信号。 * 采用了直线-直线-旋转的形式； 将旋转放在最后是考虑到其误差放大作用最大； 测量机通过11-底座固定在地面上，10-固定梁固定在底座上，9-X向运动导轨副使8-运动梁可以带动剩余部分沿X方向运动，7-Z向运动导轨架通过固定在其上的6-Z向运动导轨副使5-Z向运动架板可以带动其余部分沿Z向运动，4-密珠轴系使3-悬臂可以带动2-测头回转体和1-测头绕其回转轴旋转，定义此旋转运动为Y向 测量机的工作原理如下：1、2、3绕4回转180°，测量机即处于收拢状态。加工完毕，机床的玻璃罩打开，1、2、3绕4旋转180°，测量机恢复如图2-4所示状态，测量机进入机床加工区域进行测量。测量时，测量机在计算机的控制下，通过X向、Z向的直线运动，Y向的旋转运动，必要的时候辅以2-测头回转体的运动，带动1-测头探测被测面的表面轮廓。测完一个被测面之后，机床的回转工作台转过一个分度角，进行下一个被测面轮廓的测量，直至完成整个工件的测量。测量完毕后，测量机退出加工区域，恢复到收拢状态，加工或返修继续进行。 * * 其中h为某一点的高度，I(h)为该点的干涉信号的光强，I0为背景光强，γ是一个经验系数，lc与λ0为与光源相关的常数。 从式(2)可以看出，连续小波变换就是使用不同频率的子小波沿着时间方向扫描信号，从而计算出信号和小波之间的相干系数。当信号和小波的频率和空间特性高度匹配时，会得到一个最大的相干系数。由式(3)可以看到，干涉信号是一个受高斯函数调制的平稳信号。 * 采用了直线-直线-旋转的形式； 将旋转放在最后是考虑到其误差放大作用最大； 测量机通过11-底座固定在地面上，10-固定梁固定在底座上，9-X向运动导轨副使8-运动梁可以带动剩余部分沿X方向运动，7-Z向运动导轨架通过固定在其上的6-Z向运动导轨副使5-Z向运动架板可以带动其余部分沿Z向运动，4-密珠轴系使3-悬臂可以带动2-测头回转体和1-测头绕其回转轴旋转，定义此旋转运动为Y向 测量机的工作原理如下：1、2、3绕4回转180°，测量机即处于收拢状态。加工完毕，机床的玻璃罩打开，1、2、3绕4旋转180°，测量机恢复如图2-4所示状态，测量机进入机床加工区域进行测量