问题标签 [numerical]

我需要一些帮助，我真的需要解决这个问题。

好吧，首先感谢您的时间...

我的问题：我有一个矩阵 (826x826 double)，我想将这个矩阵与 (826x1 double) 的向量相积分，我没有这些功能。是否有命令或算法来获取矩阵相对于向量的积分？请我真的需要帮助，我是matlab的新手。

真挚地。乔治

I started to read the book "numerical recipes in C"... I try make my programs more efficients and faster... So, is the same thing, to use C or C++? I know that C++ is a super-set of C... But, are there differences between "cmath" library (c++) with "math.h" library(c)? It will intervene in the velocity of execution?... in order to that... i can use C++ without any difference with a C program?

我正在尝试使用 matlab 编写一个 3 体问题。我得到了月球在太空旋转框架中的轨迹公式。它基本上是 ydotdot, xdotdot=GM/(x^2+y^2)^3/2 公式。公式是什么，并不重要。

我面临的问题是，我应该编写一个程序，以数值方式求解月球的轨迹方程。我使用 ODE45 进行比较，因为我的目标是获得与 ODE45 相同的结果。我的最终问题是，我想以天为单位迭代时间，所以 tspan= [0 365]。最重要的是当我将引力常数转换为秒然后执行 tspace= [0 365] 我得到完全不同的结果然后如果我要执行 [0 365*3600*34] 表示一年中的秒数并且 G= 6.67e -11。看来我的单位很奇怪。

我想知道是否有人可以解释为什么在我使用 ODE 45 时会发生这种情况。为什么我不能使用 ODE45 清楚地将秒转换为天？我需要做一个额外的步骤吗？我的问题中唯一的其他变量是半径、距离和 3 个物体的质量。

太感谢了。我已经为此工作了很长时间。任何帮助将非常感激。

我正在尝试根据以下算法创建 gauss-legendre 代码：

对于 n 个点

也就是说，它创建了一个 2n 方程系统（如果我们要求对 2n-1 阶的多项式准确，

ti是n阶legendre多项式的根。legendre多项式给出：

和 wi ：

我的代码是：

我的问题是：

1) 编译器给了我 :error: 'double*' 和 'double' 类型的无效操作数到二进制 'operator*' --> 在 result=w*f(*t);

2）我不确定我是否做对了整个事情。我的意思是，如果我将所有事情组合在一起并且我是否正确实施了算法。

我正在尝试为 D 编程语言调试我的并行库。最近提交的错误报告表明，使用任务执行的某些浮点运算的低位在运行时是不确定的。（如果您阅读该报告，请注意并行化简通过以确定性方式创建任务在幕后工作。）

这似乎不是舍入模式问题，因为我尝试手动设置舍入模式。我也很确定这不是并发错误。该库经过了良好的测试（包括通过了Jinx压力测试），问题总是局限于低位，即使在单核机器上也会发生，低级内存模型问题不太成问题。浮点结果可能因调度操作的线程而异，还有哪些其他原因？

编辑：我在这里进行了一些 printf 调试，似乎各个任务的结果有时在运行中有所不同。

编辑＃2：以下代码以更简单的方式重现了此问题。它在主线程中对数组的项求和，然后启动一个新线程来执行完全相同的函数。问题绝对不是我的库中的错误，因为这段代码甚至没有使用我的库。

输出：

舍入模式：0

0.7853986633972191094

舍入模式：0

0.7853986633972437348

另一个编辑

这是我以十六进制打印时的输出：

舍入模式：0

0x1.921fc60b39f1331cp-1

舍入模式：0

0x1.921fc60b39ff1p-1

此外，这似乎只发生在 Windows 上。当我在 Linux VM 上运行此代码时，我对两个线程都得到了相同的答案。

回答： 事实证明，根本原因是浮点状态在主线程上的初始化方式与在 D 中 Windows 上的其他线程上的初始化方式不同。请参阅我刚刚提交的错误报告。

我已经从提供的源代码中复制并粘贴了 C 的数字配方用于就地 LU 矩阵分解，问题是它不起作用。

我确定我在做一些愚蠢的事情，但希望有人能够指出我正确的方向；我整天都在努力，看不出我做错了什么。

回复后更新：该项目已完成并正在运行。感谢大家的指导。

WolframAlpha说答案应该是

我越来越：

到目前为止，我已经找到了至少 3 个不同版本的“相同”算法，所以我完全糊涂了。

PS是的，我知道至少有十几个不同的库可以做到这一点，但我更感兴趣的是了解我做错了什么而不是正确答案。

PPS，因为在LU分解中，较低的结果矩阵是统一的，并且使用（我认为）实现的Crouts算法，数组索引访问仍然是安全的，L和U都可以就地相互叠加；因此为此的单个结果矩阵。

我从开发通用 .NET 应用程序的角度考虑 Sho，而过去对我来说一直很困难的事情是 .NET 平台没有标准数学库。除了商业实现之外，两个当前的 .NET 库因其深度和广度而脱颖而出：

1) Sho（封闭源代码？，较新，但它是否会被包括在核心 MS .NET 框架中）

2）Math.NET numerics（开源，虽然在“测试版”中是一个非常成熟的项目，并且与来自 MSoft 的其他人一起从事该项目）

相似之处似乎如此之多，以至于我认为他们应该真正开始合作？

出于性能考虑，我将选择限制在上述两个方面，因为它们都允许连接到 MKL、ACML、LAPACK 等的本地快速库......

还有其他人在研究类似的东西吗？

谢谢

山姆

我正在使用 C# 和 .NET，我有一个看起来像这样的正则表达式

这匹配所有非数字字符，但我不希望它匹配小数点 (.) 和负号 (-)。我怎样才能用正则表达式做到这一点？

所以我尝试了 Chris' 并对其进行了一些调整以使其工作：

（我有一个名称为“原始”的文本框）

这Original.SelectionStart = Original.TextLength;是因为每当它被替换时，它都会将选择放在开头，这对用户来说似乎有点奇怪......

我有从流行病传播中得出的微分方程。我想获得数值解。这是方程式，

t 是一个自变量，范围从[0,100]. 初始值为

起初，我打算在matlab中用ode45函数处理这些，但是我不知道如何表达系列和组合。所以我在这里寻求帮助。

**

问题是如何将方程的右边表示为 odefun，它是 ode45 函数中的一个参数。

**

我不知道这东西叫什么名字，也不知道怎么形容，所以标题可能有点误导。

第一个附加图是用 pyplot 创建的。我想画一条贯穿所有图表的直线，而不是我目前使用的三个红点。在pyplot中可能吗？第二张图片是我正在寻找的。