耗散耦合矩阵第二大特征值的先收缩后反幂算法

无权无向耦合动态网络所对应的耗散耦合矩阵的第二大特征值是判断网络同步稳定性的重要指标。试图改善目前一般直接调用Matlab中函数eig计算1 000阶以上耗散耦合矩阵第二大特征值复杂度高、时间长的问题,利用耗散耦合矩阵具有一个零特征值及其对应的特征向量为[1,1,…,1]T的特点,提出了一种先收缩后反幂算法,证明了收缩矩阵的特征值与原矩阵的非零特征值误差为零,其对应特征向量相等,并导出对于1 000阶以上矩阵,先收缩后反幂法所需乘法次数比Matlab中调用函数eig所用的QR算法大幅度减少。数值计算验证了理论分析的正确性。