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

diudiu1

Jacobicopy=Jacobi;   
   detaPQUcopy=detaPQU;
0 t2 r2 B% A/ L- ?. o0 i   detaFEcopy=detaFE;
9 g- N4 ^+ ^( q+ ~. y    for j=num:-1:2
        if flagofPV(j)==1
) V  }5 f; Z* x  ~" e            Jacobicopy(2*(j-1),:)=[];
. S5 j% ^) c2 s& j) u            Jacobicopy(:,2*(j-1))=[];  
, Q) Y, y  N6 @2 L" L) J            detaPQUcopy(2*(j-1),:)=[];
            detaFEcopy(2*(j-1),:)=[];
        end
; b( [2 e1 _2 K    end   
  r" u8 ^( s6 q& P   
detaPQUcopy=detaPQUcopy ;      
. q# Y6 k0 z% ? Jacobimodified=Jacobicopy   %Jacobi修正系数
invofJ=inv(Jacobicopy)
. p# a! J( ?. a" Y$ h detaFEcopy=invofJ*detaPQUcopy

! J7 j/ {- W$ Y; f
------------------------------------------------------------------------------------------------

# a% H. k5 r9 Y& TJacobimodified =

   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
) a6 o7 q; K# m, y& S" y         0         0         0  -31.7460         0         0         0

. u# `2 t& J$ \, S  D( h" ^3 PWarning: Matrix is singular to working precision.
4 B( h! f2 E' s5 K, C1 n> In yyg at 170

invofJ =
! |# s3 n! l6 t8 e- K
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
* o3 o/ c, X, P& j# ?* X4 M   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
( E( Q% L) H' o  e   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
0 f/ O# Q4 g0 O' m8 g   Inf   Inf   Inf   Inf   Inf   Inf   Inf