牛顿方程组方法 MATLAB 中的方程组

Question

我曾试图实现牛顿的方法方程组，但我不认为它工作正确。我应该得到[6022.5 -60]度的角，但现在我得到[995 90 -230]度。有人能看到我的代码有什么问题吗?

hold on
l = 2;
L1 = 1;
L2 = 1;
L3 = 1;
m1 = 1;
m2 = 1;
% The system of equations to solve
%l = L1cos(u1)+L2cos(u2)+L3cos(u3);
%0 = L1sin(u1)+L2sin(u2)+L3sin(u3);
%0 = m2tan(u1)-(m1+m2)tan(u2)+m1tan(u3);
u = [pi/3 pi/8 -pi/3]
for i = 1:10
    F = [L1*cos(u(1))+L2*cos(u(2))+L3*cos(u(3))-l, L1*sin(u(1))+L2*sin(u(2))+L3*sin(u(3)), m2*tan(u(1))-(m1+m2)*tan(u(2))+m1*tan(u(3))];
    DF = [-L1*sin(u(1)) -L2*sin(u(2)) -L3*sin(u(3)); L1*cos(u(1)) L2*cos(u(2)) L3*cos(u(3)); m2*sec(u(1))^2 -(m1+m2)*(sec(u(2))^2) m1*sec(u(3))^2];
    u1 = u' + inv(DF)*F';
    u = u1';
end
rad = u1;
deg = (180/pi)*u1;

Answer 1

尝试替换以下行

u1 = u' + inv(DF)*F';

u1 = u' - inv(DF)*F';

当我做这个改变时，我得到的解是[60,0，-60]，它正确地解了三个方程。