首页 文章
  • 6 votes
     answers
     views

    解决具有不同模的联系方程组

    有没有算法来解决在不同模空间中表达的方程组?例如,考虑这个方程组: (x1 + x2 ) % 2 = 0 ( x2 + x3) % 2 = 0 (x1 + x2 + x3) % 3 = 2 该系统的解决方案之一是: x1 = 0 x2 = 2 x3 = 0 我怎么能算术地找到这个解决方案(不使用暴力算法)? 谢谢
  • 17 votes
     answers
     views

    SymPy - 任意数量的符号

    我正在编写一个函数来解决任意数量的联立方程 . 方程的数量由函数的一个参数设置,每个方程由许多符号构成 - 与方程一样多的符号 . 这意味着我不能简单地对方程式进行硬编码,甚至不需要将方程组合在一起所需的符号;该函数需要能够处理任意数量的方程 . 所以,我的问题是,我如何生成符号列表? 我有一个可能的解决方案,但我的直觉告诉我,它不会非常有效 . 如果有更好的方法,请告诉我 . 我是SymPy的新...
  • 1 votes
     answers
     views

    用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 votes
     answers
     views

    在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 votes
     answers
     views

    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 votes
     answers
     views

    用Mathematica求解非线性矩阵系统

    我试图用Mathematica找到矩阵方程的非线性系统的解 . 解决方案空间太复杂,无法尝试使用Solve,因此我尝试使用FindInstance . 这是一个非线性系统的最简单的例子,我试图找到一个解决方案: T + A + C + CBA = 0 I - A - C - ABC = 0 其中A,B,C是任何矩阵,Det [T] = 1,I是维度为2x2的单位矩阵 . 我还希望所有矩阵只有整数...
  • 0 votes
     answers
     views

    用Python解决5个线性方程组

    我尝试过使用矩阵,但它失败了 . 我看过外部模块和外部程序,但没有一个有效 . 如果有人可以分享一些有用的提示或代码,谢谢 .
  • 0 votes
     answers
     views

    如何用约束求解这个线性方程?

    我有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 votes
     answers
     views

    使用sympy解决三次方程:错误的答案

    任务是解决qubic方程 . 这个方程的系数是依赖于x的函数,所以我在每个点求解方程并得到一个根数组 . 问题在于有些点,其中由sympy给出的所有三个根都是虚构的,但是具有实系数的qubic方程是不可能的 . 此外,不同的sympy版本给出了不同的答案 . import numpy as np import pylab import sympy as sp from sympy import s...
  • 1 votes
     answers
     views

    用符号表达式求解线性方程组

    嗨,我试图用mathematica解决线性方程组 . 我有18个方程和18个未知数,系数矩阵具有满秩 . 所有条目都是象征性的,因为我试图通过分析解决问题 . 不幸的是,Mathematica永远不会停止评估 . 我准备了一个最小的工作示例: n = 18 A = Table[AA[i, j], {i, 1, n}, {j, 1, n}]; A // MatrixForm x = Table[xx...
  • 5 votes
     answers
     views

    如何解决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 votes
     answers
     views

    求解非线性方程

    我求解满足以下积分方程的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 votes
     answers
     views

    如何解决极大的方程组

    我需要解决一个庞大的方程组 . 我现在需要解决的系统有30个方程和30个未知数,但我需要能够使用同一个系统来解决1,000个方程和1,000个未知数 . 我不需要知道所有未知变量的答案,我只需要知道其中一些变量的答案(30个方程式中的2个,以及1,000个方程式中的10个) . 我需要解决的方程式各自为: (c0 * x0) + (c1 * x1) + ... + (c_n * x_n) 其中 ...
  • 0 votes
     answers
     views

    求解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 votes
     answers
     views

    在MATLAB / Maple /等上有效地解决符号线性方程组

    基本上我有一个线性方程组,其中每个系数都是一个变量 . 例如,对于2x2系统: a11 * x1 a12 * x2 = f1 a21 * x1 a22 * x2 = f2 x1,x2 = ....(带有a11,a12,a21,a22,f1,f2的符号表达式) 我尝试通过Cramer的规则在MATLAB中解决它,以符号形式编写系统的矩阵并计算适当的决定因素 . 问题是,所需的时间随着系统...
  • 0 votes
     answers
     views

    Matlab:优化器/求解器的准确性差

    我很难在Matlab上找到根发现问题的准确性 . 我有一个函数, Lik(k) ,并希望找到 k 的值 Lik(k)=L0 . 基本上,问题是各种内置的Matlab求解器( fzero , fminbnd , fmincon )没有像我想要的那样接近解决方案 . Lik() 是一个用户定义的函数,它涉及广泛的编码以计算数值逆拉普拉斯变换等,因此我不包括完整的代码 . 但是,我已广泛使用此功能,它...
  • 3 votes
     answers
     views

    如何在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...
  • 1 votes
     answers
     views

    如何在MATLAB中获得该方程的所有解?

    我想解决以下等式:tan(x)= 1 / x 我做了什么: syms x eq = tan(x) == 1/x; sol = solve(eq,x) 但这只给出了解决方案的一个数值近似 . 之后我读到了以下内容: [sol, params, conds] = solve(eq, x, 'ReturnConditions', true) 但这告诉我它无法找到明确的解决方案 . 如何在某个给定范围...
  • 2 votes
     answers
     views

    在matlab中求解欠定方程组

    我有一个欠定线性方程组Ax = b(即比方程更多的未知数),我想在matlab中求解 . 我知道这通常意味着无数的解决方案,但我也知道解决方案应该是正整数且小于一定数 . 我能找到满足这些额外要求的所有解决方案吗? 这个问题来自于一个未知的矩阵,我知道每行和每列的总和 . 例如未知矩阵找 0 3 2 0 0 2 4 1 2 1 0 0 行总和已知 5 7 3 列总和知道 2 6 6 1 我已...
  • 5 votes
     answers
     views

    fsolve总是返回猜测/估计

    我第一次使用scipy的optimize.fsolve函数来找到方程的根 . 问题是我用作猜测/估计值的数字是我得到的答案(在大约8位小数内) . 当使用full_output = True时,我得到退出标志为'1',这应该意味着'解决方案收敛',这是我理解的最好意味着输出确实是方程的根 . 我知道有一些有限数量的不同根(间隔开来),因为当我绘制方程时,我可以看到它们 . 此外,当我将起始点输入到...
  • 0 votes
     answers
     views

    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 votes
     answers
     views

    SymPy - 任意数量的符号

    我正在编写一个函数来解决任意数量的联立方程 . 方程的数量由函数的一个参数设置,每个方程由许多符号构成 - 与方程一样多的符号 . 这意味着我不能简单地对方程式进行硬编码,甚至不需要将方程组合在一起所需的符号;该函数需要能够处理任意数量的方程 . 所以,我的问题是,如何生成符号列表? 我有一个可能的解决方案,但我的直觉告诉我,它不会非常有效 . 如果有更好的方法,请告诉我 . 我是SymPy的新手...
  • 0 votes
     answers
     views

    使用Mathematica根据一组参数求解方程组

    使用Mathematica,我想解决一个非线性方程组 f(x,y;m)=0 g(x,y;m)=0 其中m是方程所依赖的外部参数 . 由于方程非常难,我所做的是首先明确地修复m,然后找到m的特定值的解(通过FindRoot []) . 但这是非常缓慢和低效的(更具体地说,在我的情况下,确实我有三个参数) . 我想要的是一种自动方法,用数字方法研究m的一组离散值的解决方案,修复它而不是我不时明确地改...
  • 2 votes
     answers
     views

    求一个字符串的等式程度? [关闭]

    In JAVA:- 给出字符串中的等式: String equation = "4*x^3-19*x^2+2*x-1=0"; ,如何找到它的度数? 我想到了通过字符串查找插入符号'^'的位置,得到插入符号之后的数字,其中最高的是等式的程度 . 但是,如果数字超过一位数怎么办?实现这个就是设计一个方程解析器! 所以你能告诉我其他任何方式吗? 编辑: 我想知道它的完成方式,例如使用...
  • 2 votes
     answers
     views

    为一组三角方程找到一个以上的根: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 votes
     answers
     views

    使用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 votes
     answers
     views

    在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_...
  • 3 votes
     answers
     views

    用python求解一个超越方程组

    假设我有以下四个方程式: cos(x)/ x = a cos(y)/ y = b a b = 1 c sinc(x)= d sinc(y) 对于未知变量 x, y, a 和 b . 请注意 cos(x)/x=a 有多个解决方案 . 类似于变量 y . 我只对 x 和 y 值感兴趣,这是第一个正面的根(如果重要的话) . 您可以安全地假设 a, b, c 和 d 是已知的实常...
  • 1 votes
     answers
     views

    在MATLAB中绘制二阶方程的解

    你能帮我解决下面这个问题:我想解决一个带有两个未知数的二阶方程,并用结果绘制一个椭圆 . 这是我的功能: fun = @(x) [x(1) x(2)]*V*[x(1) x(2)]'-c V is 2x2 对称矩阵, c 是一个正常数,有两个未知数, x1 和 x2 . 如果我使用fsolve求解方程,我注意到解决方案对初始值非常敏感 fsolve(fun, [1 1]) 是否有可能在不提供精...
  • 1 votes
     answers
     views

    数学游戏问题随机生成的方程数

    您好我的同事编码器我不知道如何解决使用随机生成的数学方程式制作多选游戏的问题 . 即使在控制台内,我当前的编码也没有返回任何答案或数字 . 请不要告诉我具体怎么做,请帮帮我 . 我会喜欢一些指示和解释 . 所以我对如何做到这一点的原始角度是使用Math.random对象将两个随机生成的数字乘以10,这样我得到0-10中的整数,然后我想显示它们在框内标有带id标签的问题 . 现在我真的需要帮助才能让...

热门问题