0

我在水柱的每个深度都有一些深度测量值和相应的区域:

depth = [0,2,4,6,7,9,10];
Area = [2000,1890,1873,1654,1523,1433,1200];

figure(1);
plot(Area,depth);set(gca,'ydir','reverse');

我想计算指示的每个深度之间的水量,然后从中计算整个水量。首先,我插入了深度和面积值以提高计算的准确性:

dz = 0.1;
newD = min(depth):dz:max(depth);
newA = interp1(depth,Area,newD);

figure(2);
plot(newA,newD);set(gca,'ydir','reverse');

我从哪里开始计算每个深度之间的水量?

4

2 回答 2

1

要获得体积,您将面积乘以高度。因此,由于您的 dz 是恒定的,要获得所选深度的水量,您需要乘以

newA*dz

要获得总成交量,您可以将其相加

sum(newA*dz)

您的计算是在假设水域随深度线性变化的情况下完成的(这就是您使用 interp1 实现的)。我相信您可以设法以分析方式集成线性函数,而不是使用您当前的方法。您可以轻松获得准确的数字。

于 2012-09-17T13:16:23.163 回答
0

您实际上想将面积乘以总深度,而不是深度步长:

volume = newA.*newD;
totalVolume = sum(volume);
于 2012-09-17T16:52:05.783 回答