matlab求解含两个累加公式的方程,使用fsolve 函数求解含两个参数的多个方程问题，希望有会的人帮...

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

使用fsolve 函数求解四个方程，a，b是参数，t(1)_t(4)为变量，求解t(1)随a，b变化

function M= SHI_2( t )

a = 0:5;

b = 0:5;

y(1)=exp(20.37788-4.60517*t(1))+exp(30.14084+log(a)-4.60517*t(1))+exp(38.40712+2*log(a)-4.60517*t(1))+exp(45.2458+3*log(a)-4.60517*t(1))+exp(50.31148+4*log(a)-4.60517*t(1))+exp(48.99901+5*log(a)-4.60517*t(1))+exp(2.32561-2.30259*t(1))+exp(-15.45035)+exp(-44.18661+2.30259*t(1))+exp(-74.30442+4.60517*t(1))+exp(15.24311-4.60517*t(1))+exp(15.24311-4.60517*t(1))+exp(22.03574+3*log(a)-2.30259*t(1))+exp(-8.77285+2*log(a))+exp(-40.04195+log(a)+2.30259*t(1))-t(3);

y(2)=exp(log(a))+exp(22.56533+log(a)-2.30259*t(1))+exp(30.14084+log(a)-4.60517*t(1))+2*exp(38.40712+2*log(a)-4.60517*t(1))+3*exp(45.2458+3*log(a)-4.60517*t(1))+4*exp(50.31148+4*log(a)-4.60517*t(1))+5*exp(48.99901+5*log(a)-4.60517*t(1))+3*exp(22.03574+3*log(a)-2.30259*t(1))+2*exp(-8.77285+2*log(a)+exp(-40.04195+log(a)+2.30259*t(1))-b-2*t(4));

y(3)=exp(log(t(3)))-t(4);

y(4)=exp(22.56533+log(a)-2.30259*t(1))+exp(-2.30259*t(1))-exp(-32.65066+2.30259*t(1))+2*exp(20.37788-4.60517*t(1))-2*exp(log(t(2)))+2*exp(30.14084+log(a)-4.60517*t(1))+2*exp(38.40712+2*log(a)-4.60517*t(1))+2*exp(45.2458+3*log(a)-4.60517*t(1))+2*exp(50.31148+4*log(a)-4.60517*t(1))+2*exp(48.99901+5*log(a)-4.60517*t(1))+exp(2.32561-2.30259*t(1))-exp(-44.18661+2.30259*t(1))-2*exp(-74.30442+4.60517*t(1))+2*exp(15.24311-4.60517*t(1))+exp(22.03574+3*log(a)-2.30259*t(1))-exp(-40.04195+log(a)+2.30259*t(1));

M=[y(1); y(2);y(3);y(4)];

end

t0=[0;7;0;0];

p= fsolve(@SH_2,t0)