首先得敲个重点：你给出的三个方程左边的表达式完全一模一样，但右边却分别是2388、2633、2763——这是个矛盾到根本无解的方程组啊！同一个式子怎么可能同时等于三个不同的数？这就是为啥你用solve啥结果都拿不到的原因。

大概率是方程输入出了笔误 #

我猜你肯定是输入的时候写错了，比如原本三个方程应该是针对不同的变量取值计算的？举个例子，可能括号里的常数项不是都为5+z，而是不同的数值？比如改成这样（你得根据自己的实际需求调整）：

syms x y z
eqn1 = (x*(1+y*(5+z)))^(-1/2) == 2388;
eqn2 = (x*(1+y*(6+z)))^(-1/2) == 2633;
eqn3 = (x*(1+y*(7+z)))^(-1/2) == 2763;
solutions = solve([eqn1, eqn2, eqn3], [x, y, z], 'Real', true);

这种情况下，solve就能正常尝试求解符合条件的实数解了。

非线性方程组的数值求解备选 #

如果你的方程组是非线性的，符号求解solve可能效率低或者拿不到解析解，那可以试试MATLAB的fsolve数值求解函数，步骤给你列出来：

• 先把方程组转成残差形式（让等式两边相减等于0），定义匿名函数：

fun = @(vars) [
    (vars(1)*(1+vars(2)*(5+vars(3))))^(-1/2) - 2388;
    (vars(1)*(1+vars(2)*(6+vars(3))))^(-1/2) - 2633;
    (vars(1)*(1+vars(2)*(7+vars(3))))^(-1/2) - 2763;
];

• 给个初始猜测值（得根据你的问题合理设置，别瞎给）：

initial_guess = [1e-10, 1, 0]; % 示例值，你得自己调整

• 调用fsolve求解，还能开显示看迭代过程：

options = optimoptions('fsolve', 'Display', 'iter');
solution = fsolve(fun, initial_guess, options);

最后再提一嘴 #

先把你的方程组笔误修正了，确保三个方程左边不是完全一样的——这是当前问题的核心。改完之后，再根据需求选符号求解solve或者数值求解fsolve来拿a、b、c（对应代码里的x、y、z）的值就行。

内容的提问来源于stack exchange，提问作者Lorenzo Metilli