C语言间接平差程序.docVIP

下载本文档

37
0
约3.94千字
约 6页
2018-05-02 发布于河南
举报
版权申诉

C语言间接平差程序.doc

1、本文档共6页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

C语言间接平差程序

教材《误差理论与测量平差基础》第二版武汉大学出版社 P108页的例7-1的运行结果：源程序： #define N 5 /*N是观测值个数*/ #define T 3 /*T是必要观测数*/ #includestdio.h #includemath.h float Nbb[T][T],Nb[T][T],W[T][1],x[T][1]; main() { float D(float a[T][N],float b[N][N],float c[N][T]); float K(float a[T][N],float b[N][N],float c[N][1]); float G(float a[T][T]); float F(float ca[T-1][T-1]); float DM(float a[1][N],float b[N][N] ,float c[N][1]); int i,j,m,n; float B[N][T],BT[T][N],V[N][1],VT[1][N],P[N][N],C[N][1],Bx[N][1],f,g,h,x1; printf(请输入V的系数B[N][T]:\n); for(i=0;iN;i++) for(j=0;jT;j++) scanf(%8f,B[i][j]); printf(请输入观测值的权阵P[N][N]:\n); for(i=0;iN;i++) for(j=0;jN;j++) scanf(%8f,P[i][j]); printf(请输入常数C[N][1]:\n); for(i=0;iN;i++) for(j=0;j1;j++) scanf(%8f,C[i][j]); for(i=0;iN;i++) for(j=0;jT;j++) BT[j][i]=B[i][j]; g=D(BT, P, B); h=K(BT, P, C); f=G(Nbb); for(i=0;iT;i++) for(j=0;j1;j++) { x[i][j]=Nb[i][0]*W[0][j]; for(m=1;mT;m++) x[i][j]+=(Nb[i][m]*W[m][j]); } for(i=0;iT;i++) x[i][0]=x[i][0]/f; for(i=0;iN;i++) for(j=0;j1;j++) { Bx[i][j]=B[i][0]*x[0][j]; for(m=1;mT;m++) Bx[i][j]+=(B[i][m]*x[m][j]); } for(i=0;iN;i++) V[i][0]=(Bx[i][0]-C[i][0]); for(i=0;iN;i++) for(j=0;j1;j++) VT[j][i]=V[i][j]; x1=DM(VT,P,V); x1=x1/(N-T); printf(参数x[T][1]=\n); for(i=0;iT;i++) printf(%15f,x[i][0]); printf(\n); printf(改正数V[N][1]=\n); for(i=0;iN;i++) printf(%15f,V[i][0]); printf(\n单位权中误差x1=%15f,sqrt(x1)); printf(\n协因数阵Qxx[T][T]:\n); for(i=0;iT;i++) { for(j=0;jT;j++) printf(%15f,Nb[i][j]/f); printf(\n); } } float G(float a[T][T]) { int i,j,m,n; float c[T-1][T-1],y=0; for(i=0;iT;i++) for(j=0;jT;j++) { for(m=0;mT;m++) for(n=0;nT;n++) { if(minj) c[m][n]=a[m][n]; if(minj) c[m-1][n]=a[m][n]; if(min