R中的多项式结

Question

这个问题是关于以图形方式表示由R中的表达式here和here给出的结多项式的可能方法。

不幸的是，由于我的数学知识不太好，我没有走得太远。这就是我试过的

phi <- seq(-30, 30, len=1001) 
fx <- phi^5 - 28*phi  #these are the polynomials from the website above
fy <- phi^7 - 32*phi^3

，在使用情节时给出（fx，fy）

我想它应该给出八字结。我可以在指定的时间间隔内找到上述函数的零，这一切似乎都没问题，因为两个函数都有一对非常相似的零值，如结图中所示。有任何想法吗？

Answer 1

您的代码是正确的。你只需要放大：

phi <- seq(-2.6, 2.6, len=1e5) 
fx <- phi^5 - 28*phi  
fy <- phi^7 - 32*phi^3

plot(fx, fy, type="l")

Answer 2

尝试：

phi <- seq(-3, 3, len=1001) 
fx <- phi^5 - 28*phi  
fy <- phi^7 - 32*phi^3
plot(fx,fy,xlim=c(-50,50),ylim=c(-150,150))

您必须调整轴的极限。 HF