C语言间接平差程序(共6页).doc
精选优质文档-倾情为你奉上 教材误差理论与测量平差基础第二版武汉大学出版社P108页的例7-1的运行结果： 源程序： #define N 5 /*N是观测值个数*/#define T 3 /*T是必要观测数*/#include<stdio.h>#include<math.h>float NbbTT,NbTT,WT1,xT1;main() float D(float aTN,float bNN,float cNT); float K(float aTN,float bNN,float cN1); float G(float aTT); float F(float caT-1T-1); float DM(float a1N,float bNN ,float cN1); int i,j,m,n; float BNT,BTTN,VN1,VT1N,PNN,CN1,BxN1,f,g,h,x1; printf("请输入V的系数BNT:n"); for(i=0;i<N;i+) for(j=0;j<T;j+) scanf("%8f",&Bij); printf("请输入观测值的权阵PNN:n"); for(i=0;i<N;i+) for(j=0;j<N;j+) scanf("%8f",&Pij); printf("请输入常数CN1:n"); for(i=0;i<N;i+) for(j=0;j<1;j+) scanf("%8f",&Cij); for(i=0;i<N;i+) for(j=0;j<T;j+) BTji=Bij; g=D(BT, P, B); h=K(BT, P, C); f=G(Nbb); for(i=0;i<T;i+) for(j=0;j<1;j+) xij=Nbi0*W0j; for(m=1;m<T;m+) xij+=(Nbim*Wmj); for(i=0;i<T;i+) xi0=xi0/f; for(i=0;i<N;i+) for(j=0;j<1;j+) Bxij=Bi0*x0j; for(m=1;m<T;m+) Bxij+=(Bim*xmj); for(i=0;i<N;i+) Vi0=(Bxi0-Ci0); for(i=0;i<N;i+) for(j=0;j<1;j+) VTji=Vij; x1=DM(VT,P,V); x1=x1/(N-T); printf("参数xT1=n"); for(i=0;i<T;i+) printf("%15f",xi0); printf("n"); printf("改正数VN1=n"); for(i=0;i<N;i+) printf("%15f",Vi0); printf("n单位权中误差x1=%15f",sqrt(x1); printf("n协因数阵QxxTT:n"); for(i=0;i<T;i+) for(j=0;j<T;j+) printf("%15f",Nbij/f); printf("n"); float G(float aTT) int i,j,m,n; float cT-1T-1,y=0; for(i=0;i<T;i+) for(j=0;j<T;j+) for(m=0;m<T;m+) for(n=0;n<T;n+) if(m<i&&n<j) cmn=amn; if(m>i&&n<j) cm-1n=amn; if(m<i&&n>j) cmn-1=amn; if(m>i&&n>j) cm-1n-1=amn; if(i+j)%2=0) Nbji=F(c); else Nbji=(-1)*F(c); for(m=0;m<T;m+) y+=(a0m*Nbm0); return (y); float F(float caT-1T-1) int i,j,m,n,s,t,k=1; float f=1,c,x,sn; for (i=0,j=0;i<T-1&&j<T-1;i+,j+) if (caij=0) for (m=i;camj=0;m+); if (m=T-1) sn=0; return (sn); else for (n=j;n<T-1;n+) c=cain; cain=camn; camn=c; k*=(-1); for (s=T-2;s>i;s-) x=casj; for (t=j;t<T-1;t+) cast-=cait*(x/caij); for (i=0;i<T-1;i+) f*=caii; sn=k*f; return (sn); float D(float aTN,float bNN ,float cNT) int i,j,m; float dTN; for(i=0;i<T;i+) for(j=0;j<N;j+) dij=ai0*b0j; for(m=1;m<N;m+) dij+=(aim*bmj); for(i=0;i<T;i+) for(j=0;j<T;j+) Nbbij=di0*c0j; for(m=1;m<N;m+) Nbbij+=(dim*cmj); return (Nbb00); float K(float aTN,float bNN,float cN1) int i,j,m; float dTN; for(i=0;i<T;i+) for(j=0;j<N;j+) dij=ai0*b0j; for(m=1;m<N;m+) dij+=(aim*bmj); for(i=0;i<T;i+) for(j=0;j<1;j+) Wij=di0*c0j; for(m=1;m<N;m+) Wij+=(dim*cmj); return (W00); float DM(float a1N,float bNN ,float cN1) int i,j,m; float d1N,x; for(i=0;i<1;i+) for(j=0;j<N;j+) dij=ai0*b0j; for(m=1;m<N;m+) dij+=(aim*bmj); for(i=0;i<1;i+) for(j=0;j<1;j+) x=di0*c0j; for(m=1;m<N;m+) x+=(dim*cmj); return (x); 程序说明： 1) 用该程序前，根据具体情况输入N和T； 2) 该程序中的(N-T)自由度(即多余观测数)必须大于等于2，不然程序运行时会出错，原因在于求行列式的逆时，有语句for (s=R-2;s>i;s-)，R=1时s=-1。专心-专注-专业
