Advertisement
thanhqtran

rbc_log.mod

Jul 11th, 2025
428
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 1.05 KB | None | 0 0
  1. // RBC log
  2. var y c i n r w k A;
  3. varexo e;
  4.  
  5. // params
  6. parameters alpha beta sigma phi delta rhoa;
  7. alpha = 0.63;
  8. beta = 0.985;
  9. sigma = 2;
  10. phi = 1.5;
  11. delta = 0.025;
  12. rhoa = 0.95;
  13.  
  14. // model block
  15. model(linear);
  16.  
  17. #Rss = (1/beta)-(1-delta);
  18. #Wss = (1-alpha)*((alpha/Rss)^(alpha/(1-alpha)));
  19. #Yss = ((Rss/(Rss-delta*alpha))^(sigma/(sigma+phi))) *(((1-alpha)^(-phi))*((Wss)^(1+phi)))^(1/(sigma+phi));
  20. #Kss = alpha*(Yss/Rss);
  21. #Iss = delta*Kss;
  22. #Css = Yss - Iss;
  23. #Lss = (1-alpha)*(Yss/Wss);
  24.  
  25.  
  26. //1-Consumption-labor supply
  27. sigma*c + phi*n = w;
  28. //2-Euler equation
  29. (sigma/beta)*(c(+1)-c) = Rss*r(+1);
  30. //3-Law of motion of capital
  31. k = (1-delta)*k(-1) + delta*i;
  32. //4-production function
  33. y = A + alpha*k(-1) + (1-alpha)*n;
  34. //5-demand for capital
  35. r = y - k(-1);
  36. //6-demand for labor
  37. w = y - n;
  38. //7-Equilibrium condition
  39. Yss*y = Css*c + Iss*i;
  40. //8-Productivity shock
  41. A = rhoa*A(-1) + e;
  42. end;
  43.  
  44. steady(maxit=1000);
  45. check;
  46. model_diagnostics;
  47. model_info;
  48.  
  49. shocks;
  50. var e;
  51. stderr 0.01;
  52. end;
  53.  
  54. stoch_simul(order=1, irf=80, tex) y c n w r i k A;
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement