interp2_lq.m

来自「关于样条插值均匀与非均匀的示例」· M 代码 · 共 22 行

figure;
x = 0:400:2800;
y = 0:400:2400;
Z = [1180 1320 1450 1420 1400 1300 700 900;
    1230 1390 1500 1500 1400 900 1100 1060;
    1270 1500 1200 1100 1350 1450 1200 1150;
    1370 1500 1200 1100 1550 1600 1550 1380;
    1460 1500 1550 1600 1550 1600 1600 1600;
    1450 1480 1500 1550 1510 1430 1300 1200;
    1430 1450 1470 1320 1280 1200 1080 940]

[X Y] = meshgrid(x, y);
[X1 Y1] = meshgrid(0:100:2800, 0:100:2400);

Z1 = interp2(X, Y, Z, X1, Y1, 'spline');
mesh(X1, Y1, Z1);
xlabel('X');
ylabel('Y');
figure;
[c, h] = contour(X1, Y1, Z1, [500:100:1600]);
clabel(c, h);