differential-equations 相关问题

一个关联某个未知函数本身及其各种阶数的衍生物的等式。

PDE解决:'NoneType'没有属性'toarray'错误

我正在学习如何在Python中采用有限差分方法,以泊松方程为例。但是,当我运行时,一切都很好,直到A = ASparse.toarray()逆向代码行...

回答 1 投票 0

如何以相同的步长将Euler方法与二阶Runge-Kutta方法进行比较?

我有两种求解数值微分方程问题的算法,一种称为Euler方法,另一种称为二阶Runge Kutta(RK2)。本质上,欧拉方法和RK2近似为...

回答 1 投票 0

Numpy操作数不能与形状(3,)(2,)(3,)一起广播

我正在学习微分方程课程,并且代码是用于SIR模型的Runge Kutta方法。我不是一个出色的程序员,并且在调用中遇到dY_j + = dt * a [j,l] * ks [l]时遇到了麻烦...

回答 1 投票 0

如何求解该系统的IVP微分方程

我有一个微分方程组,例如dy / dt = f(t,y,y1,y2,....,yn,x),dy1 / dt = f(t,y,y1,y2,。 .yn,x)、. 。 。 dyn / dt = f(t,y,y1,y2,.yn,x),x(t)= f(t,y1,y2,... yn,x)并且我有值...] >

回答 1 投票 -1

我如何求解带有空气阻力的弹丸运动的二阶微分方程?

等式为:d ^ 2 r / dt ^ 2 = -c / m(dr / dt)+ g,其中r是弹丸的位置,c是阻力系数,m是弹丸的质量,g是重力引起的加速度。假设...

回答 2 投票 -1

我如何在PYTHON中求解带有空气阻力的弹丸运动的二阶微分方程?

等式为:d ^ 2 r / dt ^ 2 = -c / m(dr / dt)+ g,其中r是弹丸的位置,c是阻力系数,m是弹丸的质量,g是重力引起的加速度。假设...

回答 1 投票 -1

对于这种ODE系统,Julia中的微分方程为什么给我一个“无匹配方法”错误?

我想在Julia中解决由七个耦合ODE组成的系统。我的代码是:函数bsys!(dv,v,k,t)dp,thp,dr,thr,dx,thx,phi = va = exp(t)rp = r_y(a)rr = r_r(a)rx = r_x(a)...

回答 1 投票 0

绘制微分方程时,Bokeh JS回调

我对Python和Stackoverflow还是比较陌生,如果已经问过这个问题,我深表歉意,但是我已经搜索了很长一段时间,并且无法真正找到要解决的方法。 ...

回答 1 投票 0

使用fsolve和Euler 3BDF的SIR模型

嗨,我被要求使用MATLAB中的fsolve命令来求解SIR模型,而Euler 3向后指向。我对如何继续感到很困惑,请帮忙。这是我到目前为止所拥有的。我创建了一个函数...

回答 1 投票 1

odeint:无法根据规则“安全”将数组数据从dtype('complex128')转换为dtype('float64')

以下代码给出错误:无法根据规则'safe'从numpy.fft导入numpy,将数组数据从dtype('complex128')强制转换为dtype('float64'),从scipy.integrate导入fft。 。

回答 1 投票 0

绘制微分方程组的一个解,该微分方程的系统改变R中时间的参数方面

我一直在执行此代码来绘制仓本模型的解,我的目标是绘制一个改变耦合系数K的解,然后将所有图形放在一起以查看耦合的方式...

回答 1 投票 0

解决Python中的一阶和二阶微分方程组

我需要解决以下微分方程组:$ \ frac {dx_1} {dt} = -k_1x_1 + k_2x_2-(K_R)x_1y_1 $ $ \ frac {dx_2} {dt} = k_1x_1-k_2x_2-k_3x_2-( K_R)x_2y_2 $ $ \ frac {dx_3} {dt} = k_3x_3 $ $ \ ...

回答 1 投票 0

字典的写函数与solve_ivp一起使用

我正在尝试使用solve_ivp解决大型微分方程组。从scipy import集成def dXdt(X,t):返回np.array([dadt(X,t),dbdt(X,t),dcdt(X,t),dddt(X,t])sol =。 。

回答 1 投票 0

在GNU Octave中实施Euler方法

我正在阅读Chapra和Canale撰写的“工程师的数值方法”。在其中,他们为实现Euler方法(用于求解常微分方程)提供了伪代码。这是...

回答 2 投票 0

解决来自sympy的二阶微分方程组

我正在使用二阶拉格朗日方程组进行多自由度动力学问题。我用sympy来了解运动方程。现在这些方程在计算导数后变得相当长,...

回答 1 投票 0

如何求解和绘制R中的微分方程?

我想为指数增长求解并绘制一个微分方程,但是我不太了解如何使用deSolve库。我的方程是N = N_0 * e ^(rt),我尝试的代码是...

回答 1 投票 2

MethodError:没有方法匹配-(:: Int64,:: Array {Int64,1})

我尝试在Julia文档中尝试使用此示例。我的尝试是使细胞分裂成两个部分,每个部分的蛋白质量占一半。使用OrdinaryDiffEq constα= ...

回答 1 投票 4

Python求解线性方程组的子集

所以我知道sympy和numpy的线性代数求解器,但是我尝试做的是不同的。我没有方程式来形成可解矩阵,所以我可以使用这些库求解,所以我是...

回答 1 投票 0

OverflowError:(34,'结果太大')关于使用欧拉方法的微分方程的分辨率

我正在尝试解决一个物理问题的微分方程:高尔夫球掉落。我正在使用欧拉分辨率方法,并且我有以下代码:def F(Z,t):res = [Z [1],(-K / m)* Z [...

回答 1 投票 -1

求解边值问题集

我正在尝试解决由4个微分方程给出的一组边值问题。我在python中使用bvp_solver,并且收到错误消息,指出“除法中遇到无效的值”。我...

回答 1 投票 1

© www.soinside.com 2019 - 2024. All rights reserved.