一个关联某个未知函数本身及其各种阶数的衍生物的等式。
我制作了一个带有二次拖动的最小水平运动模拟器,可以在此处找到(通过左右箭头键控制)。作为练习,除了...
我正在尝试使用JuliaCall :: julia_eval在R脚本中的for循环中评估字符串。虽然我可以使用deSolve软件包在R中完成此操作,但是当...
我对Matlab还是相当陌生,并试图创建一个Matlab函数,该函数应用修改后的Euler方法来近似某个微分方程的解。当我在...
使用MATLAB的pdepe对(未)耦合的PDE产生奇怪的错误结果,时间加倍
[我正在尝试使用pdpe在1d中求解两个耦合的反应扩散方程,即$ \ partial_t u_1 = \ nabla ^ 2 u_1 + 2k(-u_1 ^ 2 + u_2)$ $ \ partial_t u_2 = \ nabla ^ 2 u_1 + k(u_1 ^ 2-u_2)$解决方案...
我有一个代码可以评估4个联立方程的参数值。当a + b(存储在results $ ab中)大于3000时,我特别想存储所有参数组合。...
边界条件在热方程和Crank-Nicholson的有限差分解中的应用
下面的代码求解一维热方程,该方程代表了一根棒,其端部在初始条件10 * np.sin(np.pi * x)下保持零温度。 Dirichlet边界条件(零...
我正在尝试使用Python中的odeint通过数值方式解决以下系统::我的问题是k1取决于p的值。当p <= 0.01时,k1 = 0,否则k1 = 5。我不确定如何写...
我想使用Matlab解决化学工程中的一些迭代问题。因此,我要解决的问题是计算串联绝热反应器的相互转换,这是...
我正在尝试计算方程式所描述的速度的时程:dV / dt = g −(C_d / m)* V ^ 2。 g = 9.81,m = 1.0,C_d = 1.5。为此,我需要在c ++中创建一个使用...
我的微分方程如下:%d / dt [x; y] = [m11 m12; m11 m12] [x; y] mat = @(t)sin(cos(w * t))m11 = mat(t)+ 5; m12 = 5; m21 = -m12; m22 = -m11;所以我的矩阵是...
解决散热式与隐式Euler和共轭梯度线性解法非零狄氏的BC?
许多用户询问如何解决热方程,u_t = u_xx,具有非零狄利克雷BC和用共轭梯度为内部线性解算器。这是一个常见的简化PDE问题...
我正在为数据集拟合分布。我需要根据给定的函数/公式调整分布的参数。我不知道如何解决R中的问题数据集......
使用FiPy对圆柱几何中的扩散方程的稳态解决方案与从另一软件(例如,Mathematica)获得的解决方案相当不同。等式是:$ 0 = \ ...
我想描述化学反应的动力学以及我对以下形式的微分方程中的反应模型结果(简化)的想法:y1'(t)= y1(t)+ y2(t)其中y1来自...
如何在matlab中使用interp1函数对symfun类型的函数?
我在matlab脚本中有一个微分方程组的解决方案。我将我的函数定义如下:syms V_x(t)V_y(t)X(t)Y(t); ode_V_x = diff(V_x,t)== - (B / m)* V_x; ode_V_y = diff(V_y,t)= ...
用于简单重力模拟的Python scipy.integrate.odeint失败
我正在尝试编写一个围绕原点旋转的质量的非常简单的重力模拟。我使用scipy.integrate.odeint来积分微分方程。问题是我得到了......
我正在尝试使用Python解决微分方程组。我编写了一个使用Euler方法的算法,我需要10 ^ -6 s-1的时间步长,持续100秒。那是10 ^ 8数据......
我的任务是随着时间的推移计算沿着1d杆的温度分布。我使用了这个矩阵,后面的时间等效公式就是这个。我的代码如下 - 我已经生成了...
我想制作一个微分方程,其中在一年中的特定时期存在青少年生长的峰值。在今年剩下的时间里,有一个背景出生率。我会 ...
我有一个ODE事件作为options = odeset('RelTol',1e-11,'Events',@ eventfunction); [时间,值] = ode45(@ Eq,时间,x0,选项);函数[value,isterminal,direction] = eventfunction(〜,y)value = y(1)+ ...