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

diudiu1

Jacobicopy=Jacobi;   
- t( y1 k' j- I! P   detaPQUcopy=detaPQU;
% D' Z( S8 }! J  E2 D1 a   detaFEcopy=detaFE;
/ O2 ]( a9 W! v7 V3 ?    for j=num:-1:2
. k7 X4 O9 }' k5 g: v3 K        if flagofPV(j)==1
+ v% m0 e5 Z5 R, v+ w; \            Jacobicopy(2*(j-1),:)=[];
            Jacobicopy(:,2*(j-1))=[];  
) b5 F; _  e1 x! k6 y- G            detaPQUcopy(2*(j-1),:)=[];
            detaFEcopy(2*(j-1),:)=[];
        end
3 ?) d' R) ?' }8 L- \    end   
   
detaPQUcopy=detaPQUcopy ;      
Jacobimodified=Jacobicopy   %Jacobi修正系数
  ^' h- |' P" DinvofJ=inv(Jacobicopy)
detaFEcopy=invofJ*detaPQUcopy

3 u3 D3 D% n) d. A6 k( {4 d( O
------------------------------------------------------------------------------------------------

Jacobimodified =
2 {0 k/ ~% c, M! @  t6 S0 d* B
& g/ D. R1 w  L* s( Z   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
* d$ b4 n8 X+ D5 n- p   -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
8 C6 k& E: J6 O; F& @+ r' e' ~         0         0         0  -31.7460         0         0         0
; `  Y% k4 f8 d$ \5 O9 c! M+ l' ?
  V8 t5 T$ f2 c# |7 YWarning: Matrix is singular to working precision.
> In yyg at 170

invofJ =

: Y: f; |6 F, q% i3 j2 v   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
( v0 m1 F9 E7 u' [   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
$ v' a% V0 _! f   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
% u* B/ F  u: l6 y1 m( o3 |; ]3 u   Inf   Inf   Inf   Inf   Inf   Inf   Inf