这是我在这里的第一个问题，所以如果格式化关闭了，很抱歉。import numpy as npimport matplotlib.pyplot as plt
r = odeint(model, r0, t)
plt.
浏览 22提问于2021-03-01得票数 4
1回答
我在思考数值分析中的一个基本问题。
众所周知，在离散常微分方程时，二阶方法比一阶方法更精确，因为二阶方法的截断误差为O(dx^2)，一阶方法的截断误差为O(dx^2)。例如，域是0到10000，网格大小是1000，那么dx = 10。在这种情况下，由于dx^2 = 100和dx = 10，二阶方法是否不像一阶方法那样精确？我们在处理大规模问题时会遇到这种情况，比如气候建模<
浏览 0提问于2016-11-18得票数 1
2回答
我在试着预测篮球的抛射运动。请注意，我没有考虑任何空气阻力的影响，只考虑了G的影响。为此，我使用了欧拉方法。def euler(euler_center, euler_velocity):[[331.3899068333333, 162.3366666666667]
[331.
浏览 0提问于2018-06-18得票数 1
Introduction单d是一阶导数双d是二阶导数
通过数值求解控制LaGrange方程，用Matlab对一些动态系统进行了数值模拟。基本上是一组二阶常微分方程。我从Mathworks (下面的教程链接)上找到了一个关于如何求解基本的二阶常微分方程组的很好的教程。根据本教程，我通过获得两个二阶常微分
浏览 4提问于2017-01-25得票数 0
1回答
微分方程- ODEINT
、、、
我必须用Python中的ODEINT来求解两个微分方程，它们是：z''(t) = a * (1/y(p) * y'y2z1' = z2
z2' = z''(t) = a *
浏览 0提问于2018-11-30得票数 0
http://diffeq.sciml.ai/latest/tutorials/ode_example.html 我正在尝试使用Julia (DifferentialEquations.jl)中的常微分方程求解器来求解假设系统是二维的，每个粒子的运动方程是由它相对于时间的位置的二阶常微分方程来描述
浏览 10提问于2019-10-14得票数 2
回答已采纳
如果我有一个用一阶微分方程表示的二阶微分方程，它描述了一个质量-弹簧-阻尼器系统，当我不知道一阶微分方程时，我怎么能用欧拉法来画这个方程呢？我想用MatLab来做这件事。我只想要一个简单的大纲，告诉你该怎么做。
二阶导数是d2(t+1) = (-1/m)*(c*d1+k*y)，其中c，m，k是常数，y初始是1，d1是从0开始的一阶微分，t是</e
浏览 1提问于2011-01-28得票数 1
回答已采纳
让我的RK2算法对某个二阶线性微分方程起作用，我遇到了一些问题。我已经在下面发布了我当前的代码(带有提供的参数)。由于某些原因，每次迭代时，y1的值与真实值的偏差会更大。我们将非常感谢您的任何意见。
浏览 0提问于2019-12-08得票数 0
我正在尝试对三体问题进行建模，其中我有三个广义坐标(一个径向坐标和两个角度坐标)和三个二阶(耦合)微分方程。我想看看系统在改变初始条件时是如何演变的，rho。我的脚本中有什么错误吗？python在不同的初始条件下求解常微分方程，将rho，theta，phi和V向量保存在orden中来操作它们。考虑到我想要身体的轨迹，对于这些rho和phi的值，我
浏览 2提问于2016-01-28得票数 2
我现在正在做一个项目，我已经对多体机械系统的运动进行了建模。这是一个复杂的、非光滑的系统，系统的运动并不总是由同一组方程描述，这意味着如果我用数值方法解决这个问题，我必须有时停止积分(并用新的初始条件和方程重新启动它)。现在，我想使用Matlab的事件函数在代码中实现这一点。这个概念很简单，但是像往常一样，当涉及到为实际问题实现它时，documentation有点难以描述。尽管如此，我还是写了一些代
浏览 9提问于2019-02-21得票数 0
回答已采纳
我是netlogo的新用户。我有一个反应系统(转化为常微分方程)，可以用Matlab求解。我想在netlogo中开发相同的模型(用于与matlab结果的比较)。我对time/tick感到困惑，因为netlogo使用"ticks“表示时间增量，而Matlab使用秒为单位的时间。如何将我的matlab秒转换为刻度数？有人能帮我写代码吗？模型是：
A + B ---> C (wi
浏览 2提问于2014-12-30得票数 1
我有一个二阶微分方程，我想用python来求解它。问题是，对于其中一个变量，我没有0的初始条件，只有无穷远的值。有人能告诉我应该为scipy.integrate.odeint提供哪些参数吗？它能被解决吗？Theta需要在时间上找到。它的一阶导数在t=0等于零。theta在t=0中不为人所知，但它在足够大的时间内会变为零。其余的都是已知的。因为近似
浏览 4提问于2013-01-21得票数 5
我在excel电子表格中得到了这些数据。所以导入之后，我只需要做velocity = cumtrapz(t,y)，然后做position = cumtrapz(velocity)
浏览 0提问于2014-04-25得票数 0
回答已采纳
我有一个二阶微分方程系统，我想要解决它，我有使用变量y1，dy1等的符号形式。例如：dy2 = 2*y1 + y2^3functiondydt = eqs(t,y)现在，在我手头的问题中，符号方程式太长了，甚
浏览 17提问于2018-02-03得票数 0
我有一个加速计，我想用它来检测我的宝宝在晚上的运动，有没有人试图写代码来实现类似的东西，任何指针。import adxl345accel.setRange
浏览 0提问于2014-10-11得票数 0
我目前正在学习一些拉格朗日力学，这在很大程度上进行得很好，但我决定尝试一个计算机练习，我不知道如何让Python求解这些方程。我试图解决的二阶非线性常微分方程是：Y的初始值是y0 = 0。
求解一个普通的一阶方程很容易，因为你只需创建一个函数，将它设置为某个值，然后使用odeint。没问题。但是现在我得到了一个二阶方程，我不知道如何告诉odeint或函
浏览 26提问于2019-08-05得票数 2
我已经在网上找了一段时间了，为一个物理引擎寻找集成的方法，我正试图编写一个有趣的代码(我喜欢那里的书呆子:P)。我找到了欧拉法、RK4和Verlet (以及时间校正版本)。令人惊讶的是，在我所有的研究中，我没有找到我想要做的事情的技术名称！也许描述我的具体问题会让我的问题更清晰。假设我想要模拟一个球在圆形(或球形，一旦我实现了3d)引力场中移动。这个球将遇到力向量，它可以用来计算球在该特定刻度上所处
浏览 0提问于2012-02-24得票数 1
回答已采纳
问题是罐头上写的是什么。我正在尝试用数值方法解决一个边值问题，我的朋友问我，求解器是否适用于这些条件。求解器的页面没有给出我们拥有的条件，即bc(y(a),y(b), p) = 0，但是我们问题的形式是y(0) = some constant value和y'(b) = 0，给出了我们的Neumann条件，你需要重写函数来有一个像打靶法一样的一阶降阶吗？
浏览 29提问于2019-12-09得票数 1
回答已采纳
2回答
磁学模拟
、、
假设我在一个n*m像素的2D曲面上有p个节点，我希望这些节点相互吸引，这样它们之间的距离越远，吸引力就越强。但是，如果两个节点之间的距离，比如说d(A，B)小于某个阈值，比如k，那么它们就开始相互排斥。谁能给我一些关于如何随着时间的推移更新节点坐标的代码。我的东西有点像下面的代码，这是开始做吸引，但寻找一些建议。(附注:我不能使用现有的库来完成此操作)。
浏览 0提问于2010-01-25得票数 3
回答已采纳
我应该写一个MATLAB函数，它取一阶常微分方程，形式为y‘(T)= a*y(t) +b，初始点y(t0)=y0作为输入，并计算解的前15个点。还会绘制前15个点的解决方案曲线。我们要解的方程是:y‘(T)= 4*y(t)+1，初始点y(0)=0。
对于这个函数，我写了下面的代码，但这给了我一个关于y的错误。我应该如何正确地实现euler函数？
浏览 0提问于2012-12-17得票数 0
回答已采纳
云服务器
ICP备案
云直播
对象存储
实时音视频
腾讯云开发者
