如何求解关于指定变量的渐近非线性系统
渐近非线性系统是指系统的动态行为在无穷远处趋于稳定的非线性系统。求解关于指定变量的渐近非线性系统可以通过以下步骤进行：
- 确定系统的数学模型：根据实际问题，建立描述系统行为的数学模型。这可以是基于物理原理的方程组，也可以是基于数据的统计模型。
- 线性化系统：对于非线性系统，可以通过线性化的方法来近似求解。线性化可以通过泰勒级数展开或者雅可比矩阵的计算来实现。
- 求解线性系统：对于线性化后的系统，可以使用线性系统的求解方法，如矩阵运算、特征值分析等来求解系统的稳定性和动态行为。
- 分析非线性特性：通过对线性化后的系统进行分析，可以得到系统的稳定性、振荡性、收敛性等特性。这些特性可以帮助我们理解系统的行为，并对系统进行优化和控制。
- 仿真和实验验证：通过数值仿真或者实验验证，可以验证求解得到的系统模型和特性是否与实际系统一致。这可以帮助我们进一步改进和优化系统设计。
在求解关于指定变量的渐近非线性系统过程中，可以使用腾讯云提供的相关产品来支持系统建模、仿真和实验验证等工作。例如：
- 腾讯云人工智能平台（https://cloud.tencent.com/product/ai）：提供了丰富的人工智能算法和工具，可以用于系统建模和仿真。
- 腾讯云数据库（https://cloud.tencent.com/product/cdb）：提供了高性能、可扩展的数据库服务，可以用于存储和管理系统数据。
- 腾讯云云服务器（https://cloud.tencent.com/product/cvm）：提供了弹性、可靠的云服务器，可以用于系统的部署和运行。
- 腾讯云区块链服务（https://cloud.tencent.com/product/bcs）：提供了安全、可信的区块链服务，可以用于系统的数据验证和交互。
通过利用腾讯云的相关产品，可以更高效地进行渐近非线性系统的求解和优化工作。
相关·内容
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
我想知道，这就是APM求解器背后的理论方法。
我正在解决一个问题，在这个问题中，我必须解一个没有任何约束的非线性方程组。我正在使用BPOPT解算器，它工作得很好。我只想知道解决这些系统的方法是什么。这是牛顿的方法？准牛顿法？我没有找到任何关于这件事的文件。你能告诉我在哪里能找到这样的信息吗？
是否也有可能预先知道gekko求解器需要的RAM来解一个含有(N)变量的非线性方程组？
谢谢
浏览 3提问于2020-06-17得票数 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
回答已采纳
我想开发一个用于解线性和非线性方程组的桌面应用程序。我正在考虑在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回答
下面是我用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
2回答
假设我试图优化给定一些数据的函数的参数。假设参数是一个三元素向量.使用constrOptim函数，我可以设置参数之和必须大于或等于1 .的条件(使用constrOptim函数的ui和ci参数)。
但是如何限制三元素参数向量只有2个非零元素呢？我想不出一种方法-我想指定一个4元素参数向量，其中第4个元素是向量的前3个元素的函数(即。在向量的前3个元素中计数非零的数目，然后将其限制在一个特定的数目内)，但不能使任何东西工作。有什么想法(非常感谢)吗？
#Some simple code as context:
data=cbind(rnorm(5)+1,rnorm(5)+1,rnorm(5)+1)
浏览 4提问于2016-04-13得票数 0
回答已采纳
我想知道并试图搜索它，但我不知道如何谷歌它，如何解决线性方程组，其中方程是平等的，彼此。我能够为两个方程编写一个具体的算法，但不能对N个方程进行动态计算。我将给出三个例子(我的方程大致是怎样的)：
C1，C2，C3，X -都不知道，但我不需要知道X的结果。
可以这样解释(最后一个方程C1 + C2 + C3 =1不在这里集成)：
请不要试图解决这个问题，我不确定这些方程是否有结果。我只是随机输入系数。但这就是我的方程式的样子。只有不同的系数。最好是在蟒蛇，但如果我知道如何用矩阵做数学的话，我会很高兴。我试着用两个未知数计算它，最后我得到了二次方程，所以在有三个未知数的情况下
浏览 10提问于2022-03-24得票数 1
回答已采纳
在Dymola建立了一个具有朗肯循环的太阳能电站模型，虽然初始化工作正常，但是在2840s的时候，由于非线性方程组的求解失败，导致仿真停止。如下面的截图所示，Dymola建议提供更好的起始值，但这没有任何意义，因为失败发生在模拟期间，而不是初始化过程中。有了更好的起始值，情况根本不会改善。我的问题是:我应该如何处理在模拟过程中无法求解非线性方程？ ?
浏览 71提问于2020-07-06得票数 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
1回答
在MiniZinc (windows )中，我如何解析:平板锌:错误: FD解算器后端不支持‘`var’类型的变量。我意识到我需要一个不同的解决程序，但找不到一个程序安装一个和首选项dlg似乎不起作用。另外，还不清楚哪个解算器会起作用。
浏览 8提问于2015-01-31得票数 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
回答已采纳
1回答
我需要解决一个问题:假设湖中的鱼类种群P(t)在时间t= 0被疾病攻击，导致鱼类停止繁殖(因此出生率为β= 0)，死亡率δ(每周每条鱼的死亡人数)与1/√P成正比。如果最初湖中有900条鱼，6周后剩下441条鱼，那么湖中所有的鱼需要多长时间才会死亡？我不明白如何添加条件P (6) = 441？ r=np.array([6])
p=np.array([441])
def fun (y , t, r, k):
P=y
dydt = -k*math.sqrt(P)
return dydt
y0 = np.array([900])
t = np.linspace(1, 9,
浏览 8提问于2019-08-29得票数 1
我正在尝试用MATLAB的bvp4c来调试我的“奇异雅可比”问题。我的问题是:哪些雅可比人正在被评估？我们有常微分方程y‘= f(x，y)的非线性系统。我假设f_i关于y_j的导数被放入一个矩阵中，并在y的某个值上求值。也许这个矩阵对于bvp4c算法来说应该是可逆的，尽管我发现该算法在不是这样的情况下是收敛的。
边界条件g(ya，yb) = 0的残差是什么？我假设g_i相对于ya_j的导数是取的。这给出了一个矩阵。此外，还可以求g_i相对于yb_j的导数。这也给出了一个矩阵。bvp4c会评估这些雅可比人吗？如果是这样，那么如果任何g不包含ya或yb，那么它们将始终是单数。有没有人能给我一些解释
浏览 16提问于2020-02-14得票数 0
我读过关于非线性算术和未解释函数的文章。我仍然是SMT世界的新手，所以如果我没有使用正确的词汇表，或者这是一个糟糕的问题，我很抱歉。
对于下面的代码，有一个不相关的顶级断言( (assert (> i 10)) )上面的断言放在堆栈上。然而，Z3用Reals (第一次推到第一次pop)返回unsat。我认为这与Z3试图使用Int求解器有关，因为第一个断言是在Int上，而Z3将e1分配给(/ 1.0 2.0)，这是一个没有Int表示的数字，因为有约束(assert (< e3 1)) (如果去掉这个约束，它就能工作)。使用(check-sat-using qfnra-nlsat)解决
浏览 3提问于2013-10-29得票数 1
回答已采纳
我有一个问题，可以简化为下面的陈述(实际的方程和函数要复杂得多)：
系统：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
1回答
我正在尝试解决设施选址问题。这是我的代码： set S;
param prod{i in S};
param distri{i in S};
param fixed{i in S};
param cap{i in S};
param demand{i in S};
var x{i in S, j in S}, >= 0;
var y{i in S}, binary;
minimize obj :
sum{i in S} fixed[i]*y[i] +
sum{i in S, j in S} x[i,j]*(prod[i] + distri[i])*y[i];
s.t.
浏览 64提问于2020-09-10得票数 0
回答已采纳
如何在java中实现带非线性约束的非线性优化？我目前正在使用org.apache.commons.math3.optim.nonlinear.scalar.noderiv，，并且我已经读到，默认情况下没有一个优化器(例如我目前正在使用的优化器，SimplexOptimizer)会接受约束，而是必须通过实现MultivariateFunctionPenaltyAdapter或MultivariateFunctionMappingAdapter类将受约束的参数映射到无约束的参数。然而，据我所知，即使使用这些包装器，仍然只能实现线性或“简单”约束。我想知道是否有办法包括非线性不等式约束？
例如，假设
浏览 0提问于2014-12-02得票数 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
回答已采纳
1回答
如何计算3D中倾斜的向量的上升？
、、、、
我有一个数据，它显示了3D空间中矢量的起点和终点坐标，这些矢量的方向是围绕3折螺旋轴：
x y z
0 38.522003 5.600998 129.203995 # start of v1
1 23.854996 66.576996 112.487000 # end of v1
2 4.417000 40.182999 121.309998 # start of v2
3 65.761993 27.550995 104.285004 # end of v2
4 50.272003
浏览 5提问于2018-02-23得票数 0
有人能向我展示Microsoft Solver foundation 3.0中约束非线性优化的示例或评论吗？与Matlab的fmincon相比如何？或者有没有更好的用于约束非线性优化的.net库？谢谢,
浏览 8提问于2011-06-12得票数 11
回答已采纳
1回答
SymPy符号数学库()为求解奇异方程和方程组提供了一个求解模块。下面是一个用法的例子：
import sympy as sym
x,y,z = sym.symbols('x,y,z')
c1 = sym.Symbol('c1')
eq1 = sym.Eq(2*x**2+y+z,1)
eq2 = sym.Eq(x+2*y+z,c1)
eq3 = sym.Eq(-2*x+y,-z)
result = sym.solve([eq1,eq2,eq3],(x,y,z))
print(result)
在尝试了一系列自定义方法来解决以下非线性系统后，我使用了solve()，
浏览 0提问于2020-12-04得票数 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
1回答
z3实数求幂
、、、
我遇到了一种奇怪的情况，z3py为逻辑上相同的问题生成了两个不同的答案。
版本1：
>>> import z3
>>> r, r2, q = z3.Reals('r r2 q')
>>> s = z3.Solver()
>>> s.add(r > 2, r2 == r, q == r2 ** z3.RealVal(0.5))
>>> s.check()
unknown
版本2
>>> import z3
>>> r, r2, q = z3.Real
浏览 3提问于2016-04-08得票数 0
1回答
或工具制定收入
、、、
我正在寻求关于如何在OR工具中解决一个问题的建议。 背景是，我是一个游戏所有者，在轮班时，我可以设置不同数量的游戏站点。假设我有2个站点和5个客户，那么假设分布最均匀，一个站点有3个客户，另一个站点有2个客户。 当一个车站的顾客数量不同时，游戏速度会降低。为了最大化游戏速度，如果我有5个客户，最好的解决方案是有5个工作站，这样每个客户每个都有1个工作站。但这增加了运营成本。 在这种情况下，我如何表示总游戏速度，其中涉及到客户到车站的均匀分布，并且游戏速度取决于车站的客户数量？
浏览 14提问于2021-09-21得票数 1
1回答
我正在寻找有关解决程序APOPT的文档。我在论文中使用了Gekko，因此我不得不引用并解释代码的每个部分。是否有视频/文章/网站解释解题器的工作原理？它是否使用分支和绑定来解决优化问题？APOPT总是被描述为非线性MIP的求解器？它也适用于线性MIP吗？
浏览 4提问于2020-10-13得票数 3
回答已采纳
我正在尝试学习如何使用Z3未解释的函数来显示n=2k+1 => log(m) + k*log(m*m) == n*log(m)。为了做到这一点，我使用了类似如下的方法
mylog = Function('mylog', IntSort(), IntSort())
mylog_rule1 = mylog(x*y) == mylog(x) + mylog(y)
mylog_rule2 = mylog(x**y) == y*mylog(x)
#mylog_rule3 = y*mylog(x) == mylog(x**y) #is this rule needed ?
浏览 2提问于2013-03-01得票数 0
我想知道我们如何在OpenModelica中定义求解非线性函数的迭代极限？是否可以记录每个时间步的平均迭代次数或每个时间步的迭代次数？
浏览 16提问于2021-01-12得票数 1
3回答
我完全意识到凸优化包，就像线性代数包一样，应该是你使用的东西，而不是实现的东西。然而，出于纯粹的教育目的--有没有好的资源--关于如何实现凸优化包的链接/书籍？(就像具有二次约束的二次程序一样？)
谢谢!
浏览 1提问于2010-07-08得票数 3
回答已采纳
1回答
试着解那个方程中的B
、、、
我正在试着解决那个长表达式中的变量B。它给了我一个错误，它不能转换为浮点数。我真的不确定这到底是什么意思，以及我需要做什么来修复它并获得B的值。我对Python非常陌生，所以如果有任何帮助，我将不胜感激。
代码：
import matplotlib.pyplot as plt
import numpy as np
import math
from sympy.solvers import solve
from sympy import Symbol
W = 1*(10**(-4))
mu = 300
Cox = 6.906*(10**-6)
Leff = 1*(10**(-4))
V
浏览 1提问于2018-05-25得票数 0
3回答
在不求助于渐近表示法的情况下，单调的步长计数是获得算法时间复杂度的唯一方法吗？如果不计算每一行代码的步数，我们能得到任何程序的大O表示吗？
细节:尝试找出几种数值分析算法的复杂性，以确定哪种算法最适合解决特定问题。例如，在求解方程的Regula-Falsi或Newton-Rhapson方法中，意图是评估每种方法的确切复杂性，然后决定(放置'n‘的值或任何参数)哪种方法不太复杂。
浏览 1提问于2010-07-28得票数 5
回答已采纳
我试着用solve(g)找到满足公式g = And(ForAll([i, j, k], Implies(And(k <= 0, A * i + B * j + C * k + D <= 0), k + i - j <= 0)),Not(And(A==0,B==0,C==0)))的A，B，C，D的值。这有许多可能的解决方案，一种是A=1,B=-1,C=D=0，但Z3似乎无法做到这一点(solve(g)不会终止)。
Z3能做这种非线性(但简单)的公式吗？也许我可以为此指定一些量化宽松策略或其他什么？
浏览 3提问于2013-02-19得票数 2
回答已采纳
以下代码
(declare-fun f (Int) Real)
(declare-fun g (Int) Real)
(declare-const x Int)
(declare-const y Int)
(assert (= (* (f x) (g y)) 1.0))
(check-sat)
(get-model)
返回“未知”，即使有一个明显的解决方案。消除f和g的参数(有效地使它们成为常量？)有预期任务的"sat“的结果。我想，我的问题是:没有解释函数的算术有什么特别之处？
顺便说一句，用+替换*也会导致"sat"，所以问题不是关于未解释的函数本身，而是关于它们是
浏览 0提问于2017-03-31得票数 1
回答已采纳
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
回答已采纳
我试图编写一些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
2回答
这是一个实际的组合优化问题。
我们被赋予了某一产品的一大套价值主张。价值主张有不同的类型，但每种类型都是独立的，并为整个产品增加了同等的好处。在构建产品时，我们可以包含每种类型的任何非负整数个“单元”。然而，在添加了某一类型的第一个单元后，该类型的附加单元的边际效益不断减少。事实上，新单位的边际收益是添加新单位后该类型单位的数量的倒数。我们的产品必须至少有一个某种类型的单位，并且由于这个要求，我们必须对整体价值进行一些小的修正。
假设T[]是一个数组，表示产品的特定生产运行中每种类型的数量。则整体值V由(伪代码)给出：
V = 1
For Each t in T
V = V * (t
浏览 4提问于2011-06-14得票数 10
回答已采纳
我对与中的多学科组连接的求解器感到困惑。
他们本质上构造的是什么矩阵？它们会聚到什么地方？很可能我遗漏了一些非常明显的东西。有没有这方面的文献。
浏览 0提问于2018-09-28得票数 2
3回答
具有任意约束的R中的优化
、、、、
我已经在Excel中完成了，但需要在R中运行一个适当的模拟。
我需要最小化函数F(x) (x是向量)，同时具有sum(x)=1的约束，x中的所有值都是[0,1]和另一个函数G(x) > G_0。
我已经和optim和constrOptim试过了。他们都没有给你这个选择。
浏览 3提问于2014-02-06得票数 1
回答已采纳
假设这里是非线性实算法上的一组约束
pred1 = (> (- (* (- v2_x v0_x) (- v1_y v0_y)) (* (- v2_y v0_y) (- v1_x v0_x))) 0)
pred2 = (> (- (* (- v1_x v0_x) (- v2_y v0_y)) (* (- v1_y v0_y) (- v2_x v0_x))) 0)
事实上，如果我们这么做
Z3_solver_assert(ctx,solver,pred1);
Z3_solver_assert(ctx,solver,pred2);
b = Z3_solver_check(ctx, sol
浏览 3提问于2012-11-27得票数 2
我使用的是z3 C++ API，在非线性整数算法方面存在一些问题，z3似乎对幂和乘法的处理方式不同。
x*x*x > y*y是sat，x^3 > y^2是unknown
#include <z3++.h>
#include <iostream>
using namespace std;
int main() {
z3::context c;
auto x = c.int_const("x");
auto y = c.int_const("y");
z3::solver sol(c);
浏览 3提问于2015-10-29得票数 1
回答已采纳
1回答
这是我关于stackoverflow的第一个问题。抱歉，如果它有点笨拙。 我正在尝试解决描述层状土壤中的梁的静态平衡。就作用在梁上的力而言，有一些外力，还有土对外力的反作用力。显然，行动和反应必须平衡。反作用力是来自不同土层的力的总和。土的响应是位移的非线性函数。基本上，这个特殊的平衡是一个标量，因为它就是反应=作用的总和。 另一方面，梁内的内力也必须与作用在所讨论的梁单元上的局部力相平衡。假设我们有40个这样的元素。因此，我面临的问题概括如下： 我需要找到F1(x1，x2，...x40)的根，用x1，x2...x40描述桩在土壤中的位移。我传递并初步猜测了位移，Fsolve返回了一个向量长
浏览 43提问于2020-05-03得票数 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
8回答
我对机器学习和神经网络很陌生。我知道如何建立一个非线性分类模型，但我目前的问题有一个连续的输出。我一直在寻找关于神经网络回归的信息，但我所遇到的只是关于线性回归的信息--没有关于非线性情况的信息。这很奇怪，因为为什么会有人用神经网络来解决简单的线性回归呢？这不就像用核弹杀死苍蝇吗？
所以我的问题是:是什么使神经网络非线性？(隐藏层？非线性激活函数？)或者我对“线性”这个词有一个完全错误的理解--线性回归神经网络能够准确地建模出比y=aX+b更复杂的数据集吗？“线性”一词是否与“逻辑”相反？
(我计划使用TensorFlow，但TensorFlow线性模型教程以二进制分类问题为例，因此这对我也没
浏览 12提问于2016-12-20得票数 37
回答已采纳
1回答
我刚开始使用openMDAO，所以我想这应该是一个简单的问题。 我有一个复杂的组(具有显式的组件和循环)，我可以求解一组输入变量： y1，y2，...，yi =组(x1，x2，...，xj) 我现在尝试做的是匹配两个输出(y1_target和y2_target)，改变来自组(x1，x2)的两个输入，即从组中添加两个等式，例如， y1_target - y1 =0 y2_target - y2 =0 我知道这应该通过两个隐式组件来完成，但是我如何强制只更改x1和x2呢？ 谢谢,
浏览 14提问于2021-04-05得票数 0
1回答
定理证明工具z3花了很多时间来解决一个公式，我相信它应该能够很容易地处理。为了更好地理解这一点，并可能优化我对z3的输入，我想看看z3在求解过程中生成的内部约束。从命令行使用z3时，如何打印z3为其后端求解器生成的公式？
浏览 1提问于2012-10-28得票数 7
回答已采纳
1回答
我正在尝试使用Pyomo解决一个优化问题。为此，我需要声明一个二维向量和一个三维向量。这两个向量只能存储0,1。
Sri和Xrij
R=3 V=8 1<=i，j<=V 1<=r<=R
我尝试在pyomo中使用range
model.IDXV = range(v+1)
model.IDXR = range(r+1)
model.x=Var(model.IDXR,model.IDXV,model.IDXV,within=Binary,initialize=0)
model.s=Var(model.IDXR,model.IDXV,within=Binary,initiali
浏览 83提问于2018-08-12得票数 0
回答已采纳
1回答
我们遇到了性能问题，我认为这是Z3中处理非线性算术的部分。这里有一个简单的具体的Boogie示例，当使用Z3 (版本4.1)验证时，需要很长时间(大约3分钟)才能完成。
const D: int;
function f(n: int) returns (int) { n * D }
procedure test() returns ()
{
var a, b, c: int;
var M: [int]int;
var t: int;
assume 0 < a && 1000 * a < f(1);
assume 0 < c &&
浏览 2提问于2012-09-20得票数 5
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
云直播活动推荐
腾讯云开发者
