C语言空间后方交会源代码.doc

《C语言空间后方交会源代码.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.#include#include#definen4//控制点个数#definePI3.14159265structcoordinate{doublex;//像点坐标doubley;doubleXt;//控制点坐标doubleYt;doubleZt;};//voidinverse(doublec[6][6])//矩阵求逆//{//inti,j,h,k;//doublep;//doubleq[6][12];//for(i=0;i<6;i++)//构造高斯矩阵//for(j=0;j<6;j++)//q[i][j]=c[i][j];//f

2、or(i=0;i<6;i++)//for(j=6;j<12;j++)//{//if(i+6==j)//q[i][j]=1;//else//q[i][j]=0;//}//for(h=k=0;k

3、0;k--,h--)//消去对角线以上的数据..//for(i=k-1;i>=0;i--)//{//if(q[i][h]==0)//continue;//p=q[k][h]/q[i][h];////p=q[i][h]/q[k][h];//for(j=11;j>0;j--)//{//q[i][j]*=p;//q[i][j]-=q[k][j];//}//}//for(i=0;i<6;i++)//将对角线上数据化为1//{//p=1.0/q[i][i];//for(j=0;j<12;j++)//q[i][j]*=p;//}//for(i=0;i<6;i++)//提取逆矩阵

4、//for(j=0;j

5、=dim;j

6、*dim*2+jj];}}}base=tMatrix[i*dim*2+i];if(flag==0){printf("求逆矩阵过程中被零除，无法求解!");}//exit(0);}for(intj=0;j

7、{tMatrix[i*dim*2+k]/=base;}}for(i=0;i