32得票5回答
Python中用于二维多项式的`polyfit`等效方法是什么?

我想要找到一个最小二乘的解来求解a系数z = (a0 + a1*x + a2*y + a3*x**2 + a4*x**2*y + a5*x**2*y**2 + a6*y**2 + a7*x*y**2 + a8*x*y) 给定长度为20的数组 x、y 和 z,我想要找到一个类似于 nu...

10得票2回答
使用切比雪夫多项式基函数拟合多项式。

我一直在使用Matlab中的polyfit函数,来对数据进行线性最小二乘多项式拟合。根据我的阅读,这使用标准的多项式基(单项式基)。我读到说使用Chebyshev多项式基进行拟合可以获得更高的数值稳定性,所以我想这样做。Matlab是否有这个选项?

8得票2回答
识别图表的上升趋势或下降趋势。

我正在尝试使用Python(标准线图)读取数据并将其绘制到图表上。请问如何以编程方式分类某些点在图表中是上升趋势还是下降趋势?哪种方法最优?这肯定是一个已解决的问题,存在数学方程来确定这一点吧? 以下是一些包含上升趋势和下降趋势的样本数据 x = [1,2,3,4,5,6,7,8,9,10...

15得票4回答
如何使用拉格朗日插值法计算多项式系数

我需要使用拉格朗日插值多项式计算多项式的系数,这是我的作业,我决定用Javascript完成。 以下是拉格朗日多项式(L(x))的定义: 拉格朗日基函数多项式的定义如下: 计算特定X的y值(W(x)函数)很简单,但我需要计算多项式系数([a0,a1,...,an]数组)。我需要这...

7得票2回答
我该如何创建处理多项式的函数?

我有一些关于多项式的问题,花了大约四个小时,但我就是搞不定。我是Python和编程的新手,我已经试着在纸上解决它,但我还是不知道。 编写并测试一个Python函数negate(p),用于否定由其系数列表p表示的多项式,并返回一个新的多项式(表示为列表)。换句话说,编写一个使一列数字变为负数...

18得票3回答
在Python中实现Horner方法遇到的问题

我已经编写了三种不同方法来计算多项式的代码。Horner方法应该是最快的,而朴素方法应该是最慢的,对吧?但是为什么计算时间与我的预期不符呢?有时候迭代法和朴素方法的计算时间却恰好相同。这是怎么回事呢?import numpy.random as npr import time def Hor...

8得票1回答
百斯托法的初始二次逼近

拜斯托根查找法需要非常好的初始近似值来收敛二次因子。 我尝试了各种常数、随机数、分数(通过Lin计算-a1/a2,-a0/a2)等方法,但都没有成功。 请问有没有人知道选择这些因子的好方法? 例如: 1*x^8 + 118*x^7 + 1*x^6 + 2*x^5 - 2*x^4 - 3...

15得票2回答
在Java中进行复系数多项式根查找

我正在尝试找到一种在Java中计算具有复系数的多项式根的方法(即与MATLAB中的roots()函数相当简单的等价物)。 我准备重新编码一个根查找算法,该算法构建伴随矩阵,然后使用广义特征值分解来找到根,但是为此我需要一个处理复值矩阵运算的库。 我浏览了一段时间,但似乎没有什么令人信服的选...

9得票2回答
使用Python寻找多项式根 -- 可能存在Numpy扩展Bug

我正在使用Numpy获取多项式的根。Numpy提供了一个名为“polynomial”的模块。 对于方程x^2 + 5*x + 6 = 0,我的手动计算结果是x = -2和x = -3。(简单易懂) 但是我的代码给出了错误的答案:array([-0.5, -0.33333333])(相反?)...

25得票1回答
R使用哪种正交多项式?

我正在尝试在R中匹配以下代码中的正交多项式: X <- cbind(1, poly(x = x, degree = 9)) 但是在Python中。 为了做到这一点,我实现了自己的方法来提供正交多项式: def get_hermite_poly(x,degree): #s...