捕食者能以Holling'；s IV型功能性反应，是否在R中生成？

利用R，我试图用捕食者-食饵等参线图分析捕食者-食饵相互作用，以及Holling的IV型功能反应

这是我正在使用的脚本：

library(deSolve)
library(phaseR)
model.LV <- function(t, y, parameters){
  with(as.list(parameters),{
    N<-y[1] 
    P<-y[2]
    dN <- r*N*(1-N/K)-(1.511*N/((N^2/9.2724)+N+1.483))*P*N
    dP <- -b*P+(1.511*N/((N^2/9.2724)+N+1.483))*a*P*N

    list(c(dN,dP))
  })
}

params.LV<-c(r=0.8,K=10,b=4.2,a=6)

data.LV<-as.data.frame(lsoda(c(N=3,P=3),seq(1,100,by=0.5), model.LV, params.LV))

plot(data.LV$N, data.LV$P, type="l", col="blue", main="Trajectory for Type-IV Functional Response",
     xlab="Prey density N", ylab="Predator density P", xlim=c(0,5), ylim=c(0,5))
nullclines(model.LV, x.lim=c(0,5),y.lim=c(0,5), parameters=params.LV, system="two.dim", colour=c("green","red"), add=TRUE)

库（deSolve）
图书馆（相位器）
model.LV这给出的是轨迹，而不是等高线的位置。通过（1）求解常微分方程的平衡点，（2）绘制由这些平衡点定义的曲线，可以得到等线。如果没有，则可以得到类似捕食者-食饵系统（不稳定焦点）的行为。你（可能）正确地模拟了动力学，但你问题的关键词是分析。要做到这一点，您需要做一些数学运算，然后绘制分析公式（您可以使用
曲线
绘制），而不是模拟。有关模拟此类系统的更多详细信息，请参阅@jaimedash谢谢。我确实成功地使模拟工作，并且正在跟踪您提供的链接。干杯