欢迎来到天天文库
浏览记录
ID:48623571
大小:47.00 KB
页数:11页
时间:2020-01-30
《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;j5、=dim;j6、*dim*2+jj];}}}base=tMatrix[i*dim*2+i];if(flag==0){printf("求逆矩阵过程中被零除,无法求解!");}//exit(0);}for(intj=0;j7、{tMatrix[i*dim*2+k]/=base;}}for(i=0;i
5、=dim;j6、*dim*2+jj];}}}base=tMatrix[i*dim*2+i];if(flag==0){printf("求逆矩阵过程中被零除,无法求解!");}//exit(0);}for(intj=0;j7、{tMatrix[i*dim*2+k]/=base;}}for(i=0;i
6、*dim*2+jj];}}}base=tMatrix[i*dim*2+i];if(flag==0){printf("求逆矩阵过程中被零除,无法求解!");}//exit(0);}for(intj=0;j7、{tMatrix[i*dim*2+k]/=base;}}for(i=0;i
7、{tMatrix[i*dim*2+k]/=base;}}for(i=0;i
此文档下载收益归作者所有