1、教材誤差理論與測量平差基礎(chǔ)第二版武漢大學(xué)出版社P108頁的例7-1的運行結(jié)果：Q匚;Dvcumentsand5ettings.AdministratorS.Debug9.eKe11-hl幘輸入U的系數(shù)BWIX100kl1001001-1001晴輸入觀測值的權(quán)陣pHH：E.90000037000002.5000003.3000004睛輸入常SctHlEll：11-752944-2.035289改正數(shù)U和丄=11-7529449.211767艮莊誤差“二11-752944-2.035289改正數(shù)U和丄=11-7529449.211767艮莊誤差“二24.391144協(xié)大多女陣QxxtHETJ:0-

2、2044780.0944740-094474S.1G85210-042707Q.Q7G1810-23014e整數(shù)xTjri=PressanjItejtocontinue-7.248829-2.035289-8-786460-7.2488290.0427079.97&1S10.171424源程序：#defineN5/*N是觀測值個數(shù)*/#defineT3/*T是必要觀測數(shù)*/#includevstdio.h#includevmath.hfloatNbbTT,NbTT,WT1,xT1;main()floatD(floataTN,floatbNN,floatcNT);floatK(floataTN,

3、floatbNN,floatcN1);floatG(floataTT);floatF(floatcaT-1T-1);floatDM(floata1N,floatbNN,floatcN1);inti,j,m,n;floatBNT,BTTN,VN1,VT1N,PNN,CN1,BxN1,f,g,h,x1;printf(”請輸入V的系數(shù)BNT:n);for(i=0;iN;i+)for(j=0;jT;j+)scanf(%8f,&Bij);printf(請輸入觀測值的權(quán)陣PNN:n);for(i=0;iN;i+)for(j=0;jN;j+)scanf(%8f,&Pij);printf(請輸入常數(shù)CNl:n

4、);for(i=0;iN;i+)for(j=0;jl;j+)scanf(%8f,&Cij);for(i=0;iN;i+)for(j=0;jT;j+)BTji=Bij;g=D(BT,P,B);h=K(BT,P,C);f=G(Nbb);for(i=0;iT;i+)for(j=0;jl;j+)xij=Nbi0*W0j;for(m=l;mT;m+)xij+=(Nbim*Wmj);for(i=0;iT;i+)xi0=xi0/f;for(i=0;iN;i+)for(j=0;jl;j+)Bxij=Bi0*x0j;for(m=l;mT;m+)Bxij+=(Bim*xmj);for(i=0;iN;i+)Vi0=

5、(Bxi0-Ci0);for(i=0;iN;i+)for(j=0;jl;j+)VTji=Vij;xl=DM(VT,P,V);xl=xl/(N-T);printf(參數(shù)xT1=n);for(i=0;iT;i+)printf(%15f,xi0);printf(n);printf(改正數(shù)VN1=n);for(i=0;iN;i+)printf(%15f,Vi0);printf(n單位權(quán)中誤差x1=%15f,sqrt(x1);printf(n協(xié)因數(shù)陣QxxTT:n);for(i=0;iT;i+)for(j=0;jT;j+)printf(%15f,Nbij/f);printf(n);floatG(floa

6、taTT)inti,j,m,n;floatcT-1T-1,y=0;for(i=0;iT;i+)for(j=0;jT;j+)for(m=0;mT;m+)for(n=0;nT;n+)if(mi&ni&nj)cm-1n=amn;if(mj)cmn-1=amn;if(mi&nj)cm-1n-1=amn;if(i+j)%2=0)Nbji=F(c);elseNbji=(-1)*F(c);for(m=0;mT;m+)y+=(a0m*Nbm0);return(y);floatF(floatcaT-1T-1)inti,j,m,n,s,t,k=1;floatf=1,c,x,sn;for(i=0,j=0;iT-1&

7、jT-1;i+,j+)if(caij=0)for(m=i;camj=0;m+);if(m=T-1)sn=0;return(sn);elsefor(n=j;ni;s-)x=casj;for(t=j;tT-1;t+)cast-=cait*(x/caij);for(i=0;iT-1;i+)f*=caii;sn=k*f;return(sn);floatD(floataTN,floatbNN,floatcNT)inti,j,m;floatdTN;for(i=0;iT;i+)for(j=0;jN;j+)dij=ai0*b0j;for(m=1;mN;m+)dij+=(aim*bmj);for(i=0;iT;

8、i+)for(j=0;jT;j+)Nbbij=di0*c0j;for(m=1;mN;m+)Nbbij+=(dim*cmj);return(Nbb00);floatK(floataTN,floatbNN,floatcN1)inti,j,m;floatdTN;for(i=0;iT;i+)for(j=0;jN;j+)dij=ai0*b0j;for(m=1;mN;m+)dij+=(aim*bmj);for(i=0;iT;i+)for(j=0;j1;j+)Wij=di0*c0j;for(m=1;mN;m+)Wij+=(dim*cmj);return(W00);floatDM(floata1N,floatbNN,floatcN1)inti,j,m;floatd1N,x;for