Python - 求解y的多项式
Python是一种高级编程语言，具有简洁、易读、易学的特点。它被广泛应用于各个领域的软件开发、数据分析、人工智能等方面。
求解y的多项式是指通过给定的多项式表达式，计算出对应的y值。在Python中，可以使用多种方法来求解y的多项式，下面介绍两种常见的方法：
- 使用numpy库：
Numpy是Python中常用的数值计算库，提供了丰富的数学函数和数组操作功能。可以使用numpy库中的polyval函数来求解多项式的值。具体步骤如下：
- 导入numpy库：
import numpy as np - 定义多项式系数：
coefficients = [1, 2, 3]（例如多项式为y = 1 + 2x + 3x^2） - 定义x的值：
x = 2 - 使用polyval函数求解y的值：
y = np.polyval(coefficients, x) - 打印结果：
print(y)
- 导入numpy库：
推荐的腾讯云相关产品：腾讯云服务器（https://cloud.tencent.com/product/cvm）
- 使用sympy库：
Sympy是Python中的一个符号计算库，可以进行符号运算、求解方程等。可以使用sympy库中的symbols和Poly函数来求解多项式的值。具体步骤如下：
- 导入sympy库：
import sympy as sp - 定义符号变量：
x = sp.symbols('x') - 定义多项式表达式：
polynomial = 1 + 2*x + 3*x**2（例如多项式为y = 1 + 2x + 3x^2） - 定义x的值：
x_value = 2 - 使用subs函数将x替换为具体的值，并求解y的值：
y = polynomial.subs(x, x_value) - 打印结果：
print(y)
- 导入sympy库：
推荐的腾讯云相关产品：腾讯云函数（https://cloud.tencent.com/product/scf）
以上是使用Python求解y的多项式的两种常见方法。根据具体的需求和使用场景，选择合适的方法来求解多项式的值。
相关·内容
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
回答已采纳
clear all;close all;clc
syms x y
[solx,soly] = solve(y+cos(x) == 0, (x-1)^2+(y-1)^2-4 == 0)
我可以得到x1，y1坐标是0.41864，-0.91364。但是，我找不到另一个根。有人能帮我吗？
浏览 0提问于2021-03-12得票数 0
1回答
我只想对$(x，y)$形式的椭圆曲线中的笛卡尔点进行采样。给定x或y中的一个，我希望能够计算另一个的可能性。
例如:我有EllipticCurve('11a')，我想知道它在(x，y)坐标中与x轴相交的位置。曲线在哪里与x=5相交呢？或者它在哪里与y=3相交？
有人能帮我做这件事吗？
浏览 1提问于2017-01-26得票数 0
1回答
我有一组多项式方程和实数的线性不等式。如何确定是否存在解决方案？我不需要符号解，也不需要保证答案是正确的，但是算法大部分时间都应该给出正确的答案。
例如:考虑由实数、x^2+y^2-1和y的单多项式x组成的多项式系统。线性不等式规定，我在一个小方格上寻找一个解，即：0.500<x<0.501，0.865<y<0.866。然后，该系统是：
x^2+y^2-1==0,
0.500<x<0.501,
0.865<y<0.866
这个系统有无限多的解决方案，例如x=0.5005、y=Sqrt(1-x^2)=0.865737...。
理想情况下，语言应该是
浏览 4提问于2022-05-26得票数 2
给出了由径向畸变和切向畸变定义的针孔模型的畸变系数D = k1，k2，p1，p2，k3。如何用这些系数来解构图像呢？我知道cv2.undistort()函数。有人能解释我如何实现非失真函数来执行迭代最小化吗？
# TODO: Iterative minimization required
r2 = x_ * x_ + y_ * y_
# Remove tangential distortion
x_tangential = x_ - (2 * distort_coeffs[2] * x_ * y_ + distort_coeffs[3] * (r2 + 2 *
浏览 3提问于2019-02-08得票数 1
回答已采纳
我正在研究python所关注的ode数值集成问题，特别是我发现了scipy.integrate函数solve_ipv。我尝试了中显示的示例，但是在与Lotka示例相关的代码中可能有一个错误：
def lotkavolterra(t, z, a, b, c, d):
x, y = z
return [a*x - b*x*y, -c*y + d*x*y]
sol = solve_ivp(lotkavolterra, [0, 15], [10, 5], args=(1.5, 1, 3, 1))
t = np.linspace(0, 15, 300)
z = sol.sol(t)
浏览 1提问于2021-01-12得票数 4
回答已采纳
是否有任何有效的算法可用于生成方程的所有整数解，如下面的算法？
(8+3n)m = 11 n∈{0,1}，m∈ℤ+
(5+(7+3x+2y)a+3z)b =30xx，y，z∈{0,1}，a，b∈ℤ+
理想情况下，我希望能够分别生成n，m和a，b，x，y，z的所有有效整数值的集合。至少，我想用一种方法来检验方程是否是可解的。考虑到这些方程是非线性的，我可以想象用于求解简单丢番图方程的典型方法会在这里失败。
我真的很感激我能得到的任何帮助！
浏览 0提问于2018-10-10得票数 1
1回答
我儿子正在学习如何计算抛物线的公式，在他的可汗学院课程中使用一个直角和焦点。(a，b)是焦点，k是directrix作为y=k的参数。我想向他展示一种使用Sympy检查结果的简单方法；编程对巩固内部算法有很大帮助。第一步显然是建立方程。
Parabola = Eq(sqrt((y-k)**2),sqrt((x-a)**2+(y-b)**2))
我首先为y解决了这个问题，然后展示了如何替换值并导出方程，因此：
Y = solve(Parabola,y)
这是一个合理的形式，收集了1/(2b-2k)到外面。然后，将焦点和方向性的值代入方程，得到方程y= 1/6*(x**2+16*x+49)，这是正
浏览 3提问于2016-05-29得票数 3
回答已采纳
这里有没有人能帮我解决以下问题？下面的代码计算适合给定数据集的最佳多项式，即指定次数的多项式。
不幸的是，无论数据集是什么，通常是6度或更高，MATLAB得到的是一个完全错误的拟合。通常，拟合曲线完全远离数据，以一种指数方式向下看。(参见示例: degree = 8)。
x=[1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5] % experimental x-values
y=[4.3 6.2 10.1 13.5 19.8 22.6 24.7 29.2] % experimental y-values
degree=8; % specify the degree
A = ze
浏览 0提问于2014-12-16得票数 2
1回答
我正在试着解决以下第四批BVP
Y‘=*y
我的x变量是一个有100个节点的行空间。正如你所看到的，K是同一个length=100的向量，使得方程非齐次。但是，当我按solve时，会出现以下错误：
Cell In [11], line 18, in fun(x, y)
17 def fun(x, y):
---> 18 ans = vector-np.multiply(C,y[0])
19 return np.vstack((y[1],y[2],y[3],ans))
ValueError: operands could not be broadcas
浏览 14提问于2022-11-03得票数 0
3回答
我有一组用未知方法插值的点，或者更准确地说，这个方法是已知的，但它可以是几种方法之一--可以是多项式插值，样条，简单线性……还有一条线，让我们现在假设它是以y= ax + b的简单形式给出的。
对于插值，我不知道使用了哪种方法(即函数是隐藏的)，所以我只能确定某个x的y，以及给定的y值的x。
找到两者之间的交叉点的通常方法是什么？
浏览 4提问于2010-09-25得票数 1
2回答
我有一个带有n顶点的加权无向图G。其中两个顶点是X和Y。
我需要找到一条最短的路径，它从X开始，在Y结束，经过G的所有顶点(以任何顺序)。
我怎么能做到这一点？
这不是旅行商问题:我不需要访问每个顶点一次，我也不想回到第一个顶点。
浏览 0提问于2016-04-16得票数 6
1回答
用二维已知点拟合三维曲面
、、、、
在一张覆盖未知三维物体(例如，原则上为零高斯曲率)的纸上的点的三维位置Y的最佳方法是什么，知道：
2D纸上所有点的二维位置(X)
点子集的三维位置(y in Y，for f(X)=Y)？
如果已知的三维位置没有标记(即当我们知道某些点的三维位置，但我们不知道它们对应的是哪一个二维点)时，是否可以通过添加刚性约束来近似a解(例如，希望纸张的变形最小)？
浏览 5提问于2016-02-10得票数 0
我正在尝试找到一种像样的方法(用渐近)来规范化右手边的常数的线性等式。例如：
x + 1 = y
变成：
x - y = -1
这是我的老生常谈的解决方案：
def canonical_linear(f):
""" canonicalise a linear equality """
lhs = f.lhs - f.rhs
const_l = [ x for x in lhs.as_ordered_terms() if x.is_number ]
if len(const_l) == 0:
const
浏览 2提问于2012-06-22得票数 1
回答已采纳
我正在用Z3库和Python语言学习数学，我不知道如何才能找到基于拉格朗日插值多项式的多项式？例如，
I have
f(0) = y0
f(1) = y1
f(2) = y2
....
我想找出f(X)
我使用Scipy工作得很好，但当我切换到Z3库时，它没有更多的文档可以找到:(
浏览 1提问于2014-04-29得票数 0
3回答
约简概念中的一个非常复杂的问题
、、、、
我已经研究了很多关于约简的知识，但我有一个很糟糕的问题:我从CLRS那里学到了这个：
“……通过把解决问题A归结为解决问题B，我们用B的”易“来证明A的”易“。”
我取自克里斯托斯·H·帕帕迪米特里写的“计算复杂性”：
“……如果B简化为A，那么问题A至少和问题B一样难。”
我对这两个概念感到困惑:当我们使用易用性时，我们说问题X简化为问题Y，如果我们对Y有多项式时间算法，并且简化过程在多项式时间内完成，那么问题X在多项式时间内可解，X比Y容易，或者至少不比Y难。
但是当我们使用硬度时，我们说问题X归结为问题Y，Y比X容易，或者至少不比X难。
我真的搞糊涂了，请帮帮我。特别感谢。
浏览 0提问于2011-08-16得票数 2
回答已采纳
我在寻找一种快速算法来求解3个未知变量上的N个多项式方程组。也就是说，给定3个函数，F0(x,y,z), F1(x,y,z)... FN(x,y,z)，我想找到x, y, z这样的F0(x,y,z) = F1(x,y,z) = ... = FN(x,y,z) = 0。
我试着在几个不同的地方找到解决方案，但我只能在代数几何学或密码学等主题上找到非常高级的论文。不过，我需要的是一个简单/快速的算法，它返回一个快速的数值解。有这样的算法吗？
浏览 7提问于2015-01-11得票数 3
回答已采纳
4回答
最近我读到了，上面写着：
一般图的匹配算法可以推广到加权情形，这似乎是最难的组合优化问题之一，可以在多项式时间内求解。
我立刻想到了以下问题：
你知道其他的“P-硬”问题吗？
就目前而言，我想将P-hard定义为：“在文献中发现了一个非常晚的多项式算法(1950年以后)。”(或者，如果在多项式时间内解决问题的确定性算法已经存在，如何才能更好地定义“硬”？)
浏览 8提问于2010-02-13得票数 10
回答已采纳
我需要了解scipy.integrate.LSODA函数的机制。
我已经写了一个测试脚本，它集成了一个简单的函数。根据函数的输入可以是rhs函数、t_min、初始y和t_max。另一方面，当我运行代码时，我什么也得不到。我该怎么办？
import scipy.integrate as integ
import numpy as np
def func(t,y):
return t
y0=np.array([1])
t_min=1
t_max=10
N_max=100
t_min2=np.linspace(t_min,t_max,N_max)
first_step=0.01
sol
浏览 37提问于2019-07-08得票数 2
回答已采纳
1回答
我有一些离散点，我可以用它们在颤动中绘制样条曲线(同步图表)。但是现在我必须找到曲线上的点，也就是说，通过给定x的值，我需要y的值，我被困在这里，没有任何算法可以应用。他们是如何使用离散点制作图形的？应该有一些算法可以在这里应用，并获得这些点。请帮帮我，提前谢谢！
浏览 15提问于2020-10-22得票数 0
回答已采纳
在Jean-Philippe Aumasson关于第46页“可证明的安全性”一段的“严肃密码学”中，它说：
可证明的安全性是关于证明破坏您的密码方案至少和解决另一个已知的困难问题一样困难。这种证明称为约简。We说，破解某些密码可以归结为 X if --任何解决问题的方法 X also都会产生一种破解密码的方法。
我想知道这是否是正确的方向。假设我们将破解密码P的问题简化为某些问题X，P \leq_m^p X ，如本书所建议的(如果我正确理解的话)。那么，如果我们有一个多项式时间的X算法，我们也有一个多项式算法的P.但这并不保证如果X没有多项式算法存在，就不应该有多项式-时间算法来破坏P。事实上
浏览 0提问于2020-04-04得票数 3
1回答
我正在写一个应用程序，其中我有三个数据点，需要推导出连接它们的抛物线。基本上，我需要这个结果来推断其他数据点。就像Excel中的趋势线一样。
有没有任何人知道的框架，我可以用来实现这一点，或者我必须手动编写求解器？
干杯,
克里斯
浏览 1提问于2012-04-16得票数 0
回答已采纳
1回答
我想要将多项式拟合到有噪声的数据中，以便近似的多项式始终是原始数据的>=。例如：
x = linspace (-2, 6);
y = (x-2).^2 + 1 + 2 * randn (size (x));
function ret = delta (P, x, y)
yP = polyval (P, x);
d = yP - y;
d (d < 0) *= 1000;
ret = sumsq (d);
endfunction
P0 = polyfit (x, y, 2);
f = @(P) delta (P, x, y);
[P, FVAL] = sqp (P
浏览 1提问于2016-04-18得票数 7
我有这样的价值观
x = 0 => y = 0
x = 1 => y = 1
x = 3 => y = 27
x = 4 => y = 64
我想用JAVA创建一个多项式函数来创建函数x^3。程序应该创建函数并显示它，如果我给出任何值，它应该计算插值的值。我已经创建了一个函数，它只是使用Aitken产生了值，但它不产生函数，并且很难理解如何做这个函数。因为我不知道如何将X值作为X放入java程序中。
浏览 2提问于2011-02-12得票数 2
1回答
我需要创建一个函数get_polynom，它将获取表示点的元组(x1, y1), (x2, y2), ..., (xn, yn)的列表，并找到多项式c0, c1, ..., cn的系数。 我无法理解这个任务，唯一的提示就是函数中提供的部分： import numpy as np
def get_polynom(coords):
...
return np.linalg.solve(a, b) 有没有人做过类似的事情？只要稍微解释一下我们所期望的将会是非常棒的！ 提前感谢！
浏览 365提问于2021-09-14得票数 1
回答已采纳
1回答
我是R的初学者，我正在尝试将一条曲线拟合到(例如)可能如下所示的数据集上：
(x- value) (y-value)
105 423
115 471
125 567
135 808
145 921.5
155 1040
X值表示刺激量，y值表示运动响应(以uV表示)。这些是10个受试者的平均值，其中每个受试者的x值是相同的。
有人告诉我，这个数据集通常遵循S型拟合。我试着用下面的代码来拟合它：
fit <- lm( y ~ poly(x, 3) )
但我不确定这是否是合适的方式：
到目前为止，我的代码如下所示：
p <- ggplot (data, aes(x, y)) +
浏览 3提问于2013-10-25得票数 5
2回答
最小二乘: Python
、、、
我正在尝试实现最小二乘：
我有：$y=\theta\omega$
最小二乘解是\omega=(θ^{T}\theta)^{-1}\θ^{T}y
我试着：
import numpy as np
def least_squares1(y, tx):
"""calculate the least squares solution."""
w = np.dot(np.linalg.inv(np.dot(tx.T,tx)), np.dot(tx.T,y))
return w
这个问题是这个方法很快就
浏览 9提问于2017-10-11得票数 0
回答已采纳
3回答
任意精度多项式求解
、、、
我们一直在用GSL求解多项式。然而，我们希望用任意精度来求解多项式。我查看了GMP和Boost多精度库，但是我找不到任何用浮点系数进行多项式求解的例程。
是否有免费和开放源码的库来解决任意精度的多项式或非常高精度的多项式(小数点后200位以上)？
随着数据类型的变化，GSL多项式求解程序是否有可能成为GMP的任意精度？
使用标准算法之一，使用GMP任意精确数据类型，编写多项式求解器更容易吗？
如果不清楚，请随时发表评论。
浏览 0提问于2016-01-10得票数 1
回答已采纳
在使用sympy (当前版本)求解多项式方程时(在这种情况下，polynom为d 2)：
from sympy import solve_poly_system
solve_poly_system(4*d**2*sin(a)**2*sin(b)/cos(b)**2 - d*cos(a) + 4, d**2*sin(a)**2*sin(b)/cos(b)**2 - d*cos(a) + 8, 3*d**2*sin(a)**2*sin(b)/cos(b)**2 - d*cos(a) + 3 ,d ,a, b)
我收到以下错误：
PolynomialError: cos(a)包含发电机集的一个元素
浏览 1提问于2016-02-07得票数 1
回答已采纳
我用Matlab符号工具箱及其求解函数，并试图求解一个非线性系统的4个方程，
有4个变量：
x1 y1 x2 y2
和4个参数
delta1 delta2 alpha beta
下面的图像描述了这些方程：
下面是Matlab代码：
syms x1 x2 y1 y2 alpha beta delta1 delta2
[x1,y1,x2,y2] = solve('delta1 * x1^alpha * y1^(1 - alpha) = (1 - x2)^alpha * (1 - y2)^(1-alpha)',...
'delta2 * x2^alpha *
浏览 4提问于2012-10-10得票数 1
回答已采纳
1回答
我已经试验了两种实现最小二乘拟合算法的方法，如所示。
第一个代码是简单的教科书方法，如Wolfram关于LSF的页面所述。第二个代码重新排列方程式，以最小化机器误差。这两个代码对我的数据产生类似的结果。我将这些结果与Matlab的p=polyfit(x，y,1)函数进行了比较，使用相关系数来衡量拟合的“优度”，并对3个例程中的每一个进行了比较。我观察到，虽然所有三种方法都产生了良好的结果，至少对我的数据来说，Matlab的例程是最适合的(其他两个例程的结果彼此相似)。
Matlab的p=polyfit(x，y,1)函数使用范德蒙矩阵、V (n×2矩阵)和QR分解来求解最小二乘问题。在Matla
浏览 6提问于2010-12-02得票数 2
回答已采纳
2回答
我有1x1024矩阵。所以我想估计一个多项式方程。
X= (0:1023)'
Y= acquired data. A 1024 element vector
然后我在MATLAB中试了一下：
polyfit(x,y,5)
但MATLAB发出的警告结果异常。
Warning: Polynomial is badly conditioned. Add points with distinct X values, reduce the degree of the ...
我不明白我做错了什么？
更新
我有很多这样的数字。
Y=
-0.0000000150
...
0.00001
...
浏览 0提问于2015-08-24得票数 1
回答已采纳
2回答
我有一条三次贝塞尔曲线，其中给出了第一个和最后一个点(即P0(0,0)和P3(1,1))。另外两个点的定义如下: cubic-bezier(0.25，0.1，0.25，1.0) (X1，Y1，X2，Y2，这些值也不能分别小于或大于0或1)
现在，假设只有一个X，我需要做什么才能得到给定X的Y坐标？(我知道在某些情况下可能会有多个值，但让我们把它们放在一边。我不是在这里做火箭科学，我只是想能够每秒多次得到Y来做转换)
我设法找到了这个：，但我不明白xTarget代表什么。
哦，这也不是什么家庭作业，我只是对互联网上没有关于三次bezier曲线的可理解的东西感到有点恼火。
浏览 0提问于2011-11-22得票数 22
回答已采纳
1回答
众所周知，网络流问题可以归结为线性规划问题。然而，在解决网络流量问题时，始终需要流量为整数。因此，我认为网络流应该简化为整数线性规划。由于ILP是NP-完全问题，网络流问题也应该是NP-完全问题.但这与我们所了解到的，因为网络流的运行时间是O(Cm)！我哪里错了？是不是因为网络流问题的运行时间是伪多项式时间，就像背包问题(Wn)？我现在很困惑！
浏览 4提问于2013-10-27得票数 0
回答已采纳
我正在尝试在scikit学习中训练一个逻辑回归模型，这需要很长的时间，大约2个小时。数据集的大小是21613 x 19。我刚开始学习scikit，因此我不知道我的代码是否错误，或者它只是需要很长时间才能训练。任何关于如何提高训练速度的建议都将不胜感激！
用于训练的代码如下
# get the LogisticRegression estimator
from sklearn.linear_model import LogisticRegression
# training the model
# apply algorithm to data using fit()
clf = Logis
浏览 7提问于2020-07-09得票数 0
2回答
我完全是一个编程初学者，所以请告诉我，我的问题的答案是否非常明显和明显。
一周前，我开始学习python，在学习了使用求解方程的基础知识之后，我想出了一段代码，可以为您提供一个立方方程的atleat (仅)1的解。这是我设计的代码：-
def deg3(a,b,c,d,g):
y=a*g**3+b*g**2+c*g+d
return y
def solvedeg3equation():
e=float(input("e= ")) #for ax^3+bx^2+cx+d=0, with maximum error of e
a=float(in
浏览 9提问于2014-01-11得票数 1
回答已采纳
1回答
求解多项式方程组
、、
我有一个系统，由2个多项式组成，2个变量，复系数。一般的情况是由有限对复数组成。
NSolve[{poly1==0,poly2==0},{x,y}]
在Mathematica中，工作于低次多项式，但寻找所有根所需的时间似乎是指数的，2^deg。有没有比NSolve更高效的替代方案？用其他语言吗？我们的目标是15-25度，越高越好。
浏览 0提问于2011-09-20得票数 0
回答已采纳
当我在训练数据集上的模型时，我在error.can下面，任何人请帮助我
在训练数据集上训练模型
logreg = LogisticRegression()
logreg.fit(x_train,y_train)
logreg =LogisticRegression (solver='lbfgs', max_iter=9000)
这会产生：
/Users/aswanthmanne/opt/anaconda3/lib/python3.8/site-packages/sklearn/linear_model/_logistic.py:763: ConvergenceWarning:
浏览 0提问于2021-02-03得票数 0
2回答
多年来，我一直在跟踪解决技术--关于如何将它们应用于一个特定的谜题--“交叉梯子”( )。
为了说明问题，我无意中发现了z3，并尝试将它用于特定的问题。我使用了Python绑定，并编写了以下内容：
$ cat laddersZ3.py
#!/usr/bin/env python
from z3 import *
a = Int('a')
b = Int('b')
c = Int('c')
d = Int('d')
e = Int('e')
f = Int('f')
solve(
a>
浏览 7提问于2013-07-18得票数 2
回答已采纳
我在我的优化算法课上发现了这个问题，完整的问题是:如果我们能证明所有容量限制为100的背包问题都可以在多项式时间内解决，那么所有的背包问题都属于P。这句话是对还是错？对齐。 通过我的书和一些研究，我得出了这样的结论:首先，KP是一个NP-完全问题。使用动态编程，它可以达到伪多项式时间，但这是不够的。荒谬的是，如果我们可以证明容量限制为100的KP可以在多项式时间内求解，那么我们可以假设KP属于P。 你对我的回答怎么看？我认为最后一句话中的荒谬之处不是那么正确。
浏览 19提问于2019-02-04得票数 0
回答已采纳
我有一个n次的多项式，y(x)=a0+a1*x^1+a2*x^2+a3*x^3+a4*x^4+......+an*x^n.，其中ai是一个实数。
我的问题是，matlab中有没有一个函数可以用来找出所有y(x)>0的x范围
浏览 0提问于2010-09-22得票数 1
回答已采纳
我试着用一个已知的y来求解x的值，我能够得到多项式来拟合我的数据，现在我想知道一个选定的y会落在曲线上的x值。
import numpy as np
x = [50, 25, 12.5, 6.25, 0.625, 0.0625, 0.01]
y = [0.00, 0.50, 0.68, 0.77, 0.79, 0.90, 1.00]
poly_coeffs = np.polyfit(x, y, 3)
f = np.poly1d(poly_coeffs)
我想做0.5 =f，然后求出x的值。
我可以在WolframAlpha中通过键入以下命令来解决此问题：
0.5 = -9.1e-6*x^
浏览 1提问于2013-05-30得票数 5
回答已采纳
3回答
我试图创建一个函数来计算包含复数的4次多项式的4个根。在我寻找公式的过程中，我发现了中包含的一个相当简单的公式，由Tito向页面底部描述。
现在，我相信真正的错误并不是代码中的真正错误(因为我肯定会让人讨厌校对)，而是我对所涉及的方法的理解。我的问题是，二次根是复的，我不知道在计算四次根时如何使用复数。
他建议用两个二次方程的根导出四次根。我试着用下面的代码尽可能地模仿这个公式。我的想法是，我计算两个二次根(前提是它们只是正数-我不知道其他情况如何)，然后，利用这些结果，我可以计算出qurtic根，然后将实际值和复数值分别保存到x1,x2,x3,x4中，放入r1,r2,r3,r4,c1,c2,
浏览 6提问于2016-01-19得票数 3
回答已采纳
1回答
我有一个困难的时间跟踪基数点，为这个遵循霍纳的规则使用不动点算法。
int16_t horner(int16_t q)
{
const int16_t c0 = 0x0033;
const int16_t c1 = 0x0055;
const int16_t c2 = 0x001C;
const int16_t c3 = 0x0048;
const int16_t c4 = 0x0200;
horner_rule = c0 + q * (c1 - q * (c2 + q * (c3 - q * c4)));
return horner
浏览 2提问于2017-11-01得票数 1
1回答
我在SNARKs上找到的所有例子都显示了3级或4级的多项式，但是5级以上的多项式又是怎样的呢?在计算根时，允许SNARK不受五次多项式不可解性影响的机制是什么？有没有我能读到的高次多项式的例子？
浏览 0提问于2022-07-13得票数 0
1回答
假设我正在尝试分析一个算法，我所能做的就是使用不同的输入来运行它。我可以将一组点(x，y)构造为(样本大小，运行时间)。我想将算法动态分类为复杂类(线性、二次、指数、对数等)。理想情况下，我可以给出一个或多或少近似行为的方程。我只是不确定最好的方法是什么。
对于任何次数的多项式，我可以创建回归曲线，并提出一些适合度的度量，但我真的不知道如何对任何非多项式函数这样做。这很难，因为我以前没有任何知识，我应该尝试适合什么形状。
这可能更像是一个数学问题，而不是编程问题，但它对我来说非常有趣。我不是一个数学家，所以可能有一种更简单的既定方法，可以从一组我不知道的点上得到一个合理的函数。有没有人有办法解
浏览 0提问于2010-10-14得票数 1
回答已采纳
1回答
参数在优化时间隔内。
、、、
通常，我使用Mathematica，但现在尝试转到python，所以这个问题可能是一个微不足道的问题，所以我对此感到抱歉。
总之，在python中是否有类似于Interval{min，max}的函数？链接是：
我想做的是，我有一个函数，我试图把它最小化，但是它是一个受限的最小化，我的意思是，函数的参数只允许在某个特定的区间内。
对于一个非常简单的例子，假设f(x)是一个参数为x的函数，我正在寻找使函数最小化的x的值，但是x被约束在间隔(min，max)<code>E 212</code>内。显然，实际的问题不是一维的，而是多维的，因此不同的参数可能有不同的区间。
因为这
浏览 3提问于2014-07-12得票数 2
我有以下问题要解决:如果有两个问题p1和p2，p2是NP-完全的，并且有一个从p1到p2的多项式约简，那么p1 .
( a) NP-硬但不一定是NP-完全
( b)可能在P中，即使P=NP
( c)是NP-完全
( d)上述任何一项
我认为( c)是正确的，但我不确定，我怎样才能证明是正确的？
浏览 4提问于2022-07-24得票数 -1
我正在解决一个复杂的非线性ODE系统，使用R作为我当前项目的一部分。在测试阶段，我通过两种方式(a)使用没有事件函数的事件函数(b)实现了使用lsoda方法的解决方案。这两种解决方案产生的结果略有不同。我不明白为什么会发生这种情况，也不明白如何减少错误？
我用这两种方法尝试了Lotka方程，但仍然想知道为什么解会有差异(见附图)。虽然差异很小，但我希望尽量减少可能对我的项目来说很大的错误。正如您所看到的，事件函数并没有改变因变量，但是仍然会产生一些错误。
rm(list=ls())
library(deSolve);
#####################################
浏览 0提问于2019-05-02得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
