How to get a line from two points collected using ginput

OBJECTIVE. I have a curve, I want to select two points on the curve, fit a line between two points (overlaying on the main curve) and then find the intercept and slope of the fitted line.

MY METHODOLOGY. I can get two points using ginput function, but after, I am having trouble to find the points between the points using "find" function. also, I am unable to fit or plot the line of points from find or ginput function.

[Rik]

Why do you need find? You can calculate the slope and intercept from the two points directly.

Also, what exactly have you tried when you wanted to plot the second line? Did you forget to use the hold function?

[Faisal]

I tried polyfit to fit a line after ginput.

how to find the intercept from two points?

Answer 1

[John D'Errico]

You aparently know how to use polyfit. But it seems you do not understand what it returns.

x = 1:5;
y = [2 3 5 7 11];
P1 = polyfit(x,y,1)

P1 = 1×2

2.2000 -1.0000

The result is a vector. The first element is the slope of the vector. The second element is the y-intercept.

plot(x,y,'o')
hold on
xpl = [0,5]

xpl = 1×2

0 5

ypl = polyval(P1,xpl);
plot(xpl,ypl,'r-')

As you can see, I extended the line down to x==0, where you should see it crosses the y-axis at y==-1, the y-intercept. So polyfit gives you exactly what you need. There is no need to "find" anything. You already have it.

[Faisal]

@John D'Errico Thank you!

I succeeded to use the similar script and got what I was looking for.

Thanks again!