MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   资料存档 (https://www.labfans.com/bbs/forumdisplay.php?f=72)
-   -   问题在涉及for循环的python中再现matlab代码 (https://www.labfans.com/bbs/showthread.php?t=22344)

poster 2019-12-02 22:43

问题在涉及for循环的python中再现matlab代码
 
我已经使用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出现了两次,它也是节点的第三个元素。我不知道为什么会这样,我想知道是否有人可以告诉我如何避免此错误。



[url=https://stackoverflow.com/q/59141054]更多&回答...[/url]


所有时间均为北京时间。现在的时间是 23:20

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.