球的最小二乘拟合.doc

球的最小二乘拟合 以下是我写的程序： #include Windows.h #include math.h #include stdio.h #include stdlib.h #define MAX 10 void Input(double*matrix,int m,int n,double A[4][4])； void Output(double*matrix,int m,int n,double IA[4][4])； void InputExample()； void Inverse(double*matrix1,double*matrix2,int n,double d)； double AlCo(double*matrix,int jie,int row,int column)； double Determinant(double*matrix,int n)； double Cofactor(double*matrix,int jie,int row,int column)； void PrintMTX(double A[4][4]) { for(int i=0；i 4；i++) { printf(\n)； for(int j=0；j 4；j++) printf(%10G,A[i][j])； } printf(\n)； } void main() { double*matrix1[MAX],*matrix2[MAX]； double d； int n； //printf(请输入行列式的阶数：)； //scanf(%d,n)； n=4； //InputExample()； //printf(开始输入矩阵：\n)； int i,j； double Y[4],fac[4],array[12][3],MTX[4][4],IMTX[4][4]； ZeroMemory(MTX,sizeof(MTX))； ZeroMemory(IMTX,sizeof(IMTX))； ZeroMemory(Y,sizeof(Y))； ZeroMemory(fac,sizeof(fac))； //以下的数据是我用SolidWorks画出来的球，我在球上任意取了12个点 array[0][0]=-70.08；array[0][1]=29.27；array[0][2]=6.04； array[1][0]=-72.37；array[1][1]=24.19；array[1][2]=-1.64； array[2][0]=-65.1；array[2][1]=13.6；array[2][2]=0.58； array[3][0]=-60.11；array[3][1]=13.51；array[3][2]=4.93； array[4][0]=-61.12；array[4][1]=27.69；array[4][2]=12.16； array[5][0]=-65；array[5][1]=22.32；array[5][2]=-11.14； array[6][0]=-61.44；array[6][1]=18.31；array[6][2]=-10.46； array[7][0]=-56.37；array[7][1]=15.920；array[7][2]=-7.78； array[8][0]=-53.27；array[8][1]=15.95；array[8][2]=-5.39； array[9][0]=-51.64；array[9][1]=16.53；array[9][2]=-3.82； array[10][0]=-50.99；array[10][1]=16.92；array[10][2]=-3.12； array[11][0]=-58.94；array[11][1]=15.16；array[11][2]=7.64； for(i=0；i 12；i++) { Y[0]+=pow(array[i][0],3)+pow(array[i][1],2)*array[i][0]+pow(array[i][2],2)*array[i][0]； Y[1]+=pow(array[i][1],3)+pow(array[i][0],2)*array[i][1]+pow(array[i][2],2)*array[i][1]； Y[2]+=pow(array[i][2],3)+pow(array[i][0],2)*array[i][2]+pow(array[i][1],2)*array[i][2]； Y[3]-=pow(array[i][0],2)+pow(array[i][1],2)+pow(array[i][2],2)； } for(j=0；j 12；j++) { MTX[0][0]+=pow(arra