classdef MeasurementAjustmentforTraverseNetwork < matlab.apps.AppBase

% Properties that correspond to app components

properties (Access = public)

ErrorAjustment matlab.ui.Figure

HomeMenu matlab.ui.container.Menu

MenuSave matlab.ui.container.Menu

Menu_2 matlab.ui.container.Menu

Menu_8 matlab.ui.container.Menu

Menu_9 matlab.ui.container.Menu

Menu_10 matlab.ui.container.Menu

BMenu matlab.ui.container.Menu

PMenu matlab.ui.container.Menu

lMenu matlab.ui.container.Menu

VMenu matlab.ui.container.Menu

xMenu matlab.ui.container.Menu

HelpMenu matlab.ui.container.Menu

Menu_4 matlab.ui.container.Menu

Menu_5 matlab.ui.container.Menu

AboutMenu matlab.ui.container.Menu

Menu_6 matlab.ui.container.Menu

Menu_7 matlab.ui.container.Menu

computeButton matlab.ui.control.Button

Label matlab.ui.control.Label

TextArea matlab.ui.control.TextArea

EditField_3Label matlab.ui.control.Label

A1EditField matlab.ui.control.NumericEditField

EditField_4Label matlab.ui.control.Label

A2EditField matlab.ui.control.NumericEditField

XmEditFieldLabel matlab.ui.control.Label

A1xEditField matlab.ui.control.NumericEditField

YmEditFieldLabel matlab.ui.control.Label

A1yEditField matlab.ui.control.NumericEditField

YmEditField_2Label matlab.ui.control.Label

A2yEditField matlab.ui.control.NumericEditField

XmEditField_2Label matlab.ui.control.Label

A2xEditField matlab.ui.control.NumericEditField

UIAxes matlab.ui.control.UIAxes

bianButton matlab.ui.control.Button

jiaoButton matlab.ui.control.Button

bianLamp matlab.ui.control.Lamp

jiaoLamp matlab.ui.control.Lamp

biantipEditField matlab.ui.control.EditField

jiaotipEditField matlab.ui.control.EditField

Label_2 matlab.ui.control.Label

jiaosrcTextArea matlab.ui.control.TextArea

Label_3 matlab.ui.control.Label

biansrcTextArea matlab.ui.control.TextArea

TabGroup matlab.ui.container.TabGroup

BTab matlab.ui.container.Tab

UITableB matlab.ui.control.Table

PTab matlab.ui.container.Tab

UITableP matlab.ui.control.Table

lTab matlab.ui.container.Tab

UITablel matlab.ui.control.Table

VTab matlab.ui.container.Tab

UITable matlab.ui.control.Table

xTab matlab.ui.container.Tab

UITablex matlab.ui.control.Table

Label_4 matlab.ui.control.Label

Spinner matlab.ui.control.Spinner

Label_5 matlab.ui.control.Label

EditFieldCount matlab.ui.control.NumericEditField

end

properties (Access = public)

bian % 边数据

jiao % 角数据

N_angl

N_line

N

A1

A2

A1X

A2X

A1Y

A2Y

data

count

end

methods (Access = public)

function [fwj] = fwjjs(app,a1,a2,jiao,fwj,N_angl)

for i=1:1:N_angl

if (jiao(i,1)+1)==a1 && (jiao(i,2)+1)==a2 && fwj(a2,(jiao(i,3)+1))==0

fwj(a2,jiao(i,3)+1)=fwj(a1,a2)+jiao(i,4)+pi;

fwj=fwjjs(app,a2,jiao(i,3)+1,jiao,fwj,N_angl);

%迭代计算

elseif (jiao(i,3)+1)==a1 && (jiao(i,2)+1)==a2 && fwj(a2,(jiao(i,1)+1))==0

%观测角终点为a1 && 观测角顶点为a2 && 方位角（a2,观测角起点）未计算

fwj(a2,jiao(i,1)+1)=fwj(a1,a2)-jiao(i,4)+pi;

fwj=fwjjs(app,a2,jiao(i,1)+1,jiao,fwj,N_angl);

elseif (jiao(i,1)+1)==a2 && (jiao(i,2)+1)==a1 && fwj(a1,(jiao(i,3)+1))==0

%观测角起点为a2 && 观测角顶点为a1 && 方位角（a1,观测角起点）未计算

fwj(a1,(jiao(i,3)+1))=fwj(a1,a2)+jiao(i,4);

fwj=fwjjs(app,a1,jiao(i,3)+1,jiao,fwj,N_angl);

elseif (jiao(i,2)+1)==a1 && (jiao(i,3)+1)==a2 && fwj(a1,(jiao(i,1)+1))==0

