首先，将dynare 模型存为.mod文件；其次，你的beta 可能跟Matlab自身所设定的beta函数有冲突。你先试试吧。第二条和end的错误，我觉得是因为你在.m环境下造成的。

祝你好运。

你的模型有问题

STEADY:  The Jacobian contains Inf or NaN. The problem arises from:

STEADY:  Derivative of Equation 5 with respect to Variable y  (initial value of y: 0)
STEADY:  Derivative of Equation 6 with respect to Variable y  (initial value of y: 0)
STEADY:  Derivative of Equation 8 with respect to Variable y  (initial value of y: 0)
STEADY:  Derivative of Equation 1 with respect to Variable c  (initial value of c: 0)
STEADY:  Derivative of Equation 2 with respect to Variable c  (initial value of c: 0)
STEADY:  Derivative of Equation 4 with respect to Variable k  (initial value of k: 0)
STEADY:  Derivative of Equation 6 with respect to Variable k  (initial value of k: 0)
STEADY:  Derivative of Equation 4 with respect to Variable l  (initial value of l: 0)
STEADY:  Derivative of Equation 5 with respect to Variable l  (initial value of l: 0)
STEADY:  Derivative of Equation 8 with respect to Variable l  (initial value of l: 0)
STEADY:  Derivative of Equation 2 with respect to Variable w  (initial value of w: 0)
STEADY:  Derivative of Equation 1 with respect to Variable r  (initial value of r: 0)
STEADY:  Derivative of Equation 4 with respect to Variable z  (initial value of z: 0)

STEADY:  The problem most often occurs, because a variable with
STEADY:  exponent smaller than 1 has been initialized to 0. Taking the derivative
STEADY:  and evaluating it at the steady state then results in a division by 0.
Error using dynare_solve (line 60)
An element of the Jacobian is not finite or NaN

Error in evaluate_steady_state (line 66)
            [ys,check] = dynare_solve([M.fname '_static'],...

Error in resol (line 104)
[dr.ys,M.params,info] = evaluate_steady_state(oo.steady_state,M,options,oo,0);

Error in check (line 73)
[dr,info,M,options,oo] = resol(1,M,options,oo);

Error in test (line 133)
oo_.dr.eigval = check(M_,options_,oo_);

Error in dynare (line 180)
evalin('base',fname) ;