matlab中使用[cos(b) sin(b) 0；-sin(b) cos(b) 0；0 0 1]进行旋转

            function g=rgbrotate(f,a)
            b=a*pi/180;%数字转化为角度。
            [x,y,h]=size(f);
            X=round(x*abs(cos(b))+y*abs(sin(b)));
            Y=round(y*abs(cos(b))+x*abs(sin(b)));%旋转图像后得到的新高度，“round()函数四舍五入“
            g=zeros(X,Y,h);
            M1=[1 0 -0.5*X;0 1 -0.5*Y;0 0 1 ];     %把中心平移到原点
            M2=[cos(b) sin(b) 0;-sin(b) cos(b) 0;0 0 1];  %逆时针旋转b度
            M3=[1 0 0.5*x;0 1 0.5*y;0 0 1];              %把中心从原点反平移回去。
            for i=1:X
                for j=1:Y
                    temp=M3*M2*M1*[i;j;1];
                    yy=temp(2,1);
                    xx=temp(1,1);
                    yy=round(yy);           %y四舍五入取整
                    xx=round(xx);
                    if(xx>=1&&xx<=y)&&(yy>=1&&yy<=x)  %判断的得到的(x,y)点是否在原图像上
                        g(i,j,:)=f(xx,yy,:);  %将原图像的像素点赋值给对应的旋转后图像上的点
                    end
                end
            end
            imshow(uint8(g));