% 定义起点和终点坐标 x1 = 12; y1 = 15; x2 = 170; y2 = 155; % 计算中点坐标 xm = (x1 + x2) / 2; ym = (y1 + y2) / 2; % 计算平移向量 dx = -xm + 20; dy = -ym + 35; % 定义平移矩阵 T1 = [1 0 dx; 0 1 dy; 0 0 1]; % 进行平移变换 P1 = [x1 y1 1] * T1; P2 = [x2 y2 1] * T1; % 计算旋转矩阵 theta = -45; T2 = [cosd(theta) -sind(theta) 0; sind(theta) cosd(theta) 0; 0 0 1]; % 进行旋转变换 P3 = P1 * T2; P4 = P2 * T2; % 计算反向平移向量 dx2 = xm - 20; dy2 = ym - 35; % 定义反向平移矩阵 T3 = [1 0 dx2; 0 1 dy2; 0 0 1]; % 进行反向平移变换 P5 = [P1(1) P1(2) 1] * T3; P6 = [P2(1) P2(2) 1] * T3; % 绘制原始直线 hold on plot([x1 x2], [y1 y2], 'b'); % 绘制平移后的直线 plot([P1(1) P2(1)], [P1(2) P2(2)], 'r'); % 绘制旋转后的直线 plot([P3(1) P4(1)], [P3(2) P4(2)], 'm'); % 绘制反向平移后的直线 plot([P5(1) P6(1)], [P5(2) P6(2)], 'g'); % 设置标题 title('直线变换示例'); 改成绕（20，35）旋转

要将旋转中心从原点移动到 (20, 35) 的话，需要进行如下步骤：

1. 计算线段的中点坐标。

2. 计算平移向量，使中点移动到目标位置 (20, 35)。

3. 定义平移矩阵 T1，将线段平移。

4. 计算旋转矩阵 T2，将线段绕中点旋转。

5. 计算反向平移向量，使中点移回原来的位置。

6. 定义反向平移矩阵 T3，将线段移回原来的位置。

7. 将起点和终点坐标乘以平移矩阵 T1 得到新的坐标。

8. 将新的坐标乘以旋转矩阵 T2 得到旋转后的坐标。

9. 将旋转后的坐标乘以反向平移矩阵 T3 得到反向平移后的坐标。

修改后的代码如下：

% 定义起点和终点坐标
x1 = 12; y1 = 15; x2 = 170; y2 = 155;

% 计算中点坐标
xm = (x1 + x2) / 2; ym = (y1 + y2) / 2;

% 计算平移向量
dx = 20 - xm; dy = 35 - ym;

% 定义平移矩阵
T1 = [1 0 dx; 0 1 dy; 0 0 1];

% 进行平移变换
P1 = [x1 y1 1] * T1;
P2 = [x2 y2 1] * T1;

% 计算旋转矩阵
theta = -45;
T2 = [cosd(theta) -sind(theta) 0; sind(theta) cosd(theta) 0; 0 0 1];

% 定义反向平移向量
dx2 = xm - 20; dy2 = ym - 35;

% 定义反向平移矩阵
T3 = [1 0 dx2; 0 1 dy2; 0 0 1];

% 将起点和终点坐标乘以平移矩阵 T1，再乘以旋转矩阵 T2，最后乘以反向平移矩阵 T3
P3 = P1 * T2 * T3;
P4 = P2 * T2 * T3;

% 绘制原始直线
hold on
plot([x1 x2], [y1 y2], 'b');

% 绘制平移后的直线
plot([P1(1) P2(1)], [P1(2) P2(2)], 'r');

% 绘制旋转后的直线
plot([P3(1) P4(1)], [P3(2) P4(2)], 'm');

% 绘制反向平移后的直线
plot([x1+dx2 x2+dx2], [y1+dy2 y2+dy2], 'g');

% 设置标题
title('直线变换示例');

其中，cosd和sind是MATLAB中的函数，表示求角度的余弦和正弦值。