点在任意直线上的投影点

点point (x0, y0) 投影到到直线line y = m*x+b 的 投影点 point(x1,y1).

x = linspace(-10,10,100);
m = 1;
b = 0;
y = m*x+b;
x0 = 1;
y0 = 10;
x1 = (m*y0+x0-m*b)/(m^2+1);
y1 = (m^2*y0+m*x0+b)/(m^2+1);

plot(x,y);
hold on;
scatter(x0,y0);
scatter(x1,y1);

xlim([-20 20]);
ylim([-20 20]);
axis equal