-
0 votesanswersviews
Matlab：优化器/求解器的准确性差
我很难在Matlab上找到根发现问题的准确性 . 我有一个函数， Lik(k) ，并希望找到 k 的值 Lik(k)=L0 . 基本上，问题是各种内置的Matlab求解器（ fzero ， fminbnd ， fmincon ）没有像我想要的那样接近解决方案 . Lik() 是一个用户定义的函数，它涉及广泛的编码以计算数值逆拉普拉斯变换等，因此我不包括完整的代码 . 但是，我已广泛使用此功能，它... -
0 votesanswersviews
求解非线性方程组
救命 . 我试图在MATLAB中为家庭作业分配这个非线性方程组 . 我尝试过wolfram alpha和这个online equation solver，但它们都不起作用 . 我已经尝试了我的图形计算器，它一直在说非代数变量或表达式 . 这是我的两个未知数的方程式： .75*(1100)= x*10^(6.82485-943.453/(T+239.711)) 25*1100=(1-x)*10^(... -
1 votesanswersviews
Scipy - 非线性方程组的所有解
我有一个非线性方程组，其中任何n都可以选择，因此向量x =（x1，...，xn）的长度可以不同 . 例如，系统可以是这样的： f1(x1,...,xn) = sum( xi + xi^2 ) = 0, i={1,n} f2(x1,...,xn) = sum( e^xi + xi + sin(xi*pi) ) = 0, i={1,n} 根据this example，我使用scipy库的fs... -
0 votesanswersviews
在matlab中求解大约束线性和非线性方程组的最佳方法是什么？
MATLAB中约束线性和非线性方程组的求解 我正在使用直接刚度方法在MATLAB中解决FEM问题 . 现在将问题表述为非线性方程组： KU=F or CF=U with C = K^{-1}. 问题是在3-Dd中制定的，有200多个节点 . 该模型用作变形估计器，力总是输入和变形输出 . 所以有很多线性方程 . 现在需要扩展这个问题，需要实现2种类型的约束 . 一些节点必须被约束，以便它们只能围... -
3 votesanswersviews
如何在SymPy中求解具有符号系数和非交换变量的方程？
我想解决具有符号系数的非常简单的方程： from sympy import * a, b = symbols('a b', commutative=False) x = IndexedBase('x') basis = [a, b, a * b - b * a] el = b * a - a * b coefs = [x[k] for k in range(len(basis))] eq = e... -
0 votesanswersviews
获得满足Python中给定方程的随机数量的参数
我是Python的新手，我想获得满足给定方程的1000个参数 . 更确切地说，下面是球体的等式，假设球体中心的位置 (x0, y0, z0) 及其半径都给出 . (x - x0)^2 + (y - y0)^2 + (z - z0)^2 = radius^2 我怎么能随机生成一些满足python中球形方程的 (x, y, z) ？非常感谢！ -
1 votesanswersviews
如何在MATLAB中获得该方程的所有解？
我想解决以下等式：tan（x）= 1 / x 我做了什么： syms x eq = tan(x) == 1/x; sol = solve(eq,x) 但这只给出了解决方案的一个数值近似 . 之后我读到了以下内容： [sol, params, conds] = solve(eq, x, 'ReturnConditions', true) 但这告诉我它无法找到明确的解决方案 . 如何在某个给定范围... -
9 votesanswersviews
求解具有两个未知数的两个方程的系统
求解下面两个未知数的两个方程组的系统： a1，b1，c1，a2，b2和c2由用户自己输入 . 我一直试图找到问题的数学解决方案，我似乎无法走远... 到目前为止我尝试过的是： 从第一个方程找到y . （b1y = c1-a1x，y =（c1-a1x）/ b1） 然后我在第二个等式中替换y，得到一个方程式，其中1个未知数，在这种情况下为x . 但是，我无法解决这个等式，我得到了一些奇数/... -
17 votesanswersviews
SymPy - 任意数量的符号
我正在编写一个函数来解决任意数量的联立方程 . 方程的数量由函数的一个参数设置，每个方程由许多符号构成 - 与方程一样多的符号 . 这意味着我不能简单地对方程式进行硬编码，甚至不需要将方程组合在一起所需的符号;该函数需要能够处理任意数量的方程 . 所以，我的问题是，我如何生成符号列表？ 我有一个可能的解决方案，但我的直觉告诉我，它不会非常有效 . 如果有更好的方法，请告诉我 . 我是SymPy的新... -
1 votesanswersviews
用Python和Sympy求解方程并得到数值答案
我正试图用sympy来解决方程，但我想得到一个直接的数字答案 . 我的脚本是这样的： from sympy import * A,B,V=symbols('A,B,V') eq1=Eq(630.26*(V-39.0)*V*(V+39)-A+B,0) eq2=Eq(B,1.36*10**8*(V-39)) eq3=Eq(A,5.75*10**5*V*(V+39.0)) solve([eq1,eq2,... -
1 votesanswersviews
在Mathematica中求解四个方程组（带有对数）
我试图在四个变量中解决四个方程的系统 . 我已经阅读了类似问题的一些主题，并尝试遵循这些建议 . 但是我觉得这里的日志和交叉产品有点乱 . 这是确切的系统： 7 * w =（7 * w 5 * x 2 * yz）（0.76 0.12 * Log [w] -0.08 * Log [x] -0.03 * Log [y] -0.07 * Log [7 * w 5 * x 2 * YZ]），5 * x =... -
0 votesanswersviews
Scipy - 解包含CDF的非线性方程组
我有以下方程组（简化版） . y = x*scipy.special.gdtr(a, b, f(x,y)) x = y*scipy.speical.gdtr(c, d, f(x,y)) 这里 scipy.special.gdtr 指的是伽玛分布的CDF， a ， b 是伽马CDF所采用的相应的两个参数 . 我的问题如下： 在这种情况下，是否有首选求解器？ fsolve 或 brentq 或其... -
0 votesanswersviews
用Mathematica求解非线性矩阵系统
我试图用Mathematica找到矩阵方程的非线性系统的解 . 解决方案空间太复杂，无法尝试使用Solve，因此我尝试使用FindInstance . 这是一个非线性系统的最简单的例子，我试图找到一个解决方案： T + A + C + CBA = 0 I - A - C - ABC = 0 其中A，B，C是任何矩阵，Det [T] = 1，I是维度为2x2的单位矩阵 . 我还希望所有矩阵只有整数... -
0 votesanswersviews
用Python解决5个线性方程组
我尝试过使用矩阵，但它失败了 . 我看过外部模块和外部程序，但没有一个有效 . 如果有人可以分享一些有用的提示或代码，谢谢 . -
0 votesanswersviews
如何用约束求解这个线性方程？
我有a1，a2，a3，b1，b2，b3作为已知值，所以我可以列出如下的两个方程式 . 但是，这还不足以让H解决 . 所以，还有额外的约束|| H || = 1 . matrix H = [x y z]^T [a1, a2, a3] H = 0 [b1, b2, b3] H = 0 我想问一下如何解决这个问题？ -
0 votesanswersviews
使用sympy解决三次方程：错误的答案
任务是解决qubic方程 . 这个方程的系数是依赖于x的函数，所以我在每个点求解方程并得到一个根数组 . 问题在于有些点，其中由sympy给出的所有三个根都是虚构的，但是具有实系数的qubic方程是不可能的 . 此外，不同的sympy版本给出了不同的答案 . import numpy as np import pylab import sympy as sp from sympy import s... -
1 votesanswersviews
用符号表达式求解线性方程组
嗨，我试图用mathematica解决线性方程组 . 我有18个方程和18个未知数，系数矩阵具有满秩 . 所有条目都是象征性的，因为我试图通过分析解决问题 . 不幸的是，Mathematica永远不会停止评估 . 我准备了一个最小的工作示例： n = 18 A = Table[AA[i, j], {i, 1, n}, {j, 1, n}]; A // MatrixForm x = Table[xx... -
5 votesanswersviews
如何解决Sympy中的非线性方程？
如何解决SymPy中非线性方程的形式 y = P*x + Q + sqrt(S*x + T) 我知道的地方 y(0) ， y'(0) ， y(c) ， y'(c) . 我想找 P ， Q ， S 和 T . 并将 y 表示为 x 的函数 . 我对文档感到非常困惑 . 请帮忙 . -
0 votesanswersviews
求解非线性方程
我求解满足以下积分方程的f（kx，ky），k是2d中的向量 方程： 我将f（kx，ky）映射到2d N * N网格并求解了总N ^ 2个方程组 . from scipy.optimize import newton_krylov import numpy as n L=10.0 # cut off of integral N=100 # number of devision ... -
1 votesanswersviews
如何解决极大的方程组
我需要解决一个庞大的方程组 . 我现在需要解决的系统有30个方程和30个未知数，但我需要能够使用同一个系统来解决1,000个方程和1,000个未知数 . 我不需要知道所有未知变量的答案，我只需要知道其中一些变量的答案（30个方程式中的2个，以及1,000个方程式中的10个） . 我需要解决的方程式各自为： (c0 * x0) + (c1 * x1) + ... + (c_n * x_n) 其中 ... -
0 votesanswersviews
求解python中的非线性方程
我试图找到介质波导的基本TE模式 . 我尝试解决它的方法是计算两个函数并尝试在图上找到它们的交集 . 但是，我在绘制情节时遇到了交叉点 . 我的代码： def LHS(w): theta = 2*np.pi*1.455*10*10**(-6)*np.cos(np.deg2rad(w))/(900*10**(-9)) if(theta>(np.pi/2) or theta &l... -
2 votesanswersviews
在MATLAB / Maple /等上有效地解决符号线性方程组
基本上我有一个线性方程组，其中每个系数都是一个变量 . 例如，对于2x2系统： a11 * x1 a12 * x2 = f1 a21 * x1 a22 * x2 = f2 x1，x2 = ....（带有a11，a12，a21，a22，f1，f2的符号表达式） 我尝试通过Cramer的规则在MATLAB中解决它，以符号形式编写系统的矩阵并计算适当的决定因素 . 问题是，所需的时间随着系统... -
5 votesanswersviews
fsolve总是返回猜测/估计
我第一次使用scipy的optimize.fsolve函数来找到方程的根 . 问题是我用作猜测/估计值的数字是我得到的答案（在大约8位小数内） . 当使用full_output = True时，我得到退出标志为'1'，这应该意味着'解决方案收敛'，这是我理解的最好意味着输出确实是方程的根 . 我知道有一些有限数量的不同根（间隔开来），因为当我绘制方程时，我可以看到它们 . 此外，当我将起始点输入到... -
0 votesanswersviews
9个方程的非线性系统9未知数MATLAB - 以多项式比率耦合的未知数
目标 我想解决一个带有9个未知数的9个非线性方程组，并用Matlab求解 . 所有9个未知数都作为多项式的比率耦合（参见myfun lower） Fsolve x02=[5000,5000,5000,0.4,0.4,0.4,0.4,0.4,0.4]; ctrl2=[9894+1i*0.118,9894+1i*0.118,9894+1i*0.118,0.5,0.5,0.5,0.5,0.5,0.5];... -
17 votesanswersviews
SymPy - 任意数量的符号
我正在编写一个函数来解决任意数量的联立方程 . 方程的数量由函数的一个参数设置，每个方程由许多符号构成 - 与方程一样多的符号 . 这意味着我不能简单地对方程式进行硬编码，甚至不需要将方程组合在一起所需的符号;该函数需要能够处理任意数量的方程 . 所以，我的问题是，如何生成符号列表？ 我有一个可能的解决方案，但我的直觉告诉我，它不会非常有效 . 如果有更好的方法，请告诉我 . 我是SymPy的新手... -
0 votesanswersviews
使用Mathematica根据一组参数求解方程组
使用Mathematica，我想解决一个非线性方程组 f(x,y;m)=0 g(x,y;m)=0 其中m是方程所依赖的外部参数 . 由于方程非常难，我所做的是首先明确地修复m，然后找到m的特定值的解（通过FindRoot []） . 但这是非常缓慢和低效的（更具体地说，在我的情况下，确实我有三个参数） . 我想要的是一种自动方法，用数字方法研究m的一组离散值的解决方案，修复它而不是我不时明确地改... -
2 votesanswersviews
求一个字符串的等式程度？ [关闭]
In JAVA:- 给出字符串中的等式： String equation = "4*x^3-19*x^2+2*x-1=0"; ，如何找到它的度数？ 我想到了通过字符串查找插入符号'^'的位置，得到插入符号之后的数字，其中最高的是等式的程度 . 但是，如果数字超过一位数怎么办？实现这个就是设计一个方程解析器！ 所以你能告诉我其他任何方式吗？ 编辑： 我想知道它的完成方式，例如使用... -
2 votesanswersviews
为一组三角方程找到一个以上的根：python
我使用python的fsolve使用下面的脚本找到一组方程的根： from scipy.optimize import fsolve from scipy.interpolate import spline import math import numpy as np #from mpldatacursor import datacursor %pylab from matplotlib imp... -
3 votesanswersviews
使用Sympy，python2.7求解方程组
我想解决一个方程组 . 但我希望能够精确地将值变为“获得”，并将其作为“什么”的函数 . 为了更好地理解，我从here开始，我修改了一个例子： import sympy as sp x, y, z = sp.symbols('x, y, z') rho, sigma, beta = sp.symbols('rho, sigma, beta') f1 = sigma * (y - x) f2 = x... -
1 votesanswersviews
在Python中解决符号线性方程给出了意想不到的答案
首先，请原谅即将到来的长壁测试 . 我需要解决由坐标变换引起的以下方程组： 给定v_n，v_m，x_p，z_p，y_p，s，d，解决下面的系统v_yp，v_zp： I: v_n = v_yp * (1/(s*d)) + v_xp * y_p * (1/(s*d^2)) II: v_m = v_zp * (1/(s*d)) + v_xp * z_p * (1/(s*d^2)) III: v_...