我建了一个模型，最后在求解稳态时需要解一个方程组（16个未知数、16个方程）。但是我用matlab中的fslove来解该方程组时有两个问题：（1）得到的结果却是找不到解，或者function value与0差距很大，如何能得到该方程组的解？（2）某些变量可以通过事前决定（如r_bar，可根据beta、delta确定)，但是如果放在方程组中来确定，与事前决定的值差距很大，这是为何？

求大家帮助解决以上问题，谢谢！！

方程组原式和我自己编写的matlab代码已放置附件（EquationsSolve.m中是所有变量根据外生参数求解，EquationsSolve2c.m中将某些变量先求解后作为外生参数）