Rs deSolve与Pythons odeint的区别
Rs deSolve和Python的odeint是两种不同的数值求解器，用于解决常微分方程（ODEs）的数值解。
Rs deSolve是R语言中的一个函数，用于求解ODEs。它提供了多种数值求解方法，包括基于龙格-库塔（Runge-Kutta）方法的经典求解器和基于Adams-Bashforth方法的多步法求解器。Rs deSolve在R语言生态系统中广泛使用，具有良好的稳定性和可靠性。
Python的odeint是SciPy库中的一个函数，也用于求解ODEs。它基于龙格-库塔方法，提供了一种简单而有效的求解器。odeint在Python科学计算领域非常流行，具有广泛的应用和支持。
两者的区别主要体现在以下几个方面：
- 语言：Rs deSolve是R语言中的函数，而odeint是Python中的函数。
- 生态系统：Rs deSolve是R语言生态系统中的一部分，而odeint是SciPy库中的一部分。因此，使用Rs deSolve可能需要熟悉R语言的其他功能和包，而使用odeint可能需要熟悉Python的科学计算生态系统。
- 求解方法：Rs deSolve提供了多种数值求解方法，包括经典求解器和多步法求解器。odeint基于龙格-库塔方法，提供了一种简单而有效的求解器。
- 应用领域：由于Rs deSolve和odeint分别属于R语言和Python的生态系统，它们在不同的应用领域中具有不同的使用情况。Rs deSolve在R语言生态系统中广泛应用于生物医学、生态学等领域，而odeint在Python科学计算领域广泛应用于物理学、工程学等领域。
对于Rs deSolve和odeint的具体使用方法和示例，可以参考以下链接：
- Rs deSolve：Rs deSolve官方文档
- odeint：SciPy官方文档
请注意，以上答案仅供参考，具体的使用方法和推荐的腾讯云产品可能需要根据实际情况进行调整。
相关·内容
2回答
你好，我想解决下面的第一个ODE：
dt/dr = +- cos(t)^2/cos(r)^2
我知道解是: t(r) = t(r) = arctan(tan(r)+_C1)，其中: pi/2 < t< pi/2和0< r< pi/2。
我想知道如何改进下面的代码，我的解决方案类似于图像中在t轴上趋向于+无穷大的曲线：
我的代码是：
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
"""
Eq
浏览 3提问于2018-11-06得票数 0
回答已采纳
1回答
Python中dict的定义是什么？
如果两者之间有差异，那么八度或MATLAB中的等效值是什么？
Python：
return dict(Fx = fxp, Fy = fyp, MO = MOp)
如有任何建议，我将不胜感激。
浏览 4提问于2017-01-13得票数 0
回答已采纳
1回答
我想知道，对于一个微分方程，ODEINT和solve_ivp之间有什么区别？它们之间有什么优点和缺点呢？
f1 = solve_ivp(f, [0,1], y0) #y0 is the initial point
f2 = odeint(f, y0, [0, 1], args=(a, b)) # a and b are arguments of function f
谢谢
浏览 2提问于2019-03-13得票数 5
回答已采纳
1回答
R中带函数工具的非负常数解？
、、、、
我正在尝试实现一个处理非负ODE系统的R算法。我需要像MATLAB中的ode45这样的东西来定义必须为非负的状态。 我在3年前就已经讨论过这个问题，但没有真正的解决方案。deSolve仍然不是我们要走的路。我发现了一些python代码，它们看起来很有前途。也许这在R中也是可能的。最后，我必须定义一个函数包装器，就像python中的functools一样。它的作用非常简单。以下是python包装器的代码： def wrap(f):
@wraps(f)
def wrapper(t, y, *args, **kwargs):
low =
浏览 7提问于2019-06-16得票数 0
1回答
我有一个R码，用DeSolve软件包的ode函数求解一个常微分方程。
f <- function(t, y, parms) {
with(as.list(c(y, parms)),{
rate <- M * A/(A + K) * A
dA <- -rate + D * (IN - A)
dB <- rate - D * B
Total <- A + B
Conserved<- dA + dB - D * (IN - A) + D * B
return (list(c
浏览 5提问于2022-10-20得票数 1
回答已采纳
我正在尝试将boost：：combine：：odeint与我自己的System类的实现结合起来(参见System.hpp)。
在BatchFilter类方法中使用(模板)系统对象，如下所示：
# BatchFilter.cpp
# template <typename state_type> BatchFilter class {...}
System<state_type> dynamics(this -> args,
N_true,
this -> true_dynamics_fun );
浏览 2提问于2017-12-16得票数 0
回答已采纳
我正在测试Eigen::VectorXd作为boost::odeint的state_type。我使用以下代码：
#include <Eigen/Eigen>
#include <boost/numeric/odeint.hpp>
#include <boost/numeric/odeint/external/eigen/eigen_algebra.hpp>
#include <iostream>
#include <vector>
template<class T>
struct push_back_state_and
浏览 0提问于2020-06-06得票数 0
1回答
我对Rust相当陌生，从Python的角度来看，有些事情是完全不同的。在Python中，可以通过键入.py文件从from foo import bar文件中导入单个函数，但我仍然没有在Rust中找到任何类似的函数。
我有以下文件：
.
├── main.rs
└── module.rs
内容如下：
main.rs
mod module;
fn main() {
module::hello();
}
module.rs
pub fn hello() {
println!("Hello");
}
pub fn bye() {
println!("
浏览 0提问于2018-09-18得票数 5
回答已采纳
我正在运行一个一直对我有用的代码。这次我在两个.csv文件上运行了它：“数据”(24 MB)和"data1“(475 MB)。"data“有3列，每个列约有680000个元素，而"data1”有3列，每个列有33000000个元素。当我运行代码时，经过大约5分钟的处理，我就会“被杀: 9”。如果这是一个记忆问题，如何解决呢？欢迎任何建议！
这是代码：
import csv
import numpy as np
from collections import OrderedDict # to save keys order
from numpy import genf
浏览 8提问于2015-12-14得票数 11
回答已采纳
1回答
我尝试在Debian挤压上运行boost_1_54_0中的g++4.4。
Lorenz系统工作良好，但简单的一维ode
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
/* we solve the simple ODE x' = 3/(2t^2) + x/(2t)
* with initial condition x(1) = 0.
* Analytic solut
浏览 3提问于2013-09-15得票数 3
回答已采纳
1回答
我调用了一个函数，该函数在每次通过for循环时都使用odeint (遗憾的是，我不能中断该循环中的任何内容)。但是，事情运行得比我希望的要慢得多。代码如下：
def get_STM(t_i, t_f, X_ref_i, dxdt, Amat):
"""Evaluate the state transition matrix rate of change for a given A matrix.
"""
STM_i = np.eye(X_ref_i.size).flatten()
args = (dxdt, Am
浏览 0提问于2016-03-27得票数 1
我试图使用ADO数据集向Visual数据表添加一个新记录，但没有运气。代码运行良好，没有任何异常，但是当我在没有新记录之后检查dbf时。代码片段中显示的mDataPath变量是整个数据库的.dbc文件的路径。关于底部的For循环的注意事项；我正在将收到的电子邮件正文添加到这个备忘字段中，因此我认为我需要将这个字符串添加到256个字符块中。如有任何指导，将不胜感激。
cnn1.Open("Driver={Microsoft Visual FoxPro Driver};" & _
"SourceType=DBC;" &
浏览 2提问于2010-06-16得票数 1
1回答
我想解一个矩阵微分方程，
import numpy as np
from scipy.integrate import odeint
def deriv(A, t, Ab):
return np.dot(Ab, A)
Ab = np.array([[-0.25, 0, 0],
[ 0.25, -0.2, 0],
[ 0, 0.2, -0.1]])
time = np.linspace(0, 25, 101)
A0 = np.array([10, 20, 30])
MA = odeint(d
浏览 17提问于2017-08-31得票数 1
回答已采纳
我有一个ODE，我想用从R的deSolve包调用的编译C代码来解决。所讨论的ODE是指数衰减模型(y'=-d* exp(g* time)*y)：但是在R内运行编译代码会给出R的本地deSolve的不同结果。就像在那里一样，它们被翻转180度。到底怎么回事？
C代码实现
/* file testODE.c */
#include <R.h>
static double parms[4];
#define C parms[0] /* left here on purpose */
#define d parms[1]
#define g parms[2]
/* initial
浏览 2提问于2022-03-12得票数 3
回答已采纳
我在编译这段简短的代码时遇到了问题：
更新:我有一个更简单的代码片段，现在产生了同样的错误，旧的仍然可以在帖子的末尾看到。
#include <vector>
#include <iostream>
#include <functional>
class Bar {
public:
void foo(double x, double y) {
std::cout << x << std::endl;
std::cout << y << std::endl;
}
void
浏览 2提问于2014-01-30得票数 1
我正在使用boost odeint来进化一个微分方程。时间和值是由包装类围绕一个普通数(进行高斯误差传播)给出的。
using Time = Number<double>;
using Value = Number<double>;
状态及其导数由带有声明的自定义类给出。
class State : boost::addable<State>, boost::multipliable<State, Time>, boost::dividable<State>
{
public:
void Add(Group const&
浏览 2提问于2017-06-15得票数 0
回答已采纳
我使用gprof2dot来分析我的应用程序：
./gprof2dot.py -f callgrind callgrind.out.x | dot -Tsvg -o output.svg
尽管如此，它给我提供了漂亮的图形分析，每个框中每个函数的名称都很长，而且远远超出了屏幕的大小。因为boost库对模板的使用率很高。只需查看其中一个函数名称：
std::reference_wrapper<boost::numeric::odeint::controlled_runge_kutta<boost::numeric::odeint::runge_kutta_dopri5, boost::
浏览 2提问于2015-02-10得票数 1
回答已采纳
我正在为新版本的boost而苦苦挣扎。我使用多精度的odeint。以下代码可以使用boost版本1.67.0成功编译。但是，从1.68.0版本开始，我就不能再编译了。在1.68.0版本中引入了复杂的多精度数字，我希望在我的代码中也使用这种数据类型。有人能告诉我代码中的错误在哪里吗？ #include <iostream>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>
#include <boost/multiprecision/gmp.hpp>
using na
浏览 40提问于2019-10-10得票数 2
回答已采纳
我最近偶然发现了boost.odeint库，我对它的可能性和可配置性感到惊讶。但是，由于广泛使用了scipy.integrate.odeint (它本质上是fortran中ODEPACK的包装器)，我想知道它们的性能如何比较。我知道boost.odeint还附带了并行化，这在scipy中是不可能的(据我所知)，这将大大提高性能，但我要求的是单一核心情况。然而，由于在这种情况下，我必须将boost.odeint (使用cython或boost.python)包装到python中，也许你们中已经有人这样做了？这将是一个伟大的成就，因为所有的分析可能性在python中都要高级得多。
浏览 0提问于2016-02-04得票数 5
我正在尝试将R包'deSolve‘安装在木星笔记本上，但我得到了以下错误：
> install.packages("deSolve")
Warning message in install.packages("deSolve"):
“installation of package ‘deSolve’ had non-zero exit status”Updating HTML index of packages in '.Library'
Making 'packages.html' ... done
当我尝试安装任
浏览 0提问于2020-11-01得票数 0
我正在尝试将boost库的odeint与std::vector或std::array一起使用。在这两种情况下，我对编译都没有问题，但是，在std::array的情况下，我有关于写入未知地址的运行时问题。下面给出的是MWE。 我在头文件ode.hpp中分别定义了两个类：ODE1、带有std::vector的ODE2和std::array。 #ifndef ODE_HPP_
#define ODE_HPP_
// c++ standard library headers
#include <array>
#include <iostream>
#include <
浏览 31提问于2021-09-03得票数 0
回答已采纳
让我们举个例子。在某些库中，像"scipy.integrate“这样的库调用像"odeint”这样的函数(集成函数)必须表示为"odeint(func，y0，T，...)“其中"func“是必须用两个参数(y，t0)预定义的函数的名称，y向量和t0是标量。
问题是，如果已经定义的函数"func“是用两个参数指定的，但顺序与"(t0，y)”相反，那么如何使用"odeint“。
诚挚的问候。
浏览 1提问于2017-04-25得票数 1
我目前正在使用Python和R探索Lorenz系统，并注意到ode包中的细微差别。Python和的都表示，他们使用lsoda来计算衍生品。但是，对这两个命令使用lsoda命令似乎会产生截然不同的结果。我对R中的ode函数尝试了ode45，以获得更类似于Python的东西，但我想知道为什么我不能获得完全相同的结果：
from scipy.integrate import odeint
def lorenz(x, t):
return [
10 * (x[1] - x[0]),
x[0] * (28 - x[2]) - x[1],
x[0]
浏览 49提问于2021-05-21得票数 2
回答已采纳
这是一个基本的示例方程，我试图将其拟合到一个示例数据中。 ? 我们的目标是找到最适合我的数据的k，假设数据遵循上面的等式。一种显而易见的方法是对方程进行数值积分，然后使用曲线拟合方法来最小化最小二乘，从而获得k。 使用odeint和ivp_solve进行数值集成，并在curve_fit上使用它们会产生相当大的差异。与较新的solve_ivp相比，较旧的odeint具有更好的适应性。k的最佳拟合值也有很大不同。 Best fit k using ivp = [2.74421966]
Best fit k using odeint = [161.82220545] 这背后的原因可能是什么
浏览 62提问于2021-10-10得票数 2
回答已采纳
1回答
我正在学习R来求解二阶微分方程(可能使用deSolve包)。它是我用python写成的两个一阶微分方程，下面给出。
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
def fun(X, t):
y , dy , z = X
M = np.sqrt (1./3. * (1/2. * dy **2 + 1./2. * y **2))
dz = (M*z) # dz/dt
ddy = -3.* M * dy - y # d
浏览 2提问于2017-09-06得票数 1
回答已采纳
我正在尝试使用Python来使用ODE系统来拟合曲线。然而，这两个方程都需要在任意时间点t上调用自身。如果我们的两个方程是 S'(t) = (a-b-c) * S(t) R'(t) = (a-b) * R(t) + mN(t) 我如何将这些传递给odeint？ 到目前为止，这段代码抛出了错误 UnboundLocalError: local variable 'R' referenced before assignment 我不知道如何解决这个问题。如有任何建议或帮助，我们将不胜感激。谢谢！ from scipy.optimize import minimiz
浏览 19提问于2021-08-15得票数 0
为了评估R和C++在求解ODE方面的速度差异，我在R中构建了以下ODE系统： modelsir_cpp =function(t,x){
S = x[1]
I1 = x[2]
I2 = x[3]
N=S+I1+I2
with(as.list(parm), {
dS=B*I1-mu*S-beta*(S*(I1+I2)/N)
dI1=beta*(S*(I1+I2)/N)-B*I1-lambda12*I1
dI2=lambda12*I1
res=c(dS,dI1,dI2)
return(res)
})
} 为了解决这个问题，我使用了deSolve包。 times =
浏览 42提问于2018-04-26得票数 1
回答已采纳
1回答
我是在Mac上，python是从homebrew安装的，它正在工作，我有两个定制的bash完成脚本，一个用于pip，另一个用于django。他们过去工作，但他们不再工作后，混乱的我的文件。
django完成来自，pip completion (来自)如下所示：
_pip() {
cur="${COMP_WORDS[COMP_CWORD]}"
prev="${COMP_WORDS[COMP_CWORD-1]}"
commands=$(pip --help | awk '/Commands\:/,/General Options\:
浏览 6提问于2014-05-19得票数 0
回答已采纳
我在来自deSolve的R中的ode求解器有问题。
我发现，对于应该总是等于零的变量，求解器给出了非零的中间值。在我开始尝试使用会触发的if(R>0){browser()}之类的测试进行调试之前，它看起来并不是一个问题。
我的密码在下面。提前感谢！
埃伦
library(deSolve)
simpleSIR <- function(t,states,par){
with(as.list(c(states,par)),{
S=states[1]
I=states[2]
R=states[3]
newinfections = beta*I
浏览 6提问于2016-06-06得票数 2
回答已采纳
1回答
我尝试在Mac 10.9.2 boost_1_55_0 5.1上运行[odeint复杂状态类型示例代码]。
下面的代码是解决斯图亚特-兰道振荡器的网站上的副本。
#include <iostream>
#include <complex>
#include <boost/array.hpp>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
//[ stuart_landau_system_func
浏览 0提问于2014-04-04得票数 1
回答已采纳
1回答
和在boost.odeint (boost常微分方程求解器)中有一些应用实例。
我想在不同的坐标(笛卡尔，极坐标或动作角)中使用odeint原子。
我应该对原子过载哪些操作？+，-，min，max，pow？
我可以在哪个文件中查找odeint使用的操作？
更新(1)来自的，它看起来需要"+“、"*”和abs()，max()
浏览 3提问于2012-04-11得票数 3
回答已采纳
我试图从铁锈中调用python代码(以便单元测试一些生锈代码)。所以我找到了CPython库，它允许这样做。
// src/main.rs
use cpython::{Python, PyDict, PyResult};
fn main() {
let gil = Python::acquire_gil();
}
//Cargo.toml
[package]
name = "python-rust"
version = "0.1.0"
authors = ["bashir ABDEL WAHED <bashir.abdelwahed@bi
浏览 3提问于2021-03-31得票数 0
我一直在使用WindowsSDK7.1附带的64位编译器VC++ 2010 Express。我的boost软件包版本是boost_1_55_0。
boost包odeint的示例代码可以找到，如下所示：
#include <iostream>
#include <boost/numeric/odeint.hpp>
using namespace std;
using namespace boost::numeric::odeint;
/* we solve the simple ODE x' = 3/(2t^2) + x/(2t)
* with init
浏览 1提问于2014-01-24得票数 2
回答已采纳
1回答
有没有人可以帮我解决我正在尝试使用boost odeint库进行计算的集成问题？这就是我所拥有的：
double z0 = -1.6;
double z1 = -1.0;
double dz = (z1 - z0) / 100;
double P(0.);
size_t xx = boost::numeric::odeint::integrate(my_integrand, P, z0, z1, dz);
其中my_integrand是定义为的函数(注意，目前我还没有完全定义函数计算，因为我希望它首先正确编译)：
void my_integrand
浏览 0提问于2017-03-21得票数 2
1回答
我在Debian下使用g++ 4.4，我已经安装了boost。我能够运行一个带有标题的简单boost矩阵：
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
但是我也想添加odeint头：
#include <boost/numeric/odeint.hpp>
但是，这是不可能的，因为我的/usr/include/boost/数值文件夹列表仅
conversion interval interval.hpp ublas
如何安装odeint？
浏览 5提问于2013-08-18得票数 1
回答已采纳
4回答
我想降低西皮的时间来解决一个微分方程。
为了实践，我使用了中的示例作为模板。由于odeint使用函数f作为参数，所以我将此函数编写为静态类型的Cython，并希望odeint的运行时间将显著减少。
函数f包含在名为ode.pyx的文件中，如下所示：
import numpy as np
cimport numpy as np
from libc.math cimport sin, cos
def f(y, t, params):
cdef double theta = y[0], omega = y[1]
cdef double Q = params[0], d = params[1
浏览 16提问于2017-03-16得票数 7
回答已采纳
我在编码围绕黑洞运行的光子。当光子直接进入黑洞时，由于被零除，轨迹真的很奇怪。我想忽略r< 1.5*rs的阵列中的光子，但我不知道如何忽略
我试过使用while True和if，它不起作用
h0=[t0, r2, theta, phi2, pt2, pr2, ptheta, pphi2]
T = np.linspace(0, 1000, 9000)
zz=odeint(func, h0, T, args=(rs,))
r22 = zz[:, 1]
theta22 = zz[:, 2]
phi22 = zz[:, 3]
pt22 = zz[:, 4]
pr22 = zz[:, 5]
浏览 3提问于2019-07-18得票数 1
1回答
我试图在SageMath中求解以下两个微分方程(数值)：
我的目标是获得M(r)-r.的图
我尝试了以下代码：
sage: r = var('r')
sage: M = function('M')(r)
sage: a = function('a')(r)
sage: de1 = (M*a*a*diff(M,r) + (M*M*a+6*a)*diff(a,r) + 1/(r*r) == 0)
sage: de2 = (a*r*diff(M,r) + 7*M*r*diff(a,r) + 2
浏览 3提问于2020-09-07得票数 0
回答已采纳
我有一个模板化的类，可以使用odeint的stepper类执行一些工作，我希望它是每种类型的stepper的特定(不同)工作。
/// header file ///
template<class Stepper>
class foo {
typedef typename boost::numeric::odeint::unwrap_reference< Stepper >::type::stepper_category stepper_category;
void do_work(double param);
// specific func
浏览 0提问于2013-12-13得票数 1
我正在尝试安装odeint，但是每次编译代码时，我都会得到这个错误。
c:...odeint\headmyshoulder-odeint-v2-f496df3\include\boost\numeric\odeint\config.hpp(44): fatal error C1083: Cannot open include file: 'boost/config.hpp': No such file or directory
我已经包括了包含路径。我检查了这个文件boost/config.hpp，这不是<config.hpp>的完整路径，所以我将它更改为<b
浏览 9提问于2014-10-31得票数 1
回答已采纳
2回答
我想求解并绘制指数增长的微分方程，但我不太明白如何使用deSolve库。我的方程是N= N_0 * e^(rt)，我尝试的代码是
library(deSolve)
## Time
t <- seq(0, 5, 1)
## Initial population
N0 <- 2
## Parameter values
r = 1
fn <- function(t, N0, r) with(r, list(N0 * exp(r*t)))
## Solving and ploting
out <- ode(N0, t, fn, params)
plot(out, l
浏览 4提问于2020-03-31得票数 3
回答已采纳
我正在将Python代码从2.7转换为Python3.8。
在Python2.7版本中，为了避免出现TypeError，我不得不使用降级版本的and和numpy (参见下文)。使用Python3.8，这些降级版本的fix和numpy不再可用，我得到了这个错误，我无法修复这个错误。
1.安装
先前: MacOS Catalina 10.15.7，Python2.7.16，numpy 1.9.0，Sply1.0.1
新版本: MacOS Catalina 10.15.7，Python3.8.6，numpy 1.20.3，ciply1.4.0
2.代码
当调用scipy.integrate.odei
浏览 4提问于2021-06-04得票数 0
回答已采纳
TL;DR:包安装在pip3下，但在Python3下找不到。为什么？
所有这些都发生在我的anaconda基地环境中：
因此，我一直在为tensorflow和它的版本(即将发布的另一篇文章)而奋斗，.Turns版本2.1只能在pip上使用，而不是在中。因此，在升级pip3 install --upgrade pip之后，我安装了pip3 install tensorflow==2.1.0。我打开木星-笔记本之后，发现tensorflow没有安装(运行Python3)。我先从终端检查版本，然后卸载tensorflow。它不是安装在pip下(正如预期的那样)，但它确实是在pip3下安装的。通过pi
浏览 4提问于2020-04-26得票数 1
回答已采纳
2回答
我正在使用scipy的方法integrate.odeint来求解一个二阶线性微分方程。该方法要求方程以两个一阶方程组的形式存在于两个未知数中。该方法
odeint(system_matrix,initial_conditions_matrix,time_values)
在time_values中的每个时间点输出解向量。解向量实际上是u，u‘的形式，其中u是我感兴趣的变量。所以我只想画出u。我在网上找到了一个实现这一点的方法是使用
u,u'=odeint(system_matrix,initial_conditions_matrix,time_values).T
但我不明白为什么会这样，
浏览 16提问于2017-01-14得票数 2
回答已采纳
1回答
我正在解形式的一阶初值问题：
dy/dt = f(t，y(t))，y(0)=y0
我想从给定的数值格式中获得y(n+1)，例如：
使用显式欧拉格式，我们有
y(i) = y(i-1) + f(t-1，y(t-1)) * dt
示例代码：
# Test code to evaluate different time integrators for the following equation:
# y' = (1/2) y + 2sin(3t) ; y(0) = -24/37
def dy_dt(y,t):
func = (1/2)*y + 2*np.sin(3*t)
浏览 2提问于2022-03-08得票数 1
回答已采纳
我用odeint来求解一个由4个耦合方程组成的系统，它模拟车辆行驶时的振动。我希望我的结果与我在MATLAB中得到的结果相似，但不幸的是，这种情况并没有发生。我已经多次检查我的方程，没有错误，所以这个问题必须发生在集成过程中。
我在MATLAB中编写了解决方案，以验证我从C++脚本中得到了什么。使用相同的条件，这是我从odeint获得的解决方案：
这是MATLAB中相同的解决方案：
我没想到在MATLAB中看到的微振荡会出现在odeint结果中，但大多数值甚至都不接近正确。我是否使用了错误的内存器，还是只对这个应用程序不起作用？
c++文件可以在Github上找到，名为"c
浏览 3提问于2016-12-02得票数 0
回答已采纳
1回答
我对以下使用Python (odeint)求解激光速率方程(一阶常微分方程)的小程序有问题。
当我运行程序时，总是出现一个错误：
index 2 is out of bounds for axis 0 with size 1
方程式是正确的；我不知道这个错误是什么意思。
有人能解释一下这意味着什么以及如何修复它吗？
def Rate(y,t):
D = y[0]
P = y[1]
w = 1.3
F = 0.0001
R = 90
dD_dt= w - D*P -D
dP_dt= R*D*P - R*P +F*D
dydt =
浏览 0提问于2017-06-02得票数 0
1回答
在我使用wxxima11.08.0(ubuntu12.04，Maxima版本: 5.24.0)时
我学习了一个，他的第二个例子对我没有用。
eq: y(t) = -r*c*'diff(y(t),t)+m*sin(%omega*t);
sol:desolve( eq, y(t) );
Is %omega zero or nonzero? nonzero
在我重新启动之前，Maxima就不再有反应了。
是否需要首先激活或定义maxima中的某些更改才能得到结果？
预期产出应是：
如果我用手定义方程的话，我的问题还有第二部分：
sol: y(t) = (m * sin(%omega
浏览 2提问于2013-09-29得票数 0
回答已采纳
1回答
当有阻力时，我试着表达子弹的轨迹。然而，我不能准确地表达这个图表。如何从常微分方程描述轨迹？
这是我的图表。这张图是不可信的。虽然我在设置vydot值的不同符号时遇到了困难，但它不能正常工作。
from pylab import*
from scipy.integrate import odeint
import matplotlib.pyplot as plt
import numpy as np
g=10
m=1
k=0.01
y=zeros([2])
vy0=0
vydot=200
vx0=0
vxdot=200
y[0]=vy0
y[1]=vydot
x=zeros([2
浏览 0提问于2015-04-09得票数 0
我正在使用scipy包中的Python和odeint来求解大量(~10e6)个耦合的常微分方程。方程系统可以表示为一些矩阵乘法的和，为此我使用了带有blas支持的numpy。我的问题是这需要很长的时间。当我分析代码时，我发现大部分时间都花在odeint上，而不是计算rhs。下面是来自分析器的五个最耗时的调用：
ncalls tottime percall cumtime percall filename:lineno(function)
5 1547.915 309.583 1588.170 317.634 {scipy.integrate._odepack.odein
浏览 3提问于2011-06-10得票数 5
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
