在Julia中求解初始条件为二维数组的常微分方程组
在Julia中，可以使用DifferentialEquations.jl库来求解初始条件为二维数组的常微分方程组。
首先，需要安装DifferentialEquations.jl库。在Julia的终端中，可以使用以下命令来安装：
using Pkg
Pkg.add("DifferentialEquations")安装完成后，可以使用以下代码来求解常微分方程组：
using DifferentialEquations
function myODE!(du, u, p, t)
du[1] = u[1] + u[2]
du[2] = 2u[1] - u[2]
end
u0 = [1.0, 2.0] # 初始条件
tspan = (0.0, 1.0) # 时间范围
prob = ODEProblem(myODE!, u0, tspan)
sol = solve(prob)
println(sol[end]) # 打印最终的解在上述代码中，myODE!是一个定义了常微分方程组的函数。du是微分方程组的导数，u是当前状态，p是参数，t是时间。在这个例子中，我们定义了一个简单的常微分方程组，并将其传递给ODEProblem。然后，使用solve函数来求解微分方程组，并将结果保存在sol变量中。最后，我们打印出了最终的解。
这里没有提到任何具体的云计算品牌商和产品，因为求解常微分方程组不依赖于特定的云计算平台或产品。这个方法可以在任何支持Julia语言的环境中运行，无论是本地环境还是云环境。
相关·内容
1回答
Julia是否支持边值微分代数方程？我有一个隐式的常微分方程，它的可变质量矩阵有时是奇异的，所以我必须使用DAEProblem。我的问题是x1(t)和x2(t)的两个耦合的二阶常微分方程，通过设置x1'(t) = y1(t)和x2'(t)=y2(t)，我已经将它们转换为四个一阶方程。我在域的开头和结尾有x1和x2的值，但没有y1或y2 anywhere的值，因此我需要DAE和BVP值。 This github的帖子建议这是可能的，但我担心我对机制的理解还不够深入，无法理解如何将DAEProblem与BVProblem结合起来。 我已经成功地按照数字食谱编写了多个拍摄代码来解决
浏览 30提问于2021-07-08得票数 2
我正在用pdepe解决一个二阶偏微分方程。我的问题的初始条件是一个二阶常微分方程的解，我用ode45解决了它。如何在pdepe中传递ode45解决方案作为我的初始条件？
浏览 1提问于2014-12-21得票数 0
Introduction
注:代码和指示：
单d是一阶导数双d是二阶导数
通过数值求解控制LaGrange方程，用Matlab对一些动态系统进行了数值模拟。基本上是一组二阶常微分方程。我正在使用ODE45。我从Mathworks (下面的教程链接)上找到了一个关于如何求解基本的二阶常微分方程组的很好的教程。
根据本教程，我通过获得两个二阶常微分方程(一个用于角θ，另一个用于弹簧伸长)来模拟弹性弹簧摆的运动，如下所示：
theta双素数方程：
M*thetadd*(L + del)^2 + M*g*sin(theta)*(L + del) + M*deld*thetad*(2*L + 2*de
浏览 4提问于2017-01-25得票数 0
目前，我确实使用来求解一个形式为0= F(t，y(t)，y'(t))的隐式微分方程。我想使用一个解决器，它附带了一个本机枕安装，并查看了的文档。文档中说，求解者解决了表单的一个系统：
求解一个系统的初值问题。该函数数值积分了给定初值dy / dt = f(t, y) y(t0) = y0的常微分方程组。
我确实怀念解决表单dy / dt = f(t, y, y')系统的能力。这与scipy.integrate.solve_ivp是可能的吗？
浏览 1提问于2018-06-11得票数 3
回答已采纳
引言
利用ODE45对二阶常微分方程组进行数值求解，用Matlab对一些动态系统进行了数值模拟。我从Mathworks那里找到了一个很好的教程(最后是教程链接)，介绍了如何做到这一点。
在本教程中，方程组在x和y中是显式的，如下所示：
x''=-D(y) * x' * sqrt(x'^2 + y'^2)
y''=-D(y) * y' * sqrt(x'^2 + y'^2) + g(y)
上述两个方程均为y'‘= f(x，x'，y，y')
问题
然而，我遇到的方程组中，变量不能显式地求解，如
浏览 4提问于2017-01-24得票数 1
1回答
我正在尝试使用julia提供的DifferentialEquations.jl，它工作正常，直到我尝试在二阶常微分方程上使用它。例如，考虑二阶常微分方程
x''(t) = x'(t) + 2* x(t)，带初始条件
x'(0) = 0，x(0) =1
它有一个解析解: x(t) =2/3exp(-t)+1/3exp(2t)。
为了从数值上解决这个问题，我运行以下代码：
using DifferentialEquations;
function f_simple(ddu, du, u, p, t)
ddu[1] = du[1] + 2*u[1]
end
浏览 5提问于2021-05-04得票数 3
这里有点新手。我试图解决一个由6个常微分方程组成的系统，但是有一个条件，如果(x(1)<=0.1)，那么x(6)=x(6)/2.0。我如何将其放入代码中？谢谢!
浏览 0提问于2015-07-04得票数 0
2回答
处理约束的PDE求解器
、、、、
我正在试着解一个一般形式的偏微分方程组。
F(f(x,y), f'(x,y), f''(x,y), g(x,y), g'(x,y), g''(x,y)) = 0
其中关于x和y的导数都可以取，并且f(x,y)和g(x,y)受某些约束
G(f(x,y),g(x,y)) = 0
我想知道是否有任何(最好是基于Python的)求解器(不是方法，因为我知道方法)可以处理这类问题？如果我的问题看起来很一般，我将非常感谢您的帮助和道歉。
浏览 0提问于2015-07-13得票数 1
2回答
我有一个由两个一阶常微分方程组成的系统，它们是非线性的，因此很难在封闭形式下解析求解。我想把这个常微分方程组的数值解拟合到一个数据集上。我的数据集只针对属于ODE系统的两个变量中的一个。我该怎么做呢？没有帮上忙，因为只有一个变量。
我当前导致错误的代码是：
import numpy as np
from scipy.integrate import odeint
from scipy.optimize import curve_fit
def f(y, t, a, b, g):
S, I = y # S, I are supposed to be my variables
S
浏览 0提问于2016-10-23得票数 1
1回答
我想用scipy_ivp()来解常微分方程，球在x方向以1的速度下落，y方向的速度为0。重力加速度g= 9.82，当球落地时，它的速度可能发生变化，并乘以0.9。但是，使用events参数，我发现它不能正常工作。这是我的代码，结果是：
from scipy.integrate import solve_ivp
def f_derivs_ivp(t, vars, g = 9.82):
dxdt = vars[2]
dydt = vars[3]
dvxdt = 0
dvydt = -g
return dxdt, dydt, dv
浏览 11提问于2021-12-15得票数 1
回答已采纳
2回答
我目前正在用deSolve来解一个常微分方程组，我想知道是否有任何方法可以防止微分变量值降到零以下。我看过一些关于在向量、数据框架等中将负值设置为零的帖子，但由于这是一个生物学模型(T细胞计数变为负数是没有意义的)，我需要从一开始就阻止它发生，这样这些值就不会扭曲结果，而不仅仅是替换最终输出中的负数。
浏览 5提问于2017-01-14得票数 8
回答已采纳
我有一个具有外部确定性输入(控制)和随机分量的常微分方程组。除了初始状态之外，我如何安全地(执行良好的代码风格)通过tf.contrib.integrate.odeint()将这些额外的输入参数传递给等式函数？如果有这样的方法的话。或者在外部作用域中定义它们，并从方程函数中引用它们是迄今为止唯一的方法？
浏览 8提问于2017-08-17得票数 0
1回答
当我使用牛顿方法寻找方程组的根时，我会得到以下错误消息：
警告:矩阵接近单数或标度差。结果可能不准确。RCOND = 1.625479e-33. >在新的(第16行)中
知道为什么我会收到这个错误信息吗？我知道奇异矩阵和逆矩阵，但这真的有什么用吗？如果是这样的话，是怎么做的？我能对代码做什么改变？它是第16行，它是dx = -J\f;。但我只是跟着我的数值方法教科书。一定是出了问题，但练习中写着“使用牛顿法”，所以我想这应该是可行的。我希望有人能帮我。
x = [0 0 pi/2]';
maxiter = 10;
iter = 0;
dxnorm = 1;
results
浏览 4提问于2015-11-25得票数 0
回答已采纳
我想把一个n阶常微分方程简化为一个一阶方程组。这是为数值分析做准备。我使用Sympy和Sagemath for Computer Algebra，但我没有在它们中找到任何函数来执行这种类型的约简。我不确定是否有其他人可以指出此功能是否存在于Sympy或Sagemath中。
这方面的一个例子是减少方程式：
x''' - 2x'' + x' = 0
一阶方程组：
[0 1 0
0 0 1
0 -1 2]
浏览 0提问于2014-08-31得票数 2
在MATLAB中求解约束线性和非线性方程组
我在MATLAB中用直接刚度法解决了一个有限元问题。这个问题现在被表述为一个非线性方程组：
KU=F or CF=U with C = K^{-1}.
该问题用3-Dd表示，有200多个节点。模型被用作变形估计器，力始终是变形的输入和输出。所以有很多线性方程。现在这个问题需要扩展，需要实现两种类型的约束。某些节点必须受到约束，以便它们只能围绕给定点以圆周运动方式移动。并且一些节点必须被限制在最大变形范围内。
总共有成百上千个形式为方程的方程：
a_1 X(1) + a_2 X(2) + ... - F(1) = 0
大约10个非线性形式的等式约束(也
浏览 3提问于2016-05-31得票数 0
我正在使用Matlab来尝试解决一个由三个一阶常微分方程组成的系统，但是我得到的错误信息是'syms‘需要Symbolic Toolbox。爬行器中出现错误(第1行) syms f(t) s(t) v(t) r W c h Q a k b H K e 在以前的情况下，我收到一个错误，说这个ODE系统不能显式求解(即以封闭形式)。我认为数值积分是唯一的方法。r、W、c、h等都是参数。有人能告诉我如何模拟/求解和绘制下面的ODEs吗？ syms f(t) s(t) v(t) r W c h q a k b H K e
r = 1;
W = 0.5;
c = 0.4;
h = 0.9;
q
浏览 34提问于2020-02-02得票数 0
回答已采纳
考虑以下边值问题：
y'' + e^y = 0 i.e. y(0) = y(1) = 0.
我很好奇MATLAB将如何解决这个特殊问题的有限差分法。我知道，如果我们有一个线性常微分方程，例如y'' + (e^x)y = 0，具有相同的边界条件，那么程序就相当简单。假设我们使用间隔0,1的分区为20个相等的子间隔，则以下代码将工作：
N = 19;
h = 1/N;
x = linspace(0, 1, N+1)';
A(1,1) = 1;
F(1) = 0;
for k=2:N
A(k,k-1) = -1/h^2;
A(k,k) = 2/h
浏览 6提问于2017-01-13得票数 0
4回答
我正在用c++编写程序来数值求解微分方程(常微分方程和偏微分方程组)。原因很简单--对于大规模计算，c++比matlab快得多。然而，当涉及到绘图时，matlab显然是赢家。有没有什么方法可以让matlab像在c++中绘图一样？
浏览 1提问于2014-05-26得票数 4
1回答
我想用Julia语言实现一个通用的Runge-Kutta步长函数。在Python中，我可以只传递一个函数作为这个RK4函数作为输入获取的参数之一。如果我在Julia中这样做，会有性能损失吗？
我的函数如下所示：
function uv_rk4_step(rhs,Vs,Ps)
Vs_k1 = rhs(Vs,Ps)
Vs_k1 = Ps.dt*Vs_k1
Vs_k2 = rhs((Vs+(1/2)*Vs_k1),Ps)
Vs_k2 = Ps.dt*Vs_k2
Vs_k3 = rhs((Vs+(1/2)*Vs_k2),Ps)
Vs_k3 = Ps.dt*Vs_k3
Vs_
浏览 2提问于2016-06-27得票数 2
我想使用bvpsolve来找到需要最终边界条件的常微分方程组的解。
这两个微分方程是一个动态优化问题的结果，边值来自于相关的横截性条件。
我的代码：
fun <- function(x, y, pars) {
dy1 <- -(d + n) * y[1] + y[1]^a - y[2]
dy2 <- y[2]/b * (a * y[1] ^(a-1) - d - r)
return(list(c(dy1, dy2))) }
# parameter value
n=0.01
a = 0.3
d=0.05
b = 3
r = 0.02
# initial and
浏览 0提问于2021-09-13得票数 0
我目前正在学习Matlab的常微分方程函数来解决简单的振动问题。
例如，mx''+cx'+kx=F*sin(wt)可以用以下公式求解
function dx = fun(t,x)
m=0.02; % Mass - kg
k=25.0; % Stiffness - N/m
c=0.0125; % System damping - Ns/m
f=10; % Frequency
F=5;
dx= [x(2); (F*sin(2*pi*f*t)-c*x(2)-k*x(1))/m]
然后调用ode45函数来获得位移和速度
[t,x]=ode45
浏览 6提问于2017-02-08得票数 0
回答已采纳
1回答
求解两个耦合二阶边值问题
、、、
用模solve_bvp求解了一个具有两个边界条件的单二阶微分方程。然而，现在我正试图求解二阶微分方程组；
U'‘+ a*B’=0
B'‘+ b*U’=0
在边界条件U(+/-0.5) = +/-0.01和B(+/-0.5) =0。我把它分解成一个第一常微分方程组，我试图用solve_bvp数值解它们。但是，我的解决方案只是得到满是零的数组。我相信我执行边界条件是错误的。我不清楚如何处理文档中的两个以上方程。我的尝试如下
import numpy as np
from scipy.integrate import solve_bvp
import matplotlib.pyp
浏览 5提问于2020-02-26得票数 1
回答已采纳
我如何拉链一个二维数组与“向量”行在朱莉娅？
这
X = [1 2; 3 4]
ndims(X)
Y = [-1 -2]
ndims(Y)
first(zip(X,Y))
给出(1, -1)，而我想得到([1 2], -1)。
浏览 2提问于2019-06-16得票数 3
拜托，我想用朱莉娅的NLsolve来解非线性方程组。包需要为系统的未知数设置首字母。用我的方程组(对不起，我不能在这里包括它，因为我至少需要10个名誉)，当我保留首字母0.1；1.2时，就像文档中的例子一样，我得到了“纸笔”的解。但是，如果我将首字母设置为1.1;2.2，则会收到以下错误：
DomainError:产生复杂结果的幂运算需要一个复杂的参数。将x^y替换为(x+0im)^y、复数(X)^y或类似。
请问，对于给定的方程组，我应该如何为首字母找到合适的值？
浏览 3提问于2017-08-15得票数 0
回答已采纳
1回答
有没有用傅立叶定律作为有限元的主方程并使用libmesh或deal.II库来解决二维热传导问题的例子？
浏览 3提问于2013-01-06得票数 0
回答已采纳
1回答
我正在尝试绘制一个常微分方程组的解。代码是：
tspan = [0 10];
z0 = [0.01 0.01 0.01 0.01];
[t,z] = ode45(@(t,z) odefun3(t,z), tspan, z0);
plot(z(:,3))
为什么输出绘制在区间0,60上，而不是代码中的0,10？
浏览 3提问于2018-06-10得票数 0
4回答
朱莉娅成对广播
、、、、
我想比较朱莉娅的字符串列表中的每一对字符串。一种方法就是
equal_strs = [(x == y) for x in str_list, y in str_list]
但是，如果我按以下方式使用广播：
equal_strs = broadcast(==, str_list, str_list)
它返回一个向量而不是一个二维数组。有没有办法用broadcast输出二维数组？
浏览 6提问于2022-09-02得票数 0
回答已采纳
1回答
我正试图用FiPy解决一个复杂的PDE问题。FiPY似乎不支持复杂的参数，一般的建议似乎是拆分PDE，单独解决，然后重新组合。现在还是这样吗？
其次，我看不出有任何方法来设置非常数初始条件。假设我希望我的网格在每个单元格中不具有相同的值，那么除了使用setValue之外，还有其他更简单的方法吗？
谢谢。
浏览 6提问于2014-01-23得票数 3
回答已采纳
4回答
假设idxs是定义为
3-element Vector{Vector{Int64}}:
[1, 2]
[1, 3]
[1, 4]
假设A是定义为
5-element Vector{Int64}:
6
7
8
9
10
使用Python和numpy，我可以(用伪代码)做A[idxs]，这将返回一个包含元素的二维数组。
[
[6 7]
[6 8]
[6 9]
]
如何在Julia中执行相同的索引，其中我将多次索引到一个一维向量中，并检索一个条目矩阵？
谢谢
浏览 6提问于2022-01-08得票数 3
回答已采纳
我有一个函数返回一个二维数组。如果我map这个函数到一个UnitRange，我得到一个二维数组的一维数组。有没有一种简洁的方法把这个收集到一个三维数组中？
y = rand(3,3)
f(x) = y * x
map(f, 1:9)
这将导致9-element Array{Array{Float64},2}而不是Array{Float64,3}。这也是通过理解来实现的。
浏览 0提问于2017-03-08得票数 2
回答已采纳
我想尝试使用bash shell解决线性代数方程组，但遇到了问题。输入值是一个矩阵，但是我不知道如何在bash中输入矩阵，因为它不支持二维数组。
非常感谢。
浏览 1提问于2012-10-03得票数 0
背景
我正在尝试使用ccall朱莉娅函数来使用用C编写的代码，我知道如何将数组作为参数传递给期望使用int *arg的函数。例如，尝试使用这个C函数
void sum_one(int *arr, int len)
{
for (int i=0; i<len; i++){
arr[i]++;
}
}
这个朱莉娅代码可以工作
x = collect(Cint, 1:5)
ccall((:sum_one, "/path/to/mylib.so"), Void, (Ptr{Cint}, Cint), x, 5)
问题所在
对于期望指针(int **arg)的指针被
浏览 1提问于2015-10-07得票数 8
回答已采纳
在Python中，有numpy.argmax
In [7]: a = np.random.rand(5,3)
In [8]: a
Out[8]:
array([[ 0.00108039, 0.16885304, 0.18129883],
[ 0.42661574, 0.78217538, 0.43942868],
[ 0.34321459, 0.53835544, 0.72364813],
[ 0.97914267, 0.40773394, 0.36358753],
[ 0.59639274, 0.67640815,
浏览 5提问于2016-01-13得票数 11
回答已采纳
1回答
我正在尝试用数值方法来解决下面的常微分方程 ? 使用边界条件 y(x=0)=0;y(x=20)=1 其中f(x)是我已经数值求解的另一个常微分方程的解。我正在尝试用https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.integrate.solve_bvp.html来解决我的ODE，但是我遇到的问题是f(x)被写成一个列表 coordinates x = np.arange(0,20.01,0.05) 并且我不知道如何在scipy.integrate.solve_bvp中插入列表来获得常微分方程的解。 提
浏览 11提问于2019-01-19得票数 0
回答已采纳
2回答
我试着用嵌套循环在二维数组中进行赋值。我尝试访问数组的元素，如下所示。但我弄错了。我找过了，但没找到结果。如何在二维数组中赋值Julia？
for x in 1:total
for y in 1:W
@show (x, y)
if agirliklar[x] <= y
V[x][y] = getMax(V[x-1][y], degerler[x] + V[x-1][y - agirliklar[x]])
else
print(&
浏览 1提问于2018-12-30得票数 0
我有这个二维数组：a = [ (1,2) (2,3) ; (3,4) (4,5)]，我想把它转换成： (1, 2)
(2, 3)
(3, 4)
(4, 5)
#1D Array 我怎么才能在Julia这样做呢？
浏览 52提问于2020-04-26得票数 2
回答已采纳
我在试着画出微分方程组的解，但我的解包括初始条件的时间导数，我不知道如何把点放在变量x.Does的上面，有人知道怎么做吗？谢谢。
浏览 0提问于2019-03-26得票数 1
我面临的问题如下:我想表示一个解决方案，其中给出了(t0,η)。比如说，(t0,η) = (0,1)。但是，我希望看到在整个(-1,1)域中表示的图形。
由于tspan与我的t0有着内在的联系，我如何使求解器解决给定t0的问题，而在tspan中的间隔中的初始值不同
因此，在给定的示例中，我想观察(-0.88, 0.88)中的ODE，但是要解决条件(t0, y0) = (0,1)的问题。
Ex8(c)
using DifferentialEquations
using Plots
function example2!(dy, y, t, p)
c1 = p
dy[1] = c
浏览 4提问于2020-12-20得票数 2
回答已采纳
我正在尝试使用Mathematica的NDSolve[]通过耦合的常微分方程来计算球面上的测地线：
x" - (x" . x) x = 0
问题是我只能输入x( 0 )和x'(0)的初始条件，求解器对x“=0的解很满意。问题是我在球面上的测地线有一个初始条件x"(0) = -x(0)，我不知道如何告诉mathematica。如果我将此作为条件添加，它会显示我正在将True添加到条件列表中。
下面是我的代码：
s1 = NDSolve[{x1''[t] - (x1[t] * x1''[t] + x2[t] * x2''
浏览 1提问于2010-06-17得票数 2
回答已采纳
1回答
人们可以使用Sympy来求解常微分方程。我在符号计算主题上的问题是:人们可以使用SymPy来求解一阶变系数方程吗？注意，这是ODE的一个特例。例如，如果我有一个像下面这样的等式 ? 如果可能的话，我该如何设置使用SymPy来求解这样的方程呢？
浏览 11提问于2020-11-16得票数 0
2回答
我有一道数学题，是用Maple解决的，但它必须用C++解决。问题是我找不到一个方程的数值方法(它是用内置的maple "solve“过程求解的)。我可以以某种方式将这个过程包含到我的代码中，或者找到显示Maple数值方法的完整描述的地方。这个方程是一个复杂的代数方程组。
浏览 4提问于2012-09-26得票数 1
目标
我想用Matlab求解一个9个未知数的非线性方程组.
所有9个未知数都被耦合成多个数的比例(见我的下面)
弗苏
x02=[5000,5000,5000,0.4,0.4,0.4,0.4,0.4,0.4];
ctrl2=[9894+1i*0.118,9894+1i*0.118,9894+1i*0.118,0.5,0.5,0.5,0.5,0.5,0.5];
f2 = @(x) myfun(x,Peff);
options2 = optimoptions('fsolve','Algorithm','trust-region-dogleg',
浏览 2提问于2017-01-30得票数 0
我有一个二阶微分方程系统，我想要解决它，我有使用变量y1，dy1等的符号形式。例如：
dy1 = y1 + y2
dy2 = 2*y1 + y2^3
为了使用Matlab中的常微分方程求解器，我需要一个'odefun‘，其中我以类似于以下形式表示一阶微分方程组：
function dydt = eqs(t,y)
dydt = [y(1)+y(2); 2*y(1)+y(2)^3]
现在，在我手头的问题中，符号方程式太长了，甚至不能在屏幕上正确显示：
dy1 = y1 + y2 + ... More than 25000 characters
dy2 = 2*y1 + y2^3 + .
浏览 17提问于2018-02-03得票数 0
所以我们试着解决简单的常微分方程。我们使用文本。这是本书(第18页)中的示例代码：
> with(DEtools); with(plots);
> eq1 := diff(x(t), t) = y(t); eq2 := diff(y(t), t) = -sin(x(t));
> init := [[x(0) = 1, y(0) = 0], [x(0) = 0, y(0) = 2]];
> DEplot([eq1, eq2], [x, y], t = 0 .. 110, x = -2*Pi .. 2*Pi, init, scene = [x, y], stepsize
浏览 1提问于2011-05-13得票数 0
我有一个问题，可以简化为下面的陈述(实际的方程和函数要复杂得多)：
系统：dot_X1 = -A*X1; dot_X2 = B*(X2-5) + u;
dot_表示与时间t微分，X1和X2分别是系统状态，A和B都是中间变量，u是系统控制输入，可以是自己定义的斜坡或阶跃输入。
我知道如何使用ode来解决正常的类似问题，当A和B不耦合时，例如，当A和B分别是X1，X2和u的函数时。但现在在我的问题中，A和B是耦合的，就像这样：
A = f(X1, X2, u, B), B = g(X1, X2, u, A),，f和g是非常复杂的函数。
现在，在解决常微分方程问题的每个时间间隔，例如，0，0.01s
浏览 3提问于2015-03-21得票数 1
我想在数值格式的每一步提取状态变量的值。
deSolve R包中的ode()函数，使用实现的ODE求解器之一对一个常微分方程组进行数值求解。为了做到这一点，它使用了一个动态调整的集成步骤，该步骤基于每个集成末尾的本地截断错误的值。用户基本上指定了所需的输出时间，它定义了时间步长网格，它可以有一个等于或小于输出时间的步骤。
如果我们以Lotka Volterra捕食者-猎物模型(含逻辑食饵)为例：
LVmod <- function(Time, State, Pars) {
with(as.list(c(State, Pars)), {
Ingestion <- r
浏览 2提问于2019-03-31得票数 2
回答已采纳
我试图用spicy.itegrate.odeint来集成一个微分方程组。
首先，在两个字典(x0和p)中对参数和初始条件进行采样和返回。然后创建模型并将其作为函数写入文件，大致如下(用虚拟方程)：
def model(x, t, p):
xdot = [
x['rate1'], p["a"]
x['rate2'], p["b"] * x["state1"] - p["c"] * x["state2"]
x['rate3&
浏览 2提问于2017-05-18得票数 1
回答已采纳
我正在尝试使用遗传算法来优化大约10个参数。我需要将main函数中参数的中间值传递给另一个函数，该函数基本上是线性常微分方程组的ODE求解器，参数是微分方程中的系数。然后，我会将ODE的解返回到我试图最小化的main函数中。我尝试将变量声明为全局变量，但似乎没有帮助。请建议如何最好地在两个函数之间传递变量。
浏览 1提问于2012-04-28得票数 0
我正在尝试用python求解一个模型，并将模型的未知参数与实验数据进行拟合。该模型由两个常微分方程组成，我使用scipy.integrate.solve_ivp求解它。模型的参数是未知的，因此，我想使用过多的方法来拟合它们。我的第一选择是differential_evolution，因为它可以为更复杂的模型提供非常好的结果(之前它是另一个包的一部分)。然而，问题是，当我给differential_evolution我的模型(我希望它找到计算点和实验点之间的最小二乘的全局最小值)时，它会找到模型变得不稳定的参数(由于内部步长为0，LSODA无法对其进行积分)。我试图捕捉LSODA在这些情况下抛出
浏览 4提问于2020-11-30得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
