用IDL求解方程
IDL是一种交互式数据语言（Interactive Data Language），用于科学数据分析和可视化。它提供了丰富的数学函数和工具，可以用于求解方程。
求解方程是数学中的一个重要问题，通常涉及找到使方程成立的未知数的值。IDL提供了多种方法来求解方程，包括数值方法和符号方法。
数值方法是通过迭代计算来逼近方程的解。其中一种常见的数值方法是牛顿法（Newton's method），它通过不断迭代逼近方程的根。IDL提供了用于实现牛顿法和其他数值方法的函数，如ROOTS和FSOLVE。
符号方法是通过代数运算来求解方程。IDL提供了符号计算的功能，可以对方程进行代数化简、求解和化简。其中一种常见的符号计算方法是多项式求解（Polynomial solving），它可以求解多项式方程的根。IDL提供了用于多项式求解的函数，如POLYROOT。
方程求解在科学、工程和计算领域有广泛的应用。例如，在物理学中，方程求解可以用于计算物体的运动轨迹；在工程学中，方程求解可以用于优化问题的求解；在计算机图形学中，方程求解可以用于计算曲线和曲面的交点。
对于方程求解，腾讯云提供了多种相关产品和服务。例如，腾讯云的人工智能平台（AI）提供了强大的数值计算和符号计算功能，可以用于方程求解。此外，腾讯云的大数据平台（Big Data）和云计算平台（Cloud Computing）也提供了高性能计算和分布式计算的能力，可以加速方程求解的过程。
更多关于腾讯云相关产品和服务的信息，请访问腾讯云官方网站：https://cloud.tencent.com/
相关·内容
2回答
我正在尝试对一些数据拟合一个圆圈。这需要数值求解一组三个非线性联立方程组(参见的完全最小二乘法)。
在我看来，IDL提供的NEWTON函数很适合解决这个问题。NEWTON需要一个函数的名称，该函数将为自变量的特定值计算方程系统的值：
FUNCTION newtfunction,X
RETURN, [Some function of X, Some other function of X]
END
虽然这很好用，但它要求方程系统的所有参数(在本例中是数据点的集合)都是在newtfunction中硬编码的。如果只有一个数据集需要求解，这是很好的，但是我有数千个数据集，并且不能手动为每个数据
浏览 2提问于2014-10-26得票数 0
我有，我想解决X(t)，Y(t)，Z(t)，希望用我熟悉的Octave/Matlab，但我不介意用任何其他必要的方法来解决它。
现在，Fsolve对于正则方程组是有用的，而Ode45，Lsode对于微分方程是有用的。但是，这个特殊的系统呢？请注意，底部的微分方程不仅包含Y，还包含X和Z，它们依赖于上面的两个非微分方程。
老实说，我不太确定如何使用基本的代码来解决这个系统，在花了一些时间思考之后，我决定寻求帮助。如果可能的话，我真的很感谢任何能有效地解决这个问题的指导。现在几乎所有的回复对我来说都是有用的。
浏览 9提问于2020-02-15得票数 0
回答已采纳
我知道如何用numpy来求解基本的线性矩阵方程。
但是，我有一个矩阵A，方程A^2 + xA + yI = 0，其中x和y不是向量，而是标量。I是单位矩阵，0是维数匹配A的零矩阵。
对于小型矩阵来说，这是一个非常简单的方法(当然，假设存在一个解决方案)，但是我正在练习一个编码面试，并希望用python来解决类似的问题。也许给出的矩阵很大..。
下面是一个样本矩阵A，它得到解x=-2，y=1：
np.array([[1,1,0],
[0,1,0]
[0,0,1]]
在纸上，这与解线性方程组x= -2和x+y=-1一样容易。我所面临的问题是将上述形式的方程解
浏览 4提问于2020-02-22得票数 1
回答已采纳
我想知道并试图搜索它，但我不知道如何谷歌它，如何解决线性方程组，其中方程是平等的，彼此。我能够为两个方程编写一个具体的算法，但不能对N个方程进行动态计算。我将给出三个例子(我的方程大致是怎样的)：
C1，C2，C3，X -都不知道，但我不需要知道X的结果。
可以这样解释(最后一个方程C1 + C2 + C3 =1不在这里集成)：
请不要试图解决这个问题，我不确定这些方程是否有结果。我只是随机输入系数。但这就是我的方程式的样子。只有不同的系数。最好是在蟒蛇，但如果我知道如何用矩阵做数学的话，我会很高兴。我试着用两个未知数计算它，最后我得到了二次方程，所以在有三个未知数的情况下
浏览 10提问于2022-03-24得票数 1
回答已采纳
在每一个递归循环中，我都要求解两个方程：
X=(B)*Y* inv(B)，X=X+ A‘* inv(B) * A，
我用这样的方式解决问题：
C= inv(B)_Y <=> BC = Y求解C. D= C_inv(B) <=> DB =C <=> B‘= C’
E= inv(B)*A <=> BE = A，求解E.
随着时间的推移，所有的矩阵都会发生变化，所以我不得不再做一次(或反转)每一次递归。N通常在1-10左右，可能更多，但通常是这样的。B是正定的，所以我可以用cholesky进行分解，然后求解多个右手边的方程。
这比把B反转然后用它做矩阵
浏览 7提问于2009-05-31得票数 0
回答已采纳
2回答
解/解性能差
、、、、
我用渐近解一个简单的线性方程组。
这是一个耦合的ODE，有变量的时间导数，我需要解最高导数的方程组。由于渐近不允许我解像phi_1.diff(t)这样的语句，所以我用占位符符号替换了所有的导数。
例如：
phi.diff(t).diff(t) + phi(t) =0
变成了
丁菲+ phi(t) =0
这个很好用。解决方案是正确的，我可以模拟这个系统--这是一个钟摆：
问题是，求解方程组(用linsolve)需要很长时间。
对于两个方程式，它需要2秒。对于3个方程，它仍在计算(超过10分钟后)。
编辑： @asmeurer建议我尝试解决问题。对于n=3，linsolve大约花了
浏览 5提问于2016-06-21得票数 4
回答已采纳
1回答
我试图用numpy来求解如下形式的非线性方程组：
a(y-2.7)(1-exp(-a*z)) = (x-2.7)(1-exp(-z))
b(w-2.7)(1-exp(-b*z)) = (x-2.7)(1-exp(-z))
c(w-2.7)(1-exp(-b*z)) = (y-2.7)(1-exp(-a*z)
d([y+w]/2-2.7)(1-exp(-d*z)) = (x-2.7)(1-exp(-z))
很明显，系统中的方程和未知数一样多。a，b，c，d是上述系统的常量。这是最简单的系统，在其他情况下会有更多的方程。这些方程的解具有相似的数量级，因此我知道Levenberg-Marquar
浏览 4提问于2014-06-09得票数 0
回答已采纳
我必须使用matlab分析802.11的饱和吞吐量，这是我的问题。我试图用求解函数来求解下面的参数方程(参数为m，W，a)，并得到
警告:找不到明确的解决方案
如何用matlab求解上述方程？
浏览 11提问于2013-01-15得票数 0
当我试图用Pulp在Jupyter Notebook上绘制线性规划方程时，我得到了这个错误。
AttributeError: 'NoneType' object has no attribute 'actualSolve'
上下文：
浏览 0提问于2020-07-10得票数 0
我试着用RSolve来求解数学中的两个变量递归方程，但它只是重复了我输入的内容。
能否用Mathematica求解两个变量递推方程
浏览 1提问于2011-06-09得票数 3
下面是扩散方程的一个pde。因为D= 0，它的第二阶被消除了。我是matlab的新手，我知道扩散方程与热方程有一定的相似性，但我不知道如何将这种方法应用到我的解中。
请给我你的建议。用我的真心感谢你。
另外，如果D是一个张量，而D11 neq D12 D21 neq D22，那么如何解决我的问题呢？
a = epsilon = 0.1，T0 = 5.0，b= 0.1，c= 100
扩散方程：
F(x1，x2)和F2(x1，x2)：
浏览 5提问于2011-05-14得票数 0
我经常遇到下面显示的初始化问题，上面写着some kind of non-linear solver will attempt to handle this problem。据我所知，Dymola在得到原始DAE系统的BLT形式后，用牛顿法求解方程，那么这个非线性求解器到底是什么？
浏览 2提问于2021-02-12得票数 3
回答已采纳
这是方程式。
在MATLAB中，我编写了以下代码：
mu = 305；%作为临时示例，因为它大于c，即300
syms x
eqn = ((1 + (x/(mu + 300)))^((1/2) + (150/mu)))*((1 - (x/(mu - 300)))^((1/2) - (150/mu))) - 0.2 == 0 %matlab shows the answer is -605
solve(eqn,x)
这是相同的方程式，除了我用MATLAB代替mu.MATLAB返回以下内容：
eqn = logical 0
ans = x
不知何故，我输入的方程式有误吗？这就是为什么当我不抑制方程
浏览 2提问于2017-05-02得票数 0
2回答
我想在Modelica中做以下事情:在一个特定的模型中，我有几个参数h，我想通过求解参数和其他值之间的一组隐式方程，从它们中推导出一些与时间无关的值k。由于方程本质上是隐式的，我不能简单地赋值表达式，而是必须让求解器找到解决方案。 由于参数是恒定的，我希望在系统其余部分的实际时间积分(例如，包含k作为系数的微分方程)发生之前，能够在一开始只解一次方程。 请参阅下面的示例，了解我希望能够编写的内容： model ConstantTest
parameter Real h = 2;
const Real k;
initial equation
k*k=h; // just an ex
浏览 39提问于2019-02-08得票数 5
回答已采纳
1回答
我试图画出微分方程dh/dt = I(t) - f(h)*h(t)，其中f(h)可以用f(h) =1-e^(-kh)建模，k是常数。我该怎么做呢？
浏览 4提问于2022-05-03得票数 -1
1回答
我必须写一个程序，用四种基本运算来求解方程。该程序必须如下所示：
示例输入：
SET a = 1
SET b = 3
SOLVEFOR c ( 4 * ( b - c + 2) + a ) = ( 2 - a )
输出：
c = 5.000000
我使用反向波兰符号来计算方程右侧的值(我知道那里只有已知的数字)。至于左边，那是我的问题。我知道必须有一种修改RPN来获取该变量的方法，但我不知道该如何做。有什么建议吗？
编辑:我知道这可以用二叉树来完成，但我正在尝试使用RPN。
浏览 2提问于2014-04-10得票数 1
我正在用pdepe解决一个二阶偏微分方程。我的问题的初始条件是一个二阶常微分方程的解，我用ode45解决了它。如何在pdepe中传递ode45解决方案作为我的初始条件？
浏览 1提问于2014-12-21得票数 0
6回答
我想问一下关于高效算法的问题。例如，我有一个等式：
x = y + z;
如果变量y= 1，变量z=2，那么变量x是3
但是，如果变量x=3，变量z=2，那么如何自动获得y的值呢？不创建新方程式y = x - z
我希望可以使用C#或javascript获得示例代码。
另一个例子，如果方程式是
a = (((x + y - z)/2)*10)^4
方程来自程序，用户提交3个变量的值。用户提交变量(x，y，z)或(y，z，a)或(z，a，x)或(a，x，y)
如果用户输入变量x、y和z值，则程序可以用该方程显示a的值。如果不创建a= ...
如果用户输入变量y、z和a值，则程序可以用方程显示x的
浏览 3提问于2016-03-09得票数 0
对于GAMS中的三个星号(*)，我有一些疑问，它们可能会显示在单个方程列表的末尾。我知道这是一个警告，即限制在起点是不可行的。我有一个模型，用GAMS求解后，模型状态和求解者状态为‘1’，方程似乎是可行的，但在方程的末尾显示了三个星号。
我想知道：
(一)起点是什么？
2)模型不可行吗？
我真的很感激你的帮助。
浏览 1提问于2016-10-05得票数 0
我想知道是否有任何方法可以用matlab的simulink模块得到非线性方程()的最佳逼近(dx/dt是指与t相关的x导数)。
浏览 3提问于2015-11-17得票数 0
我用MATLAB的函数'pdepe‘来解一个偏微分方程的问题，一个抛物方程。
我需要知道函数使用的是哪种数值方法，因为我必须在报告中通知这一点。
MathWorks中函数的描述是“在一个空间变量和时间内求解抛物型和椭圆型偏微分方程组的初边值问题”。是有限差分法吗？
谢谢你帮我。
浏览 4提问于2017-04-02得票数 0
回答已采纳
本文试图用牛顿-费舍尔( Newton )的(线性互补问题)(线性互补问题)的牛顿-费舍尔公式，对牛顿方程组的不同迭代子解的结果进行分析。到目前为止，我已经实现了精确的求解器-Gauss，并且用bicg 作为J*h = -p方程的子求解器(其中J是jacobian方程，p是Fischer函数的值，h是我的实现步骤)。
在代码的一部分，我实现了bicg和精确的求解器：
if itt
% use iterative solver for newton eq
while ~all(fischer(x, A*x+b) == 0) & its < max_it
浏览 4提问于2012-03-08得票数 1
回答已采纳
这是我想要解决的方程组的一个例子。
必须选择a和b，使这些积分的结果为0。但我不知道
如何与参数集成
如果我不能像Ax=b那样写出方程组，该如何求解它。如果是这样，我可以用x= A\b来求解它。
非常感谢。
浏览 6提问于2022-02-28得票数 0
首先，我不太清楚如何解释这个问题，因此不知道如何寻找。其内容如下：
我有很多方程。该方程可以求解任意给定变量，给定其他未知数，也可以同时求解多个未知数。例如，给定a和b，方程f(x, y) = a和g(x, y) = b，我们可以同时求解得到x和y。
我需要一种算法，以已知的值和方程，并返回的顺序，解决它们将得到期望的值。
例子方程：
f(a，b) = 0f(b，c) = 0
找c给定一个->用eq1找出b给a，然后用eq2找出c给定b
示例2：
f(x，y，a) = 0f(x，y，b) = 0
求x给定的a，b ->解x和y同时使用eq1和eq2
我尝试用一个图来简化这个问题，其
浏览 5提问于2022-04-21得票数 1
我在解方程组时遇到了困难。我有三个方程，每个方程都有一个已知解和三个未知数。但是，当我在MATLAB中使用solve函数时，它返回错误:我有六个方程和三个变量。
下面是我的代码片段：
syms V0 T0 X0
A=(g*X0/(2*V0^2*cos(T0)^2)-tan(T0))==a;
B=(tan(T0)-g*X0/(V0^2*cos(T0)^2))==b;
C=(-g/(2*V0^2*cos(T0)^2))==c;
soln=solve([A,B,C],[V0,T0,X0]);
我已经计算了a、b和c的标量值。g是一个常量。我不知道为什么返回我有六个方程式。
浏览 4提问于2016-01-28得票数 0
1回答
方程的确定性和随机性部分
、、、、
我正在寻找一种既能解确定性方程又能解随机方程的数值方法。在确定性的情况下，我知道四阶RK方法是一个有价值的方法，非常有效。不幸的是，还没有成功地应用于随机方程(至少据我所知)。
现在我想知道的是，是否存在一种数值方法，可以求解这两个方程(大致上，我的意思是，与解析解相比)，在这种情况下，会是什么。例如，一个解析可解的随机方程就是Black-Scholes方程。
浏览 3提问于2018-10-31得票数 0
1回答
我使用下面的代码来求解一个三次方程。
from sympy.solvers import solve
from sympy import Symbol
x = Symbol('x')
print(solve(-0.0643820896731566*x**3 + 0.334816369385245*x**2 + 1.08104426781115*x - 2.05750838005246,x))
由于它是一个具有实数系数的三次方程，所以不可能有三个不同的复根。但它给出了以下结果。
[-3.19296319480108 - 0.e-22*I, 1.43925417946882 + 0
浏览 2提问于2019-04-12得票数 4
回答已采纳
1回答
用枕氏方程求解粘度的正则化
、、、、
为了简单起见，考虑了，下面的方程是 (Burgers方程)：
让我们用一个可变时间步长的求解器(在我的例子中是scipy.integrate.ode.set_integrator("zvode", ..).integrate(T))来解决这个问题。
问题如下:如果我们在傅里叶空间中使用天真的实现，则为。
当时间步长过大时，粘度项nu * d2x(u[t])会引起超调。这可能导致解决方案中产生相当大的噪声，甚至导致(假的)发散解(即使使用刚性求解器，在这个方程的稍微复杂的版本上)。
规范此的一种方法是在步骤t+dt处求出粘度项，更新步骤变成
此解决方
浏览 10提问于2020-04-21得票数 1
1回答
我会尝试描述我的问题，但最后有一些具体的问题。
我有一些困难，解决一个涡轮机问题，在Matlab。目的是用已知的物理测量和运行条件来计算涡轮的效率。为了做到这一点，我使用了科学论文中给出的方程式：
在我的例子中的问题是，要知道效率，我需要出口气体参数(温度，压力，速度…)这又取决于效率。我主要在一个名为EES (工程方程求解器)的软件中解决了这个问题，在这个软件中，方程可以按任何顺序写入，程序对它们进行分组并自动求解。由于程序中的限制，我切换到了Matlab，这不允许我扩展这个问题。
在matlab中，我首先以一个可解的顺序将所有的方程写成一个函数，除了求解问题所需的两个退出参数之外。将两个退
浏览 1提问于2012-05-31得票数 0
我是python的新手。我有一组自治方程，试图用相空间分析来分析渐近行为。
f(x,y)=a*x*y((y**2)+a+c)
g(x,y)=a+(y**3)+((y**3)+(y**2)(x+a))
其中x和y是变量？我寻求帮助寻找临界点，雅可比和特征值，也得到相空间图？
浏览 10提问于2019-02-19得票数 0
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
回答已采纳
5回答
DLL之间的交互
、、
我有一个DLL，比如说A，它实现了求解数学方程的通用算法。我有另一个DLL，比方说B，它实现了数学方程DLL B使用DLL A来求解它的方程，但DLL B必须能够执行DLL A中的给定方程(DLL A实现一些数值近似方法，并在给定方程上“尝试”不同的值，以达到期望的数学解)现在，DLL A必须“知道”DLL B，反之亦然。
这是一个“糟糕”的设计吗？我可以把这样的代码想象成业余的cose吗？你会怎么做(请记住，DLL A实现了通用算法，这些算法被其他实现数学方程的不同DLL所使用)？
谢谢大卫
浏览 0提问于2009-07-19得票数 1
2回答
我正在开始一个项目，由几个偏微分方程的求解器(全部求解同一类型的方程)组成，用于工业目的。应用程序应该有GUI (与计算代码分离)和具有不同求解器的DLL。目的是实现可扩展性，以便创建新的更好的求解器，并将它们作为插件添加到求解器列表中。不同的求解器需要不同的参数集。我应该考虑什么样的设计模式？
浏览 1提问于2017-09-06得票数 0
我想用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
2回答
我的问题涉及Scip中的混合整数规划(MIP)：
我有以下代码：
$\min trace(X)$
subject to
$$(A+D)^TX+X(A+D)=I\\
d_i \in \left\{0,1\right\} \mbox{ for } i=1,\ldots,n$$
where A is a n*n matrix and $D=diag(d_1,\ldots,d_n)$ is a diagonal matrix.
由于矩阵约束是线性的，所以方程可以转化为线性方程组(通过Kronecker乘积和向量化操作)，但这限于小n。是否可以用Scip直接求解矩阵方程？有没
浏览 3提问于2014-03-27得票数 1
我想优化几个方程之间的共同参数，但我不知道如何同时拟合它们。
问题实质上是这样的，有四个方程要解，三个参数要优化：
a+b+c+1750=T
12=a/T*100
15=b/T*100
37=c/T*100
我想找出a，b，c的最优值，有人建议用最小二乘法吗？只有当只有一个方程需要解时，我才会熟悉。
浏览 1提问于2018-06-06得票数 0
回答已采纳
1回答
我有一个非线性方程:我在下图中画出了这个方程的左边和右边。显然，有两种解决方案。我已经尝试了scipy.optimize.fsolve，并战略性地猜测了起点，以找到两个解决方案，但这是相当不健壮的。有没有一个Python求解器，可以找到这个方程的所有解，例如，在某个有限的x值集合中？
浏览 5提问于2021-11-10得票数 0
1回答
目标优化
、、
我试图优化一个函数，使其尽可能接近于零。
其职能是：
def goal_seek_func(x: float) -> float:
lcos_list_temp = [energy_output[i] * x for i in range(life)]
npv_lcos_temp = npv(cost_capital, lcos_list_temp)
total = sum([cost_energy_capacity,
cost_power_conversion,
balance_of_pla
浏览 2提问于2022-04-10得票数 0
回答已采纳
1回答
我正在尝试将一个IDL程序翻译成Python。我必须通过以下方式解决来自SVD的结果
from scipy.linalg import svd
A = [[1,2,3],[4,5,6]]
b = [4,4,5]
u,w,v = svd(A)
这段代码运行得很好，并且可以很好地从IDL翻译过来。下一步是IDL(!)
x = svsol(u,w,v,b)
python和IDL中的u几乎相同(对于其他矩阵也是如此)。唯一的区别是维数，IDL的矩阵更大，但有很多零。从这个意义上说，Python的矩阵看起来更加压缩了。
有没有人知道类似的Python。
如果任何人需要它，这里是的手册。
浏览 0提问于2012-09-25得票数 3
回答已采纳
1回答
我有这样的方程式：
wl[n]= w[n] + Δx*v[n]
ϕl[n]= ϕ[n] + Δx*ρ[n]
ρl[n] = ρ[n] - Δt*fρ(ρ,v,w,n)
vl[n] = v[n] - Δt*fv(ρ,v,w,Δx,n)
模拟这些方程：
我可以使用ODE解算器在时间上，然后在空间中进行解算吗？或者，有没有Julia软件包可以帮助我使这些方程无条件稳定？
浏览 6提问于2021-11-19得票数 1
1回答
我试着用渐近来解决下面的三角方程，但是它给了我这个错误：
NotImplementedError: could not solve b*tan(3*atan(6/b)/2) - 9
有没有办法让我用渐近解出下面的方程？
import sympy
from sympy import *
a = sympy.symbols("a")
b = sympy.symbols("b")
eq1 = sympy.Eq(b*tan(2*a), 6)
eq2 = sympy.Eq(b*tan(3*a), 9)
result = sympy.solve([eq1, eq2]
浏览 47提问于2021-02-19得票数 1
回答已采纳
2回答
我正在寻找一个第三方自定义控件，可以处理编辑和解决简单的数学方程(与图像生成)。
我了解LaTeX并使用它的语法来生成数学方程式的图像(我在这里使用了这个实现LaTeX：的项目)，并且我找到了几个可以将方程式作为输入字符串、编译、解析和求解方程式的项目(这里有一个：)。
但是，我还没有找到任何方法，允许您使用输出的预览图像构建方程(就像LaTeX一样)，然后允许您为方程中的变量指定值，并使用输出的结果求解方程。
有没有人知道是否存在可以编辑和求解数学方程的第三方控件？我使用的是C# .NET 4.0，所以希望我可以使用.NET中的库/控件。
任何帮助或链接到类似的东西将非常感谢。
谢谢，马克
浏览 2提问于2011-04-13得票数 9
4回答
可能重复：
我需要用C#和VB编程求解线性方程组。
下面是方程的一个例子：
12.40 = a * 56.0 + b * 27.0 + tx
-53.39 = a * 12.0 + b * 59.0 + tx
14.94 = a * 53.0 + b * 41.0 + tx
我想得到a，b和tx的最佳近似。
我应该使用某种矩阵类吗？
浏览 7提问于2010-12-08得票数 6
回答已采纳
2回答
是否有一种标准方法来确定指定SAT问题的复杂性？我正在研究代数密码分析，并开始用CryptoMiniSat求解多元二次方程系统。然而，评估一个完整的密码方程系统的复杂性是很好的。它是否只取决于多项式和变量的数目？是否有一种方法来估计解方程系统所需的时间，以及我需要知道哪些参数才能找到它呢？
浏览 0提问于2012-05-13得票数 5
回答已采纳
大家好，Gekko开发人员，
我想为我的热能存储系统建立以下方程模型，
它们既是线性的，又是微分的。一些变量也是时间相关的，作为外部参数出现。从Gekko的例子中，看起来我只能添加这种格式的方程，比如mdl.Equation(T_a.dt() == mt*c_p_w*(T_b-T_a)) / d_w*c_p_w，但是对于不涉及微分的方程，我想添加为普通方程，也可以添加mdl.Equation(m[t] == mc - m[t])。
附注:我如何添加这两种样式，以及哪种模式和求解器将允许这样做。？
谢谢,
诚挚的问候,
拉胡尔
浏览 4提问于2020-11-18得票数 3
目前使用渐近来求解以下方程组：
from sympy import symbols, Eq, nsolve, exp, log, cosh
eq1 = Eq(Vmax**2/Amax * log(cosh(Amax*4.82/Vmax)),40)
eq2 = Eq(Vmax**2/Amax * log(cosh(Amax*2.8/Vmax)),20)
solution = nsolve([eq1, eq2], [Amax, Vmax], [1, 1])
我使用这个技术得到了一个解决方案。然而，我有第三个方程，我想用所有3个方程来求解这2个未知数。在Sympy中可以用3个方程解决2个未知数吗
浏览 183提问于2021-03-02得票数 2
我试图用y'' + ay' + by + c = 0 (二阶微分方程)的形式用odeint来求解python中的方程。据我理解，odeint只适用于初始条件为y(0) = y1，y'(0) = y2。而我的条件是: y'( 0 ) = 0，y'(pi/4) =0。在这种情况下，有什么方法可以使用odeint吗？下面是我的代码：
from scipy.integrate import odeint
import numpy as np
def lagrange(x, teta):
y = x[0]
dy = x[1]
xdot =
浏览 1提问于2018-02-26得票数 0
回答已采纳
1回答
当我用交感解数学问题时遇到了一些问题，我的代码是这样的，
eq1=Eq(0, -sympy.sqrt(x**2) + sympy.sqrt((x + 4)**2))
sympy.solve(eq1,x)
这个方程的答案应该是-2，但我得到了一个空的列表。
浏览 0提问于2018-09-13得票数 2
回答已采纳
1回答
我正在FORTRAN中用SOR方法求解异质结构(AlGaN系统)的泊松方程.对于特定的初始化，求解器提供NAN作为输出和停止(正如我设置了ffpe-陷阱标志)，对于不同的初始化，求解器运行良好。
方程的解是否取决于初始势的选择？
如何理解一般迭代法中NAN产生的原因？
浏览 1提问于2020-08-05得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
