计算点到直线方程的投影点坐标

点 $(x_0, y_0)$到直线 $ax + by + c = 0$的垂足的坐标 $(x_1, y_1)$ 满足以下条件:

- 点 $(x_1, y_1)$ 在直线上
- 点 $(x_1, y_1)$ 到点 $(x_0, y_0)$ 的距离是最小的

根据定义可知，点 $(x_0, y_0)$ 到直线的距离为 $$\frac{|ax_0 + by_0 + c|}{\sqrt{a^2 + b^2}}$$ 。点 $(x_1, y_1)$ 到点 $(x_0, y_0)$ 的距离为 $$\sqrt{(x_1 - x_0)^2 + (y_1 - y_0)^2}$$ 。

所以我们可以设这两个距离的平方和为 $$d(x_0, y_0) = \left(\frac{|ax_0 + by_0 + c|}{\sqrt{a^2 + b^2}}\right)^2 + (x_1 - x_0)^2 + (y_1 - y_0)^2$$ ，我们希望最小化 $$d(x_0, y_0)$$ 。

于是我们得到优化问题:

$$\min d(x_1, y_1) = \left(\frac{|ax_0 + by_0 + c|}{\sqrt{a^2 + b^2}}\right)^2 + (x_1 - x_0)^2 + (y_1 - y_0)^2$$

对 $x_0$ 和 $y_0$ 求导，得到:

$$\begin{aligned} \frac{\partial d(x_0, y_0)}{\partial x_0} &= 2\left(\frac{ax_0 + by_0 + c}{a^2 + b^2}\right)a + 2(x_1 - x_0) = 0 \\ \frac{\partial d(x_0, y_0)}{\partial y_0} &= 2\left(\frac{ax_0 + by_0 + c}{a^2 + b^2}\right)b + 2(y_1 - y_0) = 0 \end{aligned}$$

解得:

$$ x_1 = x_0 - \frac{a(ax_0 + by_0 + c)}{a^2 + b^2} $$
$$ y_1 = y_0 - \frac{b(ax_0 + by_0 + c)}{a^2 + b^2} $$

所以得到：
点到直线距离的公式：
设点 P(x0, y0)，直线方程为 Ax + By + C = 0，则点 P 到直线的距离为：
dist = |Ax0 + By0 + C| / sqrt(A^2 + B^2)
垂足的坐标为(x1, y1)，我们可以用如下的方法求出坐标：
x1 = x0 - dist * A / sqrt(A^2 + B^2)
y1 = y0 - dist * B / sqrt(A^2 + B^2)
下面是一个完整的c++实现例子：

#include <iostream>
#include <cmath>
using namespace std;
int main() {
// 点的坐标
double x0 = 1.0, y0 = 1.0;
// 直线方程参数
double A = 1.0, B = 1.0, C = 0.0;
// 计算点到直线的距离
double dist = fabs(A * x0 + B * y0 + C) / sqrt(A * A + B * B);
// 计算投影点的坐标
double x1 = x0 - dist * A / sqrt(A * A + B * B);
double y1 = y0 - dist * B / sqrt(A * A + B * B);
// 输出结果
cout << "投影点的坐标为：(" << x1 << ", " << y1 << ")" << endl;
return 0;
}