导线网平差程序源代码

程序源代码c语言程序：#include#include#include#include#include#define PI 3.1415926535898# define p2206264.806247#define MAX 50/矩阵乘法matrix mutily(double matrix 1maxmax、double matrix 2 max max、double matrix result max)Int i、j、k；双精度；for(I=0)；Id)d=fab(aIj)；ISk=I；JSk=j；If (d 1.0=1.0)return 0；If (ISk)！=k)for(j=0)；j=0；kfor(j=0)；J=i)BreakXO0=5000.0；约0=5000.0；XO15=5000.0；约15=5000.0；so0=s0；/用于查找未知点坐标近似值for(I=1)；I15I)alfI=alfI-1betaI-pi；If(alfi=(2*PI)alfI=alfI-2 * PI；XOI=XOI-1sI-1* cos(alfI-1)；约I=约I-1 s I-1 *罪(alfI-1)；/printf(X=%ft ，XoI)；for(I=1)；I15I)/寻找近似距离soI=sqrt(yoI 1-yoI)*(yoI 1-yoI)(XOI 1-XOI/寻找近似方位角，以分点讨论if(yoI 1-yoI)0(XOI 1-XOI)0)/第一个象限alfoI=atan(yoI 1-yoI)/(XOI 1-XOI)；else if(yoI 1-yoI)0(XOI 1-XOI)0)/第二象限alfoI=atan(youI 1-youI)/(XOI 1-XOI)pi；else if(yoI 1-yoI)0(XOI 1-XOI)0)/第三象限alfoI=atan(yoI 1-yoI)/(XOI 1-XOI)2 * pi；else/(yoI-yoI-1)0(XOI-XOI-1)0)/第四象限alfoI=atan(youI 1-youI)/(XOI 1-XOI)pi；/printf(alf=%ft ，alfoI)；/printf(so=%ft ，soI)；/b查找矩阵/单独计算第一个角度的系数j=0；/B0j=0.0；(约14-约0)/(so14* so14)* p/1000.0；/B0j 1=0.0；-1 *(XO14-XO0)/(so14* so14)* p/1000.0；B 0 j=-1 *(约1-约0)/(so0* so0)* p/1000.0；/b0J1=(XO1-XO0)/(so0* so0)* p/1000.0；B 0 26=(约14-约0)/(so14* so14)* p/1000.0；b027=-1 *(XO14-XO0)/(so14* so14)* p/1000.0/单独计算第二个角度的系数b1j=(yo2-yo1* so1)-(yo0-yo1)/(so0/b1J1=-1 *(XO2-XO1)/(so1* so1)-(XO)B 1 j 2=-1 *(约2-约1)/(so1* so1)* p/1000.0；b1j 3=(XO2-XO1)/(so1* so1)* p/1000.0；/寻找其他角度修正系数for(I=2；I15I)If(i14)bIj=(youI-1-youI)/(soI-1* soI-1)* p/11bIj 1=-1 *(XOI-1-XOI)/(soI-1* soI-1)bIj 2=(yoI 1-yoI)/(soI* soI)-(yoI-1-)bIj 3=-1 *(XOI 1-XOI)/(soI* soI)-(XOIbIj 4=-1 *(yoI 1-yoI)/(soI* soI)* p/1000.0bIj 5=(XOI 1-XOI)/(soI* soI)* p/1000.0；ElsebIj=(youI-1-youI)/(soI-1* soI-1)* p/1000bIj 1=-1 *(XOI-1-XOI)/(soI-1* soI-1)bIj 2=(yoI 1-yoI)/(soI* soI)-(yoI-1-)bIj 3=-1 *(XOI 1-XOI)/(soI* soI)-(XOIj=j 2；b21=0.0；/寻找侧长度修正