matlab fsolve解方程结果不对？

Question

fsolve函数求解出的结果带回原来的表达式，结果不正确。function F=root8d(x)%归一化载荷。F(1) = 12*x(2)-6*0.1046 + -0.0997*(1.2*x(2)-0.1*0.1046) + (-0.0997)^2*(-x(2)/700+0.1046/1400) - 0.0105;F(2)= -6*x(2)+4*0.1046 + -0.0997*(-0.1*x(2)+2*0.1046/15) + (-0.0997)^2*(x(2)/1400-11*0.1046/6300) - x(1);F(3)= (x(4))^2*(-0.0997/12) - (0.6*(x(2))^2-0.1*0.1046*x(2)+(0.1046)^2/15)-(-0.0997)*(-(x(2))^2/700+0.1046*x(2)/700-11*(0.1046)^2/6300)-x(3);F(4)=(1+x(3))*0.0105-(-0.0997)*x(2)+x(1)-0.125;把F（3）等号右边除了-x(3)部分带进去运算得到的结果和x(3)不同

Answer 1

首先，我们要确保fsolve函数是否收敛到正确的解。可以通过观察求解过程中的退出标志（exitflag）以及确保使用的初始值是合理的。请在调用fsolve时返回exitflag，如下所示：
options = optimoptions('fsolve','Display','iter');
[x, fval, exitflag, output] = fsolve(@root8d, x0, options);
这里，x0 是你提供给fsolve的初始值。如果退出标志（exitflag）是正数，那么求解过程应该是收敛的。如果不是，你可能需要改变初始值或者调整其他求解选项。也可以尝试将求解器显示设置为“iter”（迭代），以便在求解过程中观察详细信息。
接下来，确保你在计算F(3)的右侧时使用了正确的符号和公式。这里是将F(3)的右侧计算结果与x(3)进行比较的代码：
F3_right_side = (x(4))^2*(-0.0997/12) - (0.6*(x(2))^2-0.1*0.1046*x(2)+(0.1046)^2/15)-(-0.0997)*(-(x(2))^2/700+0.1046*x(2)/700-11*(0.1046)^2/6300);
difference = F3_right_side - x(3);
disp(['Difference between F3_right_side and x(3): ', num2str(difference)]);
如果difference值较小（接近零），则说明求解过程是正确的。如果difference值较大，可能需要检查方程是否存在误差或求解过程中是否存在问题。