Find interpolated intersection points between two curves

Question

Jorge Luis 2024-1-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2068266-find-interpolated-intersection-points-between-two-curves

编辑： Jorge Luis 2024-1-10

Hello, I have the following two curves: the blue curve (curve1.txt) and orange curve (curve2.txt).

For simplicity, I assumed the orange curve is a straight line. I would like to estimate all the interpolated intersection points based on the frequency values (x_coordinates.txt) which not necesary matches the specific value of the x_coordinates (e.g. one intersection could be 1.127 between the interpolated intersection of the two neighbouring points from both curves). I attached three .txt files for this example. I would appreciate the help. Thank you in advance for your reply.

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

请先登录，再进行评论。

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

Answer 1

Matt J 2024-1-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2068266-find-interpolated-intersection-points-between-two-curves#answer_1386511

There are a number of FEX contributions that do this, e.g.,

https://www.mathworks.com/matlabcentral/fileexchange/11837-fast-and-robust-curve-intersections?s_tid=srchtitle

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Jorge Luis 2024-1-10

编辑：Jorge Luis 2024-1-10

Thank you so much. It works perfectly on the task. It took 1min and 3 seconds to process just one components of the signal out of three components, though. I used a core i9 laptop with 2gb of ram memory. I will have to process around one hundred thoushand ground motion and each of them has three components along with other complex calculations. That means that just to estimate the intersection points it will take me 208 days.

y1=load('curve1.txt');

y2=load('curve2.txt');

x1=load(x_coordinates.txt);

x2=x1;

% finding the intersection points

[x0,y0,iout,jout] = intersections(x1,y1,x2,y2,'true');

请先登录，再进行评论。