动力系统模型阶次的确定 ORDER IDENTIFICATION OF DYNAMIC SYSTEM MODEL.pdfVIP

振动与冲击 ANDSHOCK 第27卷第1I期 JOURNAl。OFVIBRATION 动力系统模型阶次的确定 易伟建，刘翔 (湖南大学土木工程学院，长沙410082) 摘 要：为了判定状态空间方程的阶次，提出了一种新的动力系统模型定阶方法一残差期望法。该方法利用系统 的脉冲响应函数建立Hankel矩阵，对其进行奇异值分解后，得到q和鸩两个矩阵。％和旦分别为信号的主要组分和噪 声。假设系统受到白噪声干扰，并利用白噪声的零均值特性，可以准确的辨识模型的阶次。在Hankel矩阵的大小满足一 定要求的条件下，即使系统存在较大噪声，也能正确识别模型的阶次。对系统有较大噪声干扰的情况进行了数字仿真研 究，并对一钢筋混凝土框架结构的实测数据进行了对比研究。 关键词：状态空间模型；阶次辩识；残差期望；残差期望比 中图分类号：TU317+．1文献标识码：A 1985年，J．N．Juang等人提出了特征系统实现算用响应信号奇异性检测的模型阶次辨识方法。Unbe． 法(ERA)¨-31，由于它具有较高的识别精度与识别密hauen¨引则介绍了模型阶次辨识的零、极点相消法。谢 集模态的能力，很快得到发展与应用。随后出现的子 新民等¨到则利用系统数据乘积矩阵是否奇异来确定差 分方程描述多变量系统的结构参数。樊江玲¨41等提出 空间识别法(SIM)，随机子空间法(SSI)，也都具有 ERA的优点，有些方面甚至超过ERA【4J。在ERA与使用信号特征分量的能量指标来估计模型阶次，同时 SIM、SSI等算法中，都采用了奇异值分解(SVD)来求解给出一种特征频率随拟合数据变化的稳定图。 系统的特征矩阵，这样在运算过程中都会遇到确定模 根据Hankel矩阵的特性，本文提出一种新的系统 型阶次的问题。以往的时域模态识别方法及软件中， 模型阶次的判定方法，称之为残差期望法或残差期望 大多采用变化相关参数得到的稳定图来确定模型的阶 比法，它的计算量较小，而且判断准确。模拟数据与实 次¨’6J。这种方法在频谱图上标示出满足一定条件的 测数据的计算结果表明本文的方法可以有效的辨识系 稳定极点，并被认为是系统的真实极点。但这种方法 统模型的阶次。 不能完全排除噪声模态，特别是随着模型阶数的升高， 1理论分析 一些拟合模型的拟合模态往往容易趋于稳定，用稳定 图很难完全并正确确定模型阶次。理论上，对系统矩 1．1 Hankel矩阵的建立与SVD分解 阵进行SVD分解后，去掉奇异值中的零元素，所剩奇异 对于离散线性时不变系统，可用状态方程表示为： 值个数就是模型的阶次。例如，四层框架结构，按每一 ” 一’ 、7 楼层一个自由度进行测试，系统明显应该是一个四阶 【Y(k)=瓯(后) E 模型。由于奇异值共轭成对出现，这个4阶模型应有4 其中A R8。“，B∈R”。“，C∈R‘。4，龙∈R4，“∈R“，Y∈ 对8个非零奇异值，其余均为零。但由于噪声的影响， 足，A