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

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

2、c%二分點線性插值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中點樣條插值中點樣條插值clc%二分點樣條插值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