李雅普诺夫稳定性分析.ppt.pptVIP

对称矩阵的定号性(正定性)的判定(3/12) function sym_P=posit_def(P) [m,n]=size(P); if n1 for i=1:n-1 for j=i:n dia_v(j)=abs(P(j,j)); end [mindv,imin]=max(dia_v(i:n)); imin=imin+i-1; if mindv0 if imin i a=P(imin,:); P(imin,:)=P(i,:); P(i,:)=a; b=P(:,imin); P(:,imin)=P(:,i); P(:,i)=b; end for j=i+1:n x=P(i,j)/P(i,i); P(:,j)=P(:,j)-P(:,i)*x; P(j,:)=P(j,:)-P(i,:)*x; end end end end % 定义函数posit_def() % 取P矩阵的维数大小n % 若n1,则对P进行合同变换 % 对非对角线元素进行消元 % 取未消元的对角线绝对值 % 求对角线绝对值的最大者 % 将对角线绝对值最大值所在的行和列与当前行