C语言间接平差程序

1、教材误差理论与测量平差基础第二版武汉大学出版社P108页的例7-1的运行结果：Q匚;Dvcumentsand5ettings.AdministratorS.Debug9.eKe11-hl幘输入U的系数BWIX100kl1001001-1001晴输入观测值的权阵pHH：E.90000037000002.5000003.3000004睛输入常SctHlEll：11-752944-2.035289改正数U和丄=11-7529449.211767艮庄误差“二11-752944-2.035289改正数U和丄=11-7529449.211767艮庄误差“二24.391144协大多女阵QxxtHETJ:0-

2、2044780.0944740-094474S.1G85210-042707Q.Q7G1810-23014e整数xTjri=PressanjItejtocontinue-7.248829-2.035289-8-786460-7.2488290.0427079.97&1S10.171424源程序：#defineN5/*N是观测值个数*/#defineT3/*T是必要观测数*/#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的系数BNT:n);for(i=0;iN;i+)for(j=0;jT;j+)scanf(%8f,&Bij);printf(请输入观测值的权阵PNN:n);for(i=0;iN;i+)for(j=0;jN;j+)scanf(%8f,&Pij);printf(请输入常数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(参数xT1=n);for(i=0;iT;i+)printf(%15f,xi0);printf(n);printf(改正数VN1=n);for(i=0;iN;i+)printf(%15f,Vi0);printf(n单位权中误差x1=%15f,sqrt(x1);printf(n协因数阵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