C语言间接平差程序(共6页).doc

《C语言间接平差程序(共6页).doc》由会员分享，可在线阅读，更多相关《C语言间接平差程序(共6页).doc（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上 教材误差理论与测量平差基础第二版武汉大学出版社P108页的例7-1的运行结果： 源程序： #define N 5 /*N是观测值个数*/#define T 3 /*T是必要观测数*/#include#includefloat 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,flo

2、at 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;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(请输入常数CN1:n); for(i=0;iN;i+) for(j=0;j1;j+) scanf(%8f,&Cij); for(i=0;iN;i+

3、) 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;j1;j+) xij=Nbi0*W0j; for(m=1;mT;m+) xij+=(Nbim*Wmj); for(i=0;iT;i+) xi0=xi0/f; for(i=0;iN;i+) for(j=0;j1;j+) Bxij=Bi0*x0j; for(m=1;mT;m+) Bxij+=(Bim*xmj); for(i=0;iN;i+) Vi0=(Bxi0-Ci0); for(i=0;iN;i+) for(j=0;j

4、1;j+) VTji=Vij; x1=DM(VT,P,V); x1=x1/(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); float G(float aTT) int i,j,m

5、,n; float cT-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); else Nbji=(-1)*F(c); for(m=0;mT;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; fo

6、r (i=0,j=0;iT-1&jT-1;i+,j+) if (caij=0) for (m=i;camj=0;m+); if (m=T-1) sn=0; return (sn); else for (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); float D(float aTN,float bNN ,float cNT) int i,j,m; float dTN; for(i=0;iT;i+) for(j=0;jN;j+) d

7、ij=ai0*b0j; for(m=1;mN;m+) dij+=(aim*bmj); for(i=0;iT;i+) for(j=0;jT;j+) Nbbij=di0*c0j; for(m=1;mN;m+) Nbbij+=(dim*cmj); return (Nbb00); float K(float aTN,float bNN,float cN1) int i,j,m; float dTN; 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); float DM(float a1N,float bNN ,float cN1) int i,j,m; float d1N,x; for(i=0;i1;i+) for(j=0;jN;j+) dij=ai0*b0j; for(m=1;mN;m+) dij+=(aim*bmj); for(i=0;i1;i+) for(j=0;j1;j+) x=di0*c0j; for(m=1;mi;s-)，R=1时s=-1。专心-专注-专业