与Maxima循环

时间:2018-07-26 19:40:09

标签: loops for-loop maxima

for i:1 thru 3 step 1 do;
posix:arithsum(li*cos(ri(t))),1,i-1)+(li*cos(ri(t))/2);
posiy:arithsum(li*sin(ri(t))),1,i-1)+(li*sin(ri(t))/2);

我想做的是得到6个位置函数(3 x和3 y)。它应该给我类似以下的值:

pos1x:l1*cos(r1(t))/2;
pos2x:l1*cos(r1(t))+l2*cos(r2(t))/2;
pos3x:l1*cos(r1(t))+l2*cos(r2(t))+l3*cos(r3(t))/2;

那么,为什么我的代码不起作用?

1 个答案:

答案 0 :(得分:2)

这里的事物。 (1)for循环仅将一个表达式作为其循环体;通常,将多个表达式组合为(e1, e2, e3)block(e1, e2, e3)。请注意,for ... do;是不正确的语法,因为它没有循环主体-分号终止for表达式。还要注意,主体中的表达式用逗号而不是分号分隔。 (2)您可以使用下标符号来索引项目; Maxima不会自动构造符号名称,例如pos1x。而是使用下标符号:posx[1]posy[i]等。

鉴于此,这是一个解决方案。

(%i1) load (functs);
(%o1) /Applications/Maxima.app/Contents/Resources/opt/share/maxima/5.41.0/shar\
e/simplification/functs.mac
(%i2) for i:1 thru 3 step 1 do
 (posx[i]:arithsum(l[i]*cos(r[i](t)),1,i-1)+(l[i]*cos(r[i](t))/2),
  posy[i]:arithsum(l[i]*sin(r[i](t)),1,i-1)+(l[i]*sin(r[i](t))/2));
(%o2)                                done
(%i3) [posx[1], posx[2], posx[3]];
       l  cos(r (t))  3 l  cos(r (t))                          l  cos(r (t))
        1      1         2      2                         1     3      3
(%o3) [-------------, ---------------, 2 (l  cos(r (t)) + -) + -------------]
             2               2             3      3       2          2
(%i4) [posy[1], posy[2], posy[3]];
       l  sin(r (t))  3 l  sin(r (t))                          l  sin(r (t))
        1      1         2      2                         1     3      3
(%o4) [-------------, ---------------, 2 (l  sin(r (t)) + -) + -------------]
             2               2             3      3       2          2

我想l[i]r[i]也应该下标。为了解决语法问题,我更改了括号。如果您有其他需要,当然可以继续进行更改。

请注意,在此公式中,posxposy是所谓的未声明数组。未声明的数组适用于表示下标的符号变量。您可以通过listarray获取元素列表。