最小二乘复频域法PolyMa

1、最小二乘复频域法（ PolyMax）SX1201069 虞刚PolyMax模态识别方法，属于多自由度时域识别法，也称作多参考点最小二乘复频域法 （ Polyreference least squares complex frequency domain method）， 是最小二乘复频域法（LSCF）的多输入形式，是一种对极点和模态参预因子进行整体估 计的多自由度法 ,一般首先通过实验建立稳态图， 以判定真实的模态频率、阻尼和参 预因子；建立可以线性化的直交矩阵分式模型，然后基于正则方程缩减最小二乘问 题 , 得到压缩正则方程 , 于是模态参数可以通过求解最小二乘问题得到。该方法集合 了多参考

2、点法和LSCF方法的优点，可以得出非常清晰的稳态图，并且密集空间可以被 分离出来，尤其在模态较密集的系统（动力总成系统），或者FRF数据受到严重噪声污 染的情况下仍可以建立清晰的稳态图，识别出高度密集的模态，对每一个模态的频 率、阻尼和振型都有很好的识别精度 , 是国际最新发展并流行的基于传递函数的模态 分析方法。其基本思想如下：（ 1）建立频率响应函数模型多参考点最小二乘复频域识别技术（PRLSCF或 PolyMAX要以频响函数矩阵作为识别的初始数据，其数学模型采用右矩阵分式模型来描述。在频域中，系统输出o（o 1,2 ,No，其中No为输出点数）和全部输入的关系可用右矩阵分式模型（RMFD

3、来描述，右矩阵分式模型的表达式为1H o U o D 1（ 1）式中：H。ClNi理论频响函数的第o行，M是输入点数，即激励数；UoC1叫一分子多项式行向量;DoC叫Ni分母多项式矩阵。且U。 和D。 可以表示成如下形式:NUoZrBor( 01,2 ,No)( 2)r 0NDo乙 Ar( 3)r 0其中分母系数矩阵ArRN和分子系数行向量BorR1 Ni是待估计的参数。所有这些系数合并为一个矩阵(4)其中AoBo0BoNR N 1 NioRNi N 1Ni(5)Ai式（2）和式（3）中出现的多项式基函数 乙般地，有以下两种选择:式中：N 多项式阶次i .对于连续时域模型，可取为Zr(6) 式

4、中：s - 比一缩放因子，用来提高方程的数值状况。2ii.对于离散时域模型，可取为Zre j Trs（7）式中：Ts 采样周期。通常采用离散时域模型。（2）参数的线性化通过试验测量出的频率响应函数矩阵i fcNo Ni，用Ho f c1 Ni表示实测频响矩阵的第o行，o 1,2, ,N0, f 1,2, ,Nf，那么关于参数矩阵的非线性最小二乘（NLS）目标函数口表示为NoNfH丄NLSHNLS/ Onlstrof ,of ,（8）o 1 f 1式中：?H 矩阵的复共扼转置；tr ?矩阵的迹，即矩阵的主对角元素之和。通过对式（8）求极小值，便可以得到频率响应函数矩阵的右分式矩阵模型各系数的估计

5、值，即矩阵的估计值。式（8）中的加权非线性最小二乘误差函数被定义：NLS1:f, Wof Ho f, Hof Wof Uo f, o D 1f, Hof上式中Wo f是一个加权函数。一般地，为了提高估计的质量，我们采用Wovar HoHo式中：var?方差，可用相关函数求取。也可使用公式Wo1var Ho f(11)来做加权函数的。这两种加权函数都考虑了测量频响函数数据的好坏：测得频响的 方差越小，对目标函数的贡献越大。非线性误差函数可以经过一个近似的处理为一个线性的问题。实际上，通过对肿S f,右乘D f,，则可以得到一个关于参数为线性的方程，此加权线性最小这样式（12）关于参数为线性，将所

6、有频率点装配成一列,二乘（LS）方程误差：S f,为LSof ,NLSof ,Df ,1Wof Uo f, o Df ,H%fDf,（12）NH%ArWofZrfBorZrfr 0f 1,2,L ,Nf，它可用矩阵形式来表示LSoLSoXo YoJo(13)LS其中：Wo 1 zo 1,z11 ,L , zN1XoWoN fzoN fM,z1Nf ,L , zNN fC Nf N 1（14）Wo1 zo 1 , z1 1 ,L, zN1H%o 1YoMCN f Ni N 1（15）WoN fzoNf , z1N f ,L, zNN fH%o N f式中， Kronecker 积。3）缩减标准方

