MATLAB平差程序_原创精品文档.pdfVIP

MATLAB平差程序--第1页

functionBJPC2

ZB=[];ZJ=[];BC=[];

%%

formatlong

[filename,filepath]=uigetfile(*.txt,选择平差文件:);

name=[filepathfilename];

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,[2t]);

dms1=fscanf(fid1,%f,[3n+1]);

BC=fscanf(fid1,%f,[1n]);

T0=fscanf(fid1,%f,1);

T1=fscanf(fid1,%f,1);

fclose(fid1);

dms=dms1;

ZJ=dms2degrees(dms);

%%

ift==2

T0=T0;%若告诉起始方位角就直接输入

end

ift==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((Tx0)(Ty0))

T0;

end

if(((Tx0)(Ty0))||((Tx0)(Ty0)))

T0=T0+180;

end

if((Tx0)(Ty0))

T0=T0+360;

end

if((Tx==0)(Ty0))

T0=90;

end

if((Tx==0)(Ty0))

T0=270;

end

MATLAB平差程序--第1页

MATLAB平差程序--第2页

if((Ty==0)(Tx0))

T0=0;

end

if((Ty==0)(Tx0))

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);%讨论方位角