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

diudiu1

Jacobicopy=Jacobi;   
' m  X* i+ u; [3 o& Y* h" i$ l1 e   detaPQUcopy=detaPQU;
# u- j8 j* u1 a; J) v) w5 f   detaFEcopy=detaFE;
7 K6 d1 p, L5 T    for j=num:-1:2
        if flagofPV(j)==1
            Jacobicopy(2*(j-1),:)=[];
            Jacobicopy(:,2*(j-1))=[];  
            detaPQUcopy(2*(j-1),:)=[];
            detaFEcopy(2*(j-1),:)=[];
        end
    end   
   
detaPQUcopy=detaPQUcopy ;      
6 @; g5 V( v8 m9 B  z Jacobimodified=Jacobicopy   %Jacobi修正系数
; p; m% ?( D/ d4 Y8 `$ HinvofJ=inv(Jacobicopy)
+ s7 g1 E9 ^) S, ]) m detaFEcopy=invofJ*detaPQUcopy

& d1 S3 W% s. Z! Q3 h; r$ \+ v
( N, _. t- M0 s------------------------------------------------------------------------------------------------
9 m4 d3 U5 n- q5 H) S
Jacobimodified =

: U5 s) I+ u$ U! Y5 i  l$ m2 ?/ O   56.2849    4.3929    3.1120    0.8299  -63.4921         0         0
2 p( M) i" P- E1 U  z+ n% s# ~2 \/ y2 N    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
# G+ ^( J: z$ b. H, ?  -63.4921         0         0         0   63.4921         0         0
% D( @: l) d' `9 l         0         0  -31.7460         0         0   31.7460         0
- f# Q5 K$ Y6 x( {         0         0         0  -31.7460         0         0         0
3 L% }0 c1 @& e1 A8 p6 p* o- ]
* V4 P4 {" J% h/ s# }. z! y8 T& fWarning: Matrix is singular to working precision.
; q5 |& B( t) A  Z> In yyg at 170

invofJ =

   Inf   Inf   Inf   Inf   Inf   Inf   Inf
# Q- t( E4 ]( Z   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
3 G8 S; O( L' U3 u% b2 u7 G0 s   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
& g- _( N+ r" o' G3 h; M8 q! L   Inf   Inf   Inf   Inf   Inf   Inf   Inf