来自集成/uniroot 的奇怪错误
所以，我正在尝试计算隐函数的积分fx01 fx01 我得到: uniroot(fx01, interval = c(0, 10), z0 = z0) 中的错误:端点处的 f() 值不是相反符号此外: 警告消息:1:在 if (is.na(f.lower)) stop("f.lower = f(lower) is NA") 中:条件具有长度 >1 并且只使用第一个元素2:在 if (is.na
..
numerical-methods相关内容
在Matlab中构建递归函数来绘制自适应梯形正交
我想建立一个接受容差水平、上下界和函数的函数，来计算自适应梯形正交，以及绘制一个图形，比如这个: 因为我需要节点值来绘制我的图形，所以我尝试编码如下: function [node, approx] = aq(f,a,b,tol)t = (b-a)*(f(b)+f(a))/2;如果 abs((t2 - t)/3) >tol %由于 T(2)-T(1)=E(1)-E(2)=3*E(2)m =
..
将 numpy ndarray 作为 fsolve 函数的输入传递
我想对非线性方程组进行数值求解并将 numpy ndarrays 作为输入传递.考虑下面的任意代码: 将 numpy 导入为 np从 scipy.optimize 导入 fsolve定义方程(A，B，C，D):eq1 = (A - B * np.sin(C)).tolist()eq2 = [5 * B + D * np.sum(A * np.cos(C))]返回 eq1 + eq2n = 3A
..
不确定如何在 Matlab 中使用事件函数
我正在尝试绘制动态系统的状态空间图和时程图.不过有一个问题.状态空间被位于 x1 = 0 的平面分成两半.状态空间轴是 x1、x2、x3.x1 = 0 平面平行于 x2/x3 平面.x1 = 0 平面上方的状态空间由 eqx3 中的 ODE 描述，而 x1 = 0 平面下方的状态空间由 eqx4 中的 ODE 描述. 所以，平面 x1 = 0 上存在不连续性.我有一个模糊的理解，即事件函数
..
如何使用 Runge-Kutta 四阶(Matlab)执行自适应步长?
对我来说，估计的 hstep 似乎需要很长时间和长时间的迭代才能收敛.我用第一个 ODE 尝试过.基本上，您执行具有 h/2 步长的 RK4 之间的差异.请注意，要达到相同的时间步长值，您必须在 h/2 的两个时间步长后使用 y 值，以便它也达到 h. frhs=@(x,y) x.^2*y; 我的代码正确吗? clear all;close all;clcc=[];我=1;U_saved=[
..
Scipy Fsolve 在具有解的非线性方程组上失败
我有一个由 12 个非线性方程组成的系统，如下所示: 导入数学将 numpy 导入为 np# 第一:定义一些之前已经计算过的变量(为了代码工作)n_s1=[8.71557427e-02, 9.96194698e-01, 6.12323400e-17]n_s2=[5.23359562e-02, 9.98629535e-01, 6.12323400e-17]P_sens1=[ 2.00000000e
..
Runge-Kutta 4th order 使用 Python 解决二阶 ODE 系统
我正在尝试通过 runge-kutta 4 阶方法在数值上求解两个 ode 的系统.初始系统:系统解决: 而且我有非常奇怪的解决方案图...我有: 正确的图形: 我在我的 runge-kutta 中找不到问题.请帮帮我. 我的代码在这里: dt = 0.04#初始条件t.append(0)zdot.append(0)z.append(A)thetadot.append(0
..
Goldberg 的 log1p 与 gsl_log1p
我正在寻找一个简单的 log1p 可移植实现.我遇到过两种实现. 第一个在这里作为定理 4 出现http://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html, 上面的实现 double log1p(double p){易失性双 y = p;返回 ( (1 + y) == 1 ) ?y : y * ( log( 1 + y
..
C程序——泰勒系列_长公式
这个公式来自我的一个朋友——我为他修正了它.但我似乎无法弄清楚如何获得每个角度的正确正弦计算.有人可以帮我在罪恶部分获得正确的命令吗? 代码: #include#define PI 3.141592653589#define NUMBER_OF_TERMS 10双阶乘(双 x){双计数器，总数；计数器=x;总计=x；而(计数器> 1){柜台 - ;总计 = 总计 * 计数器；}总回报；}
..
数值积分:为什么我的轨道模拟会产生错误的结果?
我阅读了费曼物理讲座第 9 章，并尝试了我自己的模拟.我使用黎曼积分来计算速度和位置.虽然所有的起点都一样，但我的轨道看起来像双曲线.这是讲义:https://www.feynmanlectures.caltech.edu/I_09.html(表 9.2) 导入时间导入 matplotlib.pyplot 作为 pltx=列表()y=列表()x_in=0.5y_in=0.0x.append(x_
..
数值积分方法问题python
我正在尝试整合以下公式: 下面是我尝试使用 f(a) = sin(a) 的自制方案来执行此集成. def func(x):返回 math.sin(x)def 积分 (f, n, r, a, dtheta ):总和 = 0θ = 0而 theta f(a) = sin(a) 的一阶导数 (n=1) 是 f'(a) = cos(a).当在 a = 0 处求值时，这应该给出 cos(0) =
..
如何阻止我的 Runge-Kutta2 (Heun) 方法爆炸?
我目前正在尝试编写一些 python 代码来解决一阶 ODE 的任意系统，使用由值 alpha、gamma(两个维度为 m 的向量)和 beta(下三角矩阵)定义的通用显式 Runge-Kutta 方法由用户传入的 Butcher 表的维度 (mxm).我的代码似乎适用于单个 ODE，已经在几个不同的示例上对其进行了测试，但我正在努力将我的代码推广到向量值 ODE(即系统). 特别是，我尝试
..
找到指数总和的解
我正在尝试用 Python 对这个方程进行数值求解(numpy/scipy，一切都可用) 在这个公式中，K 是一个常数，f 和 g 是两个依赖于 E 计数器的术语(这是一个离散表示一个积分)，其中 x 是我正在寻找的变量. 例如，E 是 3 个术语: 还有 f(E) 和 g(E) 是已知的. 我从 numpy 中了解到使用“fsolve"，但我不明白如何自动生成一个函数，该
..
耦合 ODE 的 Runge-kutta
我正在用 Octave 构建一个函数，它可以解决 N 类型的耦合常微分方程: dx/dt = F(x,y,…,z,t)dy/dt = G(x,y,…,z,t)dz/dt = H(x,y,…,z,t) 使用这三种方法中的任何一种(Euler、Heun 和 Runge-Kutta-4). 以下代码对应函数: function sol = Coupled_ode(E, dfuns, step
..
使用对数来避免数值下溢的算术问题
我有两个分数列表； 说 A = [ 1/212, 5/212, 3/212, ... ] 和 B = [ 4/143, 7/143, 2/143, ... ]. 如果我们定义 A' = a[0] * a[1] * a[2] * ... 和 B' = b[0] * b[1]* b[2] * ... 我要计算A'/B'的值， 我的问题是 A 和 B 都很长，每个值都很小，
..
使用 R 求解方程
我如何在 R 中以数值方式求解以下形式的方程的根: f(r)=r*c+1-B*c-exp(-M(B-r)) 其中 M、B 和 c 是已知常数. 提前致谢. 解决方案 由于 R 无法完成此功能，您可能需要使用像 Sage 这样的超集包.Sage 包含 R 和许多其他包，可以使用网络浏览器界面执行您想要的操作.该网站是 http://www.sagemath.org/ 示例位于
..
实三次多项式的最快数值解?
R 问题:寻找最快的方法来数值求解一堆已知具有实系数和三个实根的任意三次方.据报道，R 中的 polyroot 函数对复多项式使用 Jenkins-Traub 的算法 419，但对于实多项式，作者参考了他们早期的工作.实三次或更一般的实多项式的更快选项是什么? 解决方案 充实上面 Arietta 的回答: >a 这是否比在 GSL 包中使用三次求解器更快或更慢(如上面 knguyen
..
使用对数避免数值下溢的算术问题(取 2)
我有两个分数列表； 说 A = [ 1/212, 5/212, 3/212, ... ] 和 B = [ 4/143, 7/143, 2/143, ... ]. 如果我们定义 A' = a[0] * a[1] * a[2] * ... 和 B' = b[0] * b[1]* b[2] * ... 我想计算 A' 和 B' 的归一化值 即特别是 A'/(A'+B') 和
..
Python - 实现数值方程求解器(Newton-Raphson)
我警告你，这可能会让人困惑，我写的代码更像是一个思维导图，而不是完成的代码.. 我正在尝试实施 Newton-Raphson 方法来求解方程.我想不通的是如何写这个 Python 中的方程，从上一个近似值 (xn) 计算下一个近似值 (xn+1).我必须使用循环来越来越接近真正的答案，当近似值之间的变化小于变量 h 时，循环应该终止. 如何编写方程的代码? 当近似值不再变化
..
如何获得 R 中的机器 epsilon?
在 R 中是否有一个常量可以存储机器 epsilon? 解决方案 在我的 32 位 Linux 上试试 .Machine$double.eps -- 和 .Machine机器产生这个: R>.机器$double.eps[1] 2.220e-16$double.neg.eps[1] 1.110e-16$double.xmin[1] 2.225e-308$double.xmax[1] 1.7
..