问题标签 [maxima]

我试图找出一种方法来通过算法获得在 Maxima 计算机代数系统中具有正弦项的函数的幅度和相位。这仅适用于稳态（因为 t -> infinity 和瞬态衰减）。例如，一个简单的案例是：

在这种情况下，增益为 1 / w，相位偏移为 theta，我们将忽略瞬态项 exp(-a * t)，因为我们只关心稳态增益和相位延迟，而 exp(- a * t) -> 0 作为 t -> 无穷大。我们也会忽略“+ 8”项，因为它只是一个 DC 偏移量。在我的工程课程中，我被教导这样做的方式需要大量的启发式方法和繁琐的方程式重新排列，以使它们以类似于上述的形式得到，其中答案是显而易见的。

考虑到我拥有计算机代数系统的全部功能（以及人们期望 CAS 具有的标准功能），有谁知道找到增益和相位延迟的通用算法方法。 ? 尽管我可能会在 Maxima 中实现它，但我当然会欣赏仅用数学解释的通用答案。

编辑： 我认为从我的示例中可以清楚地看出，我希望得到象征性的答案，就 w 而言。w 应该是欧米茄，代表输入的频率。我真正要问的是，是否有任何标准的数学运算可以产生增益和相位项，而无需进行大量启发式手动方程重新排列。

如何在 Macsyma、Scilab、Octave 中找到超定系统的解决方案？

我有一些复杂的方程式，我想用 Maxima 来评估。

在解决方案过程的一个步骤中，我将寄存器拆分为最重要和最不重要的组件，因此我需要能够在 Maxima 中做到这一点。

我试过用load(functs)，不，谢谢。

我正在尝试使用 为 Maxima 中的线性系统编写通用求解器linsolve(eqlist, varlist)，但不必明确指定问题的维度。

这有效，但将尺寸固定为 3：

linsolve( [ eq[0],eq[1],eq[2] ], [ a[0],a[1],a[2] ])

这不会：

关于如何让它发挥作用的任何见解？

---

问题背后的背景：这个线性系统出现在求解整数幂的有限和时，即有限多个平方、立方或一般幂的总和p。尽管有限平方和很简单，但一般解决方案却出奇地复杂：可以在此处找到讨论：递归关系的有限求和，第 2 部分。

我正在寻找将不同的“更高”数学运算与R结合使用的可能性。R 和 CAS 之间的链接或集成将是完美的解决方案。

R 和其他（数学和统计相关）系统的集成或反之亦然？它们的工作情况如何？你有什么建议？它们有多广泛（在时间、金钱和挫折方面）？

（FOSS 当然是首选 :)

我必须在 R 或者像 Maxima 这样的计算机代数系统中解决以下与卷积相关的数值积分问题。
Integral[({k(y)-l(y)}^2)dy]
其中
k(.) 是标准正态分布的 pdf
l(y)=integral[k(z)*k(z+y) dz]（标准卷积）
z和y是标量y
的域是-inf到+inf。
函数 l(.) 中的积分是不定积分。我需要在 z 上添加任何额外的假设来获得这个吗？
谢谢你。

我仍然是一个极大的新手，所以请多多包涵。我正在尝试编写自己的公式来计算矩阵的伴随（我知道 maxima 已经内置了一个，但我正在尝试自己的学习练习）。到目前为止，我有（对于 3x3 矩阵）：

我知道这可能不是最好的方法；但是，我想知道在使用 matrixmap 或 fullmapl 时是否可以访问 i 和 j 元素？

（我正在使用 wxMaxima 并且我没有很多 lisp 经验，我试图在不接触任何代码的情况下摆脱这个问题）。

我有一些功能是这样设置的：

并评估它们：

但是对于 g(1)，我得到的是符号答案而不是数字答案。使用 float() 试图获得一个数字答案，但它只是将积分中的所有项转换为浮点数。

我怎样才能得到 g(1) 作为一个数字？

我正在尝试 Maxima-fy 我的 Mathematica 框选项公式（https://github.com/barrycarter/bcapps/blob/master/box-option-value.m），但 Maxima 在一个相当简单的集成上崩溃：

使用某些值评估 upandin 崩溃：

如果没有 upandin 中的 float()，Maxima 只会以原始形式保留积分。

有人可以帮忙吗？我认为将 Mathematica 转换为 Maxima 会很容易，但现在我不确定了。

Mathematica 版本运行良好：

编辑：是否有任何类似 Mathematica 的开源程序可以在数值上近似该函数？我真的很想将开源代码发布到开源平台。

我很难找到如何做到这一点：

我希望输出是

不是