该楼层疑似违规已被系统折叠隐藏此楼查看此楼
使用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)
matlab求解含两个累加公式的方程 使用fsolve 函数求解含两个参数的多个方程问题 希望有会的人帮...