% 已知四个角点坐标
x = [0.1, 0.1, 1.1, 1.1];
y = [0.1, 1.2, 0.1, 1.2];
z = [z1, z2, z3, z4]; % 假设已知四个角点的函数值
% 定义插值点的坐标
xi = linspace(0.1, 1.1, 100); % 插值点的x坐标范围
yi = linspace(0.1, 1.2, 100); % 插值点的y坐标范围
[X, Y] = meshgrid(xi, yi); % 生成插值点的网格坐标
% 计算插值曲面片内部各点处的函数值
Z = bilinearInterpolation(x, y, z, X, Y);
% 绘制插值之后的曲面片
surf(X, Y, Z);