问题标签 [numerical-analysis]

在这里，我有一个矩阵，例如，A，其中 A=[1 0.9 0.5;0.9 1 0.9;0.5 0.9 1]，如何计算其最接近的半正定矩阵？有什么命令或算法吗？

我真的在 Modelica 规范中找不到任何答案，所以请教你们。规范指出， A tool is free to solve equations, reorder expressions and to not evaluate expressions if their values do not influence the result (e.g. short-circuit evaluation of Boolean expressions). If-statements and if-expressions guarantee that their clauses are only evaluated if the appropriate condition is true, but relational operators generating state or time events will during continuous integration have the value from the most recent event. If a numeric operation overflows the result is undefined. For literals it is recommended to automatically convert the number to another type with greater precision. 现在，我想知道，该工具可以选择在积分器步骤中多次评估表达式吗？例如（可能不是一个有效的例子，只是为了让你们知道我想知道什么:)）

这将多次打印同一时间，所以我认为正在进行某种迭代。但我真的很想得到一些消息来源的证实。

我应该在matlab代码中输入什么条件，以便使用这些公式获得二次的精确解：

直接实现这些公式在某些情况下我没有得到正确的解决方案，例如 x^2-1000001x+1

非常感谢您的帮助

我使用改进的欧拉方法编写了一个 C 代码，以定期确定振荡器的位置、速度和能量。但是，我遇到了一个问题，即振荡器的能量正在减少，尽管没有耗散项。我认为这与我更新位置和速度变量的方式特别相关，并希望得到您的帮助。我的代码如下：

我可能会错过一个非常轻微的观点，所以如果你能指出的话，我将不胜感激。我感谢您的帮助。

有谁知道寻找任何 n × n 矩阵的特征多项式的 Souriau 方法？我找到了第一个系数，很明显，但是我怎样才能找到其他系数呢？在我需要对矩阵求逆但我知道怎么做之后。

我想在我的一个文档中实现一个分位数函数LaTeX。具体来说，是 t 分布的分位数函数。是的，我说的是文档中的qt()函数的实际数值近似。R

分位数功能通常是如何实现的？无论是正则化 beta 函数，还是我在网络上任何地方都可以找到的超几何函数和LaTeX/tikZ似乎都没有计算积分或无限和/积（或者至少不是及时，我听说过大方向的选项包xint）。当我试图逼近分位数函数时，这两个函数（包括sgn(x)，但在 中很容易实现）都会不断出现。TeX

如果您知道LaTeX分位数函数的特定实现，或者您知道从哪里开始，则可以加分。

我目前正在编写一个程序，该程序可以帮助我识别数字序列中的函数，它只是写出每个数字之间的差异。

就像我给它序列 [21, 20, 18, 15, 11]，它会返回 [-1 -2 -3 -4] 帮助我识别序列中的下一个数字可能是 6。

到目前为止，这是我的代码：

这非常简单直接（我还在学习）。但是，我想知道是否有办法让代码也返回“权力”。

就像我给它 16、256、65536 一样，它应该返回模式 [1, 1]。这两个整数将表示第一项乘以 1，第二项乘以 1，等等。您知道，如 16*16 = 256 等...

我怎么能构造一个方法来做到这一点？java中是否有任何power-to功能？

提前致谢！

牛顿拉夫森法的失效分析说“对于某些函数，某些起点可能会进入一个无限循环，阻碍收敛”。我想在程序中检查它是否进入无限循环或不使用 assert 语句。如果它进入，那么程序将终止，说使用这个初始猜测是不可能收敛的。如何在程序中检测到这个循环？代码：

我正在尝试在 MATLAB 中实现 Jacobi 迭代，但无法使其收敛。我已经在网上和其他地方查看了用于比较的工作代码，但我找不到任何与我的代码相似并且仍然有效的代码。这是我所拥有的：

无论我使用什么 A 和 b，这都会爆炸。这可能是我忽略的一个小错误，但如果有人能解释什么是错误的，我将非常感激，因为这应该是正确的，但在实践中并非如此。

有没有办法用一个公式 [y = f(x)]来描述PolynomialSplineFunction或 Java 中的任何其他样条函数？