回答情况：
提问时间：
问题标签：
我不明白如何用python或python ode求解器来求解耦合PDE方程中的eta和V。(或者不需要解就可以对这些耦合方程进行数值解吗？)我花了好几天时间在这上面，但我还是不明白怎么开始！我理解下面的例子
我想画一个eta和V的时间序列，考虑到τ的不同时间序列的强制输入。
浏览 1修改于2020-04-03得票数 3
回答已采纳
1回答
这个问题问的是一个系统，它有两个质量，每个质量连接两个弹簧，看起来像这样：外部弹簧具有弹簧常数kb，而内弹簧具有常数k。我写了一些代码来寻找系统的正常模式，我得到的频率是0.160，我非常确定这是正确的。然而，当我试图用线性代数解决同样的问题时，我得到了-1.2和-1.0的特征值。我也不太明白如何编写代码，因为到目前为止，我的运动方程不适合任何python代码示例，我也不知道如何将矩阵d^2/dt^2X传递给eigh()函数。运动方程是：下面是我写的代码：
from scipy imp
浏览 2提问于2020-11-22得票数 0
1回答
.*2*r(1)*r(3));如你所见，r的每个坐标的微分方程是耦合的。我需要补充的是，每个方程都能够调用以前时间的其他坐标值，而不仅仅是它们的当前值。这样做的目的是构建延迟反馈机制。
浏览 16修改于2016-07-26得票数 1
2回答
我有一个微分方程如下：mat = @(t) sin(cos(w*t))m12 = 5;m22 = -m11 ;t = linspace(0,T,100) ; % Arbitrary 100[tf,xf] = ode45(@ddt,t,x0)
function
浏览 9提问于2017-11-10得票数 1
回答已采纳
我有一个耦合方程组:静水平衡方程、质量连续性方程和理想气体的状态方程。在数学语法中，
from scipy.constants import m_p,G,k,pi
from pylab import
浏览 10修改于2016-12-28得票数 2
回答已采纳
最初，我只是在系统确定性(无噪声项，即Zyi(t))并且得到正确的结果时，使用ODE45来求解，但当加高斯噪声时，不能再现精确的结果，即不能在不同的模拟中捕获频率灵敏度。然后在Matlab中使用SDEToolbox进行求解。我在工具箱中使用了Euler和Milstein的构建算法，但是没有用。
浏览 4提问于2015-05-27得票数 1
我一直在尝试使用odeint和solve_ivp来求解以下三个耦合微分方程系统。(这不是确切的系统。我对它进行了修改，只是为了提高问题的可读性。)稍后，我将把它推广到由数百个耦合方程组成的系统。
浏览 2修改于2020-09-28得票数 0
我正在试着解一组微分方程，但我一直很难做到这一点。我的微分方程包含一个表示从1到n的数字的"i“下标。我尝试实现一个forloop，如下所示，但我一直收到这个索引错误(错误消息如下)。我尝试更改初始条件(y0)和其他值，但似乎都不起作用。在这段代码中，我使用了solve_ivp。代码如下：import matplotlib.pyplot as pltfrom scipy.integrate import solve_ivp
X = y[0]
Y =
浏览 2提问于2021-02-25得票数 0
我对odeint有点迷惑。那么下面的表达式是指y[1] =y'和y'[1]= a*y[0]+b*y[1]吗？我有点困惑，因为表达式没有说方程的左边。这里有一个例子：from pylab import * # for p
浏览 0修改于2012-02-07得票数 4
回答已采纳
我想要解一个由七个耦合二阶微分方程组成的边值问题。有7个函数，y1(x),...y7(x)，每个函数都由一个形式的微分方程来描述。因此，我们可以把它重写为一个由14个耦合的一阶ODE组成的系统：dzi/dx = -(1/x)*zi - Li(y1,...我想用scipy.integrate.solve_bvp函数来求解这个方程组。但是，我很难理解文档()中描述的函数的确切输入参数应该是什么。
该函数需要的第一个参数是可调用的fun(x,y)。
浏览 6提问于2022-07-15得票数 1
回答已采纳
我们有一个10个ode的耦合系统。耦合出现在最后一个方程中。 我考虑使用一个10乘2的矩阵作为初始条件。我在这里也使用了一个具有相同标题的类似问题，但我仍然得到了相同的错误(“输入参数太多”)。-Is有没有其他的计算方法来表示耦合？
浏览 27修改于2019-09-10得票数 0
然而，如果我改变了公式中定义的等式的顺序，我的代码就会运行，我的图形也会发生变化。有人能告诉我这是为什么吗？现在我不知道这个系统的正确图形是什么。 def myFunction(r,t):L_L = 20 #draughtL_d = 4 #ukcw_d = 60 #width vesselb_R = 1.0b_D = 1.0M_s = 40000*10^3
w_r = r[1]u_2 =
浏览 26修改于2019-09-24得票数 0
我尝试过对非耦合方程这样做，但似乎也有一个问题。我必须准确地使用这个函数来解决这个问题，所以其他函数不是选项。
浏览 0修改于2019-05-30得票数 1
回答已采纳
这段python代码可以解决一个非耦合微分方程： import numpy as npimport numbax0 = 0h = 0.002
dxdt= x*3 + 3y 在上面的代码中，我如何将这两个系统实
浏览 46修改于2020-09-09得票数 1
回答已采纳
我有一个由2个方程和已知初始条件组成的问题。这些方程式是：dy/dt = (a1*sin(x) -a5 + a6*x + a7*sin(y-x)) / ((dx/dt)*a8*cos(y-x))我试图在MATLAB上绘制x与t和y与t的关系图，但我不确定如何用数值或解析的方法来解决这个问题。任何帮助都将不胜感激！！
浏览 0修改于2018-04-04得票数 0
用这些方程求解\phi(x，y，t)，p(x，y，t)，C_{a}(x，y，t)和C_{b}(x，y，t)，即空间和时间。我知道在Matlab中有一个函数pdepe( )用于求解1-D中抛物型椭圆偏微分方程的初边值问题。我想知道如何使用这个函数或Matlab中的其他函数来解决下面描述的问题，即二维和耦合问题。
浏览 3修改于2012-10-27得票数 4
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号