Question

Matlab求体积用mesh得出一个曲面,用什么方法求曲面和XOY面之间的体积比较精确,怎么计算?如果要曲面插值拟合和漆

Answer 1

假设10个离散点是A=[x1,y1,z1;x2,y2,z2;.;x10,y10,z10];

x=A(:,1);

y=A(:,2);

z=A(:,3);

N=100; %插值点数,自己确定,暂写100

x0=linspace(min(x),max(x),N);

y0=linspace(min(y),max(y),N);

[X,Y,Z]=griddata(x,y,z,x0',y0,'v4');%插值拟合曲面.当然,由于你没有说明xy平面上的条件,所以这里拟合的曲面Z=f(X,Y)定义在x0,和y0的矩形区间内；

figure;mesh(X,Y,Z);

%体积的近似值用下代码

L=sum(sum(abs(Z)));

S=abs((y0(2)-y0(1))*(x0(2)-x0(1)));

V=L*S; %体积