问题标签 [mupad]

我正在尝试在 m 文件中从 Mupad 调用 betaQuantile 函数，以便可以在 Matlab2014b 中运行它。我正在查看这两个似乎相关的帮助文件，但仍然无法弄清楚代码
http://au.mathworks.com/help/symbolic/mupad_ref/stats-betaquantile.html?s_tid=gn_loc_drop

http://au.mathworks.com/help/symbolic/call-built-in-mupad-functions-from-the-matlab-command-window.html

我正在尝试mtaylor从 MatLab 中的 MuPAD 引擎运行该函数，它提供了函数的多元泰勒展开。但是，它一直告诉我我正在尝试围绕一个无效点进行扩展。这是我尝试过的最小工作示例：

为什么这不起作用？

How do I easily compare coefficients of two polynomials in MuPad?

被问到似乎更清楚；但我第一次搜索时找不到它。

有谁知道如何在 Matlab 的 Mupad 会话中找到数组的大小？

它们都不适合我。

我想stats::swGOFT在 MuPAD 中使用该功能。我有一个名为 的数值向量r。我用了

错误是

使用 mupadengine/feval 时出错（第 157 行）
MuPAD 错误：错误：某些数据的类型无效。

所以我尝试了一种更直接的方法，它奏效了：

但这不起作用：

我的变量r是一个 1146×1 双向量。显然我不能手动输入所有的数字。那么，如何将向量变量传递r给 MuPAD 函数stats::swGOFT呢？

我正在尝试将附加的方程移植到 Matlab Mupad 中。

我尝试了以下 Mupad 语言代码：

我试图弄清楚正确的 Mupad 代码应该是什么来解释这个等式。最后这个选择器更令人困惑。选择器被描述为

δ()是一个选择器，即δ() 1 x = if x>0, else δ() 0 x = ;

我在 MuPAD 5.7.0 版（MATLAB R2011b）中发现了一个奇怪的行为，我想知道这是否是一个错误，如果不是，我做错了什么。理想情况下，我还想知道为什么 MuPAD 会这样做。

考虑大小为 3x3 的数组C，其中的元素具有一些示例值。我想将此数组视为数组数组，因此使用级联索引。

当两个索引都是不同嵌套范围的局部变量时，即当第一个索引的范围比第二个索引的范围宽时，显然会出现问题。如果第一个索引是常数，则没有问题。

当我输入：

我得到以下结果：

我希望代码中的第 3 行和第 5 行（输出中的第 2 行和第 4 行）产生相同的结果g+h+i：相反，第 5 行产生a+e+i，这似乎是 的对角线C。

当我对product而不是做同样的事情时sum，结果甚至更奇怪，但可能会揭示更多关于“错误”来源的信息，特别是 DOM_VAR(0,2)：

我得到：

我可能在这里走错了路，但我怀疑创建了一个闭包，试图保存周围作用域的局部变量，这些变量在创建闭包时是不确定的。此外，替换似乎在某个时候停止，它被eval().

实际问题

我试图解决的实际问题如下：

问题：计算形式的多项式

在Wikipedia上，定义了以下形式：

为了实现这一点，我们需要定义向量 alpha 的集合，其和等于 m。这些对应于长度为 m 和可能的零元素的 n 的可能组合的集合：

例如，总和

将要求这些权力组合，其中每一个都是一个向量 alpha：

此外，我们需要多项式系数。每个这样的系数取决于向量 alpha 和幂 n：

例如，第二个组合出现的次数是：

对所有组合求和产生：

幂是正确的，但系数不正确。这似乎与此问题中首先陈述的归结抽象问题有关。

我一直在运行 Matlab R2014a 的计算机上运行使用“求解”的代码。有问题的行如下所示：

符号z

如您所见，我正在尝试解决 z。给出求解的参数中列出的所有其他变量（Sum0、Sum1、Thing1 等）都是先前计算的标量。

当我在安装了 R2014a 的计算机上运行我的程序并显示 ZETA_NOW 时，可以看到一组五个真实而复杂的解决方案。真正的解决方案是我想要的解决方案，并且根据与已知结果的比较，我确信它是正确的。

当我在运行 R2014b 的计算机上运行我的程序并再次显示 ZETA_NOW 时，五个答案显示为：

此外，对于某些输入值，没有给出真正的解决方案，在另一台计算机上运行相同的输入值时，情况并非如此。

是否有一些我需要更改的“解决”选项，例如允许的迭代次数？

想象一下，我在 MuPad Notebook 中定义了两个变量：

对于产品

我得到显示：

如果我使用hold，我可以得到表达式：

但现在我想同时拥有显示的表达式和结果。对我来说似乎合乎逻辑的两个选项不起作用：

和

如何显示表达式和结果？ 如果在两行中就可以了，但我更喜欢在一行中：

z = xy = 10

我正在尝试使用以下代码评估定积分：

但它给出了以下错误：

但是，当我在没有某些功能的情况下评估积分时：

我得到的答案是Inf。任何想法？