numerical-methods 相关问题

Python-平方Equatin极限数

我已经编写了代码来解决Python中的二次方程式。但是它有很大的错误，例如“ 1e10”和“ -1e10”。有数字解决方案或python解决方案吗？在这里输入代码...

用Python解决耦合的PDE

我不理解如何使用python或python ode求解器求解耦合的PDE方程中的eta和V。 （或者有可能在没有...

从python中的非线性方程组中找到复 杂的根

我一直在测试一种已发表在文献中的算法，该算法涉及在Matlab和Python中求解一组'm'非线性方程。非线性方程组涉及...

[3d折线沿法线方向的偏移

假设3d折线（即3d点段的多边形链）具有为每个点指定的法线。是否有任何算法可以计算点位于...

累积总和与简单总和有什么区别？

假设我有一个一维数组，可以沿z方向对其进行积分。最简单的方法是用z方向上的网格间距乘以总和来代替积分，在这种情况下，我只是...

Scala的数字图书馆

我正在寻找一个库，以便在Scala（或Java，尽管可以使用scala函数的东西会更好！）中至少具有以下功能：L-BFGS最小化器...

浮点运算：派生摆动

PDE的前向和后向集成

我有以下PDE，其最终条件和边界条件如下PDE Equaton所示。我试图及时向后整合此pde，但是，解决方案似乎不合时宜。我...

数值计算阶乘和多项式的组合

我正在尝试编写一个简短的 C++ 例程来计算给定整数 j > i（通常它们位于 0 和 100 之间）和复数 z（以 |z| 为界）的以下函数 F(i,j,z) < 100),...

def randvec(n): vector = np.ndarray((1,n),dtype=float) for i in range(0,n): vector[i] = np.random.uniform(0,100) print(vector[i]) test = randvec(4) print(test)

“ W4方法的实现”（牛顿-拉普森扩展）

我可以对不可微分的函数使用自动微分吗？

我正在测试不同求解器的性能，以最小化由矩量模拟方法得出的目标函数。鉴于我的目标函数不可微，我想知道是否...

在odeint Python中使用Dfun（Jacobian）

这是一个相当普遍的问题，其中一部分可能适用于耦合ODE的任何数值模拟，其中一部分可能仅适用于scipy中的odeint方法。集成Python ...

如何减少2D连接域上集成的集成时间

我需要在简单连接的域上计算许多2D积分（大部分时间都是凸的）。我正在使用python函数scipy.integrate.nquad进行此集成。但是，时间...

3个稀疏矩阵的有效积，产生密集的中间物

我有3个矩阵，它们都是稀疏矩阵，A，B和C。我需要取AB的矩阵乘积，这将导致密集的矩阵。之后，我需要AB的元素明智乘积（元素明智*）C。C ...

Matlab欧拉方法EOM

f = @（t，y）[y（2）; cos（t）-y（2）-4 * y（1）] a = 0 b = 1 N = 5000 h = 1 y（1）= 0 y（2）= 0.1 t（1）= a对于n = 1：1：N t（n + 1）= t（n）+ hy（n + 1）= y（n）+ h * f（t （n），y（n））；我正在尝试结束图（t，y）...

快速求解数值方程

我正在使用scipy.optimize.newton_krylov解决以下问题。但是，它相当慢（这完成了数千个点）。即使我最初的猜测非常接近解决方案，也可以...

整数和数值优化（nlminb）R

我遇到了涉及包含未知变量的积分的数值估计的优化问题。数值估计积分很简单，只需使用...

Laguerre方法

我正在尝试为以下方程式实现Laguerre的方法：1 /（99 * x + 1）= 2。 （一般形式更复杂-n次多项式1 /（a * x + 1）+ ... + 1 /（z * x + 1）= res的多项式）其中a，b，... z> = 0和0 ]

代码错误，计算sqrt（a）的方法

计算sqrt（a）的一种方法，a> 0是X（n + 1）=（a +（X（n）* X（n-1））/（X（n）+ X（n -1）），n = 1，2，…，且X0 = 1且X1 = a（即，已知lim n-> Xn的无穷大= sqrt（a）...