问题标签 [equations]

我正在尝试为数学中的非线性演化方程规定自由边界条件，并且我想作为我所做的是否正确的第二意见。

边界条件已标有注释，即（自由边界条件）

我还想为固定边界条件运行它。

我正在使用 FindRoot 求解三个变量（H0D、H0S 和 H1S）中的三个非线性方程。除了感兴趣的三个变量之外，这些方程中还有四个参数我希望能够改变。我的参数和我想改变它们的范围如下：

CF∈{0,15} , CR∈{0,8} , T∈{0,0.35} , H1R∈{40,79}

问题是我的非线性系统可能对此参数范围的一部分没有任何解决方案。我基本上想问的是，是否有一种聪明的方法可以准确地找出我的参数范围的哪一部分允许真正的解决方案。

我可以在循环中运行 FindRoot，但由于非线性，FindRoot 对初始条件非常敏感，因此经常出现错误消息可能是因为初始条件不好而不是没有解决方案。

有没有办法让我找出参数空间的工作原理，而不是手动插入 10^4 个参数值组合并玩弄初始条件并希望 FindRoot 给我一个解决方案？

非常感谢，

我想知道是否有人可以帮助我使用 MatLab 解决 Lotka-Volterra 方程。我的代码似乎不起作用。我执行以下操作：

第1步 -

我创建了一个名为 pred_prey_odes.m 的文件，其中包含以下代码：

在将以下代码输入命令窗口之前，我保存了文件并确保它位于当前目录中：

但是，没有任何情节弹出。事实上，在 matlab 中什么都没有发生，我什至无法清除命令窗口。如果我输入 clc 什么都不会发生...

任何帮助，将不胜感激！

谢谢！

——斯内哈·英古瓦

我是编程新手（Python 是我的第一语言），但我喜欢设计算法。我目前正在研究一个方程组（整数），我找不到任何参考资料来解决我的特定问题。

让我解释。

我有一个方程式（如果你愿意，可以做一个测试）：

其中 a 是某个常数。

我的问题是，变量 z 的计数方式与斐波那契数列非常相似，变量 x 是 z 的步长。所以我的意思是（对于斐波那契数列），在 z 序列的第一项，x = 0，在 z 序列的第二项，x = 1。我需要求解 y。

确定z的具体过程如下

我需要扫描（跳过）z < x 的值以测试 y 的整数解的条件。

这看起来可能吗？

我对这些问题在数学和编程上都是新手。如果有人可以建议使用可以解决以下问题的 c++ 库，我将不胜感激。

给定常量：

{x_1, ..., x_n}, {y_1, ..., y_n}, {z_1, ..., z_n}, C, & variables {q_1, ..., q_n}

最大化：sum(i = 1..n} q_i*x_i

受制于：C - sum(i = 1..n){ sum(j = 1..q_i) [y_i + (j-1)*z_i ] } >= 0 AND q_i >= 0

所有常量都是大于零的整数。q_i's也是整数。

所以我试图解决{q_1, ..., q_n}

我想要一种可以解决 nxn 行列式并返回值的计算机算法。最推荐使用 C++ 语言，其中输入是 2D 数组，N，输出是值。

数学中有一种详尽的方法可以解决大小超过 3 x 3 的行列式。如果有人能找到解决行列式的代码片段，那将很有帮助。

考虑到“n”个线性方程组的系数，主要目的是求解“n”个变量。一旦 'n' 超过 4，这件事在纸上就会变得一团糟。有一种求解此类方程的“行列式”方法。虽然一旦'n'超过3，它在纸上也会变得凌乱，但如果我能找到计算机如何解决它，它是可行的。谢谢！

假设我在一些变量之间有一些已知的关系，例如a = b + c和一个起始解决方案，例如a=2, b=1, c=1。我正在编写一些代码，例如，给定更新a=3，我更新 and 的值，b以便c仍然满足关系，例如设置b=2（显然有很多可能性，我只需要一个）。

在实践中，有很多关系，它们不是线性的。该代码通过用户为每个变量提供函数作为其他变量的函数来工作，例如a(b,c)、b(a,c)和c(a,b)。然后我构建一个依赖关系图，当一个值发生变化时，我会执行广度优先搜索之类的操作来更新其他值。这似乎可行，但我不禁认为我正在重新发明轮子，这是一些众所周知的计算机科学/图论问题。也许甚至还有一个包裹？

谁能提供一些链接或一些关于我正在解决的问题的见解？

您好我正在尝试将 muParser 与用 VC++6 编写的遗留 MFC 应用程序一起使用，将代码迁移到更高版本不是一种选择。

有没有人在这种环境中使用 muParser 取得任何成功。任何示例源代码都会很棒。

我的任务是为以下公式编写 MIPS 指令代码：

我在理解公式的实际含义时遇到问题。

据我了解，我们正在将一个 int n 传递给双递归程序。

因此，对于 f(0)，等式将是：

如果 n=10，则方程为：

我知道我根本没有做到这一点，因为它不会是递归的。您可以阐明方程式的实际含义的任何信息都会很棒。一旦我理解了方程式，我应该能够编写 MIPS 代码。

我需要求解一个二维泊松方程，即 for AX=B 中的方程组，其中 A 是 n×n 矩阵，B 是 n×1 向量。作为 2D Poisson 问题的离散化矩阵，我知道只有 5 个对角线不会为空。Lapack 不提供解决此特定问题的函数，但它具有求解带状矩阵方程组的函数，即 DGBTRF（用于 LU 因式分解）和 DGBTRS。现在，5条对角线是：主对角线，主对角线上方和下方的第一条对角线以及主对角线上方和下方的两条对角线，m对角线wrt。在阅读了有关带存储的 lapack 文档后，我了解到我必须创建一个 (3*m+1)-by-n 矩阵来以带存储格式存储 A，我们将这个矩阵称为 AB。现在的问题：

1) dgbtrs 和 dgbtrs_ 有什么区别？英特尔 MKL 两者都提供，但我不明白为什么

2) dgbtrf 要求带存储矩阵是一个数组。我应该按行还是按列线性化 AB？

3）这是调用这两个函数的正确方法吗？