不同步长的欧拉法。如何更改算法的代码以考虑步长的不同值？

文章/答案/技术大牛

发布

2回答

python、for-loop、matplotlib、integration、numerical-methods

我有一个数值微分方程问题的算法，叫做欧拉法。从本质上讲，欧拉法近似于微分方程的解。我的函数是针对单个步长(值为h)工作的，但我正在尝试更改代码，以允许我循环遍历3个不同的值h(通过将h从单个值更改为可能值的</em

浏览 28提问于2020-04-25得票数 2

回答已采纳

1回答

R中Lotka-Volterra码的数值求解

r、ode

所以最初，我开始使用一个简单的ODE系统(Lotka模型)，没有噪声项。的初始值为x= 10和y= 10。α、β、delta和γ的参数值分别为1.1、0.4、0.1和0.4 (模拟此)。)X <- sde.sim(X0=c(10,10), T = 10, drift=d, sigma=s) 然而，这似乎并没有产生一个良好的循环行为的捕食者和猎物群体预期输出我使用R中的</em

浏览 0提问于2019-02-12得票数 1

2回答

对于物理引擎的数值积分，有哪些好的算法？

c++、physics-engine

我已经在网上找了一段时间了，为一个物理引擎寻找集成的方法，我正试图编写一个有趣的代码(我喜欢那里的书呆子:P)。我找到了欧拉法、RK4和Verlet (以及时间校正版本)。令人惊讶的是，在我所有的研究中，我没有找到我想要做的事情的技术名称！也许描述我的具体问题会让我的问题更清晰。假设我想要模拟一个球在圆形(或球形，一旦我实现了3d)引力场中移动。显然，我不能在C++中使用无穷小<e

浏览 0提问于2012-02-24得票数 1

回答已采纳

2回答

MATLAB:在一个绘图中绘制不同的数值逼近

matlab、plot、numerical-methods

我有下面的matlab代码，用于通过欧拉法近似微分方程：b=0.6;dt=(b-a)/Steps; x(i+1,1)=x(i,1)+dt*(x(i,1)*x(i,1)+1);plot(x) 我希望能够在一个图中绘制几个不同步长的解图

浏览 0提问于2012-11-12得票数 1

回答已采纳

1回答

离散化微分方程并求出每个时间步中的参数值

matlab、numerical-methods、differential-equations

我想通过在t上使用for循环来找到f(t)的值。控制f的方程式是：x和y在2D中是坐标。我写了以下代码：f(1,1)= a*rand(1,1); f(1,t)= cos(f(1,t-1)) g(x,y,t) dt;但我不确定这种方法是否正确有没有人能告诉我它是正确的还是错误的？

浏览 0提问于2017-05-25得票数 0

2回答

是否有办法强制ode() [deSolve包]在ode函数的每个集成步骤上提供输出？

r、ode

我想在数值格式的每一步提取状态变量的值。 deSolve R包中的ode()函数，使用实现的ODE求解器之一对一个常微分方程组进行数值求解。为了做到这一点，它使用了一个动态调整的集成步骤，该步骤基于每个集成末尾的本地截断错误的值。用户基本上指定了所需的输出时间，它定义了时间步长网格，它可以有一个等于或小于输出时间的步骤。如果我们以Lotka Volterra捕食者

浏览 2提问于2019-03-31得票数 2

回答已采纳

1回答

Simulink:消除时滞

simulink

我试图运行一个闭环系统上的simulink，如下所示。当我运行它时，它会得到这个结果。正如您所看到的，使step函数达到1有一个轻微的时间延迟，闭环输出也要到1.5秒左右才会启动。我知道这是现实世界中会发生的事情，但我想知道是否有办法消除这个时间延迟，并使输出显示“理想的”结果。步骤输入是一个标准的步骤输入，我在设置上唯一改变的是模拟时间为8秒，而求解器是一个固定步长ode1 (euler)。

浏览 0提问于2017-05-22得票数 0

回答已采纳

1回答

Python :运行时警告:在正方形、加法、乘法、减法中遇到溢出

python、runtime-error、overflow

我从Python开始，所以我可能会问一个不那么微妙的问题，但经过相当多的研究后，我无法解决这个错误。我实际上正在尝试使用Gray-Scott模型来解决一个物理问题，但我被困在代码的最末尾:结果不被视为数字，并且在平方、加法、乘法和减法中遇到溢出。谢谢！以下是我试图解决的问题的初始条件：Du, Dv, F, k = 0.00016, 0.00008, 0.035, 0.065 T

浏览 3提问于2015-12-15得票数 1

回答已采纳

1回答

在NetLogo中求解ODEs，欧拉vs R-K vs R求解器

netlogo、runge-kutta

在我的模型中，每个智能体在每个节拍上求解一个ODE系统。我采用了欧拉法(类似于NetLogo中的系统动力学模型)来求解这些一阶常微分方程。然而，为了获得稳定的解决方案，我被迫使用非常小的时间步长(dt)，这意味着使用这种方法进行模拟的速度非常慢。我很好奇，有没有人有建议可以更快地解决这些常微分方程？我正在考虑实现Runge-Kutta (具有更大的时间步长？)正如这里所做的

