N-R法求解雅可比矩阵出错 Matrix is singular to working precision. 矩阵奇异
Jacobicopy=Jacobi;detaPQUcopy=detaPQU;
detaFEcopy=detaFE;
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 ;
Jacobimodified=Jacobicopy %Jacobi修正系数
invofJ=inv(Jacobicopy)
detaFEcopy=invofJ*detaPQUcopy
------------------------------------------------------------------------------------------------
Jacobimodified =
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
-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
0 0 0-31.7460 0 0 0
Warning: Matrix is singular to working precision.
> In yyg at 170
invofJ =
Inf Inf Inf Inf Inf Inf Inf
Inf Inf Inf Inf Inf Inf Inf
Inf Inf Inf Inf Inf Inf Inf
Inf Inf Inf Inf Inf Inf Inf
Inf Inf Inf Inf Inf Inf Inf
Inf Inf Inf Inf Inf Inf Inf
Inf Inf Inf Inf Inf Inf Inf
页:
[1]
