0

我需要提交一篇 matlab 论文，有一个非常奇怪的问题。我写了一个旋转 3d 向量的函数。我称之为旋转。

rotate 得到一个向量 v、一个角度 theta、一个角度 alpha 和一个标量 r。该函数将旋转将向量与 Z 轴创建的角度增加 theta，将其围绕 Z 轴旋转（增加与 X 轴创建的角度），然后将其拉伸 r。

例如：

v=(0,0,1)
theta= pi/2
alpha=0
r=1
rotate(v,theta,alpha,r) will return (0,1,0).

我的问题是，每当我调用函数旋转时，我都会收到一条错误消息：??? 未定义的函数或变量“u”。

==> 在 51 x = u(1) 处旋转时出错；

而且在我编写的整个代码中，没有一个名为“u”的函数或变量。

这是我的代码。基本上，我将给定的向量转换为球坐标，然后添加 theta 和 alpha。我认为它应该工作不？

function [output] = rotate(v,theta,alpha,r)
if(isnumeric(v))
    [i,j]=size(v);
    if(i>j)
        for i=1:3
            sum = sum+pow(v(i,1),2);
        end
        sum=sqrt(sum);
        output(1,1)=sum*r;
        output(2,1)=acos(v(3,1)/sum)+theta;
        output(3,1)=atan(v(2,1)/v(1,1))+alpha;
        if((output(2,1)>pi)||(output(2,1)<-1*pi))
            prompt={'Invalid values for second cordinate, more than pi or less than -pi'};
            return
        end
        if((output(3,1)>2*pi)||(output(3,1)<-2*pi))
            prompt={'Invalid values for third cordinate, more than 2 pi or less than -2pi'};
            return
        end
        r=output(1,1);
        angle1=output(2,1);
        angle2=output(3,1);
        output(1,1)=r*sin(angle1)*cos(angle2);
        output(2,1)=r*sin(angle1)*sin(angle2);
        output(3,1)=r*cos(angle1);
    else
        for j=1:3
        sum=sum+pow(v(1,j),2);
        end
        sum=sqrt(sum);
        output(1,1)=sum*r;
        output(1,2)=acos(v(1,3)/sum)+theta;
        output(1,3)=atan(v(1,2/v(1,1))+alpha;
        if((output(1,2)>pi)||(output(1,2)<-1*pi))
             prompt={'Invalid values for second cordinate, more than pi or less than -pi'};
             return;
        end
        if((output(1,3)>2*pi)||(output(1,3)<-2*pi))
            prompt={'Invalid values for third cordinate, more than 2 pi or less than -2pi'};
            return
        end;
        r=output(1,1);
        angle1=output(1,2);
        angle2=output(1,3);
        output(1,1)=r*sin(angle1)*cos(angle2);
        output(1,2)=r*sin(angle1)*cos(angle2);
        output(1,3)=r*cos(angle1);
    end
else
    prompt={'not numeric'};
    return
end
end
4

1 回答 1

1

您不在正确的目录中，或者您的路径设置不正确。您正在运行内置rotate函数，它与您提供的参数有问题。

如果您键入edit rotate.m，您会注意到该功能不是您的。

正确设置路径，和/或重命名函数。

于 2013-05-29T02:24:38.110 回答