浏览 0提问于2015-04-10得票数 1

1回答

爬山算法-最佳步长

predictive-modeling、optimization、hyperparameter-tuning、hyperparameter

我正在实现一个标准爬山算法，以优化预测模型的超参数。爬山算法正在作为两阶段方法的一部分加以应用：将大值应用于超参数的网格搜索，以找到“最佳”起点。将爬山算法应用于具有大量不同随机起点的空间中，找出局部最优解。我不认为对于网格搜索中传递的所有值来说，步骤大小的一个值是合适的，因为例如，1e-04和1e-03之间

浏览 0提问于2020-05-07得票数 0

回答已采纳

3回答

Runge (RK4)游戏物理集成

math、integration、physics

上对导数进行平均化，为什么平均到4阶的导数与用较小的时间步长进行简单的欧拉积分不同？可以用简单的术语解释RK4是如何工作的吗？ RK4利用了这样一个事实，如果我们使用函数的高阶导数，而不仅仅是一阶或二阶导数，我们就可以得到更好的函数逼

浏览 5提问于2009-11-03得票数 55

回答已采纳

2回答

使用python绘制地球轨道图时出现的问题

python、python-3.x、matplotlib

我正在尝试使用辛积分器在SI中编写地球轨道的代码，我的尝试如下：import matplotlib.pyplot as plt #Set parameters

浏览 1提问于2020-05-06得票数 1

2回答

我在C++中正确地实现了欧拉方法吗？

c++、numerical-methods、ode

我被要求写一个C或C++程序来解决the given differential equation问题这必须使用欧拉法在数值上实现。用户应该能够在program.It开始时输入速度(V)、x的初始值(0)和最终时间( T )，还应该绘制时间0 我觉得我让程序运行得很好，但在方程方面，我遇到了正确实现Euler方法的问题。这是我创建的程序。任何反馈/建议都将不胜感激。如果您需要更多信息，请让我知道。

浏览 37提问于2021-10-25得票数 0

回答已采纳

1回答

如何用欧拉法逐步逼近f'(x)

python、calculus

这是俄亥俄州立大学在线微积分课程的一个问题。很简单，但是这个老家伙想写一个Python程序，而不是长篇大论地解决它。把头撞在墙上：f(−3)=0。使用重复线性逼近(也称为欧拉法)在x=−1上近似此函数的值def fprime(a):f(-3) =4.5) = -2.25 which lead

浏览 1提问于2015-11-08得票数 1

回答已采纳

2回答

OpenMDAO不想脱离默认设置

openmdao

我有一些几何图形，我已经将其参数化到网格预处理器中，该处理器为我运行的外部代码生成输入。它的输出是由其输入定义的几何图形，调用我的外部求解器，然后从CSV文件中摄取输出。我已经让它运行了，但在优化方面肯定没有最小化函数的价值。事实上，它尝试过的值都稍微偏离了起始点。我有5个影响几何形状的设计变量，它们被设置为不同的范围，让我们称它们为a，b，c，d和e。每个变量的标称起始

浏览 8提问于2021-07-08得票数 0

1回答

在Dymola使用固定步长解算器时，如何选择步长？

modelica、dymola

我想做一个实时的模拟，如果我想要在Dymola使用固定步长的求解器，有不同的步长，结果可能会有点不同，那么有什么标准的步骤来选择步长吗？或者，我需要做大量的计算来证明步长独立，就像在CFD领域，我需要证明网格独立吗？

浏览 7提问于2021-02-18得票数 1

回答已采纳

3回答

如何计算算法的精确复杂度？

algorithm、numerical-methods、numerical-analysis

在不求助于渐近表示法的情况下，单调的步长计数是获得算法时间复杂度的唯一方法吗？如果不计算每一行代码的步数，我们能得到任何程序的大O表示吗？细节:尝试找出几种数值分析算法的复杂性，以确定哪种算法最适合解决特定问题。例如，在求解方程的Regula-Falsi或Newton-Rhapson方法中，意图是评估每种方法的确切复杂性，然后决

浏览 1提问于2010-07-28得票数 5

回答已采纳

1回答

实现牛顿-拉夫森迭代法

matlab、numerical-methods、newtons-method

我正在尝试使用牛顿-拉夫森迭代来实现一个后向欧拉方案。我知道，在每次迭代中，人们都会进行初始猜测，计算残差并解决变化。在我的例子中，更改为del w。之后，我知道将该值添加到w^m，并在下一次m次迭代中获得w的更新值。我知道随着迭代的进行，要检查解的收敛性。我遇到的问题是如何将时间步长dt实现为t=0: Tmax /dt，其中Tmax是10，我对时间步长

浏览 1提问于2011-02-25得票数 1

3回答

NSV，Euler集成混淆

math、language-agnostic、physics、game-physics

代码如下： vPos.y += vVelocity.y * timeDelta; 在欧拉法中，每次反弹时球都会反弹得更高。这是由于我的代码中的错误，还是由于该方法的不准确性。我读到过，在Euler集成中，

浏览 1提问于2011-06-23得票数 0

点击加载更多