楼主如果要用MATALB,需要编2个M文件,一个用来定义函数,另一个用来估计参数。可用fminsearch命令或非线性最小二乘命令lsqnonlin()函数。
不过楼主的问题可以直接在EVIEWS中用非线性最小二乘估计,我算了一下,效果并不好。楼主的方程有点象Logistic回归,但也不是,我在EVIEWS中估计的结果如下:(只算出K1=0.032275,其他算不出来)
Dependent Variable: Y
Method: Least Squares (Gauss-Newton / Marquardt steps)
Date: 05/02/18 Time: 12:49
Sample: 1988
Included observations: 30
Failure to improve ssr (non-zero gradients) after 0 iterations
Coefficient covariance computed using outer product of gradients
WARNING: Singular covariance - coefficients are not unique
Y=C(1)*(1+C(2)*EXP(C(3)-C(4)*T))^(-1/C(2))
Coefficient Std. Error t-Statistic Prob.
C(1) 0.032275 NA NA NA
C(2) 0.000000 NA NA NA
C(3) 0.000000 NA NA NA
C(4) 0.000000 NA NA NA
R-squared -0.538757 Mean dependent var 3897.182
Adjusted R-squared -0.716306 S.D. dependent var 5400.228
S.E. of regression 7074.724 Akaike info criterion 20.69001
Sum squared resid 1.30E+09 Schwarz criterion 20.87684
Log likelihood -306.3502 Hannan-Quinn criter. 20.74978
Durbin-Watson stat 0.021762