1

我必须使用三个变量 x、y 和 z 绘制曲面，这些是三个关节角 q1、q2 和 q3 的函数。但我不知道如何绘制这个。

for q1=-170:5:170
    for q2=-65:5:140
        for q3=-180:5:70
            c1 = cosd (q1);
            c2 = cosd (q2);
            c3 = cosd (q3);
            s1 = sind (q1);
            s2 = sind (q2);
            s3 = sind (q3);
            x =(320*c1)-(975*c1*s2)-(200*c1*c2*s3)-(200*c1*c3*s2);
            y =(320*s1)-(975*s1*s2)-(200*c2*s1*s3)-(200*c3*s1*s2);
            z = (975*c2)+(200*c2*c3)-(200*s2*s3)+680;
        end
    end
end
4

1 回答 1

2

您需要考虑一些要点。

  1. 不需要for循环
  2. 你的系数向量需要相同的长度
  3. 您需要使用元素乘法运算符.*

编辑：我实际上首先误解了你的问题，现在应该是正确的。

所以你的代码看起来像：

q1x = -170:5:170;
q2x = -65:5:140;
q3x = -180:5:70;

[x1, x2, x3] = meshgrid(q1x',q2x',q3x');
Q = [x1(:),x2(:),x3(:)];

q1 = Q(:,1);
q2 = Q(:,2);
q3 = Q(:,3);

c1 = cosd (q1);
c2 = cosd (q2);
c3 = cosd (q3);
s1 = sind (q1);
s2 = sind (q2);
s3 = sind (q3); 
x  = (320.*c1) - (975.*c1.*s2) - (200.*c1.*c2.*s3) - (200.*c1.*c3.*s2);
y  = (320.*s1) - (975.*s1.*s2) - (200.*c2.*s1.*s3) - (200.*c3.*s1.*s2);
z  = (975.*c2) + (200.*c2.*c3) - (200.*s2.*s3)     +  680;

这是相当大的数据量。-x,y,z都是你表面的单点。

现在我真的不知道你的“表面”应该是什么样子。请说明您想获得什么。

scatter3会给你这个可爱的蘑菇：

scatter3(x,y,z,5,z)
colormap(jet)
view(45,10)

蘑菇

从散射数据中获取真实表面有点棘手，但您可以对“来自散射数据的表面”或其他东西进行一些研究......

于 2013-10-28T17:21:22.433 回答