U =数据(:1);%(:1)首先意味着所有列
所以你是一个列向量
(z =数据:4);%(:,4)意味着所有第四列
z是一个列向量的大小一样
trapz (z, U, 3)
请求中的值z trapazoid求和沿着三维完成,使用你的内容作为第三维坐标。然而我们知道,z是一个列向量,沿着第三维是1,所以它的长度不匹配你的长度。
如果它发生,x, y, z与价值建立常规的长方体U每一点……你没有足够的U值。
所以我们必须推断出你没有正则网格:你有散点。你不能使用trapz意义做卷集成在这种情况下。
或许你可以重塑你的x, y, z形成规律的长方体。如果是这样那你需要知道他们是订单,你需要知道边际边的长度。你将重塑,等
marglen = [17 11 16];
X =重塑(X, marglen);
Y =重塑(Y, marglen);
Z =重塑(Z, marglen);
UU =重塑(U, marglen);
一旦完成将成为有意义的开始使用trapz各种维度,
trapz(挤压(X (:, 1, 1)), trapz(挤压(Y (1: 1)), trapz(挤压(Z (1 1:)), UU, 3), 2), 1)
这里挤()部分的目的是提取边缘惟一值为每个方向。
