find first time where data crosses 1

Question

Benjamin 2019-4-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/455824-find-first-time-where-data-crosses-1

评论： Mark Sherstan 2019-4-11

I have data stored in variables r and f. At some point the variable f goes above 1. How can I return the value of r the FIRST time f goes above 1? Interpolation of the r-value to find where f goes above 1 may be needed to be more precise. Any help would be greatly appreciated!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Mark Sherstan 2019-4-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/455824-find-first-time-where-data-crosses-1#answer_370231

在 MATLAB Online 中打开

This example should help you with your data set

x = -2*pi:pi/12:2*pi;
y = 2*sin(x);
idx = find(y>1);
solution = idx(1)
fprintf("Solution at index %0.0f, x = %0.2f, y = %0.2f\n",solution,x(solution),y(solution))
figure(1)
hold on
plot(x,y)
plot(x(solution),y(solution),'*r')

2 个评论
显示无隐藏无

Benjamin 2019-4-11

编辑：Benjamin 2019-4-11

在 MATLAB Online 中打开

Does this actually interpolate though? I have this:

    ix = find(y > 1, 1, 'first');
    ix_r(count) = r(:,ix);
    ix_r = ix_r';

but the problem is that it just returns r for when f exceeds the threshold of 1. How can I interpolate to find exactly where it would cross 1? then interpolate r

Mark Sherstan 2019-4-11

You can use the built in function interp1 as you have a region of interest.

请先登录，再进行评论。

Answer 2

Star Strider 2019-4-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/455824-find-first-time-where-data-crosses-1#answer_370236

在 MATLAB Online 中打开

I created ‘zci’ to do just that:

zci = @(v) find(v(:).*circshift(v(:), [-1 0]) <= 0);                    % Returns Approximate Zero-Crossing Indices Of Argument Vector

Then either use the interp1 function to do the interpolation, or create your own linear interpolation function:

x = 0:0.2:5;                                                            % Create Data
y = x.^2 - rand(size(x))*5;                                             % Create Data
zci = @(v) find(v(:).*circshift(v(:), [-1 0]) <= 0);                    % Returns Approximate Zero-Crossing Indices Of Argument Vector
yzi = zci(y-1);
lintrp = [[x(yzi(1));x(yzi(1)+1)] ones(2,1)] \ [y(yzi(1));y(yzi(1)+1)]  % Linear Interpolation
xint = (1-lintrp(2))/lintrp(1);                                         % y=1 X-Value
figure
plot(x, y,    xint, 1, '+')
grid