7、程加权线性最小二乘估计表达式为式中：oLSH LSLStrooo1NoTTRoSootro1oSoTTo16)Ro Re XoH XoRN1N1同时，目标函数（ 16）等价于17)THl LS tr Re J J式中， J 是 Jacobian 矩阵，被如下定义X10 L0Y10JX2 LM丫2CN°Nf N 1 No Ni(18)MM OMM00LXno丫叫为使1 LS值最小，将Ils对系数矩阵o和求导，并令其为零1 LSo2 Ro o So0o 1,2,L ,No(19)1 LSN°2S： oo 1To0(20)由式（19）得到。Ro1So，把它代入式（20）得No2T

8、oS：Ro1Soo 1M0(21)其中，MNo2ToS；Ro1SorM N 1 Ni N10o 1由式（19）和（20）得到标准方程，经过整理，此标准方程的表达式为R10L0S0R2L0S21MMOMM22M2Re JH J00LRNoSNoNoNooS：SLSN0Too 10(22)式（21）即为“缩减”标准方程，其中矩阵M维数为Ni N 1 Ni N 1，比标准方程式（22）中的Re JHJ的维数N。M N 1 N。M N 1要小的多。4）求解缩减标准方程通过求解“缩减”标准方程，便可得到分母系数矩阵根据线性方程组的求解理论，先对系数矩阵 施加一个约束。假如，设定系数矩阵中的一个系数矩阵块

9、等于正则常数矩阵（例如设系数矩阵的最后一个矩阵块N1I Ni ） , 在这种前提下，缩减标准方程变为23)其中A M 1: Ni N,1: NiN系数矩阵的最小二乘估计为?LSXI NiA 1B24)一旦求得了?LS ，那么通过Ro1So就可得到所有的分子系数 ?LS,这种方法考虑了标准方程的结构特性 ，比直接求解方程 （22）要快得多 。确定了分母系数矩阵 后 通过求解 的伴随矩阵的特征值和特征向量，这样就可以得到了系统的极点和相应 的模态参与因子。方程如下0IL0000L00MMOMM V V.（25）00L0IA0TA1TLAN 2AN 1上式中，V,CNoN NoN，矩阵V的最后Ni行

10、就是模态参与因子；对角阵的角元记 录为i（i 1,2丄，N°N）由不稳定的数学极点和稳定的物理结构点两部分组成。记稳定的物 理结构极点为r e rTs，通过对这些物理结构极进行转换，便可得出结构的固有频 率r和模态阻尼比r ；关系式如下(26)（5）计算频率点和阻尼比点根据信号与系统基本理论中对系统稳定性的描述：系统的全部极点落于s域左半 平面（不包括虚轴），且满足有界输入有界输出原则，系统是稳定的。复特征矩阵 中的复特征值总是以共轭对的形式出现，同时也包含实数（虚轴上），在求解频率 点i和阻尼比点i时，对于每个共轭对只取其中一个进行分析，且不考虑实数。复特征矩阵 中的对角元i e江

11、，由式(26)， i用Re i i Im i描述，则Re i i Im i e ii %TseiTse i °%TseiTscos%Tsi sin%Ts(27)Re ii Imie iTs(28)Imi亠，、arctan%s( 29)Re i所以丄TsIn(30)% arctan m-(31)TsRe i由此可求得频率i和阻尼比ii(32)ii在求得的频率i和阻尼比i包含有结构的固有频率r和模态阻尼比r，因此，必须对所有求得的i和i进行有效的分析和选取，以确定系统真实的固有频率和阻尼 比。建立稳态图就是一种行之有效的方法。(6)建立稳态图在模态分析中，稳态图是帮助实验者分离结构物理极点和数学极点的一个有力工具，如图1所示。通过逐渐增大多项式的阶次N，且进行相应的重复性分析计算可以建立起稳态图。图1模态分析的稳