%观测角终点为a2 && 观测角顶点为a1 && 方位角（a1,观测角起点）未计算

fwj(a1,(jiao(i,1)+1))=fwj(a1,a2)-jiao(i,4);

fwj=fwjjs(app,a1,(jiao(i,1)+1),jiao,fwj,N_angl);

end

end

end

function [X,Y] = zuobiao(app,A,AX,AY,bian,fwj,cc,X,Y,N_line)

X(A,1)=AX;

Y(A,1)=AY;

%已知起算点输入

for i=1:1:N_line

if (bian(i,2)+1==A) && (cc(bian(i,3)+1,1)==0)

cc(bian(i,3)+1,1)=1;

X(bian(i,3)+1,1)=AX+bian(i,4)*cos(fwj(bian(i,2)+1,bian(i,3)+1));

Y(bian(i,3)+1,1)=AY+bian(i,4)*sin(fwj(bian(i,2)+1,bian(i,3)+1));

[X,Y]=zuobiao(app,bian(i,3)+1,X(bian(i,3)+1,1),Y(bian(i,3)+1,1),bian,fwj,cc,X,Y,N_line);

elseif(bian(i,3)+1==A) && (cc(bian(i,2)+1,1)==0)

cc(bian(i,2)+1,1)=1;

X(bian(i,2)+1,1)=AX+bian(i,4)*cos(fwj(bian(i,3)+1,bian(i,2)+1));

Y(bian(i,2)+1,1)=AY+bian(i,4)*sin(fwj(bian(i,3)+1,bian(i,2)+1));

[X,Y]=zuobiao(app,bian(i,2)+1,X(bian(i,2)+1,1),Y(bian(i,2)+1,1),bian,fwj,cc,X,Y,N_line);

end

end

end

function [B,bian,jiao,X01,v1,xxm,v,P,sgm0,sgmdianmm]= dxw(app,A1,A1X,A1Y,A2,A2X,A2Y,bian,jiao,N,N_angl,N_line)

%输入数据：*****************************************

%A1点 X Y

%A2点 X Y

%起算基础

%bian边 边号 起始点 终止点 原始数据 平差数据

%jiao角 角起点 角顶点 角终点 原始数据 平差数据

%N总点数

%N_angl总观测角数

%N_line总边数

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%计算方位角近似值

%------------------------------------------------------------------------%

fwj=zeros(N,N);%总共N个点

%初始化

%A1加1的目的：点号存储是从0开始，matlab矩阵索引从1开始

fwj(A1+1,A2+1)=atan2((A2Y-A1Y),(A2X-A1X));%计算已知相邻点的方位角

fwj=fwjjs(app,A1+1,A2+1,jiao,fwj,N_angl);

for i=1:1:N

for j=1:1:N

while fwj(i,j)>2*pi

fwj(i,j)=fwj(i,j)-2*pi;

end

while fwj(i,j)<0

fwj(i,j)=fwj(i,j)+2*pi;

end

end

end

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%计算近似坐标 均为列向量

%------------------------------------------------------------------------%

%存储2次不同起算数据计算后的坐标均值

X0=zeros(N,1);

Y0=zeros(N,1);

%

X1=zeros(N,1);

Y1=zeros(N,1);

X2=zeros(N,1);

Y2=zeros(N,1);

%初始化

global cc;

cc=zeros(N,1);

[X1,Y1]=zuobiao(app,A2+1,A2X,A2Y,bian,fwj,cc,X1,Y1,N_line);

%输入数据：*****************************************

%起算点序号 起算点X 起算点Y 边原始数据 方位角数据 cc标记数组 坐标X近似值 坐标Y近似值

cc=zeros(N,1);

[X2,Y2]=zuobiao(app,A1+1,A1X,A1Y,bian,fwj,cc,X2,Y2,N_line);

for i=1:1:N

%总计N个点 对两次不同起算基础计算后的结果 进行平均

X0(i,1)=(X1(i,1)+X2(i,1))/2;

Y0(i,1)=(Y1(i,1)+Y2(i,1))/2;

end

X0(A1+1)=A1X;

Y0(A1+1)=A1Y;

X0(A2+1)=A2X;

Y0(A2+1)=A2Y;

%还原已知点坐标

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%初始化坐标值的改正数 N个点 每个点有X,Y 总计N*2

%------------------------------------------------------------------------%

x0=zeros(N*2,1);

for i=1:1:N

x0(2*i-1,1)=X0(i,1);%奇数行存X

x0(2*i,1)=Y0(i,1); %偶数行存Y

end

x0([2*A1+1,2*A1+2,2*A2+1,2*A2+2],:)=[];%参数的近似值

%删除已知点的X,Y坐标改正数

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%计算角度观测值的系数矩阵 jiaoxishu

