MATLAB平差程序VIP免费

function BJPC2 ZB=[];ZJ=[];BC=[]; %% format long [filename,filepath]=uigetfile('*.txt','选择平差文件:'); name=[filepath filename]; fid1=fopen(name,'rt'); if(fid1==-1) disp('文件未打开，请重试！'); return; end n=fscanf(fid1,'%f',1);%输入边的个数 b=fscanf(fid1,'%f',1);%输入多余观测个数 E=fscanf(fid1,'%f',1);%输入测角中误差 t=fscanf(fid1,'%f',1);%输入坐标个数 ZB=fscanf(fid1,'%f',[2 t]); dms1=fscanf(fid1,'%f',[3 n+1]); BC=fscanf(fid1,'%f',[1 n]); T0=fscanf(fid1,'%f',1); T1=fscanf(fid1,'%f',1); fclose(fid1); dms=dms1'; ZJ=dms2degrees(dms); %% if t==2 T0=T0;%若告诉起始方位角就直接输入 end if t==4 %若没有告诉起始方位角由坐标反算 Xa=ZB(1,1);Xb=ZB(1,2); Ya=ZB(2,1);Yb=ZB(2,2); Tx=Xb-Xa; Ty=Yb-Ya; T0=atan(Ty/Tx)*180/pi;%对方位角的讨论 if ((Tx>0)&&(Ty>0)) T0; end if (((Tx<0)&&(Ty>0))||((Tx<0)&&(Ty<0))) T0=T0+180; end if ((Tx>0)&&(Ty<0)) T0=T0+360; end if ((Tx==0)&&(Ty>0)) T0=90; end if ((Tx==0)&&(Ty<0)) T0=270; end if((Ty==0)&&(Tx>0)) T0=0; end if ((Ty==0)&&(Tx<0)) T0=180; end end a0=T0; x0=ZB(1,t/2); y0=ZB(2,t/2); A=[0,0]; FWJ=[];%开始计算近似方位角 J=1; while(J<=n+1) belta=ZJ(J);%讨论方位角 if J==1 a=a0+belta; else a=a+belta; end if a>=180 a=a-180; else a=a+180; end if a>=360 a=a-360; end FWJ(J)=a; J=J+1; end FWJ; %% J=1;%开始计算近似坐标A while(J<=n) if J==1 A(1,1)=x0+BC(1)*cos(FWJ(1)*pi/180); A(2,1)=y0+BC(1)*sin(FWJ(1)*pi/180); else A(1,J)=A(1,J-1)+BC(J)*cos(FWJ(J)*pi/180); A(2,J)=A(2,J-1)+BC(J)*sin(FWJ(J)*pi/180); end J=J+1; end A; %% W=[];%开始计算闭合差 if t==2 T1=T1; end if t==4 Xa=ZB(1,3);Xb=ZB(1,4); Ya=ZB(2,3);Yb=ZB(2,4); Tx=Xb-Xa; Ty=Yb-Ya; T1=atan(Ty/Tx)*180/pi; if ((Tx>0)&&(Ty>0)) T1; end if (((Tx<0)&&(Ty>0))||((Tx<0)&&(Ty<0))) T1=T1+180; end if ((Tx>0)&&(Ty<0)) T1=T1+360; end if ((Tx==0)&&(Ty>0)) T1=90; end if ((Tx==0)&&(Ty<0)) T1=270; end if((Ty==0)&&(Tx>0)) T1=0; end if ((Ty==0)&&(Tx<0)) T1=180; end end W(1,1...