基于MATLAB的控制网平差程序设计第五章源代码.docVIP

观测数据读入程序，rddat1函数（85页） global net ed dd sd dd1 pn x0 y0 m1 m2 m3 ms pp e d sid md g f dir ni si ma s t az aa bb cc rt rr tt global pathname filename x0=[];y0=[];e=[];d=[];sid=[];g=[];f=[];dir=[];si=[];ni=[];s=[];t=[];az=[];pn=[]; [filename,pathname]=uigetfile(*.txt,请选择原始数据); fit1=fopen(strcat(pathname,filename),rt); if(fit1==-1) msgbox(Input File or Path is not correct,Warning,warn); return; end net=fscanf(fit1,%d,1); [a]=fscanf(fit1,%d,3); ed=a(1);dd=a(2);dd1=a(3); sd=ed+dd; [pn]=fscanf(fit1,%d,sd); [a]=fscanf(fit1,%f,2*ed); for i=1:ed x0(i)=a(2*i-1); y0(i)=a(2*i); end [a]=fscanf(fit1,%d,3); m1=a(1);m2=a(2);m3=a(3); isid=0; [a]=fscanf(fit1,%f,2); ms=a(1); pp=a(2); [a]=fscanf(fit1,%d %d %f,3*m1); for i=1:m1 e(i)=a(3*i-2);d(i)=a(3*i-1);sid(i)=a(3*i); end [e,i1]=chkdat(sd,pn,e); [d,i2]=chkdat(sd,pn,d); i3=0; isid=i1+i2+i3; idir=0; md=fscanf(fit1,%f,1); [a]=fscanf(fit1,%d %d %f,3*m2); for i=1:m2 n1(i)=a(3*i-2); n2(i)=a(3*i-1); unk(i)=a(3*i); end [n1,i1]=chkdat(sd,pn,n1); [n2,i2]=chkdat(sd,pn,n2); i3=0; ik=1; si(1)=1; for i=1:sd ii=0; for j=1:m2 if(n1(j)==j) ii=ii+1; g(ik)=n1(j); f(ik)=n2(j); dir(ik)=unk(j); ik=ik+1; end end ni(i)=ii; si(i+1)=si(i)+ni(i); end idir=i1+i2+i3; iaz=0; if(m30) ma=fscanf(fit1,%f,1); [a]=fscanf(fit1,%d %d %f,3*m3); for i=1:m3 s(i)=a(3*i-2); t(i)=a(3*i-1); az(i)=a(3*i); end [s,i1]=chkdat(sd,pn,s); [t,i2]=chkdat(sd,pn,t); i3=0; iaz=i1+i2+i3; end kk=isid+idir+iaz; if(kk0) msgbox(Error by function rddat1,Warning,warn); return; end fclose(all); 误差方程与法方程的组成函数-obnorm函数（90页） function obnorm global ed dd dd1 ni si e d g f s t global m1 m2 m3 ms pp md ma x0 y0 sid dir az c fit1 fit2 glob