poster
2019-12-02, 22:43
我已经使用matlab构建了以下代码:
numero=60; a=zeros(numero,1) b=zeros(numero+1,1) for i=1:numero+1 a(i)=-cos(pi*(i-1)/numero) end figure plot(a,b, '*') 它用于计算阶数为的切比雪夫多项式的节点,并将其存储在称为a的向量中。我需要在python中重现此内容。我尝试解决的办法是
from mpmath import chebyt, chop, taylor import numpy as np import sympy as sp numero=60 nodes = [] for i in range(numero+2): auxiliary=-np.cos(np.pi*(i-1)/numero) nodes.append(auxiliary) [float(i) for i in nodes] nodes.sort() print(nodes) 但是,python的输出中存在问题。第一个事实是列表的第二个数字-0.9986295347545738出现了两次,它也是节点的第三个元素。我不知道为什么会这样,我想知道是否有人可以告诉我如何避免此错误。
更多&回答... (https://stackoverflow.com/q/59141054)
numero=60; a=zeros(numero,1) b=zeros(numero+1,1) for i=1:numero+1 a(i)=-cos(pi*(i-1)/numero) end figure plot(a,b, '*') 它用于计算阶数为的切比雪夫多项式的节点,并将其存储在称为a的向量中。我需要在python中重现此内容。我尝试解决的办法是
from mpmath import chebyt, chop, taylor import numpy as np import sympy as sp numero=60 nodes = [] for i in range(numero+2): auxiliary=-np.cos(np.pi*(i-1)/numero) nodes.append(auxiliary) [float(i) for i in nodes] nodes.sort() print(nodes) 但是,python的输出中存在问题。第一个事实是列表的第二个数字-0.9986295347545738出现了两次,它也是节点的第三个元素。我不知道为什么会这样,我想知道是否有人可以告诉我如何避免此错误。
更多&回答... (https://stackoverflow.com/q/59141054)