∈>>x=-3+6*rand(200,1) y=-2+4*rand(200,1); z=(x.2-2*x)*exP(-x.2-y.2-x.*y); plot(x,y,'x ), figure, plot(x,y, z,'x,) 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
使用‘ cubic3和‘v43算法: >>[x1,y1]= meshgrid(-3:2:3,-2:2:2); z1= griddata(x,y,z,x1,y1,’ cubic’); surf(x1, y1, z1) z2=griddata(x, y, z, x1,y1,2v4); figure; surf(x1, y1, z2) 比较: >>z0=(x1.2-2*x1)*exp(-x1.-2 y1.^2-x1.*y1); surf(x1, y1, abs(z0-z1)): figure surf(x1, y1, abs(z0-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 使用 ‘cubic’ 和 ‘v4’ 算法： 比较：
例88 →给定的样本点在x-y平面分布较均匀,现 在人为剔除某些点,表明已知数据分布 不均匀,这时再进行插值分析,观察插 值效果 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 例 8.8 给定的样本点在x-y平面分布较均匀，现 在人为剔除某些点，表明已知数据分布 不均匀，这时再进行插值分析，观察插 值效果
→剔除在以(-1,-1/2)点为圆心,以0.5为半 径的圆内的点 通>x=-3+6*rand(200,1); y=-2+4*rand(200,1); z=(x.2-2*x).*exp(-x.^2-y.2-x.*y) ii=find((x+1).2+(y+0.5).2>0.5^2); X=x(11);y=y(11);z=z(11); plot(x,y,x') t=[o:1:2*pi,2*p立];x0=-1+0.5*cos(t); yO=-0.5+0.5*sin(t); line(xo, yo) 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 剔除在以 点为圆心，以0.5为半 径的圆内的点
用新的样本点拟合出曲面: >>[x1,y1]= meshgrid(-3:2:3,-2:2:2); z1=griddata(x,y, z, x1, y1,'v4,) surf(x1, y1, z1) 误差分析: >>z0=(x1.^2-2*x1)*exP(-x1.2 y1.^2-x1.*y1); surf(x1, y1, abs(z0-z1)) axis([-3,3,-2,2,0,0.1]) 误差的等高线图: > contour(x1, y1, abs(z0-z1), 30) hold on, plot(x, y,'x); line(xo, yO 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 用新的样本点拟合出曲面： 误差分析： 误差的等高线图：