例86 根据下述函数生成一些较稀疏的网格数据 z =f(, g=(a-2 c) -y-ay 进行各种插值拟合,并比较拟合结果 绘制已知数据的网格图 >>[x,y]= meshgrid(-3:6:3;-2:4:2); z=(x.2-2*x)*exP(x.2-y.2-x.*y); surf(x,y, z) axis([-3,3,-2,2,-0.7,1.5]) 2/20/2021星期六, 高等应用数学问题的 MATLAB求解 2008-9-6,13:1107 Slide 1(of 11) 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 例 8.6 根据下述函数生成一些较稀疏的网格数据 进行各种插值拟合，并比较拟合结果 绘制已知数据的网格图：
默认插值算法进行插值: >>[x1,y1]= meshgrid(-3:2:3,-2:2:2); z1=interp2(x, y, z, x1, y1) surf(x1, y1, z1) axis([-3,3,-2,2,-0.7,1.51) 立方插值和样条插值: 6 > z1=interp2(x,y, z, x1, y1, 'cubic,) z2=interp2(x, y, z, x1, y1,'spline') surf(x1,y1, z1), figure; surf(x1, y1, z2) 2/20/2021星期六, 高等应用数学问题的 MATLAB求解 2008-9-6,13:1107 Slide 1(of 11) 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 默认插值算法进行插值： 立方插值和样条插值：
误差比较: >>z=(x1.2-2*x1).* expP(-x1.2-y1.2-x1.*y1); surf (x1, y1, abs(z-z1)), figure surf(x1, y1, abs(z-z2)) 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 误差比较：
81.4二维一般分 布数据的插值问题 → griddata()函数的调用格式: z=griddata(o, y0, 20, a, y, ',) 其中, eT0,90,x0是已知的二维样本点 ·,y是期望的插值位置 z表示插值的结果,维数和,y一致 ·‘ⅵ43是 MATLAB4.0版本中提供的插值算法 2/20/2021星期六, 2008-9-6,13:1107 Slide 1(of 11) 高等应用数学问题的 MATLAB求解 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 griddata() 函数的调用格式： 其中， 是已知的二维样本点 是期望的插值位置 表示插值的结果，维数和 一致 ‘v4’ 是MATLAB 4.0版本中提供的插值算法. 8.1.4 二维一般分 布数据的插值问题
例8.7 给定z=f(x,y)=(x2-2x)e y--9 在矩形区域x∈-3,3]y∈[-2,2]内随机生成 组样本点(x1yz) 生成一组x的值,以这些值为已知数据,用 griddata()进行插值处理,并误差分析 生成并显示已知的样本点: 2/20/2021星期六, 高等应用数学问题的 MATLAB求解 2008-9-6,13:1107 Slide 1(of 11) 东北大学信息学院
高等应用数学问题的MATLAB求解 东北大学信息学院 Slide 1 (of 11) 2/20/2021星期六, 2008-9- 6, 13:11:07 例 8.7 给定 在矩形区域 内随机生成 一组样本点 生成一组 的值，以这些值为已知数据，用 griddata()进行插值处理，并误差分析。 生成并显示已知的样本点：