至于解方程,请去阅读论坛《 MATLAB最基础教程(超链接,自己点进去看) 》《 论坛常见问题归纳(超链接,自己点进去看) 》等帖子,或者入门书籍如《 邢树军译版的matlab教程(超链接,自己点进去看) 》,里面全有介绍和讲解
比如对第一个方程的数值求解,要写成
c=0.001;
k=0.001;
q=c.*(1-k)./(c.^2.*(1-k)-k);
sigma_x=0.01;
sigma_delta=0.001;
alpha_z=sigma_x./(sigma_x+sigma_delta);
sigma_u=0.05;
sigma_a=0.01;
sigma_m=0.025;
alpha_m=sigma_u./(sigma_u+sigma_m);
alpha_a=sigma_u./(sigma_u+sigma_m+sigma_a);
f = @( b ) -b + (-q.*(1-c.*q).*b.^2.*alpha_z.*sigma_x./(alpha_a.*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)).*sigma_u./(c.^2.*((alpha_a-alpha_m).*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)+b.^2.*sigma_x+(q.*(1-c.*q).*b.^2.*alpha_z.*sigma_x./(alpha_a.*sigma_u+(q.*b).^2.*alpha_z.*sigma_x)).^2.*(alpha_a.*sigma_u+q.*b.^2.*alpha_z.*sigma_x));
X = fzero( f, 1 )
f( X )