通过数值逼近（与符号计算相反）解决数学问题的算法。
我尝试过调整欧拉方法代码，并使用欧拉来计算第一个值，所以我可以在verlet中使用两个值，但是当我绘制图形时，我只得到两条垂直的直线...
Maple 2023plot() 函数对于复合函数的工作非常奇怪。我是不是做错了什么？
我正在尝试制作一个在 1 和 0 之间振荡的周期函数。我制作了一个在间隔上按照我需要的方式振荡的函数 <0,7), when i tried to make it periodical i wrote a mod...
Python 2.7 的用户输入功能存在问题，无法让这些输入与程序一起使用
我对 python 世界很陌生，尽管我已经用 php 做了很多工作......这是我的案例...... 我正在用 python 2.7 为我的小程序编写一些代码。 在该程序中，我需要接受 2 个用户输入，机器人...
问题是找到 Z 使得 epsilon(E)（方程 21）最小化。 Z 是一个 MxN 矩阵，这就是我们需要找到的。 Zx和Zy也是已知的MxN矩阵。 Dx 和 Dy 是...
Haskell 中缺乏有限元方法实现 - 有什么具体原因吗？
我很好奇为什么 Haskell 或任何函数式语言中似乎缺乏有限元方法 (FEM) 实现。鉴于 Haskell 的纯函数性质，我预计......
我目前正在尝试在Python中计算以下表达式 该表达式来自 t 分布，v 是表示自由度的标量，mu 和 sigma 是向量和矩阵
我正在尝试编写代码来在线性方程组上实现雅可比方法。 到目前为止我已经写过这个： #雅可比方法 我=0 z = np.array([[-1.,3.,2.,11.],[1.,1.,1.,6.],[0...
我正在运行梯度下降来找到非线性方程组的根，我想知道你如何检测该方法是否陷入局部最小值，因为我相信随着设置...
Python 中的 Hermite 插值，无需使用 scipy 库
def hermite_interpolation(x, y, yp, xi): n = len(x) - 1 结果 = 0.0 对于范围 (n + 1) 内的 j： 项 = y[j] 对于范围 (n + 1) 内的 i： 如果 i != j 且 x[j] != x[i]： ...
我正在实施 Dekker 的数值方法来求根，但我对更新间隔有点困惑。维基百科上提到了这一点 选择新的对位点使得 f(a_k+1) 和 f(b_k...
我正在使用 Verlet 积分（更具体地说，是 Verlet 积分维基百科页面“非常量时间差”部分中的最后一个方程），以及 C++ 和 Eigen 库...
我正在尝试实现 Dekker 的求根方法，但我似乎无法找出算法似乎无法正确收敛到解决方案并且有时会突然停止的问题......
如何在没有库的情况下数值求解隐式函数（方程）a*y^2 - b*x*y + c = 0？
我想知道如何用 x 来求解 y a*y**2 - b*x*y + c = 0 使用数值方法，例如 Newton-Raphson、Bisection 等，无需库，除非用于定义...
我有一个Matlab代码 模型 = createpde(); GeometryFromEdges（模型，@lshapeg）； applyBoundaryCondition(模型,"狄利克雷", ... “边缘”，1：模型.Geom...
我有一个函数 f(x1,…,xN)，其中变量 xi 仅采用二进制值，xi = 0 或 1。输入空间是以汉明距离为度量的空间。 我想在数值上优化这个函数...
我在 np.gradient 中遇到了相对平滑的数据的高方差。 假设我们要计算 dx/dt： 将 numpy 导入为 np 将 matplotlib.pyplot 导入为 plt x = np.array([13.11149679, 13.2141427 ,...
我在 np.gradient 中遇到了相对平滑的数据的高方差。 假设我们要计算 dx/dt： 将 numpy 导入为 np 将 matplotlib.pyplot 导入为 plt x = np.array([13.11149679, 13.2141427 ,...
我正在开发一个用 JAX 编写的微分方程求解器。我遇到的常见工作流程是这样的： 将 jax.numpy 导入为 jnp 从 jax 导入 jit # 要集成的函数。 @jit 定义
我编写了一个 Gauss Seidel 函数，它接受一个系数矩阵、一个保存常量的向量数组和一个整数 n，其中 n-1 是矩阵的长度。它返回一个数组，其中包含...
我想获得 scipy 函数 scipy.optimize.curve_fit() 背后的一些直觉。下面，我写了一个简单的版本，试图最小化最小二乘函数之和（权重均为1），