- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
matlab 常见经典平差 程序
希望本人编写的经典平差可以对 matlab 初学者以及测绘专业的学生有用
By cowboy
function void()
sprintf(请选择平差类型)
sprintf(1: 参数平差\n2: 条件平差\n3:具有条件的参数平差\n4: 具有参数的条件平差\n)
a=input(请输入对应号码\n);
switch (a)
case 1
% 参数平差 V=AX-L
sprintf(1:参数平差 V=AX-L)
n=input(请输入n=);
t=input(请输入t=);
A=input(请输入系数矩阵A=);
L=input(请输入常数项矩阵L=);
P=diag(input(请输入权阵P=))
%P=input(请输入权阵P=);
N=A*P*A;
U=A*P*L;
sprintf(开始计算)
X=inv(N)*(U);
V=A*X-L;
Qx=inv(N);
u=sqrt((V*P*V)/(n-t));
fid=fopen(data.txt,wt); %写入文件路径,文件输出
fprintf(fid,=====================参数平差: V=AX-L =====================\n);
fprintf(fid,输出n,t\n) ;
fprintf(fid,%12.5f %12.5f\n,n ,t) ;%n,t
fprintf(fid,输出矩阵A\n)
[m,n]=size(A); %A
for i=1:1:m
for j=1:1:n
if j==n
fprintf(fid,%12.5f\n,A(i,j));
else
fprintf(fid,%12.5f\t,A(i,j));
end
end
end
fprintf(fid,\n);
fprintf(fid,输出矩阵L\n) ;
[m,n]=size(L); %L
for i=1:1:m
for j=1:1:n
if j==n
fprintf(fid,%12.5f\n,L(i,j));
else
fprintf(fid,%12.5f\t,L(i,j));
end
end
end
fprintf(fid,\n);
fprintf(fid,输出矩阵P\n) ;
[m,n]=size(P); %P
for i=1:1:m
for j=1:1:n
if j==n
fprintf(fid,%12.5f\n,P(i,j));
else
fprintf(fid,%12.5f\t,P(i,j));
end
end
end
fprintf(fid,\n);
fprintf(fid,输出矩阵X\n);
[m,n]=size(X); %X
for i=1:1:m
for j=1:1:n
if j==n
fprintf(fid,%12.5f\n,X(i,j));
else
fprintf(fid,%12.5f\t,X(i,j));
end
end
end
fprintf(fid,\n);
fprintf(fid,输出矩阵V\n);
[m,n]=size(V);
您可能关注的文档
最近下载
- 实验室检测培训总结.docx VIP
- 大学学110周年庆典综艺晚会舞美灯光音响舞台工程投标文件(技术标).docx
- 《中华人民共和国安全生产法》试题.doc VIP
- 公路工程质量检验评定标准(JTG F801-2017)培训课件.pptx
- ASCO 胰腺癌领域治疗新进展.pptx
- 高三读后续写题库练习题55篇(含范文解析).pdf VIP
- 2024届湖北省七市州高三3月联考语文试题评讲课件.pptx
- 食品生产企业食品安全主要主体责任清单、每日食品安全检查记录.pdf VIP
- 急性上消化道出血急诊诊治流程--危重病课件.ppt
- 湖北省七市州2024届高三下学期3月联考二模语文试题及答案解析.docx
文档评论(0)