第三次作業(yè)-插值

1、第三次作業(yè)：插值第三次作業(yè)：插值數(shù)據(jù)文件inter2.mat文件為一曲面上的點(diǎn)（xyz）的坐標(biāo)可用load函數(shù)調(diào)入后，用mesh繪制其圖形。要求：（1）利用這些點(diǎn)的信息，分片段線型插值、分片樣條插值，來計(jì)算該曲面上位于相鄰兩插值節(jié)點(diǎn)的中點(diǎn)、三分點(diǎn)的坐標(biāo)處的坐標(biāo)（2）根據(jù)計(jì)算出的點(diǎn)和題設(shè)給定點(diǎn)的數(shù)據(jù)信息，選用相應(yīng)軟件畫出其曲線圖；（3）計(jì)算到插值后的曲面的距離為常數(shù)（比如是0.1個(gè)點(diǎn)位）的其他兩個(gè)曲面點(diǎn)的集合。中點(diǎn)線性插值中點(diǎn)線性插值cl

2、c%二分點(diǎn)線性插值loadinter2.matxyz[mn]=size(x)x0=[]y0=[]z0=[]fj=1:n1x0=[x0(x(:j)x(:j1))2]y0=[y0(y(:j)y(:j1))2]endfi=1:mfj=1:n1z0(ij)=z(ij)(z(ij1)z(ij))(y0(ij)y(ij))(y(ij1)y(ij))endend%subplot(211)figure(1)mesh(xyz)holdontitle()X

3、0=[]Y0=[]Z0=[]fi=1:mfj=1:8X0(i2j1)=x(ij)Y0(i2j1)=y(ij)Z0(i2j1)=z(ij)ifj~=8X0(i2j)=x0(ij)Y0(i2j)=y0(ij)Z0(i2j)=z0(ij)endend中點(diǎn)樣條插值中點(diǎn)樣條插值clc%二分點(diǎn)樣條插值loadinter2.matxyz[mn]=size(x)x0=[]y0=[]z0=[]fj=1:n1x0=[x0(x(:j)x(:j1))2]y0=

4、[y0(y(:j)y(:j1))2]endfi=1:mfj=1:n1z0(ij)=griddata(xyzx0(ij)y0(ij)cubic)endendfigure(1)mesh(x0y0z0)X0=[]Y0=[]Z0=[]fi=1:mfj=1:8X0(i2j1)=x(ij)Y0(i2j1)=y(ij)Z0(i2j1)=z(ij)ifj~=8X0(i2j)=x0(ij)Y0(i2j)=y0(ij)Z0(i2j)=z0(ij)enden