%---------數(shù)據(jù)準(zhǔn)備-------------------------------------------------------------
clc;clear?all;
global?d
global?jd
global?bc
flag=-9999999;wzdx=flag;
%---------打開文本文件-----------------------------------------------------------
[filename1filepath]=uigetfile（‘*.txt‘‘請選擇點坐標(biāo)數(shù)據(jù)文件‘）;
fd=fopen（strcat（filepathfilename1）‘rt‘）;
[filename1filepath]=uigetfile（‘*.txt‘‘請選擇角度值數(shù)據(jù)文件‘）;
fj=fopen（strcat（filepathfilename1）‘rt‘）;
[filename1filepath]=uigetfile（‘*.txt‘‘請選擇邊長值數(shù)據(jù)文件‘）;
fb=fopen（strcat（filepathfilename1）‘rt‘）;
%fd=fopen（‘點坐標(biāo).txt‘‘r‘）;
%fj=fopen（‘角度值.txt‘‘r‘）;
%fb=fopen（‘邊長值.txt‘‘r‘）;
%-------------------------------讀取已知點坐標(biāo)并編號--------------------------------
ds=0;dmz=[];
while?（~feof（fd））
????ds=ds+1;
????dd（ds）.dm=fscanf（fd‘%s‘1）;dmz=[dmz?dd（ds）.dm];
????dd（ds）.dx=fscanf（fd‘%f‘1）;
????dd（ds）.dy=fscanf（fd‘%f‘1）;
????dd（ds）.dh=ds;
????t0=ds;???????????????????????????????????????%已知點的個數(shù)
end
%-------------------------------讀取角度觀測值，并形成角度矩度jdz---------------------
js=0;jdz=[];
while（~feof（fj））
????js=js+1;
????jjd（js）.czm=fscanf（fj‘%s‘1）;
????jjd（js）.zzm=fscanf（fj‘%s‘1）;
????jjd（js）.yzm=fscanf（fj‘%s‘1）;
????jjd（js）.jd=fscanf（fj‘%f‘1）;
????jjd（js）.zwc=fscanf（fj‘%f‘1）;
????jdz=[jdz?jjd（js）.jd];
end
%-------------------------------讀取邊長觀測值，并形成邊長矩陣bcz-----------------------
bs=0;bcz=[];
while（~feof（fb））
????bs=bs+1;
????bbc（bs）.qdm=fscanf（fb‘%s‘1）;
????bbc（bs）.zdm=fscanf（fb‘%s‘1）;
????bbc（bs）.bc=fscanf（fb‘%f‘1）;
????bbc（bs）.zwc=fscanf（fb‘%f‘1）;
????bcz=[bcz?bbc（bs）.bc];
end
fclose（fd）;fclose（fj）;fclose（fb）;
%------------------------------統(tǒng)計未知點個數(shù)，放到t1中-------------------------------
t1=0;tt=t0;
for?i=1:js
????d1=jjd（i）.czm;d2=jjd（i）.zzm;d3=jjd（i）.yzm;
????dh1=-999;dh2=dh1;dh3=dh2;
????for?j=1:tt
????????if?d1==dd（j）.dm;
????????????dh1=dd（j）.dh;
????????????break
????????end
????end
????for?j=1:tt
????????if?d2==dd（j）.dm;
????????????dh2=dd（j）.dh;
????????????break
????????end
????end
????for?j=1:tt
????????if?d3==dd（j）.dm;
????????????dh3=dd（j）.dh;
????????????break
????????end
????end
????if?dh1==-999
????????tt=tt+1;t1=t1+1;
????????dd（tt）.dm=d1;dd（tt）.dx=flag;
????????dd（tt）.dh=tt;dd（tt）.dy=flag;
????????dh1=dd（tt）.dh;
????end
????if?dh2==-999
????????tt=tt+1;t1=t1+1;
????????dd（tt）.dm=d2;dd（tt）.dx=flag;
????????dd（tt）.dh=tt;dd（tt）.dy=flag;
????????dh2=dd（tt）.dh;
????end
????if?dh3==-999
????????tt=tt+1;t1=t1+1;
????????dd（tt）.dm=d3;dd（tt）.dx=flag;
????????dd（tt）.dh=tt;dd（tt）.dy=flag;
????????dh3=dd（tt）.dh;
????end
????jjd（i）.cdh=dh1;jjd（i）.zdh=dh2;jjd（i）.ydh=dh3;
end
%-------------------------------------------形成平差用的角度觀測值矩陣------------------
?for?i=1:js
?????%disp（strcat（num2str（jjd（i）.cdh）‘--‘num2str（jjd（i）.zdh）‘--‘num2str（jjd（i）.ydh）‘--‘num2str（jjd（i）.jd）））;
?????jd（i1）=jjd（i）.cdh;jd（i2）=jjd（i）.zdh;jd（i3）=jjd（i）.ydh;jd（i4）=jjd（i）.jd;jd（i5）=jjd（i）.zwc;
?end
?disp（strcat（‘未知點個數(shù)=‘num2str（t1）））;
?%---------------------------