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 に答える