函数逼近与曲线(面)拟合的MATLAB程序

运行后屏幕显示三元线性拟合值及其图形（略）.

例7.9.2 设节点(X,Y,Z,W)中的X,Y和Z分别是在区间[ 3,3]和[ 2.5,3.5],Y=Z上的15个随机数，W是函数w 2 xe在(X,Y,Z)的值，拟合点(xi,yi,zi)中的xi=-3:0.2:3， yi=-2.5:0.2:3.5，zi=yi， 用'linear'方法计算拟合数据的值，并作 x2 y2 z2其图形.

解 输入程序

>> x=rand(15,1); y=rand(15,1);

X1=-3+(3-(-3))*x;

Y1=-2.5+(3.5-(-2.5))*y;Z1=Y1;

[X,Y,Z] = meshgrid(X1,Y1,Z1);

W=2+X.* exp(-X.^2 - Y.^2- Z.^2);

xi=-3:0.2:3; yi=-2.5:0.2:3.5; zi=yi;

[X2,Y2,Z2]=meshgrid(xi,yi,zi);

W1=griddata3(X, Y, Z, W, X2,Y2,Z2,'linear');

slice(X2,Y2,Z2,W1,[-1 0 1.5],2,[-2 3]),

shading flat,lighting flat,

xlabel('x'), ylabel('y'), zlabel('z'),

title('被拟合函数W=2+X exp(-X^2 - Y^2- Z^2)');

hold on,colorbar('horiz'), view([-3 5])

运行后屏幕显示三元线性拟合值及其图形（略）.