问题标签 [maxima]

我想在 wxMaxima 中做一些简单的矩阵操作。我面临的问题是，Maxima 总是根据矩阵的各个元素提供答案。这是一个例子：

我想扩展 (A+B)(A+B)`，其中 A 和 B 是 2x2 矩阵。然后我想看看输出：

AA`+AB`+BA`+BB`,

但我得到的是一个像这样的矩阵：

[(b2+a2)^2+(b1+a1)^2 .......]

[.................................. .......] （我刚刚放了点而不是输入所有元素）

那么有没有办法强制千里马在不查看矩阵内部元素的情况下操纵矩阵？

谢谢

（来自讲义参考）为了使 Gauss-Seidel 和 Jacobi 方法收敛，需要检查系数矩阵是否对角占优，即对角元素应该在其列中的所有元素中具有最大值. 如果它还不是对角线占优，请使用 pivoting。对于对角占优的矩阵，应满足以下条件：（这也称为收敛）

是否有任何预定义的函数可以在maxima中使用来实现收敛，或者我应该使用交换进行循环以及应该使用哪些约束？假设矩阵的大小为 3x3，包含非零元素。

我已经看到了一些相关的问题，但答案在 matlab 中。

链接：是否有检查矩阵是否对角占优（行占优）的功能

那么，我怎样才能做到最大值呢？

我正在尝试学习如何使用 Maxima。集成出现问题：

你怎么看？

如何展平由Maxima中的其他矩阵组成的矩阵？

我的意思是，如果我将一个矩阵作为一个新矩阵中的一个元素，我将无法直接访问底层元素。外矩阵的维数没有变化！

我有以下代码：

当我尝试prob_exp使用一些具体的值调用时s，例如：

我收到以下错误消息：

有没有办法强制s作为值传递？

to_poly_solve最大值中的函数返回 a%union但我想使用列表。

既然不起作用，如何将%unions 变成列表？listify

如何定义一个以最大值返回闭包的函数？

我的具体例子：我有一个表格列表：

我想根据等号之前的部分提取一个元素。

这是我的尝试：

这如我所愿。

由于真正的列表要长得多，而且我不想复制粘贴代码，我现在想像这样使用它：

但是第一部分已经不起作用：

在 Maxima 中，我想使用函数对任意数量的变量x[1]...中的微分方程组进行数值求解 。我在变量名中使用方括号，以便我可以在 for 循环中索引变量。这是当 n = 2 时的示例调用：x[n]rkrk

不幸的是，我得到了错误：

如果我手动替换X[1]withx1和X[2]with x2，错误就会消失：

似乎rk不喜欢变量名中的方括号，但我需要它们来索引 for 循环中的变量。我该怎么办？

我愿意计算一个不完全明确的函数的 Frechet/Gateaux 导数，我的问题是：最有效的方法是什么？你会推荐我使用哪种语言？

准确地说，我的问题是我有一个函数，比如 F，它是多维函数对（即从 R^n 到 R^k）的乘积之和的欧几里得范数的平方。

AFAIK，如果我使用 Maple 或 Maxima，他们会要求我明确公式中涉及的函数，而我想让它们保持抽象。然后，我需要计算 Frechet/Gateaux 导数以保持表达式简单。确实，当我按照标准方式进行时，我开始将欧几里得范数的平方发展为平方和，并且有很多索引。我的目标是做一个三阶整数余数的泰勒展开式，在我看来，这个表达式在人类看来是不可行的（公式超过一页 A4 纸）。

因此，我更喜欢使用 Frechet/Gateaux 导数，这将允许我保留标量积而不是总和。

由于所涉及的函数与其导数有一些相似之处（由于存在指数），因此只需了解少量规则。所以我想我可以自己制作这样一个专用的计算机代数系统。

我开始学习 LISP，因为我读到它对我的问题很有效，但我现在有点迷茫，因为这种语言非常不同，我仍然习惯于用 C/Python/Perl 来思考。 ..

这是另一个问题：您是否有一些关于如何制作符号计算的代数系统的课程或文章的链接（最好是在 LISP 中）？欢迎任何建议。

非常感谢您的回答。

直到今天，我一直在 Maxima 中使用假单位，当我意识到有支持转换和其他很酷功能的单位包时。我发现 ezunits 是其家族中最完整的软件包，但我发现它存在问题。

每次我想使用它，我都必须编译它。换句话说，运行该行会load(ezunits)启动一个编译过程并将文件的位置作为输出输出ezunits.mac。但是，如果要打开一个新的 maxima 会话并load(ezunits)在其上运行，则编译将再次执行。人们会期望只编译一次，然后只加载一个编译文件。

有没有办法解决这个问题？