球的最小二乘拟合.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 A44)；void Output(double*matrix,int m,int n,double IA44)；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 A44)for(int i=0；i 4；i+)printf(n)；for(int j=0；j 4；j+)printf(%10G,Aij)；printf(n)；void main()double*matrix1MAX,*matrix2MAX；double d；int n；/printf(请输入行列式的阶数：)；/scanf(%d,&n)；n=4；/InputExample()；/printf(开始输入矩阵：n)；int i,j；double Y4,fac4,array123,MTX44,IMTX44；ZeroMemory(MTX,sizeof(MTX)；ZeroMemory(IMTX,sizeof(IMTX)；ZeroMemory(Y,sizeof(Y)；ZeroMemory(fac,sizeof(fac)；/以下的数据是我用SolidWorks画出来的球，我在球上任意取了12个点array00=-70.08；array01=29.27；array02=6.04；array10=-72.37；array11=24.19；array12=-1.64；array20=-65.1；array21=13.6；array22=0.58；array30=-60.11；array31=13.51；array32=4.93；array40=-61.12；array41=27.69；array42=12.16；array50=-65；array51=22.32；array52=-11.14；array60=-61.44；array61=18.31；array62=-10.46；array70=-56.37；array71=15.920；array72=-7.78；array80=-53.27；array81=15.95；array82=-5.39；array90=-51.64；array91=16.53；array92=-3.82；array100=-50.99；array101=16.92；array102=-3.12；array110=-58.94；array111=15.16；array112=7.64；for(i=0；i 12；i+)Y0+=pow(arrayi0,3)+pow(arrayi1,2)*arrayi0+pow(arrayi2,2)*arrayi0；Y1+=pow(arrayi1,3)+pow(arrayi0,2)*arrayi1+pow(arrayi2,2)*arrayi1；Y2+=pow(arrayi2,3)+pow(arrayi0,2)*arrayi2+pow(arrayi1,2)*arrayi2；Y3-=pow(arrayi0,2)+pow(arrayi1,2)+pow(arrayi2,2)；for(j=0；j 12；j+)MTX00+=pow(arrayj0,2)；MTX01+=arrayj0*arrayj1；MTX02+=arrayj0*arrayj2；MTX03-=arrayj0；MTX11+=pow(arrayj1,2)；MTX12+=arrayj1*arrayj2；MTX13-=arrayj1；MTX22+=pow(arrayj2,2)；MTX23-=arrayj2；MTX10=MTX01；MTX20=MTX02；MTX21=MTX12；MTX30=MTX03；MTX31=MTX13；MTX32=MTX23；MTX33=12.0；Input(matrix1,n,n,MTX)；PrintMTX(MTX)；printf(nY0=%f,Y0)；printf(nY1=%f,Y1)；printf(nY2=%f,Y2)；printf(nY3=%fn,Y3)；d=Determinant(matrix1,n)；if(d=0)printf(这个矩阵不可逆。n)；elseInverse(matrix1,matrix2,n,d)；printf(它的逆矩阵是：n)；Output(matrix2,n,n,IMTX)；for(i=0；i 4；i+)printf(n)；for(int j=0；j 4；j+)double tt=0；tt=MTXi0*IMTX0j+MTXi1*IMTX1j+MTXi2*IMTX2j+MTXi3*IMTX3j；printf(%5.2ft,tt)；printf(n)；printf(n解算出来的系数是：)；for(i=0；i 4；i+)fac0+=IMTX0i*Yi；fac1+=IMTX1i*Yi；fac2+=IMTX2i*Yi；fac3+=IMTX3i*Yi；/printf(nfac%d=%f,i,faci)；printf(na=%f,fac0/2.0)；printf(nb=%f,fac1/2.0)；printf(nc=%f,fac2/2.0)；/printf(nfac3=%f,fac3/2.0)；/for(i=0；i 4；i+)/printf(nfac%d=%f,i,faci)；/printf(n)；double R；R=pow(fac0/2.0,2)+pow(fac1/2.0,2)+pow(fac2/2.0,2)-fac3；printf(nR=%f,pow(R,0.5)；Sleep(1)；getchar()；/return 0；void Input(double*matrix,int m,int n,double A44)int i,j；for(i=0；i m；i+)matrixi=(double*)malloc(n*sizeof(double)；/printf(Please input the matrix：n)；for(i=0；i m；i+)for(j=0；j n；j+)/double tt=；*(matrixi+j)=Aij；void Output(double*matrix,int m,int n,double IA44)int i,j；for(i=0；i m；i+)for(j=0；j n；j+)IAij=*(matrixi+j)；printf(%10Gt,IAij)；printf(n)；void InputExample()printf(数据间用Tab键间隔，输完一行按回车。例如：n)；printf(12 26 48n)；printf(-1 24 10n)；printf(2-6 14n)；void Inverse(double*matrix1,double*matrix2,int n,double d)int i,j；for(i=0；i n；i+)matrix2i=(double*)malloc(n*sizeof(double)；for(i=0；i n；i+)for(j=0；j n；j+)*(matrix2j+i)=(AlCo(matrix1,n,i,j)/d)；double Determinant(double*matrix,int n)double result=0,temp；int i；if(n=1)result=(*matrix0)；elsefor(i=0；i n；i+)temp=AlCo(matrix,n,n-1,i)；result+=(*(matrixn-1+i)*temp；return result；double AlCo(double*matrix,int jie,int row,int column)double result；if(row+column)%2=0)result=Cofactor(matrix,jie,row,column)；else result=(-1)*Cofactor(matrix,jie,row,column)；return result；double Cofactor(double*matrix,int jie,int row,int column)double result；int i,j；double*smallmatrMAX-1；for(i=0；i jie-1；i+)smallmatri=(double*)malloc(sizeof(double)*(jie-1)；for(i=0；i row；i+)for(j=0；j column；j+)*(smallmatri+j)=*(matrixi+j)；for(i=row；i jie-1；i+)for(j=0；j column；j+)*(smallmatri+j)=*(matrixi+1+j)；for(i=0；i row；i+)for(j=column；j jie-1；j+)*(smallmatri+j)=*(matrixi+j+1)；for(i=