用单变量Python求解非线性方程
非线性方程是指未知数的幂次大于1的方程，例如x^2 - 2 = 0。求解非线性方程可以使用数值方法，其中一种常用的方法是使用Python编程语言中的数值计算库来实现。
在Python中，可以使用SciPy库中的optimize模块来求解非线性方程。optimize模块提供了多种数值求解器，其中包括用于求解非线性方程的函数fsolve。
下面是使用fsolve函数求解非线性方程的示例代码：
from scipy.optimize import fsolve
# 定义非线性方程
def equation(x):
return x**2 - 2
# 求解非线性方程
result = fsolve(equation, 1) # 传入方程和初始猜测值
print(result)在上述代码中，首先定义了一个名为equation的函数，该函数表示要求解的非线性方程。然后使用fsolve函数传入该方程和一个初始猜测值来求解非线性方程。最后打印出求解结果。
对于非线性方程x^2 - 2 = 0，使用初始猜测值1进行求解，得到的结果是[1.41421356]，即方程的解为1.41421356。
推荐的腾讯云相关产品和产品介绍链接地址：
- 腾讯云函数计算（SCF）：https://cloud.tencent.com/product/scf
- 腾讯云弹性容器实例（TKE）：https://cloud.tencent.com/product/tke
- 腾讯云云服务器（CVM）：https://cloud.tencent.com/product/cvm
- 腾讯云数据库（TencentDB）：https://cloud.tencent.com/product/cdb
- 腾讯云对象存储（COS）：https://cloud.tencent.com/product/cos
- 腾讯云人工智能（AI）：https://cloud.tencent.com/product/ai
- 腾讯云物联网（IoT）：https://cloud.tencent.com/product/iot
- 腾讯云移动开发（MPS）：https://cloud.tencent.com/product/mps
- 腾讯云区块链（BCS）：https://cloud.tencent.com/product/bcs
- 腾讯云元宇宙（Metaverse）：https://cloud.tencent.com/product/metaverse
请注意，以上链接仅供参考，具体产品选择应根据实际需求进行评估和决策。
相关·内容
我想知道并试图搜索它，但我不知道如何谷歌它，如何解决线性方程组，其中方程是平等的，彼此。我能够为两个方程编写一个具体的算法，但不能对N个方程进行动态计算。我将给出三个例子(我的方程大致是怎样的)：
C1，C2，C3，X -都不知道，但我不需要知道X的结果。
可以这样解释(最后一个方程C1 + C2 + C3 =1不在这里集成)：
请不要试图解决这个问题，我不确定这些方程是否有结果。我只是随机输入系数。但这就是我的方程式的样子。只有不同的系数。最好是在蟒蛇，但如果我知道如何用矩阵做数学的话，我会很高兴。我试着用两个未知数计算它，最后我得到了二次方程，所以在有三个未知数的情况下
浏览 10提问于2022-03-24得票数 1
回答已采纳
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
回答已采纳
1回答
我会尝试描述我的问题，但最后有一些具体的问题。
我有一些困难，解决一个涡轮机问题，在Matlab。目的是用已知的物理测量和运行条件来计算涡轮的效率。为了做到这一点，我使用了科学论文中给出的方程式：
在我的例子中的问题是，要知道效率，我需要出口气体参数(温度，压力，速度…)这又取决于效率。我主要在一个名为EES (工程方程求解器)的软件中解决了这个问题，在这个软件中，方程可以按任何顺序写入，程序对它们进行分组并自动求解。由于程序中的限制，我切换到了Matlab，这不允许我扩展这个问题。
在matlab中，我首先以一个可解的顺序将所有的方程写成一个函数，除了求解问题所需的两个退出参数之外。将两个退
浏览 1提问于2012-05-31得票数 0
大家好，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
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
3回答
用倍频程求解方程
、、
我有一个简单的方程想要解
num1=-2
num2=-3
x+num2=num1
x+-3=-2
x=1
我怎么才能在八度范围内做到这一点。在matlab中，我可以做y =(‘x-3= -2')，但是这在倍频程3.8.1中不起作用，这是我使用的版本。我怎样才能得到解这类方程的八度呢？
我对解决方案的数值感兴趣。
浏览 2提问于2014-09-05得票数 3
回答已采纳
我想知道，这就是APM求解器背后的理论方法。
我正在解决一个问题，在这个问题中，我必须解一个没有任何约束的非线性方程组。我正在使用BPOPT解算器，它工作得很好。我只想知道解决这些系统的方法是什么。这是牛顿的方法？准牛顿法？我没有找到任何关于这件事的文件。你能告诉我在哪里能找到这样的信息吗？
是否也有可能预先知道gekko求解器需要的RAM来解一个含有(N)变量的非线性方程组？
谢谢
浏览 3提问于2020-06-17得票数 2
回答已采纳
我想解决一个由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
我需要一个算法来解决像这样的方程：
8*p^2+8*q^2+8*p-16=0
p^2+q^2+12*p+20=0
如果有人能告诉我algorithm的名字，那就足够了。我可以遵循任何与矩阵相关的算法，但这些线性内部的二次型导致了一个问题。
浏览 2提问于2012-12-17得票数 0
回答已采纳
在Dymola建立了一个具有朗肯循环的太阳能电站模型，虽然初始化工作正常，但是在2840s的时候，由于非线性方程组的求解失败，导致仿真停止。如下面的截图所示，Dymola建议提供更好的起始值，但这没有任何意义，因为失败发生在模拟期间，而不是初始化过程中。有了更好的起始值，情况根本不会改善。我的问题是:我应该如何处理在模拟过程中无法求解非线性方程？ ?
浏览 71提问于2020-07-06得票数 1
1回答
我理解数值求解的逻辑，即matlab从未知值开始循环，直到方程收敛到一个值为止。然而，我不明白的是如何正确地进入我的方程式。我认为使用ode45函数是最好的方法。我有下面的等式U^(n+1) = U^n - (t'*3250/10)-(t'/2)(.004(v^n)^2/10。我怀疑我的方程需要不同的形式，但是我不确定让matlab求解方程的正确方法。
我试着在matlab中按原样输入方程，但是它抱怨v和n是未知变量，我不确定如何处理它们。求解这个方程的最终目标是找到v的值。
C=@(t,v) u^n-(3250*t'/10)-(t'/2)*((.004*(v^
浏览 1提问于2014-09-16得票数 0
我试图用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
回答已采纳
我有一套4套政策制订和评价报告：
du/dt + A(u) * du/dx = Q(u)
其中，u是一个矩阵，包含：
u=[u1;u2;u3;u4]
A是一个4×4矩阵。A和q是u=u1;u2;u3;u4.的函数Q为4*1。
但我的问题是：
如何在MATLAB中求解上述方程？
如果用Matlab的PDE函数进行求解，能不能将其转换成一个简单的函数，而不是从Matlab的现成函数中得到的？
是否有任何方法可以显式地计算A、和Q。我的意思是，在每一个时间步骤中，我从前一个时间步骤的数据中计算出A和Q，并在使程序运行速度更快的方程中加入新的值。
浏览 8提问于2014-01-25得票数 1
2回答
我是Freefem++的新手，我想要解决的问题是费舍尔方程：
du/dt = d^2u/dx^2 + d^2u/dy^2 + k * u * (1-u)
du/dn = 0 - border condition
我试图用弱形式重新描述这个问题，但是Freefem在公式上出现了一个错误：
problem Fisher(uh, vh) = int2d(Th)(uh * vh/dt + Grad(uh)' * Grad(vh) ) - int2d(Th)(k * uh * vh) + int2d(Th)(uh0 * vh/dt) - int2d(Th)(k * uh * vh * uh);
浏览 4提问于2014-10-10得票数 0
回答已采纳
1回答
本文尝试在MATLAB中求解一个基于解矢量的节点热源四面体有限元的热扩散问题。非线性方程组如下所示：
B_U‘+ A_U = q(T)
B为热吸附矩阵，A为电导率矩阵，Q为源项，U为温度。我使用亚当斯-巴什福思/梯形规则预测器-校正方案与Picard迭代，然后一个时间步长控制。源项的温度是在最后一步的温度和预测器的温度之间精确地计算出来的。以下是预测校正器代码的简化版本。源的计算是一个函数。
% predictor
K0 = t(n)-t(n-1);
Upre(dirichlet) = u_d_t(coordinates(dirichlet,
浏览 1提问于2014-12-05得票数 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
我正在尝试建立一个丝杠传动的理论模型，将其作为力元素参与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
2回答
Sympy“不相等”方程
、、、
有什么方法可以解决包含≠的代数方程吗？
例如: 2-x≠3
所以，也许看起来是这样的：
print(solve(not_equal(2-x, 3)))
浏览 3提问于2022-07-03得票数 0
回答已采纳
1回答
有没有人建议使用类似于iOS中的GRG型非线性求解器的等同GRG？
它可以是一个c/c++的解决方案，可以很容易地集成到iOS中，即可以求解类似于EXCEL GRG非线性求解器的非线性方程组。
敬请指教。
浏览 0提问于2015-06-25得票数 4
我试图编写一些R代码来执行约束优化，其中约束不是通常的(g(x) < c)类型。
我试图解决的问题类型的一个例子是：
考虑到每年的假设价格，一家公司正试图确定未来5年的年产量。它每年可以产生最多X，但在任何给定的3年窗口，产量不能超过2X。此外，在任何特定年份，除非产量为0倍，否则产量不得低于0.5倍。其他约束也同样不可能写成一个集合方程(例如，如果X> a，Y< c)。
我很乐意在Excel中解决这类问题，因为在这种限制下非常灵活是相对容易的(例如，每年都有一个if语句，它查看当前、过去和明年产生的数量，如果小于2X，返回0，在单元格中求和if语句，并将该变量的约束设置为0
浏览 3提问于2022-07-24得票数 0
3回答
我有个大问题。我需要用一个C++函数或类来解决一个由3个变量的3个方程组成的非线性系统。我想过用牛顿-拉夫森方法来求解。不幸的是，我没有找到可以为我做这件事的源代码。会有人知道这样的程序吗？我快要决定自己构建它了。谢谢
浏览 0提问于2013-06-01得票数 1
回答已采纳
我试图找到参数值，以便用给定的一组点最好地拟合下面的模型函数。
Model:
x
y = ------------------------
1- (a1*x)^2 + (a2*x)^4
where parameters are: a1, a2
(1)我用这个方程找到了最优参数值，但它给出了任意x(红色图)的y=~0值。
(3)令人惊讶的是，使用以下只包含一个参数的模型得到了更好的拟合。
Model:
x
y = -------------
1- (a1*x)^2
where para
浏览 7提问于2022-08-09得票数 1
2回答
我有一道数学题，是用Maple解决的，但它必须用C++解决。问题是我找不到一个方程的数值方法(它是用内置的maple "solve“过程求解的)。我可以以某种方式将这个过程包含到我的代码中，或者找到显示Maple数值方法的完整描述的地方。这个方程是一个复杂的代数方程组。
浏览 4提问于2012-09-26得票数 1
我想开发一个用于解线性和非线性方程组的桌面应用程序。我正在考虑在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
我有个问题
我想知道如何用matlab求出这样的二阶微分方程：
x''(2nd order) = A + B sin(x)
x = ?
A,B are parameters like 1,2
浏览 1提问于2014-06-13得票数 0
回答已采纳
1回答
当我试图用RSymPy库求解一个方程时，我得到了NotImplementedError:多元方程误差。
我不知道我错了什么地方，或者我应该使用其他的多元方程库，或者有一个解决办法来使它工作？
此外，如果不对其进行硬编码，我也无法将变量形式的值传递给方程参数(ymax、ymin和ec50)。
library(rSymPy)
ymax <- Var("ymax")
ymin <- Var("ymin")
ec50 <- Var("ec50")
sympy("var('x y')")
sympy(
浏览 4提问于2017-02-24得票数 0
回答已采纳
Wolfram系统建模器有一个名为的函数，youtube视频( 20:58)介绍了“方程浏览器”，我想知道Dymola是否具有相同的功能，这将有助于调试非线性系统。
另外，它还欢迎其他一些解决方案，比如使用python来可视化方程结构。
这是“等式浏览器”的截图。
浏览 1提问于2020-05-08得票数 6
回答已采纳
1回答
我有一个线性方程Ax = b，我试图用Python中的纸浆来最小化(Ax)^2。A的尺寸是(1000,500)，b是(1000，)。
到目前为止，我已经尝试过：
import pulp
mse = pulp.LpProblem("Example_Problem", LpMinimize)
vars = pulp.LpVariable.dicts('x', range(len(A[0])), lowBound=0,upBound =1,
cat='Integer')
for row, rh
浏览 4提问于2018-07-26得票数 0
我有一个非线性系统，它有二次输入(操纵变量，MV)。为了使用GEKKO对该系统进行MPC模拟，我定义了以下函数来处理二次输入： def NTVin(T, u):
n = T.shape[0]
Tntv = zeros(n)
u2 = zeros_like(u, dtype = float64)
for i in range(len(u2)):
u2[i] = copy(u[i].VALUE)
for i in range(n):
Tntv[i] = dot(dot(u2, T[i,:,:]), u2)
retu
浏览 8提问于2019-10-02得票数 3
回答已采纳
3回答
如何优化非线性方程组的解？
、、、、
我有一些非线性方程，例如：
Y= f1(X)
Y= f2(X)
..。
Y= fn(X)
一般来说，它们没有精确的解，所以我使用来解决它们。方法是基于迭代的，我正在寻找优化计算的方法。减少计算时间的方法是什么？避免计算平方根或其他数学函数？也许我应该在C++代码中使用程序集(解决方案是用C++编写的)？
浏览 2提问于2010-12-08得票数 2
1回答
二维点间距离参数化的有效方法
、、、、
我们一直在开发一个小的简单的"CAD“解决方案，使我们可以参数化一些特定的，简单的形状的宽度和长度。
例如，考虑形成三角形的下列顶点集。其中任意两个点形成一条线。因此，改变点之间的距离就是改变线的宽度。
我们已经对如何处理这个问题进行了严格的讨论。
我们讨论过的问题是：
维护所有顶点之间所有关系的方程列表。假设我们有点A，B，C。让W是某个用户定义的参数。该形状的约束方程是Bx = Ax + W、By = Ay和Cx = Bx等。
它的复杂性是巨大的，但它起作用。
也许把每个顶点建模成图中的一个节点.？
在这一领域中广泛使用的合适方法是什么？
浏览 3提问于2017-08-18得票数 0
回答已采纳
在MATLAB中求解约束线性和非线性方程组
我在MATLAB中用直接刚度法解决了一个有限元问题。这个问题现在被表述为一个非线性方程组：
KU=F or CF=U with C = K^{-1}.
该问题用3-Dd表示，有200多个节点。模型被用作变形估计器，力始终是变形的输入和输出。所以有很多线性方程。现在这个问题需要扩展，需要实现两种类型的约束。某些节点必须受到约束，以便它们只能围绕给定点以圆周运动方式移动。并且一些节点必须被限制在最大变形范围内。
总共有成百上千个形式为方程的方程：
a_1 X(1) + a_2 X(2) + ... - F(1) = 0
大约10个非线性形式的等式约束(也
浏览 3提问于2016-05-31得票数 0
1回答
这是我关于stackoverflow的第一个问题。抱歉，如果它有点笨拙。 我正在尝试解决描述层状土壤中的梁的静态平衡。就作用在梁上的力而言，有一些外力，还有土对外力的反作用力。显然，行动和反应必须平衡。反作用力是来自不同土层的力的总和。土的响应是位移的非线性函数。基本上，这个特殊的平衡是一个标量，因为它就是反应=作用的总和。 另一方面，梁内的内力也必须与作用在所讨论的梁单元上的局部力相平衡。假设我们有40个这样的元素。因此，我面临的问题概括如下： 我需要找到F1(x1，x2，...x40)的根，用x1，x2...x40描述桩在土壤中的位移。我传递并初步猜测了位移，Fsolve返回了一个向量长
浏览 43提问于2020-05-03得票数 0
回答已采纳
1回答
下面是我用scipy.optimize.leastsq解决问题的一个最小例子
from scipy.optimize import leastsq
from numpy import array, exp, sin, cos
def MatrixFun(x, *par):
a, b, c, d = par
m11 = a*sin(x[0])+b*cos(x[1])
m12 = c*cos(x[0])*sin(x[1])
m21 = c*sin(x[0])/cos(x[1])
m22 = d*exp(x[0]*x[1])
M = array([[
浏览 1提问于2016-02-29得票数 0
我用Python构造了一个线性规划问题，并使用了CPLEX求解器。变量被表示为像docplex.mp.LinearExpr(.)这样的复合表达式。我们的目标是构造一个目标，在这个目标中，我们使用类似于幂-1或3的表达式对变量进行优化：
选择1/变量
不幸的是，我收到以下错误消息；docplex.mp.utils.DOcplexException:无法引发“docplex.mp.LinearExpr(.)”到电源-1。变量的指数必须是0，1或2。
有没有人知道是否有办法解决这个问题？那么，是否有办法将某个数字除以相当于docplex.mp.LinearExpr(.)的值？
浏览 12提问于2022-08-18得票数 0
5回答
我有六个参数方程，使用了18个不同的变量(实际上不是26个)，其中6个是未知的。
我可以坐下来几页纸，计算出每个未知数的方程式是什么，但有没有一个简单的程序解(我正在考虑在Matlab中)，会吐出我正在寻找的六个方程式？
编辑：耻辱这已经关闭了，但我想我能理解为什么。如果有人还感兴趣，方程(我相信)是非线性的：
r11^2 = (l_x1*s_x + m_x)^2 + (l_y1*s_y + m_y)^2
r12^2 = (l_x2*s_x + m_x)^2 + (l_y2*s_y + m_y)^2
r13^2 = (l_x3*s_x + m_x)^2 + (l_y3*s_y + m_y)^2
浏览 3提问于2009-12-11得票数 3
回答已采纳
我想用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
3回答
我试图绘制一个包含两个联立方程的图，但我不需要解它们，我只是想从代换中得到多个结果，比如当x为1时，或者y为0时。
我的方程式是5x + 2y = 20，y = 2x + 1
我找到的所有解都是为了解这个方程，而不是用值替换来得到多个结果。
最好是用numpy或sympy函数，我正在试着学习那些哈哈
浏览 1提问于2021-05-29得票数 1
我正试图使用Apache-Commons-数学求解器来解决"V“的以下数学方程。
方程式：
V= log(X/ (V-1))
以下是MWE；
import org.apache.commons.math3.analysis.UnivariateFunction;
import org.apache.commons.math3.analysis.solvers.*;
public class Test {
public static void main(String[] args) {
UnivariateFunction function = v ->
浏览 6提问于2017-09-29得票数 0
回答已采纳
我将尽量以更好的方式解释我的怀疑：
我试图用PDE工具箱(Matlab)求解反应扩散方程，得到解(u)的语法是：
抛物-解抛物型PDE问题
该函数给出了标量PDE问题的有限元解: u1 =抛物线(u0，tlist，b，p，e，t，c，a，f，d)
c，a，f，d是抛物型方程的系数：
d(∂u/∂t)−∇⋅(c∇u)+au=f,
但我想要这个：
u=parabolic(u0,tlist,b,p,e,t,c,@coeffunction,f,d);
u0是初始解，tlist是一个向量数组，其中我记录了解，b对应于边界条件，(p，e，t)是网格的元素，@ coefficient是对应于系数&#
浏览 20提问于2016-06-16得票数 2
我有一个应用程序，其中我需要使用迭代技术来求解变量。基本上，我需要猜测a的值，使K= 1.2；以下是我拥有的参数和方程：
b = .02;
M = 3;
Sg = 4.5;
以下是方程式：
X = pi*alpha/2;
F = sqrt((1/X)*tan(X))*(.923 + .199*(1 - sin(X))^4)/cos(X);
K = F*Sg*sqrt(pi*a)
当我猜测.010125的值时，我得到了K的适当值1.2时，我如何使用带有初始猜测值的迭代技术来获得a的.010125的适当值？
下面是完整的代码：
%ABS Plastic Under Moment Load
clc
浏览 0提问于2017-02-20得票数 0
1回答
我很难理解如何将不等式方程加到函数中。 例如： 下面是这些包： import numpy as np
from scipy.optimize import fsolve 下面是我想要使用的公式： x1 >= 0.4 and x1 <= 0.7
x2 >= 0.2 and x2 <= 0.4
5x2**2 + 2x1**3 = 2 这就是我要创建的函数： myFunc(z):
x1 = z[0]
x2 = z[1]
F = np.empty((3))
F[0] = x1 >= 0.4 and x1 <= 0.7 #
浏览 151提问于2021-04-10得票数 0
我在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回答
具有多重方程的枕方程组
、、、
你好，当我有一个方程组时，我在使用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
回答已采纳
1回答
所以我知道有很多关于这个话题的问题/答案，但是我还没有找到我的尝试中到底出了什么问题。我有两个非线性函数f(x,y)和g(x,y)，我正在试图求解这个系统。
f(x,y) - g(x,y) = 0
f(x,y) + g(x,y) = c
其中c是一些正常数。我一直在使用这个问题的答案：中描述的片段，但是我面临着一些问题。如果我为我的代码运行这段代码段，它将返回x和y值，以便只满足系统中的第二个方程--即，它返回x和y，从而得到f(x,y) + g(x,y) = c，而对于另一个方程，则返回f(x,y) - g(x,y) != 0。在使用scipy.optimize.root函数时，我会遇到完全
浏览 4提问于2020-06-04得票数 0
回答已采纳
1回答
所以我试图解一个有3个未知因子的方程。我决定用nlm。
我定义了函数F，它接受三个参数，放在一个向量中，我想要找到的是vetor X，它验证了以下方程：
F(X)-F(X1)-F(X2)-F(X3)=0
所以我向LHS申请了nlm。但我得到了一些奇怪的结果，而不是有一个使LHS接近于零的解决方案，而是给出了使LHS收敛到-infinite的解决方案。
谁能给我指明正确的方向。
(预先谢谢大家:)
rm(list=ls())
Ta <- 30 #commun parameter
c <- 0.09 #commun parameter
Delta_T <- c( 10,
浏览 0提问于2018-07-02得票数 0
回答已采纳
2回答
我正在编写一个程序来处理各种方程。我有一个方程列表，还有一个函数，它可以解出一个未知值的给定方程。
对于任何变量，a+b+c=0 ->都可以求解，给出另外的2个变量。
现在我需要编写一个算法，它接受关系列表、已知值列表和目标，然后找到一系列必须采取的步骤来查找未知信息。
例如：(a = b，b= c) a= 1，查找c：
结果显示，您首先使用"a = b“来转换-> b，然后将b -> c与"b = c”进行转换。
对于我可以研究的算法或搜索术语，我们将不胜感激。
浏览 4提问于2022-05-23得票数 0
回答已采纳
2回答
我有一个简单的方程式要解决，我想在matlab中做，然而，有一些年，我没有使用，并有点忘记它。所以这就是
({(387=mod(324*a+b，601))，(491=mod(381*a+b,601))}，{a，b})，
没有给我任何结果，因为我找不到一个错误。最初的方程组是：
(324a+b)mod601=387 (381a+b)mod601=491，
显然是仿射密码。谢谢!
浏览 3提问于2014-11-30得票数 0
我在寻找一种快速算法来求解3个未知变量上的N个多项式方程组。也就是说，给定3个函数，F0(x,y,z), F1(x,y,z)... FN(x,y,z)，我想找到x, y, z这样的F0(x,y,z) = F1(x,y,z) = ... = FN(x,y,z) = 0。
我试着在几个不同的地方找到解决方案，但我只能在代数几何学或密码学等主题上找到非常高级的论文。不过，我需要的是一个简单/快速的算法，它返回一个快速的数值解。有这样的算法吗？
浏览 7提问于2015-01-11得票数 3
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
