积分返回预定义区域所需的y值

Question

鉴于以下情况：

>

具有X(time[s])和Y的时间序列（此处放电[mΩ/s])

值v1(y)（与x的整数相同）。在这种情况下，一个小体积[m].

v1必须放置在峰值周围。

我认为这必须是一个迭代过程，其中拟合标准（和精确程度）也必须由用户设置。

直到现在，我还没有找到开始的方法；除了纯粹的整合。

Answer 1

您可以通过减少一些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)] )

产出：