海底地形测量图的插值模型MATLAB代码.doc
版权申诉
66 浏览量
2022-07-05
14:15:19
上传
评论 1
收藏 175KB DOC 举报
%数据输入
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)
资源评论
