从jacbian maxima创建功能

Question

我想使用Maxima作为r：

的函数进行线性稳定性分析

f(x):=rx + x^3 - x^5
A:solve(f(x)=0,x)
J:jacobian([f(x)],[x])

现在对于A中的每个元素，我想检查J的符号作为r的函数。一般来说，我想要一个来自r的函数，它告诉我是否存在任何具有正实部的J的特征值。

Answer 1

也许你已经知道了，但是：Maxima中的乘法用星号表示。所以你必须写：

f(x):=r*x + x^3 - x^5;

到目前为止，我没有看到您的方法有任何问题。雅可比矩阵是一个1乘1矩阵，因此计算特征值是微不足道的。然后将x的值替换为x，并查看实部（函数realpart）。