欢迎来到天天文库
浏览记录
ID:45035100
大小:84.50 KB
页数:21页
时间:2019-11-08
《C语言空间后方交会源代码》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
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
2、];实用文档//for(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;k3、or(h=k=5;k>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;4、i<6;i++)//提取逆矩阵//for(j=0;j5、0;i6、x[i*dim*2+jj]+=tMatrix[ii*dim*2+jj];}}}base=tMatrix[i*dim*2+i];if(flag==0){printf("求逆矩阵过程中被零除,无法求解!");}//exit(0);}for(intj=0;j7、*2+k];}}for(intk=0;k8、9.41,25273.3
3、or(h=k=5;k>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;
4、i<6;i++)//提取逆矩阵//for(j=0;j5、0;i6、x[i*dim*2+jj]+=tMatrix[ii*dim*2+jj];}}}base=tMatrix[i*dim*2+i];if(flag==0){printf("求逆矩阵过程中被零除,无法求解!");}//exit(0);}for(intj=0;j7、*2+k];}}for(intk=0;k8、9.41,25273.3
5、0;i6、x[i*dim*2+jj]+=tMatrix[ii*dim*2+jj];}}}base=tMatrix[i*dim*2+i];if(flag==0){printf("求逆矩阵过程中被零除,无法求解!");}//exit(0);}for(intj=0;j7、*2+k];}}for(intk=0;k8、9.41,25273.3
6、x[i*dim*2+jj]+=tMatrix[ii*dim*2+jj];}}}base=tMatrix[i*dim*2+i];if(flag==0){printf("求逆矩阵过程中被零除,无法求解!");}//exit(0);}for(intj=0;j7、*2+k];}}for(intk=0;k8、9.41,25273.3
7、*2+k];}}for(intk=0;k8、9.41,25273.3
8、9.41,25273.3
此文档下载收益归作者所有