r - 将 MathCAD 代码转换为 R 时如何绘制极坐标图

Question

我有一个MathLAD 代码的例子。此代码的结果之一是极坐标图。

我试图将代码重写为 R（代码如下）。我的极地图是：

编辑。在迈克尔伯德评论之后，我修复了代码。

问题。有人能给出一个想法如何在整个圆圈上扩展极坐标图吗？

    K<-8
N<-6

x <- seq(-10, 10, 0.01)

f1 <- function(x){
result <- 1/N*sin(pi*x)/sin(pi/N*x)
return(result)
}

f2 <- function(x){
result <- 1-abs(f1(x-N/2))
return(result)
}

par(mfrow=c(1,2))
plot(f2(x), type="l")

phi <- seq(0, 360)
xphi <- phi/2/pi*N*K 

polar.plot(f2(xphi), 
           polar.pos = seq(0,360, length.out=length(phi)), 
           labels="", rp.type = "polygon", main="f2")