在圆上旋转直线- Matlab

% 定义圆的参数
x0 = 0;  % 圆心x坐标
y0 = 0;  % 圆心y坐标
r = 5;   % 圆的半径

% 定义直线的参数
x1 = -3;  % 直线上第一个点的x坐标
y1 = 3;   % 直线上第一个点的y坐标
x2 = 3;   % 直线上第二个点的x坐标
y2 = -3;  % 直线上第二个点的y坐标

% 计算直线与圆的交点
syms x y
eq1 = (x - x1)*(y2 - y1) - (y - y1)*(x2 - x1);  % 直线的方程
eq2 = (x - x0)^2 + (y - y0)^2 - r^2;            % 圆的方程
sol = solve(eq1, eq2, x, y);                    % 求解交点坐标

% 绘制旋转后的直线
theta = linspace(0, 2*pi, 100);  % 旋转角度范围
x_circle = x0 + r*cos(theta);    % 圆上点的x坐标
y_circle = y0 + r*sin(theta);    % 圆上点的y坐标

figure;
hold on;
plot(x_circle, y_circle, 'b');   % 绘制圆
plot([sol.x(1), sol.x(2)], [sol.y(1), sol.y(2)], 'r');  % 绘制旋转后的直线
axis equal;
xlabel('x');
ylabel('y');
title('Rotation of Line on Circle');
legend('Circle', 'Rotated Line');