基于仿射矩阵的倍频程/ matlab旋转/仿射矩阵函数

Question

我可以创建一个蓝色的图形，这是正确的，但当我试图旋转它在90度左右，点(0,0)，图的红色是创建。它应该绕点(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，并试图针对t而不是x进行绘图：

siga=[t;y;z];

应该是

siga=[x;y;z];

和

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

应该是

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