ode

3热度

1回答

我想寻求特定的ODE y'' - y' - 2y = 4x^2 我做了如下脚本y： syms x A0 A1 A2 ypa = A2*x^2+A1*x+A0; % y_p assume cyp = diff(ypa,2) - diff(ypa) - 2*ypa % according to ODE P1 = 4*x^2; P2 = cyp ; % Equating P1 and P2 C

1热度

1回答

Scilab中ODE的奇行为：方程dy/dx = A没有正确解决

我仍在学习Scilab（5.5.2），因此我正在编写并运行测试代码以熟悉本软件。要测试的数值微分方程解算器，我开始容易从方程DY/DX = A，其具有作为溶液Y =斧+ C（线方程）。这是我写的代码： // Function y = A*x+1 function ydot=fn(x, A) ydot=A endfunction A=2; //Initial condition

0热度

1回答

对每个网格单元应用一组常微分方程

我正在开发一个基于代理的模型来模拟传染病在由栖息地多边形（或连接单元的团块）组成的异质景观中传播疾病的传播。为了简化模型，我考虑了包含每个单元的多边形ID的栖息地网格（或栅格）。另外，我还有与每个多边形ID相关的流行病学参数。在每个时间步，参数值在多边形中变化。因此，在每个时间步骤更新数据帧landscape（见下文）。这里是在t = 0的例子： landscape <- data.frame(p

2热度

1回答

用scipy解决一个BVP solve_bvp

我有一个3个微分方程组（从我认为的代码中可以看出）具有3个边界条件。我设法在MATLAB中用一个循环来解决这个问题，如果它要返回一个错误，那么一点一点地改变最初的猜测而不终止程序。不过，在scipy的solve_bvp，我总是得到一些的答案，虽然是错的。所以我一直在改变我的猜测（它不断改变答案），并且对我从实际解决方案中得到的数据给出了非常接近的数字，但它仍然不起作用。也许还有其他一些代码问题，由

0热度

1回答

使用推力的ODE求解器的CUDA编程

我想用CUDA中的6个变量来解决使用推力的ODE。我的程序在这里。 #include <iostream> #include <cmath> #include <utility> #include <cstdlib> #include <thrust/device_vector.h> #include <thrust/reduce.h> #include <thrust/funct

2热度

1回答

使用scipy.integrate.odeint

当我试图解决的颂歌基本上类似于这一个，但多一个弹簧和阻尼器==>http://scipy-cookbook.readthedocs.io/items/CoupledSpringMassSystem.html 我有一个制度如何正确实施时间因变量但是，由于我想要实现的参数之一是时间依赖性的，所以存在一些小问题。我第一次尝试是下列之一： import scipy as sci import numpy

0热度

1回答

通过scipy文档试图解决颂歌的示例代码中的错误是什么？

我想了解如何解决python和scipy中的ODE，这似乎是一个很好的起点。我已经习惯使用odeint，现在正在尝试学习如何使用ode。我试着运行scipy docs中的示例代码，但它返回一个错误。我已将以下代码复制到错误旁边。 CODE from scipy.integrate import ode y0, t0 = [1.0j, 2.0], 0 def f(t, y, arg1):

0热度

1回答

具有小时间求解DDE for循环

步骤在我想在MATLAB使用for循环来解决以下DDE： xdot(t) = Ax(t) + BKx(t-h) 其中： A = [0 1 ; -1 0.1]; B = [0 ; 1]; h = 0.2; K = [-0.0469 -1.7663]; t = [0 5] 与常规方法解决，这是简单和结果是可以接受的。 sol = dde23(ddefun,lags,history,ts

0热度

1回答

用bvp4c Matlab求解边界值

我试图解决一个系统的边界值问题q' = f(q(t), a(t))与输入a在Matlab中使用bvp4c。其中q = [q1, q2, q1_dot, q2_dot]' 我的Matlab代码无法正常工作。有谁知道如何解决这个问题？ img: System Equation of the One Pendulum a是一个输入功能。初始状态：q1(0) = pi, q2(0) = 0, q1_do

-1热度

1回答

数据必须是单个矩阵Y或对的列表X，Y

我想使两条线图Ckkk和Ckk出现在同一图上。上述错误不断发生。我检查过类似的问题，没有任何修复程序有帮助。的代码如下 function cancer() clear all; clc; t0=0;tend=60;nt=tend; dt=(tend-t0)/(nt-1); t=t0:dt:tend; mambda=0.75;nambda=0.50;lambda=0.25;phi=0.