%------------------------------------------------------------------------%

%总计N_angl个观测角

jiaoxishu=zeros(N_angl,N*2);

%计算jiaoxishu

for i=1:1:N_angl

j=jiao(i,2)+1;%观测角 顶点j

k=jiao(i,3)+1;%观测角 终点k

h=jiao(i,1)+1;%观测角 起点h

rou=180*3600/pi;

%秒与弧度转换参数ρ

up1=rou*sin(fwj(j,k))/1000;

up2=-rou*cos(fwj(j,k))/1000;

%单位为 ″/mm 所以除以1000

under1=sqrt((X0(k,1)-X0(j,1))^2+(Y0(k,1)-Y0(j,1))^2);

%S0 观测边长

a1=up1/under1;

b1=up2/under1;

%ajk 与 bjk

up3=rou*sin(fwj(j,h))/1000;

up4=-rou*cos(fwj(j,h))/1000;

under2=sqrt((X0(h,1)-X0(j,1))^2+(Y0(h,1)-Y0(j,1))^2);

a2=up3/under2;

b2=up4/under2;

%i为方程行号

jiaoxishu(i,2*h-1)=a2; %h点x系数

jiaoxishu(i,2*h)=b2; %h点Y系数

jiaoxishu(i,2*j-1)=a1-a2;

jiaoxishu(i,2*j)=b1-b2;

jiaoxishu(i,2*k-1)=-a1;

jiaoxishu(i,2*k)=-b1;

end

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%计算边系数矩阵 bianxishu

%------------------------------------------------------------------------%

%总计N_line条边

bian0=zeros(N_line,1);%计算边长近似值bian0

for i=1:1:N_line

j1=bian(i,2);

%边起点

j2=bian(i,3);

%边终点

bian0(i)=sqrt((X0(j1+1)-X0(j2+1))*(X0(j1+1)-X0(j2+1))+(Y0(j1+1)-Y0(j2+1))*(Y0(j1+1)-Y0(j2+1)));

end

bianxishu=zeros(N_line,N*2);%计算边长观测值的系数

for i=1:1:N_line

j=bian(i,2);%观测边起点

k=bian(i,3);%观测边终点

up1=X0(k+1)-X0(j+1); %ΔXjk

up2=Y0(k+1)-Y0(j+1); %ΔYjk

under=bian0(i); %边的观测值

m=up1/under; %ΔX

n=up2/under; %ΔY

bianxishu(i,2*(j+1)-1)=-m;

bianxishu(i,2*(j+1))=-n;

bianxishu(i,2*(k+1)-1)=m;

bianxishu(i,2*(k+1))=n;

end

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%构造系数矩阵B

%------------------------------------------------------------------------%

B=[jiaoxishu;bianxishu];

B(N_angl+N_line,:)=[];

%删除第N_angl+N_ling行已知边长

B(:,[2*A1+1,2*A1+2,2*A2+1,2*A2+2])=[];

%删除系数矩阵中已知点的系数

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%计算P矩阵

%------------------------------------------------------------------------%

P=zeros(N_angl+N_line,N_angl+N_line);

%N_angl+N_line 总观测量

for i=1:1:N_angl

P(i,i)=1;

%角度观测的权为1

end

j=1;

for i=N_angl+1:1:N_angl+N_line

P(i,i)=4*12*12/bian(j,5)/bian(j,5);

%边长观测权

j=j+1;

end

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%计算l矩阵

%------------------------------------------------------------------------%

l1=zeros(N_line,1);

%边改正数l1初始化

for i=1:N_line

%bian0 边近似值

%bian(i,5) 边观测值

l1(i,1)=(bian(i,5)-bian0(i))*1000;%化单位为毫米

if l1(i,1)<0.0000001 && l1(i,1)>-0.0000001

l1(i,1)=0;

%低于要求精度 取为0

end

end

l1(N_line,:)=[];

%删除已知点对应边

jiao01=zeros(N_angl,1);

jiao02=zeros(N_angl,1);

jiao0=zeros(N_angl,1);

for i=1:1:N_angl %由近似坐标计算近似观测角度

jiao01(i,1)=atan2((Y0(jiao(i,3)+1,1)-Y0(jiao(i,2)+1,1)),(X0(jiao(i,3)+1,1)-X0(jiao(i,2)+1,1)));

%jk方位角

jiao02(i,1)=atan2((Y0(jiao(i,1)+1,1)-Y0(jiao(i,2)+1,1)),(X0(jiao(i,1)+1,1)-X0(jiao(i,2)+1,1)));

%jh方位角

%atan2 返回 Y 和 X 的四象限反正切

jiao0(i,1)=jiao01(i,1)-jiao02(i,1);

