在Python中解决由6个非线性方程组成的系统的问题
可以使用数值求解方法来逼近方程的解。以下是一个可能的解决方案：
- 数值求解方法：可以使用牛顿法或改进的牛顿法来解决非线性方程组的问题。这些方法基于迭代过程，通过不断逼近方程组的解。
- 牛顿法：牛顿法通过使用方程组的雅可比矩阵来进行迭代。它在每一步都通过线性近似来求解非线性方程组。牛顿法的优点是收敛速度较快。
- 改进的牛顿法：改进的牛顿法通过添加阻尼系数来防止迭代过程中出现发散的情况。阻尼牛顿法在收敛性方面更加稳定，但可能会降低迭代速度。
- Python库：SciPy是一个常用的科学计算库，其中包含了用于数值求解的方法和工具。可以使用其中的牛顿法或其它相关方法来解决非线性方程组的问题。
- 应用场景：非线性方程组的问题在科学计算、物理学、工程学等领域中广泛存在。例如，在电路设计、化学反应动力学、机械结构分析等方面都会遇到非线性方程组的求解问题。
- 腾讯云相关产品和产品介绍链接地址：由于要求不提及具体的云计算品牌商，这里无法给出腾讯云相关产品和产品介绍链接地址。
总结：在Python中解决由6个非线性方程组成的系统的问题，可以使用数值求解方法，如牛顿法或改进的牛顿法。通过使用SciPy等库来实现求解过程，并且可以将这种方法应用于各种应用场景中。
相关·内容
我想解决一个由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
1回答
具有多重方程的枕方程组
、、、
你好，当我有一个方程组时，我在使用fsolve函数时遇到了问题。
我的代码只让我有和我的系统中的方程一样多的变量，但是我们都知道，在实践中，你可以拥有比变量更多的方程，特别是当解不是一个具体的数字而是一系列的数字时，一个附加的方程可以帮助缩小“干草堆”的范围。
如何插入比变量更多的方程是我的问题？
假设我有以下非线性系统：
A/B=0.4583(3)
A/C=0.25
A/D=0.72(2)
B/C=0.54(54)
因此，我有以下代码：
from scipy.optimize import *
from numpy import *
def FUNC(arg):
A,B,C,D
浏览 0提问于2018-11-19得票数 1
回答已采纳
在Dymola中，我经常会遇到一个非线性系统初始化失败，或者是一个在大型热流体系统中很难解决的刚性系统，但是对于一个简单的系统来说，不会出现这样的问题。我的问题是：
，，所以我想知道，求解非线性系统模型的最大能力是多少？例如，我最多可以在我的模型中包含多少个非线性方程？，在Dymola中是否有允许提高求解非线性系统能力的设置？，如何在不损害模型精度的情况下减少模型中的非线性方程的数量？
浏览 3提问于2020-04-24得票数 2
回答已采纳
1回答
我有一个非线性系统，由一阶非线性差分方程组成：
x(k) = f(x(k-1)) + u(k-1)
y(k) = h(x(k)) + v(k)
在u(k-1)和v(k)是独立的情况下，零均值高斯噪声过程分别具有协方差Q和R .
我想使用matlab中的kalman函数来估计给定y的x。但是，matlab控制系统工具箱中的接口是：
[K,L,P]= kalman(sys, Q, R, N)
其中sys是状态空间模型。我的问题是:如何在matlab中为我的一阶非线性差分方程系统定义sys？
浏览 4提问于2016-08-16得票数 2
回答已采纳
我正在尝试建立一个丝杠传动的理论模型，将其作为力元素参与mbs模拟。这意味着
F=blackbox(u)
在上图中，我将系统简化为一个有代表性的模型。phi_1(t)和u_1(t)是螺母的给定运动。应计算引线的运动和返回的反作用力。它们是phi_2(t)和u_2(t)。
我写下了静态情况下的所有方程。但是我不确定我应该如何选择方程的运动来数值求解它们。
在运动学中，您可以看到有一个内部自由度。你可以选择它作为u_2，所以根据给定的phi_1和u_1，phi_2有一种可能的位置。
现在你可以根据u_2写出轴向力F和扭矩M，这样你就可以计算出影响螺纹的切向力和法向力。
在这一点上我遇到了
浏览 0提问于2013-08-23得票数 0
在Dymola建立了一个具有朗肯循环的太阳能电站模型，虽然初始化工作正常，但是在2840s的时候，由于非线性方程组的求解失败，导致仿真停止。如下面的截图所示，Dymola建议提供更好的起始值，但这没有任何意义，因为失败发生在模拟期间，而不是初始化过程中。有了更好的起始值，情况根本不会改善。我的问题是:我应该如何处理在模拟过程中无法求解非线性方程？ ?
浏览 71提问于2020-07-06得票数 1
我正在尝试实现apache commons数学优化包。我正在寻找一个"hello world“的例子来解决一个有约束的两个方程的非线性系统。我的系统是可微的，由两个独立变量组成。apache commons文档和彻底的google搜索没有给我提供任何可供我构建的代码示例--有什么建议吗？
浏览 0提问于2011-11-17得票数 3
我想知道并试图搜索它，但我不知道如何谷歌它，如何解决线性方程组，其中方程是平等的，彼此。我能够为两个方程编写一个具体的算法，但不能对N个方程进行动态计算。我将给出三个例子(我的方程大致是怎样的)：
C1，C2，C3，X -都不知道，但我不需要知道X的结果。
可以这样解释(最后一个方程C1 + C2 + C3 =1不在这里集成)：
请不要试图解决这个问题，我不确定这些方程是否有结果。我只是随机输入系数。但这就是我的方程式的样子。只有不同的系数。最好是在蟒蛇，但如果我知道如何用矩阵做数学的话，我会很高兴。我试着用两个未知数计算它，最后我得到了二次方程，所以在有三个未知数的情况下
浏览 10提问于2022-03-24得票数 1
回答已采纳
2回答
我正试图求解python中的非线性方程组。我试过使用Sympy的求解器，但它似乎不适用于for循环语句。我要为解决一个输入范围内的变量x问题。
我已经把我的代码附在下面
import numpy as np
import matplotlib.pyplot as plt
from sympy import *
f_curve_coefficients = [-7.14285714e-02, 1.96333333e+01, 6.85130952e+03]
S = [0.2122, 0, 0]
a2 = f_curve_coefficients[0]
a1 = f_curve_coeffici
浏览 5提问于2020-09-28得票数 0
回答已采纳
3回答
我有个大问题。我需要用一个C++函数或类来解决一个由3个变量的3个方程组成的非线性系统。我想过用牛顿-拉夫森方法来求解。不幸的是，我没有找到可以为我做这件事的源代码。会有人知道这样的程序吗？我快要决定自己构建它了。谢谢
浏览 0提问于2013-06-01得票数 1
回答已采纳
2回答
我有一道数学题，是用Maple解决的，但它必须用C++解决。问题是我找不到一个方程的数值方法(它是用内置的maple "solve“过程求解的)。我可以以某种方式将这个过程包含到我的代码中，或者找到显示Maple数值方法的完整描述的地方。这个方程是一个复杂的代数方程组。
浏览 4提问于2012-09-26得票数 1
我需要解决联立线性方程(5个7个未知数的方程，即一个欠定问题)，其中变量在0-1,00000的大范围内变化。有人能建议我应该使用哪种健身功能吗？
浏览 0提问于2010-09-03得票数 3
5回答
我正在尝试使用SymPy在python中求解一个方程。我有一个生成的方程式(类似于function = y(8.0-(y**3.0))，我在SymPy中使用它来创建一个新的方程式，如下所示:输出y(8.0-(y**3.0)) == 2的eq = sympy.Eq(function, 2)。但是sympy.solve(eq)似乎不起作用。
>>> from sympy import Eq, Symbol as sym, solve
>>> y = sym('y')
>>> eqa = Eq(y(8.0-(y**3.0)), 8
浏览 1提问于2010-12-15得票数 7
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
回答已采纳
如何计算一个值在数学公式中出现的次数？说明:假设在一个圆的方程式中，我想计算y=5的值出现了多少次。或者在抛物线中，方程式的值是5的多少次；或者sin(x)在给定范围内有多少次的值为1？
是否有可能以数学方式获得计数，而不是迭代或递归计算
浏览 0提问于2012-02-15得票数 0
2回答
我正在尝试对一些数据拟合一个圆圈。这需要数值求解一组三个非线性联立方程组(参见的完全最小二乘法)。
在我看来，IDL提供的NEWTON函数很适合解决这个问题。NEWTON需要一个函数的名称，该函数将为自变量的特定值计算方程系统的值：
FUNCTION newtfunction,X
RETURN, [Some function of X, Some other function of X]
END
虽然这很好用，但它要求方程系统的所有参数(在本例中是数据点的集合)都是在newtfunction中硬编码的。如果只有一个数据集需要求解，这是很好的，但是我有数千个数据集，并且不能手动为每个数据
浏览 2提问于2014-10-26得票数 0
1回答
假设我有一个复杂的连续系统模型，它由以下几个微分方程组成：
f = 1/(ura235 + 2.43*plu239 + 2.679*plu241)
d/dt ura236 = (.016*ura235 - 0.0012*uran236)*f
d/dt uran235 = -0.106*f*ura235
.
.
.
如何在Matlab中确定这一系列方程的特征值，而不必手动形成雅可比矩阵？
任何帮助都将不胜感激。
浏览 0提问于2011-01-22得票数 1
2回答
在Dymola 2021 X中显示方程关联图有一个新的特点，但关联图中的非线性方程组比平移日志中的非线性方程组多，因此我不知道如何比较这两种统计数据。如何解读这些信息？
--------------------------------------------------------------------
更新时间为12.23.2020:以下是使用Dymola2021x-HotFIX-1的结果
浏览 3提问于2020-12-02得票数 2
回答已采纳
我在Dymola (基于modelica的商业软件)中运行一个系统模型。模型研究了在强迫对流作用下多孔物质的传热传质问题。我问这个问题的问题是关于模型模拟过程中产生的误差。
具有模型当前设置的在给定的模拟时间(7200秒)内成功运行。但是它在日志文件中生成以下提示(在整个运行时间歇性地)：
ERROR: Failed to solve non-linear system using Newton solver.
To get more information: Turn on Simulation/Setup/Debug/Nonlinear solver diagnostics/Detail
浏览 3提问于2016-04-11得票数 4
回答已采纳
1回答
我会尝试描述我的问题，但最后有一些具体的问题。
我有一些困难，解决一个涡轮机问题，在Matlab。目的是用已知的物理测量和运行条件来计算涡轮的效率。为了做到这一点，我使用了科学论文中给出的方程式：
在我的例子中的问题是，要知道效率，我需要出口气体参数(温度，压力，速度…)这又取决于效率。我主要在一个名为EES (工程方程求解器)的软件中解决了这个问题，在这个软件中，方程可以按任何顺序写入，程序对它们进行分组并自动求解。由于程序中的限制，我切换到了Matlab，这不允许我扩展这个问题。
在matlab中，我首先以一个可解的顺序将所有的方程写成一个函数，除了求解问题所需的两个退出参数之外。将两个退
浏览 1提问于2012-05-31得票数 0
1回答
我想知道在下面的场景中我最好的方法是什么。我有8个未知数，然而实际上是无限数量的非线性方程，这使得系统超定。
未知数：
U M V N J S W N
方程式：
U*M + V * Catime1 – V*M – Mgtime1 = 0
J*M + W * Catime1 – W*M – Srtime1 = 0
U*N + V * Catime2 – V*N – Mgtime2 = 0
J*N + W * Catime2 – W*N – Srtime2 = 0
U*S + V * Catime3 – V*S - Mgtime3 = 0
J*S + W * Catime3 – W*S - Srt
浏览 0提问于2013-05-28得票数 0
据我所知，多元二次方程组存在多项式时间算法，即(eXtended线性化)。但我不知道对于一般的多元三次方程组是否存在多项式时间算法。有人能给我举个例子吗？非常感谢！
浏览 2提问于2014-01-05得票数 4
回答已采纳
我想开发一个用于解线性和非线性方程组的桌面应用程序。我正在考虑在Matlab中使用C#。我以前从未使用过matlab，但我正在努力为我的项目学习。
function F = myfun(x)
F = [2*x(1) - x(2) - exp(-x(1));
-x(1) + 2*x(2) - exp(-x(2))];
**x0 = [-5; -5]; % Make a starting guess at the solution**
options=optimset('Display','iter'); % Option to
浏览 1提问于2011-10-12得票数 1
1回答
Z3Py:从方程组生成抽象公式
、、、、
我的例子:方程组
伪码约束基
a = b+c
∧ e = a*c
∧ a = +2 ; some replaceable concrete values
∧ c = +18
溶液
b = -16
∧ e = -32
我想要的信息
在方程组中，我想得到以下知识：
抽象公式，我可以用它从给定的值(在约束基中)计算变量值(解决方案)。
(就像在高中，老师不仅希望看到结果，还希望看到这样一个经过改造的抽象公式。)
公式如.
b = a-c ; is an equivalent transformation from `a = b+c`
∧ e = (a-c)*c ; is
浏览 3提问于2013-02-28得票数 1
我意识到S盒可以用AES实现非线性变换.然而，我不确定这如何使AES安全。例如，如果我们没有S盒，那么就可以从一组线性方程组中计算密钥：
C^1=Ax+k
C^2=AC^1+k
..。
y=AC^n+k
其中A为线性变换，k为密钥，C为中间密文，n为加密轮数，x为输入，y为最终输出。但是，如果我们加上一个S框，那么它所执行的替换就不可能表示为x，f(x)的函数，所以现在我们有：
C^1=Af(x)+k
C^2=Af(C^1)+k
..。
y=Af(C^n)+k
在我看来，这似乎也是高斯消除的牺牲品(通过将每个方程替换成下一个方程)，尽管这样一个函数用于S盒中的替换可能非
浏览 0提问于2021-06-19得票数 1
1回答
我在试着理解神经网络的隐藏层。输入层部分介绍了在将信息传递到出现关注点的隐藏层之前使用的步骤。
输入层：
根据我的理解，神经网络的第一步是通过选择“最佳”线性函数来加权输入。
Methods用于寻找最优权重(可能与问题无关)：
对于求解二次极小化问题的最优权问题，我通常是这样做的:通过求凸的全局极小(这种情况下，$A$是基矩阵，$b$是理想向量)函数(因为每个临界点都是全局最小值)或正交投影方程(因为误差正交于基的列空间，它们的点积给出了方程- $A^T(Ax-b$)，它等于$x$：$x =(A^TA){-1}A^Tb。
一旦找到了最佳的线性系数(比如$m$作为斜率，$c$作为偏差)，下一步就
浏览 0提问于2018-07-07得票数 0
回答已采纳
我在Dymola建立一个媒体库，类似于Helmholtz-Media，但是对于Ammonia+Water，一种混合物。你得到了很多不明确的可解方程。
由于Modelica介质和流体库的结构，我需要从p，h和x得到热力学状态，状态向量由d，T和x组成。
这是一个简单的示例，如何获得状态向量：
model getState_phX
parameter AbsolutePressure p = 500000 "pressure";
parameter SpecificEnthalpy h = 2500000 "enthalpy";
parameter SI
浏览 3提问于2014-10-08得票数 2
回答已采纳
1回答
我正在用R中的5个方程来估计一个非线性的看似不相关的回归(SUR)模型，并且我正在systemfit软件包上工作。一切都很顺利，直到它需要对我的方程式设置一些限制。使用软件包systemfit，函数nlsystemfit() it适用于非线性方程。但选项/参数restrict.matrix不允许用于nlsystemfit() (它适用于函数systemfit()中的线性方程)。
一个简化的例子是(我认为在这里显示数据是不相关的)：
EQ_1 <- Y1 ~ (c1 - x)*Q + c11*G11 + c12*G12 + c13*G13
EQ_2 <- Y2 ~ (c2 - x)*
浏览 33提问于2020-02-10得票数 3
回答已采纳
2回答
计算四次方差
、、
我正在使用Modelica来求解一个热传递问题的方程组，其中一个方程是辐射，它写为
Ta^4-Tb^4
有没有人可以说，如果方程写成：
(Ta-Tb)(Ta+Tb)(Ta^2+Tb^2)
浏览 2提问于2011-12-15得票数 1
回答已采纳
1回答
我正在尝试实现一个程序来拟合我的实验中的电泳数据。
目的是根据数据集中的峰面积，从电图中峰的相对面积中推导出生物分子相互作用的动力学参数。
由于所有相关的微分方程都是已知的，而且方程组有一个解析解，如下所述：
我开始输入相关的方程式(6，8，13，.从引用的手稿)在matlab。
这样创建的函数可以工作，我可以用它来模拟相互作用物种的电图。
显然，我现在想使用这个函数来拟合实验数据并检索参数(总共8，Va，Vc，MUa，MUc，k，A0，C0，基线噪声)。
其中一些显然是相互关联的。示例值可能是(给出其大小的一个概念)：
params0 = [ ...
8.44E-02
浏览 3提问于2012-09-25得票数 0
Wolfram系统建模器有一个名为的函数，youtube视频( 20:58)介绍了“方程浏览器”，我想知道Dymola是否具有相同的功能，这将有助于调试非线性系统。
另外，它还欢迎其他一些解决方案，比如使用python来可视化方程结构。
这是“等式浏览器”的截图。
浏览 1提问于2020-05-08得票数 6
回答已采纳
我想知道，这就是APM求解器背后的理论方法。
我正在解决一个问题，在这个问题中，我必须解一个没有任何约束的非线性方程组。我正在使用BPOPT解算器，它工作得很好。我只想知道解决这些系统的方法是什么。这是牛顿的方法？准牛顿法？我没有找到任何关于这件事的文件。你能告诉我在哪里能找到这样的信息吗？
是否也有可能预先知道gekko求解器需要的RAM来解一个含有(N)变量的非线性方程组？
谢谢
浏览 3提问于2020-06-17得票数 2
回答已采纳
我想用python来解决下面的方程。问题是因变量'y‘也存在于方程的右侧。第一个问题，这样的方程在数学上是如何命名的？
如果我跳过RHS中的“y”，我就能解决这个问题。但不知道如何解决如何保持它在适当的位置。
我使用了以下代码:将numpy作为np从matplotlib导入作为plt导入pyplot
A=2
B=1.3
C=0.25
D=1.25
def func(x,A,B,C,D):
y=A*np.sinh(((x)/B-C)/D) #I skipped (x-y) here
return y
x=np.linspace(-3,3,200)
y=func(x,
浏览 12提问于2019-06-04得票数 2
15回答
我想要解一组方程，线性的，或者有时是二次的。我没有具体的问题，但我经常遇到这种情况。
使用来解决这些问题是很简单的，它相当于数学的web版本。但这并不能提供iPython shell的舒适性和便利性。
有没有一个简单的库来处理python shell中的线性和二次方程？
就我个人而言，我发现使用Casio991MS科学计算器非常方便。我知道如何设置变量，如何求解方程，以及如何做很多事情。我想要一个这样的工具，最好是可以从ipython shell中使用。我很惊讶竟然没有找到任何东西。我对sage的印象还不够深刻；也许我遗漏了什么。
浏览 1提问于2009-10-29得票数 25
1回答
我使用的是倍频程3.8.1，它的工作原理类似matlab。
我有一个数千个值的数组，我只包含了三个分组，如下所示：
(amp1=0.2；freq1=3；phase1=1；是一个分组的例子)
t=0;
amp1=0.2; freq1=3; phase1=1; %1st grouping
amp2=1.4; freq2=2; phase2=1.7; %2nd grouping
amp3=0.8; freq3=5; phase3=1.5; %3rd grouping
下面的Octave / Matlab代码解决了Y的问题，所以我可以将它插入方程中，以检查值和计算值(而不是位于数组中)。
clea
浏览 6提问于2015-08-08得票数 2
回答已采纳
1回答
我想用Sympy解决简单的2D运动学逆问题。我知道x和y位置的正运动学方程。
x = l1*cos(theta1) + l2*cos(theta1+theta2)
y = l1*sin(theta1) + l2*sin(theta1+theta2)
如果我知道theta1和theta2这两个方程，如何用Sympy求解这两个方程？
浏览 2提问于2019-03-29得票数 1
大家好，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
3回答
我收集了一个有约束条件的非线性方程组，其数量可能会有所变化，我想用一些数值方法来求解。
我已经能够用Solver在Excel中解决一个简单的(一个等式)情况，但是以前没有在Python中把类似的东西放在一起，所以我希望能对方法提出一些建议。
在做了一些调查之后，它似乎是解决这些系统的一种流行方法。对于一个简单的两种方程情况，我的问题采用以下形式，为了清晰起见，分为几个部分：
第二种方程的形式是一样的。
A是常数，变量Z，S和x是每个实体i的常数，唯一独立的值是指数a和b；两个方程，两个未知数，所以应该有一个唯一的解。
如前所述，我在Excel中建立了简单的一个方程，并使用Solve
浏览 0提问于2018-10-31得票数 1
回答已采纳
1回答
用FiPy求解柱面几何中扩散方程的稳态解与用其它软件(如Mathematica )得到的解有很大不同。 方程式是： $0 = \frac{1}{r}\frac{d}{dr}\left(\frac{r}{T^{1/2}}\frac{dT}{dr}\right) + cte*T^{3/2}$ 这意味着，通过使用CylindricalGrid1D网格，我们可以将方程写成： mesh = CylindricalGrid1D(nr=100, dr=0.01, origin=0.0)
T = CellVariable(name='temperature', mesh=mesh, hasO
浏览 25提问于2019-02-13得票数 0
回答已采纳
9回答
使用Python. (Numpy，Scipy或Sympy)解决对非线性方程组的(最好的)方法是什么？
例：
x+y^2 =4
e^x+ xy =3
解决上述两个问题的代码片段将是很棒的。
浏览 12提问于2012-01-05得票数 85
回答已采纳
我试图用R中的Ryacas软件包求解一个3组方程。
library(Ryacas)
yacas("OldSolve({(f*Reg/(k + Reg)) - r*BB - ab * BB - a1*f1*BB*AA*P - fp*BB*P==0,(ab * BB + a1*f1*BB*AA*P)/e1 - m1 * AA == 0, (fp * BB * P)/ep - f1*AA*P - mp * P ==0 },{BB,AA,P})")
但我在R中得到了下面的答案
expression(list(list(BB == BB, AA == (f * Reg/(k + Reg
浏览 0提问于2018-12-14得票数 1
回答已采纳
嗨，我正在尝试实现一个简单的塑性有限元代码。我不清楚如何设置通过拉格朗日乘数应用的位移控制。在线性问题的情况下，我执行了以下操作来应用显示控制，并且它起作用了。如果能量平衡是线性的，那么将得到以下方程：
然而，在非线性问题中，K是K(u)并且是非线性的。因此能量平衡方程及其导数如下所示
感觉不对劲。我该如何评估它？如何在非线性情况下修改方程的矩阵形式？那么矩阵方程是如何修改的呢？我使用添加了用latex输入的图像。如果不清楚，请让我知道。
浏览 9提问于2021-06-24得票数 0
2回答
我需要知道如何解一个非线性方程组，但是改变一个参数，这样每次你改变那个参数就会抛给我那个系统的结果(需要所有的结果)，我想了一个for，它改变了参数，求解方程，每个结果都存储在一个电子表格中，问题是，由于你不能求解这个系统，所以我抛出了非符号数值，它们给出了一个必须求解的系统的例子：
0 = 125 +100 * cos (x) -25 * cos (a) -175 * cos (y)
0 = 100 * sin (x) -25 * sin (a) -175 * sin (y)
在要更改的参数中是a，并且将在电子表格中保留x和y的相应值。
浏览 1提问于2012-09-30得票数 0
我试图在python中使用FiPy来解决pde问题。它显示的是RuntimeError: Factor正好是奇异的。这不是在这些方程中的情况，而是在我计算的大多数方程中。
(1+αv)uv/(1+h(1+αv)u) +D∇^2 u
(1+αv)uv/(1+h(1+αv)u) -v+D∇^2v
from fipy import *
nx=ny=100
dx=dy=0.25
L=dx*nx
dt=0.01
sigma=10.0
h=0.1
C=0.8
alp=0.57 #alpha
mesh =Grid2D(dx=dx,dy=dy,nx=nx,ny=ny)
u=CellVariable(nam
浏览 14提问于2020-07-31得票数 0
回答已采纳
1回答
这是我关于stackoverflow的第一个问题。抱歉，如果它有点笨拙。 我正在尝试解决描述层状土壤中的梁的静态平衡。就作用在梁上的力而言，有一些外力，还有土对外力的反作用力。显然，行动和反应必须平衡。反作用力是来自不同土层的力的总和。土的响应是位移的非线性函数。基本上，这个特殊的平衡是一个标量，因为它就是反应=作用的总和。 另一方面，梁内的内力也必须与作用在所讨论的梁单元上的局部力相平衡。假设我们有40个这样的元素。因此，我面临的问题概括如下： 我需要找到F1(x1，x2，...x40)的根，用x1，x2...x40描述桩在土壤中的位移。我传递并初步猜测了位移，Fsolve返回了一个向量长
浏览 43提问于2020-05-03得票数 0
回答已采纳
我有一个问题，可以简化为下面的陈述(实际的方程和函数要复杂得多)：
系统：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
我在试着解非线性方程组。 为此，我成功地使用渐近在线解决了一个由三个方程组成的系统。但对于7个等式，它不再起作用，返回给我一条"Operation timed out“消息。所以，我尝试使用标准的fsolve来解决这个问题。但是，虽然我只用一组简单的3个方程就可以得到很好的渐近结果，但fsolve给了我很差的结果。 这是我写的代码(感谢我在这个论坛上找到的)： def equations(p):
w1, w2, m = p
f0 = w1+w2-1
f1 = w1*m+w2*(1-2*m)-1/3
f2 = w1*m**2+w2*(1-2*m)**2-
浏览 27提问于2019-04-29得票数 1
回答已采纳
如何使用LFSR产生二进制输出序列，但增加了线性复杂度？
我试过搜索，但在方程的细节上弄糊涂了。所以，有谁能帮我解释一下如何增加线性复杂度，比如一些图表。谢谢
浏览 0提问于2019-04-06得票数 2
我改变了我的模型中的参数，但有时它会导致模拟停止，因为病态雅可比，我不知道这个错误意味着什么。
我的问题是：
对于如何在Modelica模型中求解DAE方程，有什么解释吗?这样我就可以理解Jacobian矩阵的使用了吗？
我怎么知道哪个方程或参数导致这个问题在Dymola？
浏览 3提问于2020-07-12得票数 4
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
云直播活动推荐
腾讯云开发者
