y      t
22500  1
25152  2
30685  3
34975  4
41538  5
45685  6
51780  7
55620  8
57510  9

y.fit<-nls(y~k*(a^(b*t)),start=list(a=1.26209,b=0.31729,k=3.69407),lower=c(0,0,58000),upper=c(2,2,60000),algorithm="port")

请教各位大神，龚帕兹曲线，想用r做nls求精确初值，三和值法求出了a、b、k的初始值，但是但进去执行显示Error in nlsModel(formula, mf, start, wts, upper) :
  参数初始估计值里有奇异梯度矩阵
代码是参考的别人的，请问各位大神，是哪里出了问题