基于离散数据点的曲面插值Matlab示例代码

共包含100个离散数据点，其数值范围限定在-0.05至0.15之间，每个数据点的幅值(E)及其对应的坐标(xd,yd)均被存储于结构体变量data中，相关数据文件以matlab.mat格式进行保存。各数据点在空间中的分布情况如以下图表所示：

数据导入操作如下：data=importdata('E:/matlab.mat')；

变量pnum的计算方式为：pnum=size(data,2)-1；

初始化三个空数组xdata、ydata以及edata；

通过循环结构依次处理每个数据点，具体实现为：

for i=1:pnum

xdata=[xdata,data(1,i).xd];

ydata=[ydata,data(1,i).yd];

edata=[edata,data(1,i).E];

end

随后绘制散点图，命令为：plot(xdata,ydata,'ro')；

新建图形窗口，定义ti变量为0到100之间的偶数序列，即ti = 0:2:100；

利用meshgrid函数生成网格坐标(cx,cy)，并采用三次插值方法对原始数据进行插值运算，得到cz=griddata(xd