How to find intersection points between plot and straight line

Question

Vaultec 2015-3-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/184494-how-to-find-intersection-points-between-plot-and-straight-line

评论： Star Strider 2015-3-23

I tried to find the intersection points of the attached figure using the following code provided to me which worked for a different figure.

h1c = get(gca, 'Children');
Xdc = get(h1c, 'XData');
Ydc = get(h1c, 'YData');
maxlen = max(cell2mat(cellfun(@max, cellfun(@size, Xdc, 'Uni',0),'Uni',0)))
Xd2 = cell2mat(Xdc(2));
Yd2 = cell2mat(Ydc(2));
Xd = Xd2;
Yd1 = cell2mat(Ydc(1));
Yd = [Yd1(1)*ones(size(Xd2)); Yd2];
% CALCULATIONS: 
Ydn = diff(Yd, [], 1);              % Subtract line from curve to create zero-crossings
Zx = circshift(Ydn, [0 1]) .* Ydn;  % Use circshift to detect them
Zxi = find(Zx < 0);                 % Their indices
for k1 = 1:length(Zxi)              % Use interp1(Y,X,0) to get line intercepts as Xzx
    Xzx(k1) = interp1([Ydn(Zxi(k1)-1) Ydn(Zxi(k1))], [Xd(1,Zxi(k1)-1) Xd(1,Zxi(k1))], 0);
end
% PLOT ZERO-CROSSINGS ON FIGURE TO CHECK: 
hold on
plot(Xzx, repmat(Yd(1,1),1,length(Xzx)), '*r')
hold off

It returns with the error

Error using cellfun
Input #2 expected to be a cell array, was double instead.

And Im not quite sure how to fix it

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2015-3-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/184494-how-to-find-intersection-points-between-plot-and-straight-line#answer_172289

That is my code, so I sort of remember what I was doing when I wrote it here nine months ago (in R2014a). What version of MATLAB are you using?

I’ll do my best to help, but I can’t figure out what you want to do. I don’t see a reference line in ‘Curve figure.fig’.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Star Strider 2015-3-22

编辑：Star Strider 2015-3-22

My pleasure.

I opened R2013a to see if there was a problem or some sort of incompatibility with R2015a, since I couldn’t see the horizontal (black) line in R2015a, but could in R2013a. The code worked with your previous figure (‘Figure(7).fig’) because I got two values each for both 'XData' and 'YData' with it, one for the line and one for the curve.

With the figure in this (current) post (‘Curve figure.fig’), I’m only getting the horizontal line. I’m not getting any data for the blue curve, either in R2013a or R2015a. For some reason, the blue curve just isn’t present in the figure data. I checked both ‘Curve figure.fig’ you posted to see if the second one was an update of the first, and to see if I could recover the blue curve from it, but they both appear to be the same.

Please check it to see if you’re getting the same result. (Compare it with the figures from your previous post. I linked to it, so you can click on that link to open it.) Then re-plot the current figure, being sure the data for both the line and the curve are present in it, and re-post it. At the very least, I need the data for the blue curve. I can create a constant line at 1.7.

Vaultec 2015-3-23

Thank you very much,

So the reason your previous code didnt work was it because of a version issue?

Star Strider 2015-3-23

My pleasure.

I don’t believe it was a version issue, but that for some reason the blue curve didn’t appear as 'XData' or 'YData' when I tried to get them from the figure. (The half-power line may have overwritten them somehow, without erasing the blue curve.) I didn’t explore that in my latest code, since I had all the data in the spreadsheet file. You may want to, because a .fig file is actually a very good way to transport both the data and the concept you want to (literally) illustrate.

请先登录，再进行评论。