问题标签 [symbolic-computation]

我是单数的新手。我刚刚下载了 Singular4-0-2_64.dmg，安装了图像，右键单击，显示包内容，然后将内容文件夹移动到 Applications 目录。双击 Contents/MacOS/Singular（或 Esingular） 我在终端中运行 Singular（或通过终端中的 emacs）。但是当我在运行的 emacs 中尝试 Mx 单数时，它会说“无法打开加载文件：没有这样的文件或目录，单数”。

然后我将 .emacs-general 和 .emacs-singular 复制到主目录。将以下行添加到 .emacs

(setq load-path (cons "" load-path)) (autoload 'singular "singular" "使用默认值启动 Singular。" t) (autoload 'singular-other "singular" "请求参数并启动 Singular。" t )

将singular-emacs-home-directory 更改为/Applications/Singular.app/Contents/share/singular/emacs/ 它仍然不起作用。此外，当我用 .emacs-singular 替换 .emacs 时，替换上面的目录，它会给出错误：警告（初始化）：加载 `/Users/Me/.emacs' 时发生错误：

符号作为变量的值是无效的：/Applications/Singular.app/Contents/share/singular/emacs

我没有找到正确的奇异 emacs-home-directroty 还是有任何其他问题？

我感谢您的帮助。

另外我必须提到，我是 Macaulay2 的 emacs 用户。为此，我更改了 .emacs 和 .emacs-Macaulay2 以教 emacs 在哪里寻找 M2。Singular 有什么类似的事情要做吗？我应该在某处有一个文件 .emacs-singular 吗？

谢谢你的帮助。

我正在尝试为圆图 f(x) = 2x + a + b*sin(2*pi*x)/pi mod 1 绘制 Arnold 舌头的边界（存在周期性解的区域）。这些是当f^n(x)=x且d/dx(f^n(x)) = 1时定义，其中f^n(x)表示函数迭代n次，即f^2(x) = f(f (x))，n为周期点的周期。

我希望能够采用这两个方程并用 b 为 Arnold 舌头的边界写一个方程，所以我会得到满足方程的 x = g(b) 和 a = h(b)。然后我想针对 b 绘制 a 。

从分析上讲，我可以通过重新排列 d/dx(f(x)) = 1 for x 以这种方式解决周期 1 边界的问题，x 以 b 形式给出 x，然后将此值代入 f(x) = x 以给出 a in b. 条款 我还设法通过以下方式在 MATLAB 中使用符号方程来做到这一点。

问题是，我有办法解决第 2 期或更高的边界吗？我试过以下，

但是，我收到“找不到明确的解决方案”警告。我认为也许这个方程太复杂了，MATLAB 无法以符号方式求解。有没有办法使用符号方程让它工作？如果没有，是否有合适的数值方法来执行上述操作？

非常感谢任何帮助，在此先感谢。

我有一个函数F ，它将向量a作为输入。函数的输出和a都是长度为N的向量，其中 N 是任意的。每个分量F _n的形式为g(a(n),a(nk))，其中g对于每个分量都是相同的。

我想在 matlab 中使用它的符号功能实现这个函数并计算它的雅可比（然后使用 .m 将函数和它的雅可比存储为常规的 .m 文件matlabFunction）。我知道如何为每个输入都是可以手动处理的标量的函数执行此操作。但是在这里，我想要一个能够为任何 N 生成这些文件的脚本。有没有很好的方法来做到这一点？

我想出的一个解决方案是生成一个字符串数组 "a0","a1", ..., "aN" 并使用eval. 但这很混乱，我想知道是否有更好的方法。

谢谢！

[编辑]

这是我当前解决方案的最小工作示例：

最近我遇到了以下悖论（至少在我看来）。

我试图在 Matlab 中评估符号矩阵的等级（不确定rank()在符号情况下是否按预期工作，但无论如何），但为了简化计算产生的系数，我vpa()在矩阵中使用了。

我注意到的是，使用vpa命令后，矩阵突然变成了满秩矩阵。

有没有人有过类似我描述的事情的经验？我应该相信哪一个？