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

diudiu1

Jacobicopy=Jacobi;   
   detaPQUcopy=detaPQU;
   detaFEcopy=detaFE;
, a- {! K% l/ p! q  G1 Y4 G    for j=num:-1:2
        if flagofPV(j)==1
3 g+ B$ {) n+ n( ^7 c" O$ [8 w            Jacobicopy(2*(j-1),:)=[];
            Jacobicopy(:,2*(j-1))=[];  
            detaPQUcopy(2*(j-1),:)=[];
6 S! x7 U) A& y7 K6 ?7 z            detaFEcopy(2*(j-1),:)=[];
/ V" f4 B' x- O5 d; n7 ]( n        end
5 r* p0 V& y8 Y: n5 ]0 `    end   
6 \+ W+ r% z9 }4 u   
& i( x: H  d8 k0 g: d5 b6 I2 Z detaPQUcopy=detaPQUcopy ;      
Jacobimodified=Jacobicopy   %Jacobi修正系数
: g* b: |! Y7 u; b  j" m4 pinvofJ=inv(Jacobicopy)
detaFEcopy=invofJ*detaPQUcopy

- X, S% [8 K2 X  s8 J. m
------------------------------------------------------------------------------------------------
4 j8 `7 ]' U$ \2 k+ {3 D
Jacobimodified =

   56.2849    4.3929    3.1120    0.8299  -63.4921         0         0
1 r7 I2 \. r/ s% d+ v; l    1.4851   77.6767   -0.8299    3.1120         0         0         0
. r' L5 l0 \% T" E) P  q$ x! N' G    3.1120    0.8299   25.8604    4.7915         0  -31.7460         0
  Q0 |# q% e9 q- X0 r( _   -0.8299    3.1120    1.6223   45.6152         0         0         0
: s  u4 \( t# g  -63.4921         0         0         0   63.4921         0         0
* M: j8 U# w# F) x1 U& F. A- Q         0         0  -31.7460         0         0   31.7460         0
# S3 q- s8 T7 @: W% ]4 w, ?         0         0         0  -31.7460         0         0         0

& C9 u/ J7 j* Y$ o  m- AWarning: Matrix is singular to working precision.
  F" i& I! |3 l$ _3 V( A/ a/ _+ S> In yyg at 170
7 ?8 ]) l3 P. i5 ^
) P" T. p4 f9 r9 e& ~invofJ =
. ?0 F, M8 r* l/ ]1 O# w" @: S' }2 u
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
' w, x( p. r7 g0 V   Inf   Inf   Inf   Inf   Inf   Inf   Inf
   Inf   Inf   Inf   Inf   Inf   Inf   Inf
- O# ?$ Y/ D% B   Inf   Inf   Inf   Inf   Inf   Inf   Inf
- B. v# ~  O5 F! L% B" ?   Inf   Inf   Inf   Inf   Inf   Inf   Inf
* [4 v5 L# {  @* g9 S4 J0 I, \2 b# c   Inf   Inf   Inf   Inf   Inf   Inf   Inf
" A8 n# D/ X) z& D2 q- @1 G   Inf   Inf   Inf   Inf   Inf   Inf   Inf