Question

我有以下代码：

for n:1 thru 11 do for j:1 thru 21 do v[n,j]:1/sqrt(dp)* 
(sum(eigenfunctionsort[n,j]*exp(%i*2*%pi*m*x/dp),m,-10,10));

如果先前定义了eigenfunctionsort，x是我将在稍后整合的变量，我在m上求和。

当我打印说v [1,1]时，我得到了一个很长的令人讨厌的等式。我怎么能让Maxima把它煮成一些有意义的东西，这样我才能检查我的结果。

最佳，本

Answer 1

尝试'trigsimp'功能或者可能是map（trigsimp，your_expression）。不确定它会有所帮助，但值得一试。另请参阅'demoivre'。

Answer 2

我不知道你的矢量特征函数或你的“大长讨厌方程”是什么样的，但我经常从千里马得到复杂的特征值和特征向量，即使我知道它们应该简单而真实。

例如，

(%i1) A : matrix([1, 4, 1], [4, 1, 9], [1, 9, 1]);
(%i2) eigenvalues(A);

弄得一团糟。通过应用rectform将输出转换为笛卡尔形式，然后trigreduce减少结果的虚部，可以简化它。最后你可以转换结果浮点数：

(%i3) rectform(%)$
(%i4) trigreduce(%)$
(%i5) float(%);