利用响应数据识别模态参数的子空间在线递推算法.pdf

第22卷第1期 振 动 工 程 学 报 V01．22No．1 ofVibration Feb．2009 2009年2月 Journal Engineering 利用响应数据识别模态参数的子空间在线递推算法 张家滨，陈国平 (南京航空航天大学振动工程研究所，江苏南京210016) 摘要：基于随机子空间算法，利用矩阵空间的性质和QR分解将行空问到过去行空间的投影展开为一种用于跟踪 的修改递推模式。通过子空间跟踪算法，不断跟踪计算投影的左奇异值向量，再利用最／b-乘法求出系统的模态参 数．最后用悬臂梁作为实验模型，通过不断改变系统激励的频带范围，验证方法的有效性及稳定性。结果表明，只要 选取适当的衰减系数，该方法就可以既保证一定的识别精度，又具有良好的跟踪特性． 关键词：随机子空间lPAST；子空间跟踪；模态识别 中图分类号：0324，TU311．3文献标识码：A 文章编号：1004—4523(2009)01—0026—05 为离散的状态向量；W(五)∈R。为过程噪声，主要由 引 言 噪声干扰和模型的不准确引起的，还包括未知的输 入信息；V(正)∈R“为测量噪声，主要由传感器测量 Ho和Kalman于1965年提出了参数辨识的子误差引起的，还包括未知的输入信息；k为时刻；尸 空间方法[1】，在此基础之上，很多学者研究了多种基 为输出的个数；行为系统的阶数；A，C分别为系统状 于子空间的递推算法，并将其应用于时变结构的模 态矩阵和系统输出矩阵。 态参数识别和在线模态识别。文[2，3]将输入输出之 假设过程噪声和测量噪声均为零均值的白噪 间投影关系表示为递推格式，然后利用URV跟踪 声，其协方差矩阵为 方法进行在线模态识别，文[4]在其基础之上，用投 。 影逼近子空间跟踪算法(PAST)替换URV跟踪计 E[阱％H，]一滢汁㈣ 算时变模态参数，但由于必须用到系统的激励信息， 式中 E表示求数学期望；“为Kronecker函数。 方法的实用性受到较大的限制。 环境激励是一种非稳态的激励，其统计特性随 阵(Hankel矩阵是交叉对角线上具有元素相同的矩 时间不断变化。主响应模态随激励的变化而不断改 阵)。把Hankel矩阵的行空间分成“过去”行空间和 变。以假设激励为白噪声的传统算法，识别精度很难 “将来”行空间 达到工程要求。本文通过改进随机子空间识别算法， o‘Y0 Yx Yz … 。 利用空间跟踪技术，不断代入新的响应数据，跟踪由 Yl y2 Ys … 于激励变化而导致的响应空间的变化，然后利用空 I··· … ··· ··· 牡力一 间不变性求解模态参数。该方法既适用于在线模态 吒 参数识别，也能够实现非稳态激励下的结构模态参 肛万1 I Y∞l-1‰Yi∞Yi++■1 数识别，具有较重要的工程应用价值。 I … l∞+lYi+z弘+3 ； l… … … … 蜘蜘肌～ 1 随机子空间算法[5] LYzi一1Yu y2／+1…