写一个matlab代码，计算两个圆的交点

以下是一个简单的 MATLAB 代码，用于计算两个圆的交点：

% 已知参数
c1 = [0, 0];    % 第一个圆心坐标
r1 = 2;         % 第一个圆的半径
c2 = [2, 0];    % 第二个圆心坐标
r2 = 1.5;       % 第二个圆的半径

% 计算两个圆的交点
d = norm(c2 - c1);  % 两个圆心的距离
a = (r1^2 - r2^2 + d^2)/(2*d);
h = sqrt(r1^2 - a^2);
p = c1 + a*(c2 - c1)/d;
x1 = p(1) + h*(c2(2) - c1(2))/d;
y1 = p(2) - h*(c2(1) - c1(1))/d;
x2 = p(1) - h*(c2(2) - c1(2))/d;
y2 = p(2) + h*(c2(1) - c1(1))/d;

% 输出结果
disp(['两个圆的交点为：(', num2str(x1), ', ', num2str(y1), ') 和 (', num2str(x2), ', ', num2str(y2), ')']);

解释一下代码：

首先，我们定义了已知参数，即两个圆的圆心坐标和半径。然后，我们使用勾股定理计算两个圆心的距离d，并使用已知参数计算出两个圆的交点。具体地，我们首先计算出一个点p，它是第一个圆心到交点的向量与第二个圆心到交点的向量之间的夹角平分线上的点。然后，我们使用勾股定理和三角函数计算出交点的坐标。

请注意，这个代码只计算两个圆的交点，如果两个圆没有交点，或者有无限个交点，则需要进行适当的调整来解决这个问题。