0

在我使用 wxmaxima 11.08.0（ubuntu 12.04，Maxima 版本：5.24.0）进行练习期间，
我遵循了P.Lutus 的示例，而他的第二个示例对我不起作用。

eq: y(t) = -r*c*'diff(y(t),t)+m*sin(%omega*t);
sol:desolve( eq, y(t) );
Is  %omega  zero or nonzero? nonzero

然后 Maxima 不再有反应，直到我重新启动它。
我需要先激活或定义最大值以获得结果吗？

预期的输出应该是：
中间结果

如果我手动定义方程，我的问题还有第二部分：

sol: y(t) = (m * sin(%omega*t))              / (%omega^2*c^2*r^2 + 1) - 
            (%omega*c*m*r*cos(%omega*t))     / (%omega^2*c^2*r^2 + 1) +
            (%omega*c*m*r*%e^-((1*t)/(c*r))) / (%omega^2*c^2*r^2 + 1);

连续过程的初始条件：

init_val:-(c*m*r*(%e^-(t/r*c))*%omega)/(c^2*r^2*%omega^2+1); 
atvalue(y(t),t=0, init_val);
try2 : desolve(sol,y(t));
"Is  "%omega"  zero or nonzero?" nonzero;

上个学期还在

这里最后一个任期仍然存在。这些问题是基于三角函数的使用吗？

最好的问候，马库斯


我通过 PPA 更新到 wxMaxima 13.04.0 和 Maxima 5.29.1。现在desolve完成了，但是最后一个词似乎很复杂。 最大值 5.29.1 解决结果

使用负数最后一项和 desolve 命令执行 init_val 仍然会在等式中留下 %e^(..)*... 。

4

1 回答 1

1

您可能会在 Maxima 邮件列表上对这个问题更感兴趣。见： http: //maxima.sourceforge.net/maximalist.html

对于 Lutus 示例 2 的第一个版本，我得到：

y(t) = m*sin(%omega*t)/(%omega^2*c^2*r^2+1)
 -%omega*c*m*r*cos(%omega*t)/(%omega^2*c^2*r^2+1)
 +(y(0)*%omega^2*c^3*r^3+%omega*c^2*m*r^2+y(0)*c*r)*%e^-(t/(c*r))
  /(c*r*(%omega^2*c^2*r^2+1))$

如果 y(0) = 0，这与预期结果相同。但是，我看不出这是在哪里假设的。

之后atvalue(y(t),t=0,init_val)，我得到与 Lutus 相同的结果，即：

y(t) = m*sin(%omega*t)/(%omega^2*c^2*r^2+1)
 -%omega*c*m*r*cos(%omega*t)/(%omega^2*c^2*r^2+1)$

我正在 Linux 上使用由 Clisp 构建的 Maxima 5.31.1。

于 2013-10-01T21:13:46.960 回答