海底地形测量图的插值模型MATLAB代码整理.pdf

下载文档 降价啦

95
0
约5.27千字
约 5页
2021-12-03 发布于江苏
举报
版权申诉
保障服务

海底地形测量图的插值模型MATLAB代码整理.pdf

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

%数据输入 x_Yd=[129.0 140.0 108.5 88.0 185.5 195.5 105.5 157.5 107.5 77.0 81.0 162.0 162.0 117.5]; y_Yd=[ 7.5 141.5 28.0 147.0 22.5 137.5 85.5 -6.5 -81.0 3.0 56.5 -66.5 84.0 -38.5]; z_Ft=[ 4 8 6 8 6 8 8 9 9 8 8 9 4 9]; % 由已知 14 点计算出 14*14 点（其中有两点重复） %计算 QiQj 的值 Qi_Qj Qi_Qj=zeros(14,14); for i=1:14 for j=1:14 Qi_Qj(i,j)=sqrt((x_Yd(i)-x_Yd(j)).^2+(y_Yd(i)-y_Yd(j)).^2); end end [G_x,G_y]=meshgrid(x_Yd,y_Yd); %计算所求 196 个 G 点的坐标 %求出每个 G 点对应的深度 Z_g for m=1:14 for k=1:14 if ((m==k)|((m==12)(k==13))|((m==13)(k==12))) %排除 Qi,Qj 为同一点和重复点的情况 Z_g(m,k)=z_Ft(m); continue; else %计算 GQij=min{GQi,GQj} G_Qi=zeros(1,14); G_Qj=zeros(1,14); for i=1:14 G_Qi(i)=sqrt((G_y(m,k)-y_Yd(i)).^2+(G_x(m,k)-x_Yd(i)).^2); end for j=1:14 G_Qj(j)=sqrt((G_x(m,k)-x_Yd(j)).^2+(G_y(m,k)-y_Yd(j)).^2); end G_Qij=zeros(14,14); [G_Qi,G_Qj]=meshgrid(G_Qi,G_Qj); G_Qij=min(G_Qi,G_Qj); %计算 G 点到 Qi,Qj 两点连线的距离 G_Pij for i=1:14 for j=1:14 if sqrt((y_Yd(i)-y_Yd(j)).^2+(x_Yd(i)-x_Yd(j)).^2) G_Pij(i,j)=abs((G_x(m,k)-x_Yd(i)).*(y_Yd(i)-y_Yd(j))-(G_y(m,k)-y_Yd(i)).*(x_Yd(i)-x_Yd(j))). /sqrt((y_Yd(i)-y_Yd(j)).^2+(x_Yd(i)-x_Yd(j)).^2); end end end %计算 PQij=min{PQi,PQj} P_Qi=sqrt(G_Qi.^2-G_Pij.^2); P_Qj=sqrt(G_Qj.^2-G_Pij.^2); P_Qij=min(P_Qi,P_Qj); %