N-R法求解雅可比矩阵出错 Matrix is singular to working precision. 矩阵奇异

diudiu1

Jacobicopy=Jacobi;   
/ g$ P5 K; V( Y1 ^' [   detaPQUcopy=detaPQU;
   detaFEcopy=detaFE;
    for j=num:-1:2
( l9 {1 s5 F1 ~* D# w) \2 ^$ t8 M        if flagofPV(j)==1
            Jacobicopy(2*(j-1),:)=[];
9 ^, c& }% a" L6 ^5 E3 O            Jacobicopy(:,2*(j-1))=[];  
            detaPQUcopy(2*(j-1),:)=[];
            detaFEcopy(2*(j-1),:)=[];
        end
: |6 m. K" y4 ~$ \3 }; {    end   
" V+ u) D0 ?# d0 V. R6 {   
detaPQUcopy=detaPQUcopy ;      
Jacobimodified=Jacobicopy   %Jacobi修正系数
invofJ=inv(Jacobicopy)
detaFEcopy=invofJ*detaPQUcopy
0 O1 y3 t* h$ I
0 Q1 P* U& v% l, O
------------------------------------------------------------------------------------------------
6 j8 X1 j* d8 M$ n: ]
Jacobimodified =
5 C8 t: `5 C$ ^* R6 a# R) w1 R
6 x0 N6 Q1 g1 `  a/ I' G; p   56.2849    4.3929    3.1120    0.8299  -63.4921         0         0
    1.4851   77.6767   -0.8299    3.1120         0         0         0
    3.1120    0.8299   25.8604    4.7915         0  -31.7460         0
   -0.8299    3.1120    1.6223   45.6152         0         0         0
  -63.4921         0         0         0   63.4921         0         0
         0         0  -31.7460         0         0   31.7460         0
         0         0         0  -31.7460         0         0         0

: n' ?) b8 M+ F* ~# }8 G! aWarning: Matrix is singular to working precision.
> In yyg at 170
( [' K9 v3 p! C8 f  \' c/ d# B
5 \. B3 @% [: V4 A) O! RinvofJ =

   Inf   Inf   Inf   Inf   Inf   Inf   Inf
  _" D6 D) x$ i; c. G1 G* _   Inf   Inf   Inf   Inf   Inf   Inf   Inf
4 l, ?# w6 x. T+ g) R2 B( l# p4 m   Inf   Inf   Inf   Inf   Inf   Inf   Inf
) G% F3 B! A  w' H. W+ |& y   Inf   Inf   Inf   Inf   Inf   Inf   Inf
! h) \, r* m7 Q' R   Inf   Inf   Inf   Inf   Inf   Inf   Inf
: N' M/ I0 U5 i/ m+ s& F   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf