原创力文档- 1、本文档共14页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
MATLAB测量数据处理5.doc
观测数据读入程序，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
global a q1 pa3
文档评论（0）