0

我正在使用枫树进行差异化评估。我有一个问题。我想表达d/dt(alpha(t))这个方程中保持不变的变量(例如一部分):

-2*(diff(alpha(t), t))*sin(beta(t))*(diff(beta(t), t))*cos(psi(t))*
cos(theta(t))-2*(diff(alpha(t), t))*cos(beta(t))*sin(psi(t))*(diff(psi(t),t))*
cos(theta(t))-2*(diff(alpha(t), t))*cos(beta(t))*cos(psi(t))*sin(theta(t))*
(diff(theta(t), t))-2*(diff(beta(t), t))*sin(alpha(t))^2*(diff(alpha(t),t))*
cos(beta(t))*sin(psi(t)). 

任何帮助都是适当的。谢谢

4

1 回答 1

0

您可以通过替换来做到这一点。例如,假设涉及衍生产品的大量输出是通过运行一些我将缩写为“mycode;”的代码产生的。然后你可以这样做:

 output := mycode;
 new_output := subs(diff(alpha(t), t) = v,output);

然后,在 new_output 中,符号 diff(alpha(t), t) 的实例将被符号 v 替换,然后您可以使用 coeff 之类的函数来去除 v 的系数。这样,​​您可以找出输出的三角多项式表示是。

于 2011-06-02T22:44:17.003 回答