旋转实验： a、利用imrotate，实现图像的旋转 b、不使用imrotate，定义文件myrotate.m，其中直接定义一个旋转矩阵，逐像素的实现图像的旋转。给出代码

a) 利用imrotate实现图像的旋转：

img = imread('image.png');
rotated_img = imrotate(img, angle, 'nearest', 'crop');
imshow(rotated_img);

其中，angle是旋转角度。

b) 不使用imrotate，定义文件myrotate.m实现图像的旋转：

function rotated_img = myrotate(img, angle)
% 输入：img - 输入图像
%      angle - 旋转角度（逆时针为正方向）
% 输出：rotated_img - 旋转后的图像

% 将角度转换为弧度
theta = angle * pi / 180;

% 定义旋转矩阵
R = [cos(theta) -sin(theta);
     sin(theta) cos(theta)];

% 获取输入图像的尺寸
[h, w, ~] = size(img);

% 计算旋转后的图像尺寸
[new_w, new_h] = deal(round(abs(h * sin(theta)) + abs(w * cos(theta))), ...
                      round(abs(h * cos(theta)) + abs(w * sin(theta))));
% 定义旋转后的图像
rotated_img = uint8(zeros(new_h, new_w, 3));

% 计算旋转后的原点位置
ori = [round((new_w - w) / 2), round((new_h - h) / 2)];
ori = ori + [round(w / 2), round(h / 2)] * R;

% 逐像素进行旋转
for x = 1:new_w
    for y = 1:new_h
        coord = [x, y] - ori;  % 坐标变换
        rot_coord = round(coord * R);  % 坐标旋转
        rot_coord = rot_coord + ori;  % 坐标还原
        if rot_coord(1) >= 1 && rot_coord(2) >= 1 && rot_coord(1) <= w && rot_coord(2) <= h
            rotated_img(y, x, :) = img(rot_coord(2), rot_coord(1), :);
        end
    end
end

使用方式：

img = imread('image.png');
rotated_img = myrotate(img, angle);
imshow(rotated_img);