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

diudiu1

Jacobicopy=Jacobi;   
7 z9 w1 g5 G' q) [   detaPQUcopy=detaPQU;
! k! c+ N0 r( h( A/ `   detaFEcopy=detaFE;
    for j=num:-1:2
; a) t1 y! \& H3 i0 D- X  O        if flagofPV(j)==1
9 G( G0 P8 c$ P( C. l            Jacobicopy(2*(j-1),:)=[];
7 [% v8 B2 @4 b. ^% ~/ s5 G0 a            Jacobicopy(:,2*(j-1))=[];  
            detaPQUcopy(2*(j-1),:)=[];
1 w9 ^* m2 t! o% }' P            detaFEcopy(2*(j-1),:)=[];
3 i* {5 f  g7 ]% f        end
    end   
   
& w5 l; i- V1 r& L7 I1 } detaPQUcopy=detaPQUcopy ;      
- `0 h. b$ K. P+ i, R3 B6 m Jacobimodified=Jacobicopy   %Jacobi修正系数
invofJ=inv(Jacobicopy)
detaFEcopy=invofJ*detaPQUcopy

8 n! m+ k7 D1 q! k/ g
------------------------------------------------------------------------------------------------
1 A$ Q: J+ b% v$ F" Q( J7 i# ]8 y
5 {1 u/ z2 {* X6 P, v  h- VJacobimodified =

2 X* X5 P2 @! r  B" c& H* d   56.2849    4.3929    3.1120    0.8299  -63.4921         0         0
    1.4851   77.6767   -0.8299    3.1120         0         0         0
: C0 y8 g! Y  r8 V3 N7 Q    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
* g* W1 r5 Q$ b1 N5 f% j         0         0         0  -31.7460         0         0         0
: H8 p+ t4 s: t
4 j. V/ l/ L4 N& B, k7 QWarning: Matrix is singular to working precision.
- E$ M& t+ Y* N: W5 T" P( y> In yyg at 170
8 [' Y7 l8 |5 L& I) T4 X
% D9 U" @7 b7 O; ~invofJ =

* O. I% r' s" q' M% o2 K, u   Inf   Inf   Inf   Inf   Inf   Inf   Inf
3 h. m, c9 s+ K. k% V$ Q) {   Inf   Inf   Inf   Inf   Inf   Inf   Inf
" a& b: t7 p0 r- t- t% ^   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
* d. |" n1 G. ]$ R! g( e3 w2 a   Inf   Inf   Inf   Inf   Inf   Inf   Inf
; c! E) N2 h# r8 o% W3 _% Z$ y! O   Inf   Inf   Inf   Inf   Inf   Inf   Inf