Về cơ bản, tôi cần sao chép 'Hướng dẫn sử dụng để giải quyết các mô hình chu kỳ kinh doanh thực tế' của Hartley ( http://www.econ.ucdavis.edu/facemony/kdsalyer/LECTOUND/Ecn235a/Linearization/ugfinal.pdf ). Cụ thể, tôi muốn mô phỏng hệ thống động lực được ngụ ý bởi mô hình được chỉ định như sau:
Trong đó là tiêu dùng, h là cung lao động, k là vốn, z là quy trình công nghệ tự phát, y là đầu ra và i là đầu tư.
Tôi mô phỏng nó bằng logic sau: giả sử tại thời điểm , mọi thứ đều ở trạng thái ổn định và tất cả các giá trị là 0, từ đó chúng ta có k t + 1 . Sau đó, tại bằng cách gây sốc cho hệ thống thông qua , tôi giải quyết cho và (vì tôi đã 'sốc' và đã thu được trước đó . Sau đó, tôi cắm hai cái đó để lấy phần còn lại, cụ thể là - và lặp lại quy trình.ε c t + 1 h t + 1 z t + 1 k t + 1 y t + 1 , i t + 1 , k t + 2
Thật không may, tôi nhận được một quy trình bùng nổ không có ý nghĩa:
Tôi cũng bao gồm mã R được sử dụng để mô phỏng điều này:
n<-300
data.simulated <- data.table(t = 0, zval = 0, cval = 0, hval = 0, kval = 0, yval = 0, ival = 0)
data.simulated <- rbind(data.simulated, data.table(t = 1, kval = 0), fill = TRUE)
for (ii in 1:n){
##initial shocks
eps <- rnorm(1, mean = 0, sd = 0.007)
zt1 <- data.simulated[t == ii - 1, zval]*0.95 + eps
kt1 <- data.simulated[t == ii, kval]
##solve for ct, ht
lmat <- matrix(c(1, -0.54, 2.78, 1), byrow = T, ncol = 2)
rmat <- matrix(c(0.02 * kt1 + 0.44 * zt1, kt1 + 2.78 * zt1), ncol = 1)
solution <- solve(lmat, rmat)
ct1 <- solution[1, ]
ht1 <- solution[2, ]
##now solve for yt1 and kt2 and it1
yt1 <- zt1 + 0.36 * kt1 + 0.64 * ht1
kt2 <- -0.07 * ct1 + 1.01 * kt1 + 0.06 * ht1 + 0.1 * zt1
it1 <- 3.92 * yt1 - 2.92 * ct1
##add to the data.table the results
data.simulated[t == ii, c("zval", "cval", "hval", "yval", "ival") := list(zt1, ct1, ht1, yt1, it1)]
data.simulated <- rbind(data.simulated, data.table(t = ii + 1, kval = kt2), fill = TRUE)
}
a <- data.simulated[, list(t, cval, ival, yval)]
a <- data.table:::melt.data.table(a, id.vars = "t")
ggplot(data = a, aes(x = t, y = value, col = variable)) + geom_line()
Sy câu hỏi của tôi rất đơn giản - hệ thống được chỉ định trong bài viết vốn không ổn định và cho kết quả, hoặc tôi đã mắc lỗi ở đâu đó?