誤差理論與測量平差課程設計--水準網(wǎng)嚴密平差及精度評定(編輯修改稿)
2025-02-09 18:33 本頁面
【文章內(nèi)容簡介】 uble a[3][7]=BTP；static double c[3][1],b[7][1]=l；Matrixmul(a,b,3,7,1,c)。printf(“\n”)。for(i=0。i=2。i++){for(j=0。j=0。j++) printf(“%\t”,c[i][j]。 printf(“\n”)。}printf(“\n”)。return0。}二、矩陣的逆計算函數(shù)（求NBB1）include define M 3 void main() { float MAT[M][2*M]。 float MAT1[M][M]。 float t。 int i,j,k,l。 /***********************************************/ /*對矩陣進行初始化*/ for(i=0。iM。i++) for(j=0。j2*M。j++) MAT1[j]=39。\039。 /*對MAT1矩陣賦初值 */ for(i=0。iM。i++) for (j=0。jM。j++) scanf(%f,amp。MAT1[j])。 /*打印目標矩陣?*/ printf(原矩陣為：\n)。 for (i=0。iM。i++) { for (j=0。jM。j++) printf(%,MAT1[j])。 printf(\n)。 } /********************************************/ /*對MAT1矩陣進行擴展,MAT1矩陣添加單位陣，由M*M變成2M*2M矩陣 */ for(i=0。iM。i++) for(j=0。j2*M。j++) if (jM) MAT[j]=MAT1[j]。 else if (j==M+i) MAT[j]=1。 else MAT[j]=0。 /*對M矩陣進行變換，使得前半部分矩陣成為單位陣，則 */ /*后半部分矩陣即為所求矩陣逆陣 */ for(i=0。iM。i++) { /*對第i行進行歸一化 */ for (j=0。j2*M。j++) for(k=i+1。kM。k++) MAT[j]=MAT[j]+MAT[k][j]。 t=MAT。 for(j=i。j2*M。j++) MAT[j]=MAT[j]/t。 /*對矩陣進行行變換，使得第i 列只有一個元素不為零，且為1*/ for(k=0。kM。k++) if(k!=i) { t=MAT[k]。 for (l=i。l2*M。l++) MAT[k][l]=MAT[k][l]MAT[l]*t。 } } /*將后半部分矩陣即所求矩陣逆陣存入MAT2矩陣。*/ for(i=0。iM。i++) { for(j=0。jM。j++) MAT1[j]=MAT[j+M]。 printf(\n)。 } /*********************************************/ /*輸出所求的逆陣*/ printf(逆陣為：\n)。 for(i=0。iM。i++) { for(j=0。jM。j++) printf(%,MAT1[j])。 printf(\n)。 } }，即改正數(shù)xmain(){int i,j。
公司管理相關推薦
鄂ICP備17016276號-1