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

diudiu1

Jacobicopy=Jacobi;   
   detaPQUcopy=detaPQU;
   detaFEcopy=detaFE;
: q0 H# z+ y# g& v    for j=num:-1:2
        if flagofPV(j)==1
% j. v5 u5 r+ C- }9 s            Jacobicopy(2*(j-1),:)=[];
            Jacobicopy(:,2*(j-1))=[];  
            detaPQUcopy(2*(j-1),:)=[];
            detaFEcopy(2*(j-1),:)=[];
6 r4 W' w) G1 [! N) ]        end
; @# G* A' l) c# y' r5 u    end   
2 y, ?' H! E$ A, Q3 C8 d% H& V   
detaPQUcopy=detaPQUcopy ;      
Jacobimodified=Jacobicopy   %Jacobi修正系数
* C5 ?% X# l( y, `1 g  ZinvofJ=inv(Jacobicopy)
2 S4 J9 X) t8 Y; q. P; Z detaFEcopy=invofJ*detaPQUcopy


+ s' R+ `9 C' \1 X  W7 g& i; [------------------------------------------------------------------------------------------------

  M; C6 ^6 _% ~4 ~+ t1 AJacobimodified =
+ V& M% c2 x* b' C' t' B  j' I
6 ]( G. p5 n4 k   56.2849    4.3929    3.1120    0.8299  -63.4921         0         0
6 z) d3 l: R5 _* N% S+ j    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
, d7 l' A; V# H! f' W         0         0  -31.7460         0         0   31.7460         0
7 Y4 q. V% \+ Q1 Y! b% e/ S         0         0         0  -31.7460         0         0         0

! ]* s) I8 N7 PWarning: Matrix is singular to working precision.
> In yyg at 170

invofJ =

+ i& `7 |5 s" e+ F   Inf   Inf   Inf   Inf   Inf   Inf   Inf
9 N4 u" [" ~. k8 l4 D   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
( D# o% o! _5 R8 X+ |   Inf   Inf   Inf   Inf   Inf   Inf   Inf
, ]& g; m7 v& I- ]. l# b. v! [1 v: f   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
7 j; ~. c0 a+ D& `) `   Inf   Inf   Inf   Inf   Inf   Inf   Inf