首页 > 其他 > 详细

matlab-画地形图

时间:2019-08-20 14:32:43      阅读:190      评论:0      收藏:0      [点我收藏+]

1.画三维图

之前画曲面的三维图,运用z=x2+y2 算出z和Z,如果是给出数据的地形则没办法用公式算,为此,引入插值自动造出地形的坐标。

拟合和插值的区别:插值是必须要过点,曲线可以不光滑;拟合则是可以不过点,曲线要光滑。

例如依旧画出z=x2+y,但只给出固定点,插值出Z矩阵。

clc;
clear all;

x=[-3 -2 -1 0 1 2 3];
y=[-3 -2 -1 0 1 2 3];
z=[18 8 2 0 2 8 18];
xminn=min(x)
xmaxx=max(x)
yminn=min(y);
ymaxx=max(x);
zminn=min(z);
zmaxx=max(z);

[X,Y,Z]=griddata(x,y,z,linspace(xmaxx,xminn)‘,linspace(yminn,ymaxx),‘v4‘);
%v4是一个插值,理解成自动造出Z矩阵即可,第4第5个需要其中一个转置
surf(X,Y,Z)

技术分享图片

 

再用公式法求Z看看画出的图对比是否相同

 

clear all;
clc;
x=-3:0.05:3;
y=-3:0.05:3;
z=x.^2 + y.^2;
[X,Y]=meshgrid(x,y);
Z=X.^2 + Y.^2;
surf(X,Y,Z);

 

技术分享图片

显然,画出的图形状几乎一样,因为点与点的间距使得图形有些不同,插值有很多种,v4插值适用,就用v4吧,可以尝试其他插值是否得出想要的结果。

 

2.画四维图

除了长宽高,怎么可能画得出其他维度,往往需要体现地形上的温度、浓度等指标,这类“第四维度”一般用颜色体现出来。

clc
%手动导入表中数据后,对x、y、z、r赋值

x=one(:,1);
y=one(:,2);
z=one(:,3);%高度
r=one(:,11);%浓度

[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x))‘,linspace(min(y),max(y)),‘v4‘);%插值,造出Z

variables=one(:,1:3);
F=TriScatteredInterp(variables,r,‘linear‘);
R=F(X,Y,Z);
%以上3句,理解为造出R,这个R和前三个东西都有关系,插值用的是linear,此处不能用v4
surf(X,Y,Z,R)
title(‘4D‘);
colormapeditor %颜色编辑器,颜色较深可调色,嫌麻烦就随便弄1,工具-标准颜色图-hsv

技术分享图片技术分享图片

为了体现“第四维度”,一般需要调出颜色条。

 

matlab-画地形图

原文:https://www.cnblogs.com/shoulinniao/p/11382305.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!