Error using solve command of symbolic toolbox

조회 수: 2 (최근 30일)
Victor
Victor 2014년 4월 5일
댓글: Victor 2014년 4월 8일
I wrote a code to get a solution for the DSGE model I am working on but the S function was empty at first (when I trying to get the solution using sectors by sector) then later I had this error (when I wrote the whole equation at a go): Error using sym>checkindex (line 2429) Index must be a positive integer or logical.
Error in sym>privformatscalar (line 2376) checkindex(x);
Error in sym>privformat (line 2360) s = privformatscalar(x);
Error in sym/subsref (line 1578) [inds{k},refs{k}] = privformat(inds{k});
Please I need a genuine help. Thank you :)

채택된 답변

Walter Roberson
Walter Roberson 2014년 4월 5일
Your expression includes k(-1) and ppb(-1) which are attempts to access arrays at index negative one.
  댓글 수: 6
Walter Roberson
Walter Roberson 2014년 4월 6일
편집: Walter Roberson 2014년 4월 7일
Sometimes there just isn't a closed form solution to a problem.
When you use solve() in the form you are using, each of your numeric constants is being converted from (inherently in-exact) floating point into (exact) rational values, and then an exact solution is sought. If you want to approximate, get a numeric solution, then you need to set it up slightly differently.
Are you aware that one of your left hand sides is s-s (which would of course be 0) ?
You are attempting to solve 49 equations in 62 variables, unless some of those variables (other than coef) are being loaded with the "load" statement. You will almost never find a solution when you have more variables than equations.
Victor
Victor 2014년 4월 8일
Walter, Thank you for the help. I hope you won't give up on me. I missed out on the 50th equation but I included it. I have 62 equations with if I added the equation for the (12) shocks 'epsa' etc. The equations for the shocks are AR(1) so this may not be right ( I tried but got some wrong values). The aim is to get a balanced growth path for the equation when I write them in steady state that is why I didn't add them. Please see the attached file, I made some adjustments and commented out the last 12 equations.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Mathematics에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by