如何修改这个微分方程求解器来求解大量的变量？
要修改微分方程求解器来求解大量的变量，可以考虑以下几个方面的改进：
- 优化算法：针对大量变量的求解，可以尝试优化求解算法，提高求解效率。例如，可以采用并行计算技术，将求解任务分配给多个处理器或计算节点进行并行计算，从而加快求解速度。
- 内存管理：大量变量的求解可能会导致内存占用过高，甚至超出系统的内存限制。因此，需要对内存进行合理管理和优化。可以考虑使用稀疏矩阵等数据结构来降低内存占用，或者采用分块计算的方式来减少内存压力。
- 缓存机制：对于重复计算的变量，可以引入缓存机制，将已经计算过的结果保存起来，避免重复计算，提高求解效率。
- 并行计算：如果系统具备多核处理器或者分布式计算环境，可以考虑将求解任务分解为多个子任务，并行计算。这样可以充分利用系统资源，加快求解速度。
- 数据分析与优化：对于大量变量的求解，可以通过数据分析来发现其中的规律和特点，从而优化求解过程。例如，可以通过特征选择、降维等技术来减少变量数量，或者通过数据预处理来提高求解效率。
- 腾讯云相关产品推荐：腾讯云提供了一系列云计算相关的产品和服务，可以帮助进行大规模计算和数据处理。例如，可以使用腾讯云的弹性计算服务（Elastic Compute Service，ECS）来部署求解器，并根据需求选择合适的实例类型和规模。此外，腾讯云还提供了云原生数据库TDSQL、云存储COS等产品，可以满足大规模数据存储和管理的需求。
请注意，以上答案仅供参考，具体的修改方法和推荐产品需要根据实际情况进行评估和选择。
相关·内容
我有，我想解决X(t)，Y(t)，Z(t)，希望用我熟悉的Octave/Matlab，但我不介意用任何其他必要的方法来解决它。
现在，Fsolve对于正则方程组是有用的，而Ode45，Lsode对于微分方程是有用的。但是，这个特殊的系统呢？请注意，底部的微分方程不仅包含Y，还包含X和Z，它们依赖于上面的两个非微分方程。
老实说，我不太确定如何使用基本的代码来解决这个系统，在花了一些时间思考之后，我决定寻求帮助。如果可能的话，我真的很感谢任何能有效地解决这个问题的指导。现在几乎所有的回复对我来说都是有用的。
浏览 9提问于2020-02-15得票数 0
回答已采纳
1回答
这个问题不是关于如何使用matlab，而是试图找出当matab求解微分方程时发生了什么。当使用不同的数值方法时，结果通常是不同的。我想知道matlab中使用的是哪种数值方法。
浏览 5提问于2020-07-11得票数 0
我需要对以下形状的两个微分方程进行交互式模拟：
der(x)= A * x + c1 * y
der(y)= c2 * x + c3 * y
其中der(x)是x的时间导数
但我需要能够在模拟实时运行时更改A。
我是matlab的新手，但我稍微摸索了一下，所以我并不完全是新手
浏览 1提问于2013-11-11得票数 0
我有一个微分方程A*dx/dt + B(y-y0) =0，其中x是y的一个非常复杂的函数。
我如何用Mathematica重新排列y，得到一个函数x来解决这个问题？
谢谢
浏览 4提问于2014-05-29得票数 2
1回答
问题
我现在有一个由四个方程组成的系统。两个是二阶微分方程，两个是一阶微分方程：
最初的条件是：
x = 0 |
y = 0.3 |
f(x) = 2.05 |
f(y) = 0.55 |
所有的角度都是度的。
我已经尝试过的
我曾尝试使用，并与SciPy和NumPy合作。不幸的是，我不知道如何编程，因为这些方程是非线性的。有人能给我一些其他模块的提示吗？
浏览 4提问于2020-10-15得票数 1
回答已采纳
1回答
我试图画出微分方程dh/dt = I(t) - f(h)*h(t)，其中f(h)可以用f(h) =1-e^(-kh)建模，k是常数。我该怎么做呢？
浏览 4提问于2022-05-03得票数 -1
1回答
NumPy中的变量
、、、
在NumPy中，我尝试表示以下形式的微分方程：y' = p(t)y + g(t)，其中p(t)是n×n矩阵，g(t)是nx1矩阵。类似于：
y' = [[1,5], [2,1]] + [[e^t], [1]]
我知道如何在NumPy中表示矩阵，但是如何表示包含变量(例如，2t或e^t)的矩阵呢？
浏览 4提问于2013-04-30得票数 2
回答已采纳
我有一个问题，可以简化为下面的陈述(实际的方程和函数要复杂得多)：
系统：dot_X1 = -A*X1; dot_X2 = B*(X2-5) + u;
dot_表示与时间t微分，X1和X2分别是系统状态，A和B都是中间变量，u是系统控制输入，可以是自己定义的斜坡或阶跃输入。
我知道如何使用ode来解决正常的类似问题，当A和B不耦合时，例如，当A和B分别是X1，X2和u的函数时。但现在在我的问题中，A和B是耦合的，就像这样：
A = f(X1, X2, u, B), B = g(X1, X2, u, A),，f和g是非常复杂的函数。
现在，在解决常微分方程问题的每个时间间隔，例如，0，0.01s
浏览 3提问于2015-03-21得票数 1
给出了带a(t) = k*y'^2的微分方程a(t) = k*y'^2，其中y是t (时间)的函数。我的初始条件是y(0) = 600;和y'(0) = 0;
在MATLAB中，我知道如何用
ydd = diff(y,t,2) == -g + a(t)/m;
但是我不知道这是一个“嵌套的”非线性微分方程，我不知道如何定义它，更不用说，用MATLAB来解决它了。
浏览 5提问于2016-12-15得票数 0
回答已采纳
有没有人尝试用Modelica实现Navier Stokes偏微分方程(PDE)？我发现了空间基函数(SBF)的方法，它通过数值修改得到Dymola可以处理的常微分方程(ODE)。
致以敬意，
维克多
浏览 2提问于2015-06-17得票数 4
1回答
📷
在pbrt v3中，这本书给出了光束透过率的描述，但是我不知道如何像它所说的那样求解微分方程来得到Tr，有人能告诉我如何求解这个微分方程吗？非常感谢。
浏览 0提问于2021-11-11得票数 1
回答已采纳
1回答
我正在使用Comsol，想解决一些等离子问题，但默认的等离子模块不适合我。
我需要在一维轴对称几何中求解PDE方程。但!方程的参数是从玻尔兹曼方程(这是另一个组成部分)导出的。我试图在第一阶段计算Boltzmann方程，然后导出所需的参数，然后将它们导入到我的新模型中，但我没有成功。
有没有办法将两个组件连接起来？
我想让Comsol来求解我的偏微分方程，当需要的时候，解玻耳兹曼方程，然后把导出的参数放入偏微分方程，以此类推。
谢谢!
浏览 28提问于2017-07-10得票数 0
3回答
我这里有一个齿轮系统，我需要用MATLAB和Simulink进行仿真。
要做到这一点，首先我需要得到这个旋转(齿轮)系统的微分方程。我认为成功地做到了，如下所示
iT1 = J * theta2''+D * theta2'+K * theta2,
其中J = i^2 * Jn1 + Jn2 + J2 (同时，i是齿轮比；i = theta1/theta2 = N2/N1 = Tn2 / Tn1 = T2/T1)。
现在我需要用MATLAB中的ode23来求解这个二阶微分方程。问题是，我不知道任何初始条件，也不知道如何从这个特定的系统中得到它们。例如，y(0)
浏览 6提问于2021-02-01得票数 0
我的任务是为某个物理问题建模，并使用matlab来求解它的微分方程。我做了这个模型，但它似乎比我到目前为止学到的要复杂得多，所以我不知道如何解决这个问题。
黑色表示它是一个常量
浏览 0提问于2019-08-19得票数 0
1回答
我正在使用GEKKO来估计微分方程的参数，并且我将其中一个变量限制在0到1之间。但是，当我求解常微分方程时，我得到的值超出了这个变量的界限，所以我想知道是否有人知道GEKKO是如何找到解的，因为这可能会帮助我解决问题。下面是我用来拟合数据的代码。这给出了一个解x和u，其中x在0和1之间。
然而，之后，我尝试用scipy.integrate.solve_ivp求解这个常微分方程，用我得到的u的初始值，我得到的u的解不在这个范围内。由于它应该是唯一的，我想知道GEKKO寻找解决方案的过程是什么(它是将值投影到边界还是如何处理？)任何意见都是非常感谢的。
这是一个MVCE。如果你运行它，你可以看到，
浏览 7提问于2020-11-04得票数 2
目前，我确实使用来求解一个形式为0= F(t，y(t)，y'(t))的隐式微分方程。我想使用一个解决器，它附带了一个本机枕安装，并查看了的文档。文档中说，求解者解决了表单的一个系统：
求解一个系统的初值问题。该函数数值积分了给定初值dy / dt = f(t, y) y(t0) = y0的常微分方程组。
我确实怀念解决表单dy / dt = f(t, y, y')系统的能力。这与scipy.integrate.solve_ivp是可能的吗？
浏览 1提问于2018-06-11得票数 3
回答已采纳
1回答
我正在实现一个基于论文的卡尔曼滤波器。状态空间向量有6个变量，因为状态变量随时间的变化而变化，文中还给出了这些变量随时间(dt)的变化的微分方程。我的问题是，在python中实现这个概念时，我应该如何使用微分方程和dt，我想知道简单地将方程与dt相乘，但我认为这不是解决方案，因为它是一个微分方程。请在这方面给我指点。另外，我想问一下，由于这些方程将估计新的状态，所以它们应该被添加到更新步骤中吗？谢谢!
浏览 25提问于2020-12-16得票数 0
回答已采纳
1回答
我想找到一种优雅的方法来解决以下微分方程：
from sympy import *
init_printing()
M, phi, t, r = symbols('M phi t r')
eq = Eq(-M * phi(t).diff(t), Rational(3, 2) * m * r**2 * phi(t).diff(t) * phi(t).diff(t,t))
假设phi(t).diff(t)不是零。因此，左边和右侧都缩短了。
我就是这样找到解决办法的：
# I assume d/dt(phi(t)) != 0
theta = symbols('
浏览 3提问于2015-10-26得票数 2
回答已采纳
2回答
我试图用MatLab和ode45函数来求解微分方程。问题是如何更改最小步长？我想把它做大一点，避免太小的步子。我试着这样做：
tspan = [3, 4]; (boundaries of time line)
[t, q] = ode45('dqdt2', tspan, q0);
T-时间
Q-解查找
dqdt2 -我的文件有函数
q0 -初始条件
浏览 5提问于2017-05-03得票数 1
回答已采纳
我想用ode45求解100个微分方程。我的整个100 ode方程具有相同的模式，即以xdot= A_x + B_y + C_u的形式表示，其中xdot、A、B和C具有n_1顺序。第一个方程为x点(1)=A(1)*x(1)+ B(1)*y(1) +C(1)*u(1)，第n方程为xdot(n)= A(n)*x(n) + B(n)*y(n) +C(n)*u(n)。在ode的函数脚本文件中编写这100条方程是很费劲的。我的问题是，如何在MATLAB函数脚本中使用for循环或任何其他方法来编写这100个方程并求解ode方程？
浏览 1提问于2022-02-28得票数 0
我仍然在学习Scilab (5.5.2)，所以我正在编写和运行测试代码，以使自己熟悉这个软件。
为了检验数值微分方程求解器，我从方程dy/dx = A开始，其解为y= Ax+c (线方程)。
这就是我写的代码：
// Function y = A*x+1
function ydot=fn(x, A)
ydot=A
endfunction
A=2;
//Initial conditions
x0=0;
y0=A*x0+1;
//Numerical Solution
x=[0:5];
y= ode(y0,x0,x,fn);
//Analytical solution
y2 = A*x+1;
浏览 4提问于2017-06-22得票数 1
回答已采纳
2回答
我正在开始一个项目，由几个偏微分方程的求解器(全部求解同一类型的方程)组成，用于工业目的。应用程序应该有GUI (与计算代码分离)和具有不同求解器的DLL。目的是实现可扩展性，以便创建新的更好的求解器，并将它们作为插件添加到求解器列表中。不同的求解器需要不同的参数集。我应该考虑什么样的设计模式？
浏览 1提问于2017-09-06得票数 0
1回答
我正在尝试使用julia提供的DifferentialEquations.jl，它工作正常，直到我尝试在二阶常微分方程上使用它。例如，考虑二阶常微分方程
x''(t) = x'(t) + 2* x(t)，带初始条件
x'(0) = 0，x(0) =1
它有一个解析解: x(t) =2/3exp(-t)+1/3exp(2t)。
为了从数值上解决这个问题，我运行以下代码：
using DifferentialEquations;
function f_simple(ddu, du, u, p, t)
ddu[1] = du[1] + 2*u[1]
end
浏览 5提问于2021-05-04得票数 3
我用MATLAB的函数'pdepe‘来解一个偏微分方程的问题，一个抛物方程。
我需要知道函数使用的是哪种数值方法，因为我必须在报告中通知这一点。
MathWorks中函数的描述是“在一个空间变量和时间内求解抛物型和椭圆型偏微分方程组的初边值问题”。是有限差分法吗？
谢谢你帮我。
浏览 4提问于2017-04-02得票数 0
回答已采纳
我有个问题
我想知道如何用matlab求出这样的二阶微分方程：
x''(2nd order) = A + B sin(x)
x = ?
A,B are parameters like 1,2
浏览 1提问于2014-06-13得票数 0
回答已采纳
1回答
如果我知道dN/dt=rN(1-(N/K)) (逻辑微分方程)的话，那么摘自萨尔·汗学院的讲座。
如何求解N，并用R绘制N=f(t)？
谢谢
浏览 3提问于2014-07-28得票数 3
回答已采纳
我正在用pdepe解决一个二阶偏微分方程。我的问题的初始条件是一个二阶常微分方程的解，我用ode45解决了它。如何在pdepe中传递ode45解决方案作为我的初始条件？
浏览 1提问于2014-12-21得票数 0
我正在试着解一组微分方程，但我一直很难做到这一点。我的微分方程包含一个表示从1到n的数字的"i“下标。我尝试实现一个forloop，如下所示，但我一直收到这个索引错误(错误消息如下)。我尝试更改初始条件(y0)和其他值，但似乎都不起作用。在这段代码中，我使用了solve_ivp。代码如下：
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
from scipy.integrate import solve_ivp
def testmodel(t, y):
X = y[0]
Y =
浏览 2提问于2021-02-25得票数 0
例如，如果我有下面的代码: y''+y=0
在Matlab中，我可以假设y(1)=y和y(2)=y'，然后
dy(1)=y(2)
dy(2)=-y(1)
现在，如果我有y''+ f(x) *y=0形式的常微分方程，其中f(x)是一个函数，我可以在任何地方计算它的数值，但你没有解析形式，不能将f(X)拟合到某个函数。我的问题是，如果我想把f(x)的值放入常微分方程，那么如何在MATLAB中数值求解这个常微分方程，例如，
f(x)=0.1，当x<0.5和
当x>=0.5时，f(x)=1.0
这首颂歌可能不切实际，但它给了我一个关于我所面临的问题的
浏览 1提问于2015-06-02得票数 0
1回答
我已经得到了基本的SIR模型以下代码。我需要在没有ode()函数或任何其他函数的情况下求解微分方程。你能帮我在不使用ode()函数的情况下计算这个微分方程吗？
library(deSolve)
library(shiny)
ui <- fluidPage(
titlePanel("Choose amounts"),
fluidRow(column(width = 4,
sliderInput(inputId = "time_values", label = "Days", value = 10, min = 1, m
浏览 2提问于2020-02-18得票数 0
如果我有一个用一阶微分方程表示的二阶微分方程，它描述了一个质量-弹簧-阻尼器系统，当我不知道一阶微分方程时，我怎么能用欧拉法来画这个方程呢？我想用MatLab来做这件事。这是一个家庭作业问题，所以我还没有发布任何代码。我只想要一个简单的大纲，告诉你该怎么做。
二阶导数是d2(t+1) = (-1/m)*(c*d1+k*y)，其中c，m，k是常数，y初始是1，d1是从0开始的一阶微分，t是时间。
有什么想法吗？
谢谢:)。
浏览 1提问于2011-01-28得票数 1
回答已采纳
我正在使用scikits.odes.odeint来解决一个相对较大的常微分方程系统(40个变量，100个动力学参数)。我想通过向求解器提供一个雅可比来加速这个过程。但是，我不知道如何为此目的最好地定义雅可比矩阵并将其提供给求解器。
浏览 24提问于2019-07-05得票数 1
回答已采纳
我想解决一个由4个耦合的微分方程组成的系统，其中包含python (sympy)：
eqs = [Eq(cP1(t).diff(t), k1*cE1(t)**3), Eq(cE1(t).diff(t), -k1 * cE1(t)**3 + k6 * cE3(t)**2), Eq(cE2(t).diff(t), -k8 * cE2(t)), Eq(cE3(t).diff(t), k8 * cE2(t) - k6 * cE3(t)**2)]
当我尝试用“dsolve_system”来解决系统问题时：
solution = dsolve_system(eqs, ics={cP1(0): 0, cE1
浏览 7提问于2020-12-21得票数 0
我有一个类似下面的matlab格式的微分方程。
syms x y m g r l J
% x,y are variables, the others are constant
1: 0.5*m*(r^2*x^2+l^2*(Dx-Dy)^2+2*r*l*Dx*(Dx-Dy)*cos(y))+0.5*J*(Dx-Dy)^2=m*g*
(l*sin(x-y)-r*(1-cos(x)));
2: J*(D2x-D2y)+l^2*(D2x-D2y)-r*l*(Dx)^2*sin(y)+r*l*D2x*cos(y)-m*g*l*cos(x-y)=0,
3: x(0)=pi/2,y(0)=pi/2,
浏览 8提问于2016-09-02得票数 0
1回答
时滞方程的标准近似格式？
、、、
我想为时滞微分方程编写一种近似方法，这在某种程度上是一个利基主题，但我想尝试一下。然而，标准的Euler方法或其他Runge-Kutta方法并不一定符合这一点。怎样才能精确地近似求解时滞方程？
浏览 21提问于2020-12-31得票数 1
1回答
我试着把一个二阶微分方程输入到matlab中，求解x=0到x =1，我不知道怎么做。下面是方程式：
y'‘=1+ 0.1 \sqrt{1+(y')^2}
初始条件为零。
浏览 0提问于2011-07-14得票数 2
回答已采纳
这里有点新手。我试图解决一个由6个常微分方程组成的系统，但是有一个条件，如果(x(1)<=0.1)，那么x(6)=x(6)/2.0。我如何将其放入代码中？谢谢!
浏览 0提问于2015-07-04得票数 0
1回答
正如标题所说，我想知道如何从Matlab常微分方程中检索二阶导数。我试图模拟的系统用12个一阶微分方程来描述。
让我们来解释一下：
...
[T,Y] = ode113(@sixdofsolver,time,Y0,options,settings);
T= 1700x1向量
Y= 1700x12矩阵
现在，如果我这样做了：
[dY] = sixdofsolver(T,Y,settings)
dY = 12x1向量
我原本期望有一个和Y一样大小的矩阵。
我做错了什么？
浏览 0提问于2013-10-05得票数 0
1回答
求解两个耦合二阶边值问题
、、、
用模solve_bvp求解了一个具有两个边界条件的单二阶微分方程。然而，现在我正试图求解二阶微分方程组；
U'‘+ a*B’=0
B'‘+ b*U’=0
在边界条件U(+/-0.5) = +/-0.01和B(+/-0.5) =0。我把它分解成一个第一常微分方程组，我试图用solve_bvp数值解它们。但是，我的解决方案只是得到满是零的数组。我相信我执行边界条件是错误的。我不清楚如何处理文档中的两个以上方程。我的尝试如下
import numpy as np
from scipy.integrate import solve_bvp
import matplotlib.pyp
浏览 5提问于2020-02-26得票数 1
回答已采纳
4回答
用C#求解偏微分方程
、、、
我正在做一个项目(C#和.NET框架)，它需要我解决一些偏微分方程。有没有基于.NET框架的特定库可以让我看到，让我的工作变得更简单？
我曾经使用过MATLAb，在那里求解偏微分方程非常简单。我该如何解决这个问题？
浏览 4提问于2009-09-20得票数 5
回答已采纳
我想使用bvpsolve来找到需要最终边界条件的常微分方程组的解。
这两个微分方程是一个动态优化问题的结果，边值来自于相关的横截性条件。
我的代码：
fun <- function(x, y, pars) {
dy1 <- -(d + n) * y[1] + y[1]^a - y[2]
dy2 <- y[2]/b * (a * y[1] ^(a-1) - d - r)
return(list(c(dy1, dy2))) }
# parameter value
n=0.01
a = 0.3
d=0.05
b = 3
r = 0.02
# initial and
浏览 0提问于2021-09-13得票数 0
1回答
你好，亲爱的社区，
我现在正试图解一个二阶微分方程，时间很长了。该公式如下：
∂P/∂t =D*∂^2P/∂x^2
我应该运行几个P；数据帧的P1= df$col1，P2=df$col2，P3=df$col3和D是df$col1 4，x是df$col5。
我正在尝试使用"deSolve“来实现这一点，但目前我并没有取得多大的成功。对于如何用R来求解二阶微分方程，我将不胜感激。
浏览 1提问于2021-02-03得票数 1
回答已采纳
在Julia中，我正在积分一个向量微分方程(分量的数目是数量级为50)。当我绘制这个解决方案时，如果我绘制所有组件，显然这个图会非常拥挤，所以我只想绘制前10个组件。我找不到该怎么做。另外，我不想为每个组件绘制输出，我想绘制的不是每个时间步骤的积分输出，而是这个输出的双曲切线。我也找不到如何做这个操作。我认为关键的问题是，我不明白/找到一个微分方程的解在julia中是什么样的对象。
下面是一个极小的工作示例，它集成了微分方程并绘制了整个解。
using DifferentialEquations
using Plots
using LinearAlgebra
N=50
浏览 0提问于2019-01-08得票数 0
回答已采纳
2回答
处理约束的PDE求解器
、、、、
我正在试着解一个一般形式的偏微分方程组。
F(f(x,y), f'(x,y), f''(x,y), g(x,y), g'(x,y), g''(x,y)) = 0
其中关于x和y的导数都可以取，并且f(x,y)和g(x,y)受某些约束
G(f(x,y),g(x,y)) = 0
我想知道是否有任何(最好是基于Python的)求解器(不是方法，因为我知道方法)可以处理这类问题？如果我的问题看起来很一般，我将非常感谢您的帮助和道歉。
浏览 0提问于2015-07-13得票数 1
1回答
我正在尝试用数值方法来解决下面的常微分方程 ? 使用边界条件 y(x=0)=0;y(x=20)=1 其中f(x)是我已经数值求解的另一个常微分方程的解。我正在尝试用https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.integrate.solve_bvp.html来解决我的ODE，但是我遇到的问题是f(x)被写成一个列表 coordinates x = np.arange(0,20.01,0.05) 并且我不知道如何在scipy.integrate.solve_bvp中插入列表来获得常微分方程的解。 提
浏览 11提问于2019-01-19得票数 0
回答已采纳
2回答
我是Python初学者。我正在尝试切换matlab中的一些程序。我需要解一个刚性常微分方程，它的输入都是矩阵。在matlab中，我使用
[ttT,uT] = ode23s('SST',t,fT);
浏览 0提问于2012-01-05得票数 4
回答已采纳
1回答
我想知道，对于一个微分方程，ODEINT和solve_ivp之间有什么区别？它们之间有什么优点和缺点呢？
f1 = solve_ivp(f, [0,1], y0) #y0 is the initial point
f2 = odeint(f, y0, [0, 1], args=(a, b)) # a and b are arguments of function f
谢谢
浏览 2提问于2019-03-13得票数 5
回答已采纳
1回答
我正在尝试用Julia数值求解一个抛物型偏微分方程，但我找不到任何可访问的文档来帮助我。
下面是一个例子: t，x是一维实数。我想求解u(t，x)=u1(t，x) u2(t，x)；u满足偏微分方程
du1/dt = d^2u1/dx^2 + a11(x，u) du1/dx + a12(x，u) du2/dx + c1(x，u)
du2/dt = d^2u2/dx^2 + a21(x，u) du1/dx + a22(x，u) du2/dx + c2(x，u)
在Julia中这样做是可能的吗？这类问题可以在Matlab中使用pdepe来解决。
浏览 0提问于2018-01-18得票数 7
回答已采纳
2回答
我现在有一条振荡曲线，它是一组非线性常微分方程的解的一部分。随着时间的推移，我需要测试这条曲线的稳定性/收敛性。如何在Matlab中实现？
该图如下所示：
浏览 0提问于2013-04-26得票数 0
回答已采纳
1回答
Julia是否支持边值微分代数方程？我有一个隐式的常微分方程，它的可变质量矩阵有时是奇异的，所以我必须使用DAEProblem。我的问题是x1(t)和x2(t)的两个耦合的二阶常微分方程，通过设置x1'(t) = y1(t)和x2'(t)=y2(t)，我已经将它们转换为四个一阶方程。我在域的开头和结尾有x1和x2的值，但没有y1或y2 anywhere的值，因此我需要DAE和BVP值。 This github的帖子建议这是可能的，但我担心我对机制的理解还不够深入，无法理解如何将DAEProblem与BVProblem结合起来。 我已经成功地按照数字食谱编写了多个拍摄代码来解决
浏览 30提问于2021-07-08得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
