GM(1,1)

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

x <- c(1618,1674,1728,1753,1775,1785,1797,1815)
n<-length(x)
x1<-numeric(n);
for(i in 1:n)   ##一次累加
{
  x1[i]<-sum(x[1:i]);
}
z = diff(x1) #
m = n-3
v_m = c()
for(i in 1:m){
  y=c()
  for(j in 1:(n-i-1)){
    y = c(y,z[j+i]-z[j])
  }
  for(k in 1:(n-i-2)){
    v_m = c(v_m,log(y[k+1]/y[k]))
  }
}
v = sum(v_m)/length(v_m)
v
A = matrix(c(exp(seq(1,n)*v),seq(1,n),rep(1,n)),nrow = n, ncol = 3)
c<-solve(t(A)%*%A)%*%t(A)%*%x1;
c

# 拟合
x_1 = c()
for(i in 1:n){
  x_1 = c(x_1,c[1]*exp(v*i)+c[2]*i+c[3])
}
x_0 = c(x_1[1],diff(x_1))
x_0

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

这是啥呀