dot_product = nx * (x2-x1) + ny * (y2-y1); rx = x2 - 2 * dot_product * nx; ry = y2 - 2 * dot_product * ny; 错误使用 * 内部矩阵维度必须一致。

非常抱歉，我的回答有误。这里 nx、ny、x2、y2 都是向量，无法直接进行矩阵乘法。正确的计算方式是使用点积（dot product）来计算反射光线的方向，代码如下：

theta = atan2(y, x);
nx = cos(theta);
ny = sin(theta);

dot_product = nx .* (x2-x1) + ny .* (y2-y1);
rx = x2 - 2 .* dot_product .* nx;
ry = y2 - 2 .* dot_product .* ny;

在这里，我们使用了点乘符号 `.*` 来表示对应元素相乘，而不是矩阵乘法的乘号 `*`。这样，代码就可以正确地计算出反射光线的方向了。非常抱歉给您带来困扰。