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

《2022年2022年海底地形测量图的插值模型MATLAB代码 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年海底地形测量图的插值模型MATLAB代码 .pdf（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、%数据输入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)=sqr

2、t(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=minGQi,GQj G_Qi=zeros(1,14);G_Qj=zeros(1,14);for i=1:14 G_Qi(i)=sqrt(G_y(

3、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 名师资料总结-精品资料欢迎下载-名师精心整理-第 1 页，共 5 页 -if sqrt(y_Yd(i)-y_Yd(j).2+(x_Yd(i)-x

4、_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=minPQi,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);%Z_g_wang=zeros(14,14);%根据平面几何关系计算出由Qi,Qj 点的深度线性外推的Pij 的深

5、度 Z_p for i=1:14;for j=1:14;if(Qi_Qj(i,j)=0)|(z_Ft(i)=z_Ft(j)Z_p(i,j)=z_Ft(i);else if(P_Qi(i,j)=P_Qj(i,j)&(Qi_Qj(i,j)z_Ft(j)Z_p(i,j)=z_Ft(j)-(abs(z_Ft(i)-z_Ft(j).*P_Qij(i,j)./Qi_Qj(i,j);else if(P_Qi(i,j)P_Qj(i,j)&(Qi_Qj(i,j)z_Ft(i)Z_p(i,j)=abs(z_Ft(i)-(abs(z_Ft(j)-z_Ft(i).*P_Qij(i,j)./Qi_Qj(i,j);els

6、e if(P_Qj(i,j)=P_Qi(i,j)&(Qi_Qj(i,j)z_Ft(i)Z_p(i,j)=z_Ft(j)+(abs(z_Ft(j)-z_Ft(i).*P_Qij(i,j)./Qi_Qj(i,j);else if(P_Qi(i,j)=P_Qj(i,j)&(Qi_Qj(i,j)z_Ft(j)Z_p(i,j)=z_Ft(i)+(abs(z_Ft(i)-z_Ft(j).*P_Qij(i,j)./Qi_Qj(i,j);else if(P_Qi(i,j)=Qi_Qj(i,j)&(P_Qj(i,j)z_Ft(j)Z_p(i,j)=z_Ft(j)+(abs(z_Ft(i)-z_Ft(j).*P_

7、Qj(i,j)./Qi_Qj(i,j);else if(P_Qi(i,j)=Qi_Qj(i,j)&(P_Qj(i,j)z_Ft(i)Z_p(i,j)=z_Ft(i)+(abs(z_Ft(j)-z_Ft(i).*P_Qi(i,j)./Qi_Qj(i,j);else Z_P(i,j)=1;end 名师资料总结-精品资料欢迎下载-名师精心整理-第 2 页，共 5 页 -end end end end end end end end%计算 G 点深度 Z_g Z_g_1=ones(14,14)./(G_Pij.2+G_Qij.2+Qi_Qj.2);Z_g_2=Z_p./(G_Pij.2+G_Qij.2

8、+Qi_Qj.2);for i=1:14 for j=1:14 Z_g_3(i,j)=sum(sum(Z_g_1);end end Z_g(m,k)=sum(sum(Z_g_2./Z_g_3);end end end%程序一：插值并作海底曲面图x1=75:1:200;y1=-50:1:150;x1,y1=meshgrid(x1,y1);z1=-griddata(G_x,G_y,Z_g,x1,y1,v4);meshc(x1,y1,z1)xlabel(x),ylabel(y),zlabel(z)figure;%程序二：插值并作出水深小于5 的海域范围。x1=75:1:200;y1=-50:1:150;x1,y1=meshgrid(x1,y1);z1=-griddata(G_x,G_y,Z_g,x1,y1,v4);%插值z1(z1=-5)=nan;%将水深大于5 的置为 nan，这样绘图就不会显示出来meshc(x1,y1,z1)xlabel(x),ylabel(y),zlabel(z)%产生的 G 点中有重复点，运行时有警告名师资料总结-精品资料欢迎下载-名师精心整理-第 3 页，共 5 页 -名师资料总结-精品资料欢迎下载-名师精心整理-第 4 页，共 5 页 -名师资料总结-精品资料欢迎下载-名师精心整理-第 5 页，共 5 页 -