使用带不等式方程的fsolve求解方程
是一种数值计算方法，用于求解包含不等式约束的方程。fsolve是MATLAB中的一个函数，可以通过迭代的方式求解非线性方程组。
在使用fsolve求解带不等式方程的过程中，首先需要定义一个函数，该函数的输入参数是待求解的变量，输出是方程的值。然后，通过调用fsolve函数，将定义的函数和初始猜测值作为参数传入，即可得到方程的解。
在实际应用中，带不等式方程的求解可以应用于各种领域，例如优化问题、经济学模型、物理学模型等。通过求解方程，可以得到满足约束条件的最优解或近似解。
对于腾讯云相关产品和产品介绍链接地址，由于要求不能提及具体品牌商，无法给出具体的产品推荐。但是，腾讯云作为一家知名的云计算服务提供商，提供了丰富的云计算产品和解决方案，涵盖了计算、存储、网络、人工智能等多个领域。可以通过访问腾讯云官方网站，了解更多关于腾讯云的产品和服务信息。
总结：使用带不等式方程的fsolve求解方程是一种数值计算方法，适用于求解包含不等式约束的方程。腾讯云作为一家知名的云计算服务提供商，提供了丰富的云计算产品和解决方案，可以满足各种应用场景的需求。
相关·内容
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
1回答
我有一个非线性方程:我在下图中画出了这个方程的左边和右边。显然，有两种解决方案。我已经尝试了scipy.optimize.fsolve，并战略性地猜测了起点，以找到两个解决方案，但这是相当不健壮的。有没有一个Python求解器，可以找到这个方程的所有解，例如，在某个有限的x值集合中？
浏览 5提问于2021-11-10得票数 0
1回答
我想要找到一个方程的数值解，用Python表示为无限和。使用fsolve()和sympy的简化示例如下：
from scipy.optimize import fsolve
import math
from sympy import *
i = symbols('i', integer=True)
def f(x):
return Sum(x**i, (i,0, oo)).evalf(10)-1
print fsolve(f, 0.5)
我收到了一个错误：
AttributeError: 'list' object has no attri
浏览 4提问于2016-03-29得票数 1
回答已采纳
我想知道fsolve in scipy是如何工作的。到目前为止，我所发现的是它是一个数值求解器，可以找到非线性方程组的根。但是我找不到它使用什么方法来找到根。有人知道树根是怎么找到的吗？
浏览 6提问于2019-07-12得票数 0
1回答
我正在尝试使用优化子模块中的SciPy最小化函数来优化x-y平面中距原点0.15半径范围内的曲面方程。 ans = minimize(minimum_objective_3d, (.10,.10), args=(args), bounds=[(-.150,.150),(-.150,.150)]) 我知道我正在使用的边界参数指定了操作的矩形边界。我的问题是。有没有办法将这些边界指定为以原点为中心的圆形半径？ 提前感谢你的帮助。
浏览 4提问于2021-03-28得票数 0
2回答
MATLAB有两种求解非线性方程的方法：
：求解一个单一的非线性方程
：解一个非线性方程组
因此，可以使用以下方法求解n非线性独立方程组：
用一个循环用fzero分别求解方程
用一个循环用fsolve分别求解方程
使用fsolve一起解决它们
我的直觉是：
当复杂度(梯度计算)为0(n^2)时，循环方法比单个大n系统更快。
对于小型n，循环可能比较慢，因为一个循环在MATLAB中有很高的开销，并且可能有一些恒定的启动时间。
fzero比fsolve更快，因为它是专为一个非线性方程而做的。
问题：，对于一个n非线性独立方程组，最快的方法是什
浏览 24提问于2017-07-10得票数 3
回答已采纳
1回答
解复矩阵方程
、、、
我有一个涉及矩阵的复杂方程式：
R = expm(X)*A + (expm(X)-I)*inv(X)*B*U;
其中R、B和U是已知矩阵。
I是一个单位矩阵。
我需要解决X的问题。有没有办法在MATLAB中解决这个问题？
浏览 0提问于2019-04-28得票数 1
我有一个应用程序，其中我需要使用迭代技术来求解变量。基本上，我需要猜测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
2回答
我有一个方程'a* x +logx-b=0，(a和b是常量)‘，我想要解x，问题是我有许多常量a(相应地有很多b)。如何使用python求解此方程？
浏览 2提问于2014-11-04得票数 0
我有，我想解决X(t)，Y(t)，Z(t)，希望用我熟悉的Octave/Matlab，但我不介意用任何其他必要的方法来解决它。
现在，Fsolve对于正则方程组是有用的，而Ode45，Lsode对于微分方程是有用的。但是，这个特殊的系统呢？请注意，底部的微分方程不仅包含Y，还包含X和Z，它们依赖于上面的两个非微分方程。
老实说，我不太确定如何使用基本的代码来解决这个系统，在花了一些时间思考之后，我决定寻求帮助。如果可能的话，我真的很感谢任何能有效地解决这个问题的指导。现在几乎所有的回复对我来说都是有用的。
浏览 9提问于2020-02-15得票数 0
回答已采纳
我想开发一个用于解线性和非线性方程组的桌面应用程序。我正在考虑在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回答
具有多重方程的枕方程组
、、、
你好，当我有一个方程组时，我在使用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回答
用起点解决问题
、、、
我正在使用fsolve来解决一个非线性方程。我的问题是，根据起点的不同，解决方案会发生变化，我不确定我找到的解决方案是最合理的。这是代码
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import fsolve, brentq,newton
A = np.arange(0.05,0.95,0.01)
PHI = np.deg2rad(np.arange(0,90,1))
def f(b):
return np.angle((1+3*a**4-3*a**2)+(a**4-a**6)*
浏览 1提问于2014-02-07得票数 0
2回答
我有一个由三个方程组成的系统，我想通过MATLAB来解决，我有点困惑如何去做。
我有三个方程式：
A = R*P1
B = R*P2
C = R*P3
A, B, C、P1, P2和P3是3x1矩阵，而R是3x3矩阵。对于所有三个方程，R都是相同的。
我需要找到R，并得到A, B, C和P's。
我想使用fsolve，但是当变量是矩阵时，fsolve似乎不能工作。您建议使用的替代方法是什么？
只是编了几个要处理的数字：
P1 = [1;1;1];
P2 = [2;3;4];
P3 = [5;4;3];
R = [2 4 5; 1 5 4; 1 2 3];
这意味着：
A = [11;
浏览 3提问于2014-12-14得票数 1
回答已采纳
我在matlab中见过几个例子，但似乎找不到任何显示如何以矩阵形式传递参数的例子。
这是我的密码。
[A,b] = equationsToMatrix(eq1,eq2)
X0 = [0 0]
fsolve([A,b], X0)
这是输出
eq1 = - sx - sy/2 == 5
eq2 = - (3*sx)/2 - (3*sy)/2 == 9
A =
[ -1, -1/2]
[ -3/2, -3/2]
b =
5
9
X0 =
0 0
Error using lsqfcnchk (line 109)
If FUN is a MATLAB obj
浏览 7提问于2013-10-18得票数 1
回答已采纳
给出了这个例子：
from sympy.core.symbol import symbols
from sympy.solvers.solveset import nonlinsolve
x, y, z = symbols('x, y, z', real=True)
nonlinsolve([x*y - 1, 4*x**2 + y**2 - 5], [x, y])
{(-1, -1), (-1/2, -2), (1/2, 2), (1, 1)}
但即使是在他们网站上的实时外壳中，这也会导致一个错误：
>>> from sympy.solvers.solvese
浏览 2提问于2016-11-24得票数 8
回答已采纳
我用Matlab试图用两个变量来求解两个方程。
定义了两个函数，f2(n_1，n_2)，f3(n_1，n_2)，它们都依赖于f1(n_1，n_2)，然后定义了包含这两个函数的矢量函数G(n_1，n_2)。
后来，我定义了一个想要的陈述点，并试图解决。但是，当运行代码时，它会引发一个错误，我还不完全理解。
上述内容显示在下面的代码中：
守则：
clear, close all; clc
%Const
N0=25;
G1=1;G2=1;
a1=6;a2=3;
k1=1;k2=4;
%main
syms n_1 n_2
X_0=[-5;5];
f1=N0-a1.*n_1-a2.*n_2;
f2=f1
浏览 0提问于2018-11-14得票数 0
我有一套方程式..。
dF(a,b,c)/da = 0;
dF(a,b,c)/db = 0;
dF(a,b,c)/dc = 0;
其中a、b、c是未知变量常量，dF/d*是变量的匿名函数。我必须在一个优化问题中求解a、b和c。当系统简化为一个方程时，我使用Matlab的fzero对变量进行求解，并且工作正常。例如
var_a = fzero(@(a) dF(a)/da,0);
在注意到fzero和fsolve在某些情况下给出了截然不同的答案后，我进行了一些搜索。根据我的，fzero只适用于一个单一变量的单一方程？因此，我想选择一个最合适的方法。我以前使用过Matlab的solve，但我相信它只
浏览 1提问于2016-04-23得票数 1
我正在使用python2.7和同情。我想要解一个方程，但考虑到变量的值是受约束的(并不是整数变量的所有值)，例如-25 <= a <= 25。
我该怎么做？
我的代码
from sympy import *
from sympy import sin
from sympy.solvers import solve
a, b, c, d = symbols('a b c d')
print solve(cos(a)*sin(b) - sin(a)*cos(b)*cos(d - c), a, b, c, d)
浏览 4提问于2015-07-03得票数 0
我正在用ode15在MATLAB中求解一个刚性偏微分方程，它通常会根据初始条件而冻结。我实际上从来没有收到错误，只是它不会结束，即使在10个小时后，它应该需要大约30秒才能运行。我正在尝试不同的空间和时间节点间隔，但这很难，因为我得不到反馈。
对于fsolve，是否存在某种等价于诊断的功能？stats没有什么用处，因为它只在fsolve完成后才显示输出。
浏览 0提问于2013-06-13得票数 0
1回答
我正试图用fsolve来解一个非线性方程组，比如说
F(x; lambda ) = 0，其中lambda是参数的向量，x是我想要求解的向量。
我使用的是Matlab的解决方案。
我有两个参数lambda的值，我想解决这个系统。对于lambda的一个值，我得到了一个解决方案，这似乎没问题。
对于lambda的另一个值，我又得到了一个解(matlab退出时标记为1。但是我知道这不是一个实际的解，例如，我知道x的某些维数必须是相等的，这不是我从For得到的解的情况。)
我已经尝试过信任域和levenberg-marquardt算法，但是我没有得到更好的结果。(显式地强制这些x是相同的，似乎仍然给出了
浏览 4提问于2013-03-19得票数 2
1回答
我试图在Python中求解一个二次方程，这是一个矩阵方程。我正在寻找一个矩阵Ax，它是2x2，这样它就能满足
M_1 Ax^2 - M_2 A_x - M_3 = 0
其中M_1、M_2和M_3是已知的2x2 numpy常量数组，Ax^2本身就是Ax的矩阵乘法。
到目前为止我的代码是：
import numpy as np
from scipy.optimize import fsolve
# Calibration
gamma = 0.36
beta = 0.98
b_bar = 0.1
psi = 0.001
g = 0.66
rho_x = 0.11
rho_z = 0.95
c_
浏览 4提问于2021-01-26得票数 1
1回答
我试图用数值方法求出X^2+X+C=0的解，其中C是矩阵C=[-6,-5;0,-6]和0=[0,0;0,0]，这是一个二次方程，变量是2x2矩阵。
所以我编写了下面的matlab命令
C=[-6,-5;0,-6]
[X1,F,e_flag]=fsolve('X^2+X+C',[1,1;1,1])
其中[1,1;1,1]是我们的初步估计，或X0。
我得到了以下错误
“在内联表达式==>、X^2+X+C、未定义函数或变量'X‘中使用内联表达式(第15行)出错。
内联/feval中的错误(第34行) INLINE_OUT_ = inlineeval(INLINE_INP
浏览 1提问于2014-12-26得票数 1
回答已采纳
1回答
假设我们有三个方程：
eq1 = x1 + (x1 - x2) * t - X == 0;
eq2 = z1 + (z1 - z2) * t - Z == 0;
eq3 = ((X-x1)/a)^2 + ((Z-z1)/b)^2 - 1 == 0;
而已知的六个变量是：
a = 42 ;
b = 12 ;
x1 = 316190;
z1 = 234070;
x2 = 316190;
z2 = 234070;
因此，我们正在寻找三个未知的变量：
X , Z and t
我写了两种方法来解决它。但是，由于我需要对570万个数据运行这些代码
浏览 3提问于2016-06-13得票数 1
在Matlab中，我很难用函数来求解一组5个方程。
以下是5个方程式：
Y=a+d+e
Y+x=c+d+ 2e
2x = 4a + 2b + 2c
k1 = (d * b^3 / (a * c) )* ((P/Pref)/(a+b+c+d+e))^2
k2 = b_e/(d_c)
y、x、k1、k2、P、Pref都是我设置的参数，但希望将它们保留在函数中，以便在代码中快速更改它们以找到新的答案。a，b，c，d，e是我想要解的变量(它们是反应平衡方程的组成)。
我试图对函数中的参数进行硬编码，但这是行不通的。我只是不知道该怎么办。我所做的每一件事都会产生新的错误。最常见的是数据类型必须是"
浏览 5提问于2015-11-25得票数 0
1回答
当我使用fsolve时，程序会给我发送一个错误，但是我不知道为什么。我应该去找x(2) = 2.96e+13和x(1)= 2.72e+22。我试着用这个来解我的非线性方程:我有两个方程，两个不知道：
function y = f(x)
global x_rayon_droite
global z1
global r_droite
global x_rayon_gauche
global r_gauche
global y_rayon_droite
global y_rayon_gauche
x_rayon_droite=406
z1=2.72e+22
r_droite = 556.77
x
浏览 4提问于2016-12-02得票数 0
回答已采纳
1回答
我正在尝试用Matlab求解一个由12个方程组成的系统。因为我对变量的最小值和最大值有约束，所以我使用lsqnonlin而不是fsolve。但是，我希望优化器在输出(与每个方程所在点的平方偏差之和)足够接近于零时停止。有没有办法指定这样的停止标准？ 标准的停止条件是比较输出值与前一次迭代相比的变化，但这与我的相关性较小。
浏览 32提问于2021-01-05得票数 0
1回答
我想用fsolve在Matlab中求解一个非线性方程组，但我也要区分两个变量的函数。以下是步骤中的问题：
步骤1:我定义了非线性函数系统F。
function F = root2d(x)
F(1) = (exp(-x(1)+2)/(1+exp(-x(1)+2)+exp(-x(2)+1)));
F(2) = (exp(-x(2)+1)/(1+exp(-x(1)+2)+exp(-x(2)+1)));
end
步骤2:我想找出F(1)的导数与x(1)的导数和F(2)的导数与x(2)的导数，以便：
function F = root2d(x)
F(1) = (exp(-x(1)+2)/(1+
浏览 12提问于2022-09-30得票数 0
我有两个简单的方程式：
k = 2.013*h^0.4917 and h = 3.57*k^0.4917
当k= 5.77，h= 8.47时，可以解析地求解这些方程。我试图在Python中使用from来解决这个问题，我遵循了以下方法：
下面是我的代码：
from scipy.optimize import fsolve
def equations(x):
return[x[0] - 2.013*x[1]**0.4917,x[1] - 3.57*x[0]**0.4917]
root = fsolve(equations, [1, 1])
结果是
<ipython-input-
浏览 4提问于2021-01-15得票数 0
回答已采纳
1回答
我很难理解如何求解大O-表示法方程，比如这个：
f(n) = 10n + 5 and g(n) = n
Show that f(n) is O(g(n))
从这里拍摄：
上面写着To show f(n) is O(g(n)) we must show constants c and k such that f(n) <= cg(n) for all n >=k
我得到:从k开始，f(n)的增长速度不会超过g(n)乘以任意常数。
我不明白的是：
We are allowed to choose c and k to be integers we want as long as th
浏览 4提问于2015-07-03得票数 1
回答已采纳
我想做一个功能，以计算真正的Elo性能评级与游戏的国际象棋。
它将纳入：
game_results向量，例如表示玩家赢得第一局、输了第二局和第三局平局的向量。由此我们可以得到另一个变量的分数，这里是1.5，和另一个变量num_games，这里是3。
opponent_ratings的一个矢量，例如12001401350，它指示第一个对手被评级为1200，第二个140分和第三个1350。
我知道，真正的性能等级是x，其方程式如下：
for i = 1:num_games
score = score + (1/(1+10^((opp_ratings(i)-x
浏览 2提问于2018-10-03得票数 2
回答已采纳
我试着用解析法和狗腿法来求解一个二阶非线性方程组。我的目标函数和它的雅可比矩阵是这样的
function [F jacF]= objective(x)
F(:,1) = ((((x(:,2)./10).*k).*(x(:,1)./100)).^2).*(rZ - Rs) +(( Cmax .* ( x(:,1)./100 ) ).^2).*( w.^2.*(rZ - Rs) ) - (((x(:,2)./10).*k).*(x(:,1)./100));
F(:,2) = (x(:,2).*k).^2.*(iZ - w.*Ls) + (x(:,2).*k).^2.*x(:,
浏览 0提问于2019-09-12得票数 0
1回答
R二次规划
、、、
我有一个问题想在R中解决。
我看到，我可以使用包lsei中的函数limSolve来最小化以矩阵形式写成的线性方程组( Ax=b )，但条件是等式约束Ex=f和不等式约束Gx>=h。
然而，不是线性方程组，我现在有一个二次方程组，它可以写成t(x)Ax=b。
我看到了二次型情况下的包quadprog，但是它似乎不允许一组二次方程，只有一个方程。
有人知道我能用什么来最小化等式和不等式约束下的二次方程组吗？
这是我的例子。我试图结合三个概率P(A)，P(B)，P(C) -这就创建了7个片段v1到v7，其中v1是P(A solus)等等。v4是P((A和B)不是C)等。v7为P(A、B、C)。
浏览 2提问于2016-01-13得票数 0
我试图解一个线性方程组(11个不同的方程/限制上的5个变量)。
我的问题是，我正在寻找的解决方案不可能仅仅是，因为我试图优化预算，它必须是>=或<=，但是我不知道如何用numpy来写它。
方程式如下(只需简化几个)：
X/.07 >= 1728Y/.02 >= 4431Z/.05 >= 5409A/.005 >= 7138<code>H 113</code>B/.01 >= 4897<code>H 214</code><code>H 115<//code>A+B+X+Y+Z >
浏览 4提问于2019-12-18得票数 0
我试着用西佩的解来找出两个非线性方程组的答案。这两个方程式是：
f1 = math.log(x) + 1. - ((1. + (m - 1)*x) / m) + chi * (1 - x)**2
f2 = math.log(1 - x) - (m - 1)*x + chi*m*x**2
在这种情况下，m和chi是常量。其基本目标是寻找同时满足x和f2(x) = f2(y)的f1(x) = f1(y)、f2(x) = f2(y)。我知道x的最初猜测，y是0.3和0.99。下面是我的密码。
from scipy.optimize import fsolve
import math
# some
浏览 2提问于2020-01-06得票数 0
回答已采纳
我有一个在MATLAB中生成向量的代码，例如下面是一个三分量向量(n=3)：
a1_1 - sin((17*a1_1)/60 + a2_1/8 + a3_1/40 - 0.153233)
(15*a1_1)/16 + a2_1/4 + a3_1/32 - sin((17*a1_1)/60 + a2_1/8 + a3_1/40 - 0.0282326)
(3*a1_1)/4 + a2_1/2 + a3_1/8 - sin((17*a1_1)/60 + a2_1/8 + a3_1/40 + 0.846767)
如你所见，每个分量都是一个非线性方程。向量的三个分量形成了一个由三个非线性方程组成的系统
浏览 0提问于2015-08-03得票数 0
1回答
在python中，我希望找到以下形式的方程的根：
-x*log(x) + (1- x) *log(n) - (1-x)*log(1 -x)-k=0
其中n和k是将要指定的参数。
对根的一个额外的约束是x >= (1-x)/n，所以为了它的价值，我将过滤掉不满足这个条件的根。
我的第一个尝试是使用scipy.optimize.fsolve (请注意，我只是将k和n分别设置为0和1)：
def f(x):
return -x*log(x) +
浏览 4提问于2013-12-03得票数 1
我正试图用Python来求解一个非线性方程组。方程的形式如下：
(1) x^2 + y^2 =a
(2)(x-b)^2 + y^2 =c
其中x和y是变量，a，b，c是参数。我想要一个函数，它可以把参数a，b，c传递给它，它会返回x和y的值。
我现在拥有的是
from scipy.optimize import fsolve
def equation(var, *data):
a,b,c = data
x,y = var
eq1 = x**2 + y**2 - a**2
eq2 = (x - b)**2 + y**2 - c**2
return [e
浏览 51提问于2022-06-08得票数 0
回答已采纳
我用MATLAB的fsolve函数来求解非线性方程组。我有两个双变量(x,y);的非线性方程
我试图为这两个变量找到所有可能的根。我注意到fsolve只给出一个根。如何才能得到方程的全根？
我的代码如下：
function F = fun(guess)
x = guess(1);
y = guess(2);
F = [2*x -y - exp(-x));
-x + 2*y - exp(-y) ];
end
调用该函数：
guess = [-5 -5]
fsolve(@fun,guess);
浏览 7提问于2014-01-13得票数 1
我在试着解两个耦合的代数方程
f1(x,y) = 0;
f2(x,y) = 0;
函数f1和f2的典型数量级是10^42。我运行了matlab代码，但它显示没有找到解决方案。我认为问题是因为涉及的规模非常大。重新调整整个等式是相当乏味的。我想在delta(f)/f < epsilon(say 1e-6)时停止根查找函数(fsolve)。如何在matlab中实现此条件？对缩放问题的任何替代解决方案也是受欢迎的。
浏览 4提问于2015-04-05得票数 0
3回答
? 我想解一个基于洛伦兹方程的微分方程： def f(xyz, t, rho, sigma, beta):
x, y, z = xyz
return [sigma * (y - x),
x * (rho - z) - y,
x * y - beta * z] 我写了以下内容： def f(xyz, t, rho, sigma, beta):
x, y, z = xyz
return [sigma * y(t).diff(t) + sigma * x + beta * y -77,
浏览 60提问于2021-07-27得票数 0
回答已采纳
所以我有一个复杂的方程，我需要解。我认为最终x应该是1E22。但是这个代码的问题是它使我的整个系统崩溃了。有解决办法吗？我尝试了scipy.optimize.root，但它并没有真正解决这个数量级的任何问题(它给出了最终的答案，作为最初的猜测，无需任何迭代)。
from scipy.optimize import fsolve
import math
import mpmath
import scipy
import sympy
from sympy.solvers import solve
from sympy import Symbol
from sympy import sqrt,exp
浏览 0提问于2018-10-13得票数 1
回答已采纳
我正在尝试用Haskell重写一个相当简单的MATLAB脚本。我正在寻找一个等同于MATLAB的fsolve (数值求根程序)和ode15s (一种常微分方程求解器)的程序。在Hackage上，似乎有一个MATLAB接口包。这是要走的路吗(在速度、简单性等方面)或者有没有其他更合适的包呢？
谢谢，本
浏览 0提问于2010-07-19得票数 5
回答已采纳
嗨，我被要求使用MATLAB中的fsolve命令来解决SIR模型，并向后返回Euler 3点。我真的不知道该怎么做，请帮助我。这就是我到目前为止所拥有的。我为3BDF格式创建了一个函数，但我不确定如何继续求解非线性常微分方程组。SIR模型如下所示
并且3BDF方案被表示为
clc
clear all
gamma=1/7;
beta=1/3;
ode1= @(R,S,I) -(beta*I*S)/(S+I+R);
ode2= @(R,S,I) (beta*I*S)/(S+I+R)-I*gamma;
ode3= @(I) gamma*I;
f(t,[S,I,R]) = [-(beta*I
浏览 2提问于2020-04-15得票数 1
在讨论问题之前，我使用相同的问题搜索了堆栈溢出上的相关线程：
据我所知，关于这个错误，
提高TypeError(味精)
TypeError: the：“func”参数'fsolve_function‘的输入和输出形状不匹配
问题是输入和输出的形状是不一样的。
在下面的代码示例中，我有以下内容：
输入，initialGuess ( fsolve函数在scipy.optimize中使用的起始估计)。输入，initialGuess有坐标x，y和z的3个起始估计，因此我期望我的起始输入估计总是有三个输入。
out (非线性联立方程)。在这个例子中
浏览 1提问于2015-06-04得票数 6
回答已采纳
1回答
我正在努力学习如何解决线性规划问题，我想使用ojAlgo LinearSolver。
I使用SolverTask求解等式约束，得到一个单一的可行点。很好！INFEASIBLE.Then！，然后我用它失败的等式约束解决了一个相同的问题，说问题是LinearSolver --我添加了不等式约束，目的是不影响解，它也说它不可行。但是现在我也得到了比预期更大的结果向量。
我的问题：
如何用仅带等式约束的SolverTask?How ()重新构造问题以解决问题，同时获得与SolverTask?What相同的解决方案
我使用ojAlgo 49.2.1，在50.0.1中得到了相同的结果。
public vo
浏览 8提问于2022-04-01得票数 1
1回答
我有一个由8个方程和8个未知数组成的系统。在这组方程中没有任何数字。方程式中只有参数；其中8个是我的未知变量(以下称为未知变量)，我正在求解这些变量，还有一堆其他参数，它们不是未知值，例如alpha，beta，gamma (以下称为参数)。
我需要解方程组，这样我就可以得到参数方面的8个未知变量。
我在Maple中使用了fsolve(方程，变量)。但是，它返回一个错误，显示tau，alpha，beta，...都在方程式中，并且没有解出来。我如何才能让Maple理解这些不是我喜欢求解的变量，而是其他未知变量应该用来求解的参数呢？
下面是我的一个方程式的例子：
e1 := (L__af+L__mf
浏览 0提问于2015-09-12得票数 0
1回答
我想找到向量p (即p(1),p(2),p(3),...)的值，它使函数A(p)最大化。
我正在使用MATLAB做这件事，我找到了fsolve，我认为它可以帮助我。因此，我创建了函数A：
function A = myfun(p)
R = 0.1;
u1 = 500;
u2 = 400;
u3 = 300;
A = ( (p(1)+p(2)+p(3)) * (1/u1+1/u2+1/u3)) * ...
(1 + R*(p(1)^2+p(2)^2+p(3)^2) * (1/u1+1/u2+1/u3) );
然后，我需要求解一个方程
浏览 1提问于2014-03-20得票数 0
回答已采纳
我正在用Python模拟MATLAB代码。问题是Matlab vpasolve的确切实现在Python中不可用。我已经使用fsolve来求解方程，但在MATLAB和Python中的解决方案是不同的。
我想知道在Python中是否有其他可以像vpasolve一样工作的求解器。
浏览 119提问于2020-02-16得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
