使用八度/ matlab中的仿射矩阵的旋转/仿射矩阵函数

Question

我可以创建一个以蓝色显示的图表，这是正确的，但是当我尝试将其围绕点（0,0）旋转90度时，会创建红色图。它应该绕点（0,0）旋转大约90度。任何人都知道我的函数/代码的哪个部分导致了问题？

请参阅以下示例代码中的函数：

    %affine matrix rotation about a point rtaffine[theta,rpx,rpy,sigx,sigy,sigz]
%crx cry rpx rpy represent center of rotation
function [rotsig,theta,crpx,crpy,sigx,sigy,sigz] = rtaffine(theta,rpx,rpy,sigx,sigy,sigz) 
    rotsig=[];
    %affinematrix=[];
    siga=[sigx;sigy;sigz];

    r00 = cosd(theta); r01 = -sind(theta); r10 = sind(theta); r11 = cosd(theta);
    affinematrix = [r00, r01, rpx(1,1) - r00*rpx(1,1) - r01*rpy(1,1);...
    r10, r11, rpy(1,1) - r10*rpx(1,1) - r11*rpy(1,1);...
    0, 0, 1];   
    rotsig=affinematrix*siga; %new affine matrix
end


%radial arms
t = linspace(0,2*pi,500);
r=e^0.3063489*t;
x = r.*cos(t);
y = r.*sin(t);
plot(x,y)

hold on

%rotation
theta=90;
z = ones(size(y));
siga=[t;y;z];


rotsig=rtaffine(theta,0,0,siga(1,:),siga(2,:),siga(3,:));

plot(t(1,:),rotsig(1,:),'r-')

Answer 1

您正在传递t，并试图在{2}中使用t而不是x进行投标：

siga=[t;y;z];

应该是

siga=[x;y;z];

和

plot(t(1,:),rotsig(1,:),'r-')

应该是

plot(rotsig(1,:),rotsig(2,:),'r-')