动力系统模型阶次的确定.pdfVIP

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