differential-equations

6热度

2回答

我真的不知道从哪里开始解决这个问题，因为我没有太多的经验，但需要解决这个项目的这个部分一台电脑。我有一个2阶ODE是： m = 1220 k = 35600 g = 17.5 a = 450000 和b是1000和10000之间用500 x(0)= 0 x'(0)= 5 m*x''(t) + b*x'(t) + k*x(t)+a*(x(t))^3 = -m*g 增量

8热度

2回答

与分析解决方案相比，ODE45和Runge-Kutta方法的绝对误差

如果有人能够帮助解决以下问题，我将不胜感激。我有以下ODE： dr/dt = 4*exp(0.8*t) - 0.5*r ,r(0)=2, t[0,1] (1) 我在两种不同的方式解决（1）。通过Runge-Kutta方法（4阶），并通过ode45在Matlab中。我比较与分析解决方案，它由下式给出了两个结果： r(t) = 4/1.3 (exp(0.8*t) - exp(-0.5*t))

3热度

1回答

Routh-Hurwitz有用，当我可以计算特征值？

这是以下形式的线性常微分方程的N维系统的自学： dx/dt=Ax 其中A是该系统的系数矩阵。我知道你可以通过确定A的所有特征值的实部是否定的来检查稳定性。如果存在任何纯粹虚构的特征值A，则可以检查振荡。如果存在任何纯粹虚构的特征值A，则我在阅读的书中的作者引入了用于检测系统的稳定性和振荡的Routh-Hurwitz准则。这似乎是比计算特征值更有效的计算捷径。如果您现在可以快速找到特征值，那

4热度

1回答

确实提升了odeint是否有跳过算法？

我正在使用boost :: odeint，到目前为止我使用的是runge_kutta4步进器。现在我想切换到跳过方法，例如我的迭代步骤应该是这样的： f(t+dt) = f(t-dt) - p * f(t) 所以我需要一个多步骤的方法，但我有点失去了与文档，并希望得到一些帮助。

1热度

1回答

ode45没有足够的初始条件

我已经使用matlabFunction定义了一个函数。这里是代码： matlabFunction([a16;-((1+x16^2)/(2*x16))*a16],'vars',{x16,[a16]},'file','DE_19') 这似乎工作。但是，当我尝试使用ode45解决由matlabFunction定义的微分方程时，出现错误。下面是代码： [x,y] = ode45(@(x16,Y) DE

1热度

1回答

可以单独从Phabricator安装差速器吗？

我打算使用Phabricator进行团队代码审查，是否可以只安装和使用差异？我们需要的仅仅是代码审查工具。谢谢！

0热度

1回答

斜面情节衍生物

假设Mathematica中此代码： w[w1_] := NDSolve[{y''[x] + y[x] == 2, y'[0] == 0, y[0] == w1},y[x], {x, 0, 30}]; Manipulate[Plot[Evaluate[{y[x], y'[x]} /. w[w1]], {x, 0, 30}], {w1, -2, 3}] 的代码以这种方式工作：一个滑杆（在操纵）

1热度

1回答

Python：odeint解决可变系数（QHO）ODE

我试图通过使用odeint数值求解方程y'' + (epsilon-x^2)y = 0。我知道解决方案（QHO的wavefunctions），但odeint的输出与它没有明显的关系。我可以用常系数来解决ODE问题，但只要我选择了可变参数，我就无法解决任何我尝试过的问题。这里是我的代码： #!/usr/bin/python2 from __future__ import division i

2热度

1回答

没有足够的输入参数

我正在尝试创建一个SIR模型（求解带有两个参数的3个微分方程），并在MATLAB上绘制解决方案，并且不断收到错误'没有足够的输入参数'。我真的不知道自己出错的地方。下面是代码： function dx = sir(t,x) dx=[0; 0; 0]; beta = .003; delta = 1; dx(1)= -beta *x(1)*x(2); dx(2)=beta*x(1)*x(2)

3热度

2回答

numpy：评估矩阵中的函数，使用之前的数组作为参数计算下一个

我有一个m x n数组：a，其中整数m > 1E6和n <= 5。我有功能˚F和ģ，其由这样的：˚F（ü，ģ（ü，t））的。 ü是一个1 x n数组，t是一个标量，并且F和G返回1 x n数组。我需要˚F评估的a每个row，并使用先前估计的行为ü -array为接下来的评测。我需要做出这样的评估m。这必须非常快。我之前对整个数组的评估印象深刻，但是这个问题需要使用先前计算的数组作为计算下一个