急！求MATLAB解多元非线性方程组，最好给出程序，求大神

Question

[w,x,y,z]=solve('exp((1193*y)/20 - (543*y^2)/5 - (3139*y^(1/2))/100 + 2867/1000) - w=0', 'exp((1193*z)/20 - (543*z^2)/5 - (3139*z^(1/2))/100 + 2867/1000) - x =0', '- (w - 751/5000)/(x - 751/5000) - 8994/625 =0', '- (y - 57/625)/(z - 57/625) - 8994/625=0')

Answer 1

该多元非线性方程组不适宜用solve()函数命令来求，可以用fsolve()函数命令来求其在某区间内的数值解。求解过程如下：

func=@(x)[exp((1193*x(3))/20-(543*x(3)^2)/5-(3139*x(3)^(1/2))/100+2867/1000)-x(1),exp((1193*x(4))/20-(543*x(4)^2)/5-(3139*x(4)^(1/2))/100+2867/1000)-x(2),-(x(1)-751/5000)/(x(2)-751/5000)-8994/625,-(x(3)-57/625)/(x(4)-57/625)-8994/625]; %自定义函数

x0=[0.40181 0.075967 0.23992 0.12332]; %初值

x= fsolve(@(x)func(x),x0); %方程组的解

w=x(1),y=x(3),z=x(4),x=x(2) 

运行结果