Question

如何用matlab解下面的方程组

y=2.8761*(x^(-5.988)-1)^(0.334)
x=0.91937-0.65156*e^(-715.49*y)
用solve
[x,y]=solve('y=2.8761*(x^(-5.988)-1)^0.334','x=0.91937-0.65156*e^(-715.49*y)')
最后解的x =

[ empty sym ]

y =

[]
如何用matlab求出真正的解。
非常急。！！

Answer 1

e^()应该为exp()，不过，还是没有解析解，用数值方法。

function hh
ezplot('2.8761*(x^(-5.988)-1)^0.334-y')
hold on
ezplot('0.91937-0.65156*exp(-715.49*y)-x')
grid
X=fsolve(@fun,[1 1])
function f=fun(X)
x=X(1);y=X(2);
f(1)=2.8761*(x^(-5.988)-1)^0.334-y;
f(2)=0.91937-0.65156*exp(-715.49*y)-x;

结果：
Optimization terminated: first-order optimality is less than options.TolFun.

X =

0.91937 2.4962