插值法画三维曲线图和等高图

三维曲线

% 清理命令和变量，以免内存不够
clc;clear;

% 读数据
% 数据格式：x,y,z
data = load('hv_re.txt');
x1 = data(:, 1);
y1 = data(:, 2);
z1 = data(:, 3);

% 抽稀，以免内存不够
count    = 1;     % 新变量计数器
interval = 1000; % 抽稀间隔
for i = 1 : interval : length(x1)
    x(count, 1) = x1(i, 1);
    y(count, 1) = y1(i, 1);
    z(count, 1) = z1(i, 1);
    count = count + 1;
end

%确定网格坐标（x和y方向的步长均取0.1）
[X,Y]=meshgrid(min(x):0.1:max(x),min(y):0.1:max(y)); 
%在网格点位置插值求Z，注意：不同的插值方法得到的曲线光滑度不同
Z=griddata(x,y,z,X,Y,'v4');
%绘制曲面
figure(1)
surf(X,Y,Z);
shading interp;
colormap(jet);
% view(0, 90);
colorbar;
print(gcf, '-djpeg', 'xyz.jpg'); % save picture
 ———————————————— 
版权声明：本文为CSDN博主「itsc」的原创文章，遵循CC 4.0 by-sa版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/shanchuan2012/article/details/72772191

等高图

%等高图
Z=griddata(x,y,price,linesapce(0,28654),linespace(0,18449),'v4');
contourf(X,Y,Z);
colorbar;