%jk-jh=hjk

end

for i=1:1:N_angl %将观测角置于 0-2*pi

while jiao0(i,1)>2*pi

jiao0(i,1)=jiao0(i,1)-2*pi;

end

while jiao0(i,1)<0

jiao0(i,1)=jiao0(i,1)+2*pi;

end

end

l2=zeros(N_angl,1);

for i=1:1:N_angl %计算角度观测值的l，并将单位化为秒

l2(i,1)=(jiao(i,5)-jiao0(i))*3600*180/pi;

if l2(i,1)<0.00000001 && l2(i,1)>-0.00000001

l2(i,1)=0;

end

end

l=[l2;l1];%构成矩阵l

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%列法方程解算

%------------------------------------------------------------------------%

P(N_angl+N_line,:)=[];

P(:,N_angl+N_line)=[];

%清除权阵中多余已知点

NBB=B'*P*B;

W=B'*P*l;

xx=NBB\W;%单位mm

v=B*xx-l;

%观测值改正数v 单位为mm

xxm=xx.*0.001;

%将坐标改正值单位转化为 m

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%精度评定

%------------------------------------------------------------------------%

sgm0=sqrt((v'*P*v)/(N_angl+N_line-1-(N*2-4)));

%单位权中误差 n-t=n-u

X01=x0+xxm;%参数平差值

v1=zeros(N_line+N_angl,1);%边平差值

for i=N_angl+1:1:N_angl+N_line-1

v1(i,1)=bian(i-N_angl,5)+v(i,1)*0.001;

end

v1(N_angl+N_line,1)=bian(N_line,5);

for i=1:1:N_angl

v1(i,1)=jiao(i,5)+v(i,1)/3600*pi/180;

end

Qxx=inv(NBB);%协因数阵

Qx=zeros(N-2,1);%对x坐标求协因数

for i=1:1:N-2

Qx(i,1)=Qxx(2*i-1,2*i-1);

end

Qy=zeros(N-2,1);%对y坐标求协因数

for i=1:1:N-2

Qy(i,1)=Qxx(2*i,2*i);

end

sgmdian=zeros(N-2,1);

%每个点的点位误差 单位为mm

for i=1:1:N-2

sgmdian(i,1)=sgm0*sqrt(Qx(i,1)+Qy(i,1));

end

sgmdianmm=zeros(N-2,1);

%每个点的点位误差 单位为m

for i=1:1:N-2

sgmdianmm(i,1)=sgmdian(i,1);

end

%************************************************************************%

%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%

%将改正结果赋值给原始数据矩阵进行迭代

%------------------------------------------------------------------------%

for i=1:1:N_line-1

bian(i,4)=bian(i,5)+v(i+N_angl,1)*0.001;

end

for i=1:1:N_angl

jiao(i,4)=v(i,1)/3600*pi/180+jiao(i,5);

end

end

end

% Callbacks that handle component events

methods (Access = private)

% Code that executes after component creation

function startupFcn(app)

app.count=0;

app.bianLamp.Color="red";

app.jiaoLamp.Color="red";

app.biantipEditField.Value="请粘贴边数据！";

app.jiaotipEditField.Value="请粘贴角数据！";

app.A1EditField.Value=0;

app.A2EditField.Value=0;

app.A1xEditField.Value=0;

app.A2xEditField.Value=0;

app.A1yEditField.Value=0;

app.A2yEditField.Value=0;

app.N=0;

app.bian=[];

app.jiao=[];

app.N_angl=0;

app.N_line=0;

end

% Button pushed function: computeButton

function computeButtonPushed(app, event)

if app.bianLamp.Color(1,1)==1 || app.jiaoLamp.Color(1,1)==1

app.TextArea.Value="未读入数据或数据格式不正确！";

return;

end

%数据未成功输入 报错并停止继续运行

try

app.A1=app.A1EditField.Value;%A1点号

app.A2=app.A2EditField.Value;%A2点号

%A1 A2坐标

app.A1X=app.A1xEditField.Value;

app.A1Y=app.A1yEditField.Value;

app.A2X=app.A2xEditField.Value;

app.A2Y=app.A2yEditField.Value;

app.N=max(app.jiao(:,1:3),[],'all','omitnan')+1;

%N总点数

app.N_angl=length(app.jiao(:,1));

%总观测角数

knownLine=sqrt((app.A1X-app.A2X)^2+(app.A1Y-app.A2Y)^2);

app.N_line=length(app.bian(:,1));

app.bian=[[1:1:app.N_line]',app.bian];

%计算已知边 并加入数据

app.bian=[app.bian;[app.N_line+1, app.A1, app.A2, knownLine, knownLine]];

%添加已知边

app.N_line=app.N_line+1;

%总边数

i=1;

app.count=1;

while i==1%进行迭代运算

%dxw 导线网平差函数

[~,app.bian,app.jiao,~,~,xxm1,~,~,~,~]=dxw(app,app.A1,app.A1X,app.A1Y,app.A2,app.A2X,app.A2Y,app.bian,app.jiao,app.N,app.N_angl,app.N_line);

%输入数据：

%A1点 X Y

%A2点 X Y

%bian边： 边号 起始点 终止点 平差迭代数据 原始数据备份

%jiao角： 角起点 角顶点 角终点 平差迭代数据 原始数据备份

%N 总点数

%N_angl 观测角数量

%N_line 边数量

%输出数据：

%bian

%jiao

%X参数平差值

%L 观测边与角的平差值

%xxm 坐标改正数 以m为单位

%V 改正数

%P 权阵

%sgm0 单位权中误差

%sgmdianmi 每个点的点位误差 单位为m

[B,app.bian,app.jiao,X,L,xxm2,V,P,sgm0,sgmdianmm]=dxw(app,app.A1,app.A1X,app.A1Y,app.A2,app.A2X,app.A2Y,app.bian,app.jiao,app.N,app.N_angl,app.N_line);

app.data{app.count,1}=B;

app.data{app.count,2}=P;

l=B*xxm2-V;

app.data{app.count,3}=l;

app.data{app.count,4}=V;

app.data{app.count,5}=xxm2;

app.count=app.count+1;

if max(abs(xxm2-xxm1))>0.00000001

i=1;

else

i=0;

end

end

Nbb=B'*P*B;

Qll=B*inv(Nbb)*B';

Dll=sgm0^2*Qll;

Dll2=diag(Dll);

sgmll=sqrt(Dll2);

[sgmllmax,index]=max(sgmll(app.N_angl+1:app.N_angl+app.N_line-1,1));

%求最弱边中误差 单位为mm

sgmllmax=sgmllmax/1000/app.bian(index,4);

%求最弱边相对中误差

L(1:app.N_angl,1)=rad2deg(L(1:app.N_angl,1));

%观测角平差值弧度转换为度

dushu1=L(1:app.N_angl,1);

dushu2=string(fix(dushu1))+"°";

dushu1=(dushu1-fix(dushu1))*60;

dushu2=dushu2+string(fix(dushu1))+"′";

dushu1=(dushu1-fix(dushu1))*60;

dushu1=round(dushu1,1);

%取位精度为秒小数点后1位

dushu(:,2)=dushu2+string(dushu1)+"″";

dushu(:,1)="∠"+string(app.jiao(:,1))+string(app.jiao(:,2))+string(app.jiao(:,3))+"=";

%dushu存储角度平差值 单位为度分秒

L(app.N_angl+1:app.N_angl+app.N_line,1)=round(L(app.N_angl+1:app.N_angl+app.N_line,1),3);

%L(N_angl+1:N_angl+N_line,1)储存边观测平差值取小数点后3位，单位为m

V(1:app.N_angl)=round(V(1:app.N_angl),1);

dushu(:,3)=string(V(1:app.N_angl,1));

%dushu(:,3)观测角改正数，取小数点后1位，单位为″

V(1+app.N_angl:app.N_angl+app.N_line-1,1)=round(V(1+app.N_angl:app.N_angl+app.N_line-1,1)/1000,3);

%观测边改正数，取小数点后3位，单位为m

bianchang(:,1)=string(app.bian(:,2));

bianchang(:,2)=string(app.bian(:,3));

bianchang(:,3)=string(L(app.N_angl+1:app.N_angl+app.N_line,1));

bianchang(1:app.N_line-1,4)=string(V(app.N_angl+1:app.N_angl+app.N_line-1,1));

%bianchang储存边编号与边平差值 改正数

X=round(X,3);

%坐标改正数取小数点后位

tempX=zeros(app.N-2,1);

tempY=zeros(app.N-2,1);

cla(app.UIAxes);

%坐标区清零

hold (app.UIAxes,"on");

for i=1:app.N-2

tempX(i)=round(X(2*i-1,1),3);

tempY(i)=round(X(2*i,1),3);

end

coor=zeros(app.N,2);

coor(app.A1+1,1)=app.A1X;

coor(app.A1+1,2)=app.A1Y;

coor(app.A2+1,1)=app.A2X;

coor(app.A2+1,2)=app.A2Y;

j=1;

for i=0:app.N-3

if j-1==app.A1

j=j+1;

%跳过已知点

end

if j-1==app.A2

j=j+1;

%跳过已知点

end

coor(j,1)=X((i+1)*2-1,1);

coor(j,2)=X((i+1)*2,1);

j=j+1;

end

%coor 记录所有点的坐标 按照行号从小到大排列

plot(app.UIAxes,tempY,tempX,'bo');

plot(app.UIAxes,app.A1Y,app.A1X,'r^');

plot(app.UIAxes,app.A2Y,app.A2X,'r^');

%画点

for i=1:app.N_line

plot(app.UIAxes,[coor(app.bian(i,2)+1,2),coor(app.bian(i,3)+1,2)],[coor(app.bian(i,2)+1,1),coor(app.bian(i,3)+1,1)],'-c');

end

%连边

for i=1:app.N

text(app.UIAxes,coor(i,2)+5,coor(i,1)+5,string(i-1));

end

%标记点号

%成果数据显示：

app.TextArea.Value="*************角度平差值*************";

for i=1:app.N_angl

app.TextArea.Value=[app.TextArea.Value;dushu(i,1)+dushu(i,2)];

end

app.TextArea.Value=[app.TextArea.Value;"*************角度改正数*************"];

for i=1:app.N_angl

app.TextArea.Value=[app.TextArea.Value;dushu(i,3)];

end

app.TextArea.Value=[app.TextArea.Value;"*************边长平差值*************"];

for i=1:app.N_line-1

app.TextArea.Value=[app.TextArea.Value;bianchang(i,1)+"--"+bianchang(i,2)+" "+bianchang(i,3)];

end

app.TextArea.Value=[app.TextArea.Value;"*************边长改正数*************"];

for i=1:app.N_line-1

app.TextArea.Value=[app.TextArea.Value;bianchang(i,1)+"--"+bianchang(i,2)+" "+bianchang(i,4)];

end

app.TextArea.Value=[app.TextArea.Value;"*************坐标平差值*************"];

for i=1:(app.N)

app.TextArea.Value=[app.TextArea.Value;string(i-1)+" X:"+string(coor(i,1))+" Y:"+string(coor(i,2))];

end

app.TextArea.Value=[app.TextArea.Value;"*************点位中误差*************"];

for i=1:(app.N-2)

app.TextArea.Value=[app.TextArea.Value;string(round(sgmdianmm(i,1),3))];

end

app.TextArea.Value=[app.TextArea.Value;"*************平差后最弱边相对中误差*************"];

app.TextArea.Value=[app.TextArea.Value;string(round(sgmllmax,5))+" (1/"+string(round(1/sgmllmax,0))+")"];

app.TextArea.Value=[app.TextArea.Value;"*************平差后最弱点点位中误差*************"];

app.TextArea.Value=[app.TextArea.Value;string(round(max(sgmdianmm,[],'all','omitnan'),3))];

app.TextArea.Value=[app.TextArea.Value;"*************平差后单位权中误差*************"];

app.TextArea.Value=[app.TextArea.Value;string(round(sgm0,3))];

app.Spinner.Limits=[1,app.count-1];

app.EditFieldCount.Value=app.count-1;

value=1;

app.UITableB.Data=app.data{value,1};

app.UITableP.Data=app.data{value,2};

app.UITablel.Data=app.data{value,3};

app.UITable.Data=app.data{value,4};

app.UITablex.Data=app.data{value,5};

%数据重新初始化 为下次计算做准备

app.bianLamp.Color="red";

app.jiaoLamp.Color="red";

app.biantipEditField.Value="若需重新计算请重新读入边数据！";

app.jiaotipEditField.Value="若需重新计算请重新读入角数据！";

app.A1EditField.Value=0;

app.A2EditField.Value=0;

app.A1xEditField.Value=0;

app.A2xEditField.Value=0;

app.A1yEditField.Value=0;

app.A2yEditField.Value=0;

app.N=0;

app.N_angl=0;

app.N_line=0;

app.bian=[];

app.jiao=[];

catch ME

app.bianLamp.Color="red";

app.jiaoLamp.Color="red";

app.biantipEditField.Value="若需重新计算请重新读入边数据！";

app.jiaotipEditField.Value="若需重新计算请重新读入角数据！";

app.A1EditField.Value=0;

app.A2EditField.Value=0;

app.A1xEditField.Value=0;

app.A2xEditField.Value=0;

app.A1yEditField.Value=0;

app.A2yEditField.Value=0;

app.N=0;

app.bian=[];

app.jiao=[];

app.N_angl=0;

app.N_line=0;

throw(ME);

end

%app.UITableB.Data=dushu;

%app.TextArea.Value = uitextarea('Value',temp);

% app.TextArea.Value=app.TextArea.Value+mat2str(dushu);

% app.TextArea.Value=app.TextArea.Value+"\n";

% app.TextArea.Value=app.TextArea.Value+mat2str(bianchang);

end

% Callback function

function bianEditFieldValueChanging(app, event)

end

% Callback function

function jiaoEditFieldValueChanging(app, event)

end

% Menu selected function: MenuSave

function MenuSaveSelected(app, event)

A = app.TextArea.Value;

[file,path] = uiputfile('*.xlsx');

if path==0

return;

end

filename = fullfile(path,file);

xlswrite(filename,A);

end

% Menu selected function: Menu_2

function Menu_2Selected(app, event)

h=get(app.UIAxes,'children');

figure('visible','off');

copyobj(h,gca);

[file,path] = uiputfile('*.png');

if path==0

return

end

filename = fullfile(path,file);

saveas(gcf,filename);

end

% Menu selected function: Menu_6

function Menu_6Selected(app, event)

app.TextArea.Value=["仅限于平差科研学习！";"请勿将本程序用于商业用途";"©Li Jintao, School of Geodesy and Geomatics, Wuhan University"];

end

% Menu selected function: Menu_7

function Menu_7Selected(app, event)

app.TextArea.Value=["Email:";"lijintao.alex@qq.com";"Github:";"https://github.com/AkexStar"];

end

% Menu selected function: Menu_5

function Menu_5Selected(app, event)

app.TextArea.Value=["请将遇到的问题以邮件方式发送给lijintao.alex@qq.com";"收到后第一时间回复您"];

end

% Menu selected function: Menu_4

function Menu_4Selected(app, event)

app.TextArea.Value=[

"本程序配套于武汉大学测绘学院测量平差仿真辅助教学系统V2.0";...

"本程序中“边数据”指观测边边长数据";...

"本程序中“角数据”指观测角角度数据";...

"";...

"使用流程如下：";...

"》1、将角数据与边数据粘贴进指定文本框";...

"》2、点击读入角数据与边数据按钮";...

"》3、若数据读入成功 提示灯变绿且自动显示2个已知点点号";...

"》4、填入已知点坐标数据";...

"》5、点击平差计算按钮";...

"";...

"边数据请粘贴进上方指定区域内";...

"内容包含边对应点号与边长(m),下面为样例数据格式示例：";...

"0 -- 1 97.711";...

"1 -- 2 69.81";...

"2 -- 3 78.536";...

"3 -- 4 72.12";...

"4 -- 5 67.714";...

"5 -- 6 已知值";...

"6 -- 0 65.505";...

"3 -- 7 77.436";...

"7 -- 8 40.1";...

"8 -- 6 68.199";...

"";...

"角数据请粘贴进上方指定区域内";...

"内容包含对应角号与角度(°′″),下面为样例数据格式示例:";...

"∠0 1 2 236°00′33.5″";...

"∠1 0 6 130°33′18.9″";...

"∠1 2 3 225°51′8.5″";...

"∠2 3 4 231°45′21″";...

"∠2 3 7 281°25′11.9″";...

"∠3 4 5 213°46′6.1″";...

"∠4 3 7 49°40′8.4″";...

"∠4 5 6 246°24′48″";...

"∠5 6 0 236°45′44.4″";...

"∠5 6 8 323°59′33.9″";...

"∠0 6 8 87°13′30.7″";...

"∠3 7 8 174°53′58.8″";...

"∠7 8 6 199°35′58.9″";...

"";...

"已知点数据需手动填入，下面为样例数据：";...

"5 x=164.668(m), y=112.313(m)";...

"6 x=274.722(m), y=136.706(m)";...

"";...

"本程序默认测边误差：1/2000";...

"默认测角误差：12″";...

];

end

% Callback function

function SwitchValueChanged(app, event)

end

% Button pushed function: bianButton

function bianButtonPushed(app, event)

if app.biansrcTextArea.Value==""

app.bianLamp.Color="red";

app.biantipEditField.Value="边数据为空！请将边数据粘贴入指定位置";

return;

end

value=strings;

for i=1:length(app.biansrcTextArea.Value)

value=value+" "+app.biansrcTextArea.Value{i,1};

end

bian_temp1 = textscan(value,'%f %s %f %s');

[~,n]=size(bian_temp1);

if n<4

app.bianLamp.Color="red";

app.biantipEditField.Value="边数据格式错误！请重试或参考Help->程序说明";

return;

end

if isempty(bian_temp1{1,2})

app.bianLamp.Color="red";

app.biantipEditField.Value="边数据读入失败！请重试或参考Help->程序说明";

return;

end

bian_temp2=str2double(bian_temp1{1,4});

app.bian=[];

app.bian=[bian_temp1{1,1},bian_temp1{1,3},bian_temp2];

index=isnan(bian_temp2);

if length(index)<2

app.bianLamp.Color="red";

app.biantipEditField.Value="边数据格式错误或已知点少于两个！请重试或参考Help->程序说明";

return;

end

app.A1EditField.Value=app.bian(index,1);

app.A2EditField.Value=app.bian(index,2);

app.bian(index,:)=[];

app.bian=[app.bian,app.bian(:,3)];

app.bianLamp.Color="green";

app.biantipEditField.Value="边数据读入成功！";

end

% Button pushed function: jiaoButton

function jiaoButtonPushed(app, event)

if app.jiaosrcTextArea.Value==""

app.jiaoLamp.Color="red";

app.jiaotipEditField.Value="角数据为空！请将角数据粘贴入指定位置";

return;

end

value=strings;

for i=1:length(app.jiaosrcTextArea.Value)

value=value+app.jiaosrcTextArea.Value{i,1};

end

jiao_temp1 = textscan(value,'%c %f %f %f %f %c %f %c %f %c');

[~,n]=size(jiao_temp1);

if n<10

app.jiaoLamp.Color="red";

app.jiaotipEditField.Value="角数据格式错误！请重试或参考Help->程序说明";

return;

end

if isempty(jiao_temp1{1,2})

app.jiaoLamp.Color="red";

app.jiaotipEditField.Value="角数据读入失败！请重试或参考Help->程序说明";

return;

end

jiao_temp2=(jiao_temp1{:,5}+jiao_temp1{:,7}./60+jiao_temp1{:,9}./3600)/180*pi;

app.jiao=[];

app.jiao=[jiao_temp1{1,2},jiao_temp1{1,3},jiao_temp1{1,4},jiao_temp2,jiao_temp2];

app.jiaoLamp.Color="green";

app.jiaotipEditField.Value="角数据读入成功！";

end

% Value changed function: jiaosrcTextArea

function jiaosrcTextAreaValueChanged(app, event)

app.jiaosrcTextArea.Value = app.jiaosrcTextArea.Value;

end

% Menu selected function: Menu_8

function Menu_8Selected(app, event)

A = app.biansrcTextArea.Value;

[file,path] = uiputfile('*.xlsx');

if path==0

return;

end

filename = fullfile(path,file);

xlswrite(filename,A);

app.biantipEditField.Value="边数据已保存在："+filename;

end

% Menu selected function: Menu_9

function Menu_9Selected(app, event)

A = app.jiaosrcTextArea.Value;

[file,path] = uiputfile('*.xlsx');

if path==0

return;

end

filename = fullfile(path,file);

xlswrite(filename,A);

app.jiaotipEditField.Value="角数据已保存在："+filename;

end

% Value changed function: Spinner

function SpinnerValueChanged(app, event)

value = app.Spinner.Value;

if value==0

return;

end

if value>=app.count

app.Spinner.Value=app.count;

return;

end

app.UITableB.Data=app.data{value,1};

app.UITableP.Data=app.data{value,2};

app.UITablel.Data=app.data{value,3};

app.UITable.Data=app.data{value,4};

app.UITablex.Data=app.data{value,5};

%app.UITableP.Data=P;

%app.UITablel.Data=l;

end

% Menu selected function: BMenu

function BMenuSelected(app, event)

A = app.UITableB.Data;

[file,path] = uiputfile('*.xlsx');

if path==0

return;

end

filename = fullfile(path,file);

xlswrite(filename,A);

app.TextArea.Value="数据已保存在："+filename;

end

% Menu selected function: PMenu

function PMenuSelected(app, event)

A = app.UITableP.Data;

[file,path] = uiputfile('*.xlsx');

if path==0

return;

end

filename = fullfile(path,file);

xlswrite(filename,A);

app.TextArea.Value="数据已保存在："+filename;

end

% Menu selected function: lMenu

function lMenuSelected(app, event)

A = app.UITablel.Data;

[file,path] = uiputfile('*.xlsx');

if path==0

return;

end

filename = fullfile(path,file);

xlswrite(filename,A);

app.TextArea.Value="数据已保存在："+filename;

end

% Menu selected function: VMenu

function VMenuSelected(app, event)

A = app.UITable.Data;

[file,path] = uiputfile('*.xlsx');

if path==0

return;

end

filename = fullfile(path,file);

xlswrite(filename,A);

app.TextArea.Value="数据已保存在："+filename;

end

% Menu selected function: xMenu

function xMenuSelected(app, event)

A = app.UITablex.Data;

[file,path] = uiputfile('*.xlsx');

if path==0

return;

end

filename = fullfile(path,file);

xlswrite(filename,A);

app.TextArea.Value="数据已保存在："+filename;

end

end

% Component initialization

methods (Access = private)

% Create UIFigure and components

function createComponents(app)

% Create ErrorAjustment and hide until all components are created