matlab数值分析程序--高等数学,数值代数的matlab实现-文字版, matlab电子书, 和matlab 有关的电子书:

6.8.3 三元插值及其MATLAB程序

6.8.3  三元插值及其MATLAB程序

 

例6.8.7  设节点的坐标为,计算函数e在插值点处的三元线性插值,并作其图形.

解  输入程序

>> x=[-4,0,1,12];y=[-1,0,3,15];z=y;

[X,Y,Z]= meshgrid(x,y,z);

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

[xi,yi,zi] = meshgrid(-3:.25:10,-3:.25:3,-3:.25:13);

vi = interp3(X,Y,Z,V,xi,yi,zi),

slice(xi,yi,zi,vi,[-1 6 9.5],9,[-2 .2 9]),

shading flat,lighting flat

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

title('V=2+ x exp(-x^2 - y^2-z^2) 的三元线性插值图形')

hold on,colorbar('horiz'), view([-30  45])

运行后屏幕显示三元线性插值及其图形(略).

例6.8.8  取 n=10 ,作函数flow在插值点 处的三元三次样条插值及其图形.

解  输入程序

>> [x,y,z,v] = flow(10);

[xi,yi,zi] = meshgrid(.1:.25:10,-3:.25:3,-3:.25:3);

vi = interp3(x,y,z,v,xi,yi,zi,'spline');

slice(xi,yi,zi,vi,[2.5 7.5],0.1,[-1.2 1.5]),

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

title(' flow 的三元三次样条插值图形')

hold on,colorbar('horiz')

运行后屏幕显示三元三次样条插值及其图形(略).

欢迎转载,转载请注明来自一手册:http://yishouce.com/book/3/3060803.html
友情链接It题库(ittiku.com)| 版权归yishouce.com所有| 友链等可联系 admin#yishouce.com|粤ICP备16001685号-1