SymPy无法求解三角方程组
SymPy是一个Python库，用于进行符号计算和数学表达式的操作。它提供了一套强大的工具，可以用于解决各种数学问题，包括代数方程组的求解。
然而，SymPy在求解三角方程组方面可能会遇到一些限制。三角方程组是由三角函数组成的方程组，例如sin(x) + cos(y) = 2，其中x和y是未知数。由于三角函数的特殊性质，这类方程组的求解可能比较复杂，甚至无法找到解析解。
对于无法求解的三角方程组，可以考虑使用数值方法进行近似求解。数值方法通过迭代计算逼近方程组的解，通常可以得到一个接近真实解的数值解。在Python中，可以使用数值计算库如NumPy或SciPy来实现这些数值方法。
对于符号计算和数学表达式操作，SymPy是一个非常强大的工具。它可以用于求解各种代数方程组、微分方程、积分等数学问题。此外，SymPy还提供了绘图功能，可以用于可视化数学函数和方程。
在腾讯云的产品中，与数学计算和符号计算相关的产品包括云函数（Serverless Cloud Function）和人工智能开发平台（AI Development Platform）。云函数可以用于部署和运行自定义的数学计算函数，而人工智能开发平台则提供了一系列用于数学计算和机器学习的工具和算法。
总结起来，SymPy是一个强大的Python库，用于进行符号计算和数学表达式操作。虽然它可能无法直接求解三角方程组，但可以通过数值方法进行近似求解。在腾讯云的产品中，可以使用云函数和人工智能开发平台来进行数学计算和符号计算相关的任务。
相关·内容
给出了这个例子：
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
回答已采纳
1回答
我正在尝试用Linsolve解决一个方程组，但显然一定是误解了什么，因为我总是得到意想不到的结果。假设我想要解决以下两个方程： a + b = 0
a - b + c = 0 我期望得到的结果是： b = 0.5*c 相反，Sympy返回空集。对于非线性求解，我得到了(-a)，这也没有多大意义： >>> import sympy
>>> a, b, c = sympy.symbols('a b c')
>>> Eqns = [a + b, a - b + c]
>>>sympy.linsolve(Eqns
浏览 9提问于2020-12-15得票数 0
回答已采纳
2回答
我试图让SymPy解决一个方程系统，但它给我一个错误，如下所示：
NotImplementedError: could not solve 3*sin(3*t0/2)*tan(t0) + 2*cos(3*t0/2) - 4
有没有别的方法可以让我解方程组：
sin(x)+(y-x)cos(x) = 0
-1.5(y-x)sin(1.5x)+cos(1.5x) = 2
我使用：
from sympy import *
solve([sin(x)+(y-x)cos(x), -1.5(y-x)sin(1.5x)+cos(1.5x)-2], x, y)
浏览 19提问于2018-01-12得票数 3
回答已采纳
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
回答已采纳
我想解决一个由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
2回答
如何在python中解出如下所示的简单方程组？
X= (2/3)*y + (1/3)*
Y= (2/3)*1 + (1/3)*x
我试过SymPy，但我想不通。
解了方程式部分
from sympy import *
from sympy.solvers.solveset import linsolve
x, y = symbols('x, y')
linsolve([x - 2/3*y, y - 2/3 - 1/3*x ], (x, y))
输出：{(0.571428571428571，0.857142857142857)}
类型为'sympy.sets.sets.F
浏览 2提问于2019-02-16得票数 1
1回答
我有一个由四变量$a，b，c，d$和两个参数$i，h$组成的线性方程组，其中方程的形式大致为$$a h^3i^3+b^2i^2 +c +d=0$$。
我想得到$a，b，c，d$，关于$i，h$。
这在SymPy中是可能的吗？如果没有的话，有人能不能用其他的方式来推荐它呢？
浏览 2提问于2016-10-26得票数 1
回答已采纳
我正在做一个CFD项目，我正在使用新的CUDA 5库"cusparse“来解线性方程组。我测试了示例代码"conjugateGradientPrecond“。结果表明，使用ILU的预条件梯度算法比不使用预条件算法的共轭梯度算法需要更多的时间来获得最终结果。前者的迭代次数较少，但在"cusparseScsrsv_solve“上耗时较长，因此整体时间较长。
这是我的问题，有没有其他预处理共轭梯度可以大大减少迭代，同时不包括任何耗时的函数，如"cusparseScsrsv_solve"？
浏览 1提问于2012-12-02得票数 2
目前使用渐近来求解以下方程组：
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
1回答
考虑这样的情况，给定一个MxM矩阵A和一个向量b，我想要解某种形式的inv(A @ A.T) @ b (其中我知道A是可逆的)。据我所知，使用solve_*而不是inv总是更快。也有更有效的解决私营部门司矩阵的变体(必须是A @ A.T )，使用Cholesky因式分解。
我的问题--因为我正在构造矩阵A @ A.T，只是为了马上扔掉它--是否有一个更专门的程序来解决线性方程组与A的克矩阵，而不必构造它呢？
浏览 2提问于2021-01-01得票数 0
回答已采纳
1回答
我试图在python中求解一个方程组(带有渐近性)，但我得到了"TypeError:无效NaN比较“。这些方程是“几乎线性的”(线性方程组，但也包含一个abs部分)。
对于函数f，我试图找到两点x0，x1，使得f( x0 ) = x1和f(x1) =x0(因此这两个方程形成一个系统)。
我有：
from sympy import *
c = 2
x = Symbol('x', real=True)
f1 = c * (1/2) * (Abs(2*x-1)) # define a function with variable 'x', also Abs
浏览 6提问于2022-10-27得票数 0
回答已采纳
假设我有以下方程组：
x = 1/y + b
y = 2*x
因此，b = x-1/(2*x)
不过，如果我试图在SymPy中定义这些表达式，
from sympy import *
b = symbols('b')
x = 1/y + b
y = 2*x
它给了我
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
<ip
浏览 0提问于2015-03-24得票数 3
回答已采纳
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回答
我正在试图求解大型的任意方程组(对于流程流板应用程序来说，物理组件模型是用简单的数学编写的)。
在处理非整数时，我对渐近的solve_poly_system函数的行为感到困惑，并且无法通过文档找到它。这个简单的代码片段：
sympy.solve_poly_system([x-5,y-x])
给我(5,5)如我所期望的。但是，当我将数字更改为浮点数时：
sympy.solve_poly_system([x-5.1,y-x])
我得到了错误:引发DomainError(“不能在%s之上计算Groebner基”% opt.domain) sympy.polys.polyerrors.DomainEr
浏览 2提问于2013-12-10得票数 3
回答已采纳
我想把一个n阶常微分方程简化为一个一阶方程组。这是为数值分析做准备。我使用Sympy和Sagemath for Computer Algebra，但我没有在它们中找到任何函数来执行这种类型的约简。我不确定是否有其他人可以指出此功能是否存在于Sympy或Sagemath中。
这方面的一个例子是减少方程式：
x''' - 2x'' + x' = 0
一阶方程组：
[0 1 0
0 0 1
0 -1 2]
浏览 0提问于2014-08-31得票数 2
2回答
我需要解线性代数方程组A.X = B
矩阵A是双精度的，大小约为33000x33000，当我试图分配它时，我会得到一个错误：
无法在计算数组大小时分配数组溢出。
由于我在英特尔MKL库中使用LAPACK dposv，所以我想知道是否有一种方法可以将一个较小的矩阵传递给库函数？(因为只需要解决矩阵数组的一半)
dposv函数只需要A. 的上三角矩阵或下三角矩阵
更新:请注意，A矩阵是N，但它以lda: INTEGER作为The leading dimension of a; lda ≥ max(1, n)。所以，也许有一种方法可以将A解析为一维数组？
浏览 7提问于2016-12-31得票数 1
回答已采纳
2回答
我一直在想这个问题，但我似乎不能把它当回事。
我想要解一个矩阵，它有三个未知数x，y，z，所以它们都等于相同的数。
让我说我的方程式是：
x + 3 = A
y(2y - 2) = 2A
z(4z - 1) = A
因此，我可以构造一个矩阵，如下所示：
[(X + 3) , 0 , 0] [0] [A]
[ 0 ,(2y - 2), 0] [y] = [2A]
[ 0 , , 0, (4z -1)] [z] [A]
我知道，但那是只有当答案(A)已经知道。
我的问题是，我需要构造一个循环来强迫(A)的答案，还是有更多的仿生方
浏览 4提问于2015-02-20得票数 3
回答已采纳
2回答
我有两个耦合的非线性微分方程组。
为此，我编写了以下代码，以使用Sympy获取两个函数的值：
from sympy import *
t = symbols('t')
c, b, B, alp, mu = symbols('c b B alp mu', integer=True)
f, g = symbols('f g', cls=Function)
print(solve([Eq(f(t).diff(t)+c*f(t)+0.5*B*f(t)**2-b*g(t), 0), Eq(g(t).diff(t)+b*g(t)-c*f(t), (1-
浏览 1提问于2017-01-14得票数 0
我正在尝试使用SymPy，这是我以前从未做过的，来解决一个非线性方程组，但我不能让它工作。
以下是我尝试过的方法：
from sympy import *
a, b, c, d, e, f, x, y, z, g0, g1, g2 = symbols('a:f x:z g:3')
system = [(a * (x - y*z))/(b * (1 - z**2)) - g1,
(a * (y - x*z))/(c * (1 - z**2)) - g2,
f - d * (a * (x - y*z))/(b * (1 - z**2)) - e * (a * (y - x*z))
浏览 11提问于2019-10-22得票数 0
3回答
我必须在MATLAB中求解一个线性方程组A*x=B，其中A是对称的，它的元素依赖于指数的差异：Aij=f(i-j)。
我使用迭代求解器，因为A的大小是40000x40000。迭代求解器需要确定产品A*x，其中x是测试解决方案。该乘积的计算结果是一个卷积，因此可以用快速傅里叶变换(cputime ~ Nlog(N)而不是N^2)的方法来实现。关于这个问题，我有以下几个问题：
这是卷积圆吗？因为如果它是圆形的，我认为我必须为新的矩阵使用一个特定的索引来获得fft。是那么回事吗？
我发现很难为fft的例程编程，因为我无法理解我应该使用的索引。有没有现成的程序，我可以用快速傅立叶变换直接评估乘
浏览 11提问于2013-11-18得票数 7
1回答
给出几个SymPy矩阵M
M = Matrix([
[0.000111334436666596, 0.00114870370895408, -0.000328330524152990, 5.61388353859808e-6, -0.000464532588930332, -0.000969955779635878, 1.70579589853818e-5, -5.77891177019884e-6, -0.000186812539472235, -2.37115911398055e-5],
[-0.00105346453420510, 0.000165063406707273, -0.0
浏览 0提问于2018-12-15得票数 4
回答已采纳
我对使用Sympy真的很陌生，我正在尝试解决以下方程系统： a*cos(theta1) + b*cos(theta2) + c*cos(theta3) + d*cos(theta4) =0 a*sin(theta1) + b*sin(theta2) + c*sin(theta3) + d*sin(theta4) =0 其中theta1是输入变量，它会随时间变化(类似于theta1 = 2*t)，而除theta2和theta3之外的所有其他变量都是常量。 我希望Sympy根据输入变量为我提供theta2和theta3的表达式。 我尝试的解决方案如下： from sympy import *
浏览 37提问于2021-04-29得票数 1
回答已采纳
1回答
有没有办法用SymPy得到三个球体相交(三边测量)的解？sympy.geometry没有球体对象，所以直接的方法是不可行的。如Trilateration and the Intersection of Three Spheres所示，SymPy能否解决非线性方程组
浏览 20提问于2021-11-16得票数 1
回答已采纳
我想用数值方法解决这个方程
beta.ppf(x，a，b)-c=0
其中a，b，c是已知的常数。当我试着
from sympy.solvers import solve
from sympy import Symbol
from scipy.stats import beta
x=Symbol('x')
solve(beta.ppf(x,a,b)-c,x)
它还给了我
TypeError:无法确定关系的真值
我怎么才能修好它？
浏览 2提问于2016-08-28得票数 0
我试图得到弗洛伊登斯坦方程(Psi)的一个角度：
k1 * cos(psi) - k2 * cos(fi) + k3 - cos(psi - fi) =0
我有k1，k2，k3和fi的值。我尝试了以下几种方法：
from sympy import Symbol, solve, cos
x = Symbol('x')
realPsi = solve(k1 * cos(x) - k2 * cos(fi) + k3 - cos(x - fi), x)
我收到这个错误：
File "/usr/lib/python2.7/dist-packages/sympy/solvers/
浏览 1提问于2013-05-28得票数 1
回答已采纳
1回答
我用同情来定义符号，
a, b, c = sympy.symbols(['a', 'b', 'c'])
然后，当我试图解出下列方程组时，
sympy.solve([sympy.Eq(b - a**2.552 - c), sympy.Eq(a, 2)])
我知道解决办法了
[{b: c + 5.86446702875684, a: 2.00000000000000}]
但是，当我试图解决时，
sympy.solve([sympy.Eq(b - a**2.552 - c), sympy.Eq(b, 2)])
它似乎一直在运行(大约4小时)，没有解决办
浏览 2提问于2016-10-05得票数 4
回答已采纳
假设我有一个方程系统，如下所示，但我有一个要求，在我的方程系统中，我需要2个语句中的1个为真。
示例：
sympy.solve(eq1或eq2，eq3或eq4，eq5或eq6，...，vars)
当我放入or语句时，它不能像我预期的那样工作(第一个eq总是被计算出来)。我该如何用或按我想要的方式来表达这个方程组呢？
谢谢。
浏览 0提问于2015-10-13得票数 1
1回答
我想得到微分方程组的符号解。给出一个"NotImplementedError“。diff(t，1)效果很好。获得解决方案的最佳方法是什么？或者在符号解不可能的情况下，如何将其转换为一阶的数值解。任何帮助都将不胜感激。
import sympy as sp
import mpmath
from sympy import Function, symbols, dsolve, Derivative, erf, sqrt, power, diff
a, t, rx, ry = symbols('a t rx ry')
rx = Function("rx")(t
浏览 2提问于2019-12-22得票数 0
回答已采纳
我需要将我的MATLAB代码转换成C++，其中包括形式为xA=0的线性方程组。
我知道特征可以处理线性方程组Ax=b，我的问题是:是否有一种方法来求解线性方程组xA=b，使用C++ (Visual 2010)，A是一个稀疏矩阵？如果没有，我可以使用什么库？
谢谢你的帮助。
浏览 0提问于2013-10-07得票数 0
回答已采纳
1回答
我想要像A = ax + by和B = cx + dy这样象征性地解一个方程组，对于x和y，显式地表示。
我尝试了作为solve([A, B], [x, y])的解函数，但没有起作用。它正在返回一个空列表，[]。
我怎样才能用同情来解决它呢？
这是我想要解决的实际方程式：
from sympy import*
i,j,phi, p, e_phi, e_rho = symbols(r'\hat{i} \hat{j} \phi \rho e_\phi e_\rho')
e_rho = cos(phi)*i + sin(phi)*j
e_phi = -p*sin(ph
浏览 1提问于2022-05-29得票数 0
回答已采纳
我想知道，这就是APM求解器背后的理论方法。
我正在解决一个问题，在这个问题中，我必须解一个没有任何约束的非线性方程组。我正在使用BPOPT解算器，它工作得很好。我只想知道解决这些系统的方法是什么。这是牛顿的方法？准牛顿法？我没有找到任何关于这件事的文件。你能告诉我在哪里能找到这样的信息吗？
是否也有可能预先知道gekko求解器需要的RAM来解一个含有(N)变量的非线性方程组？
谢谢
浏览 3提问于2020-06-17得票数 2
回答已采纳
在每一个递归循环中，我都要求解两个方程：
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
回答已采纳
9回答
使用Python. (Numpy，Scipy或Sympy)解决对非线性方程组的(最好的)方法是什么？
例：
x+y^2 =4
e^x+ xy =3
解决上述两个问题的代码片段将是很棒的。
浏览 12提问于2012-01-05得票数 85
回答已采纳
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
我要用“n-解数”函数来解这个方程。不幸的是，出现了此错误：
ValueError: Could not find root within given tolerance. (435239733.760000060718 > 2.16840434497100886801e-19)
Try another starting point or tweak arguments.
守则是：
import sympy
d=[0.3, 32.6, 33.4, 241.7, 396.2, 444.4, 480.8, 588.9, 1043.9, 1136.1, 1288.1, 1408.1, 143
浏览 6提问于2022-06-07得票数 0
如何用SymPy编写方程组并输出等效的Latex？latex函数一次只接受一个表达式。
import sympy as sp
x, y, z = sp.symbols('x, y, z')
eq1 = sp.Eq(x + y + z, 1)
eq2 = sp.Eq(x + y + 2 * z, 3)
output = sp.latex() # Do something here?
浏览 8提问于2022-11-22得票数 0
我正在尝试用渐近自动求解$x**-d - sympy.exp(100)$。
import sympy
from sympy import solve, Symbol
x = Symbol('x')
d = Symbol('d', integer=True)
print(solve(x**-d - sympy.exp(100), x))
不返回解决方案[]，但解决方案应为$x=e^{-100/d}$ (例如，参见)
我做错什么了？删除integer=True时的结果相同。
然而，
print(solve(x**-d - sympy.exp(100), x))
返
浏览 8提问于2020-03-03得票数 0
我正在使用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
1回答
我想通过使用r重新排列变量的下面的等式。
P = C * ((1-(1+r)**(-n)))/r + fv*(1+r)**(-n)
to
r = blabla...
我明白，同情与这样的重新安排任务有关。所以我写了下面的代码。
# Solve the equation for r
import sympy
from sympy import symbols
P, C, r, n, fv = sympy.symbols('P C r n fv')
eq = sympy.Eq(P, C * ((1-1/(1+r)**n))/r + fv/(1+r)**n)
sympy.solve
浏览 2提问于2020-04-13得票数 1
回答已采纳
比如说，我在试图解方程组：
对于a和ɸ，其中b、c和d是复数，而LHS上的矩阵是酉矩阵。
我手头的确实成功地完成了这项工作，但也有一些边缘的情况，它错过了。
from sympy import *
def get_angles(a, b, c, d):
theta, phi, lamb = symbols('\\theta \\phi \\lambda', real=True)
a_eq = Eq(cos(theta / 2), a)
b_eq = Eq(exp(I * phi) * sin(theta / 2), b)
c_eq =
浏览 0提问于2019-05-26得票数 1
回答已采纳
4回答
当我试图用python3中的以下代码求解一个由三个方程组成的方程组时，我总是会出错：
import sympy
from sympy import Symbol, solve, nsolve
x = Symbol('x')
y = Symbol('y')
z = Symbol('z')
eq1 = x - y + 3
eq2 = x + y
eq3 = z - y
print(nsolve( (eq1, eq2, eq3), (x,y,z), (-50,50)))
以下是错误消息：
追溯(最近一次调用)：文件"/usr/lib
浏览 1提问于2017-01-17得票数 3
1回答
在SymPy中，是否可以对不定积分应用限制并对其求值？
import sympy
from sympy.abc import theta
y = sympy.sin(theta)
Y_indef = sympy.Integral(y)
Y_def = sympy.Integral(y, (theta, 0, sympy.pi / 2))
Y_def.evalf()生成一个数字。
我正在寻找像Y_indef.evalf((theta, 0, sympy.pi/2))这样的东西来获得相同的答案。
浏览 4提问于2013-02-08得票数 4
回答已采纳
2回答
我目前正在C++中进行流体模拟，算法的一部分是求解稀疏线性方程组。为此，人们建议使用艾根图书馆。我决定用我写的这个简短的程序来测试它：
#include <Eigen/SparseCholesky>
#include <vector>
#include <iostream>
int main() {
std::vector<Eigen::Triplet<double>> triplets;
triplets.push_back(Eigen::Triplet<double>(0, 0, 1));
浏览 1提问于2018-09-04得票数 2
回答已采纳
1回答
当我用交感解数学问题时遇到了一些问题，我的代码是这样的，
eq1=Eq(0, -sympy.sqrt(x**2) + sympy.sqrt((x + 4)**2))
sympy.solve(eq1,x)
这个方程的答案应该是-2，但我得到了一个空的列表。
浏览 0提问于2018-09-13得票数 2
回答已采纳
1回答
--我正试图用数值方法来求解一个方程，它可以包含在Python中的截断。我想知道最好的图书馆和方法是什么？以下是有关此问题的更多详细信息：
方程式每次都变了。从人的角度来看，方程应该非常简单；它们包括诸如+、-、*、/之类的公共运算符，它们有时还具有截断函数(截断为整数)或限制函数(将括号中的值限制在两个提供的边界之间)或(很少)多个变量。有几个例子(这些是单独的例子，而不是方程组)将是：
TRUNCATE(VAR_1 + 300) - 50.4 = 200(VAR_2 + VAR_3)*3 = 35LIMIT(3,5)(VAR_4) = 8VAR_5 = 34
(这并不完全是方程的样子，因为
浏览 7提问于2021-07-15得票数 0
回答已采纳
我试图用函数来解未知x的线性方程组A*x=b。下面是一个很好的例子：
import numpy as np
import scipy.linalg as linalg
A = np.array([[ 0.18666667, 0.06222222, -0.01777778],
[ 0.01777778, 0.18666667, 0.01777778],
[-0.01777778, 0.06222222, 0.18666667]])
b = np.array([0.26666667, -0.26666667, -0.4])
x = li
浏览 4提问于2020-01-10得票数 2
在Dymola建立了一个具有朗肯循环的太阳能电站模型，虽然初始化工作正常，但是在2840s的时候，由于非线性方程组的求解失败，导致仿真停止。如下面的截图所示，Dymola建议提供更好的起始值，但这没有任何意义，因为失败发生在模拟期间，而不是初始化过程中。有了更好的起始值，情况根本不会改善。我的问题是:我应该如何处理在模拟过程中无法求解非线性方程？ ?
浏览 71提问于2020-07-06得票数 1
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回答
本质上，我想做的是取一个复数的绝对值，它在实部和虚部都有一个参数，并导出这个参数的绝对值。然后，我想知道那个参数的值，导数是零。所以本质上，我想知道带参数的复数绝对值的Min/Max。
我在使用同情时遇到了一个问题，我不太知道如何解决它/如果这是他们的问题的话。在我把它报告为窃听器之前，我想问一下。好的，这基本上就是我想要做的(实现比较复杂，但是错误是一样的)
0 import sympy
1 import mpmath
2 sympy.init_printing() # enable pretty printing
3 a,b = sympy.symbols("a,b") #
浏览 0提问于2018-06-16得票数 1
回答已采纳
1回答
我试图找到一个4x4矩阵R的特征值和特征向量。为此，我使用了R.eigenvals()和R.eigenvects()命令，但是python返回空括号。我知道特征值是1和e^(+-2ix)。我做错了什么吗？
import pprint
from sympy import *
R=Matrix([[(cos(x))**2, -cos(x)*sin(x), -cos(x)*sin(x), (sin(x))**2],[cos(x)*sin(x), (cos(x))**2, -(sin(x))**2, -cos(x)*sin(x)], [cos(x)*sin(x), -(si
浏览 4提问于2017-05-02得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
