鉴于以下内容:
x
(时间[s])和y
(此处放电[m³/ s])V1
(相同单位积分y
),小于所有x
的积分。在这种情况下,体积小[m³]。我想计算一下:
y
值y_V1
使得y = y_V1
线和y
曲线之间的积分等于V1
。
下图显示了这一点,橙色区域是V1
,我想要在y
轴上的圆圈值:V1
必须放在峰值附近。我认为这必须是一个迭代过程,其中必须由用户设置拟合标准(和精确度)。
到现在为止,我还没有办法开始;除了纯粹的整合。
想法是指定一个区域。应计算包围该区域的峰值左右y值。
编辑
如果应用了接受的答案,则结果如此。
您可以通过减少一些y
值来实现此目的,直到达到您的区域目标。有关详细信息,请参阅以下评论
% Input data
x = 0:0.01:pi;
y = sin(x);
target = 1; % Target area
yi = max( y ); % Initialise yi to be max possible y
dy = 0.001; % Step change in yi
Ai = 0; % Area each iteration
thresh = 0; % Threshold for stopping loop
while target - Ai > thresh && yi >= min(y)
yi = yi - dy;
ix = y >= yi;
% Approximate integral above the line
Ai = trapz( x(ix), y(ix) - yi );
end
% Plot
figure(1); clf; hold on
plot( x, y );
patch( x(ix), y(ix), [1,0.5,0.5], 'facealpha', 0.5 );
plot( x, ones(size(x))*yi, '--', 'linewidth', 2 )
xlim( [min(x),max(x)] )
输出: