最小二乘复频域法(PolyMax)

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

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

3、关系可用右矩阵分式模型（ RMFD）来描述，右矩阵分式模型的表达式为H oU o1（1）D式中： H oC l Ni 理论频响函数的第 o 行， N i是输入点数，即激励数；U oC l Ni 分子多项式行向量；D oC N iNi 分母多项式矩阵。且 U o和 D o可以表示成如下形式：NU oZ rBor（ o 1,2 , N 0 ）（2）r 0N（3）D oZrArr 01欢迎下载。精品文档式中： N 多项式阶次其中分母系数矩阵Ar RNi Ni和分子系数行向量 Bor Rl Ni是待估计的参数。所有这些系数合并为一个矩阵。TTTT（4）1N其中Bo0A0oBo1R N 1 Ni ，oA

4、1R Ni N 1 N i（5）BoNAN式（ 2）和式（ 3）中出现的多项式基函数Z r，一般地，有以下两种选择： . 对于连续时域模型，可取为ijZr（6）s式中： s1Ni 缩放因子，用来提高方程的数值状况。2. 对于离散时域模型，可取为ZrejT r s（7）式中： Ts 采样周期。通常采用离散时域模型。（ 2）参数的线性化通过试验测量出的频率响应函数矩阵CNoNiC1 Ni表Hf，用 H o f示实测频响矩阵的第 o 行， o 1,2, N o ,f1,2, N f，那么关于参数矩阵的非线性最小二乘（ NLS）目标函数可表示为N oN fNLSHNLS（8）trf ,nlsof ,o

5、o 1f 1式中： ?H 矩阵的复共扼转置；2欢迎下载。精品文档tr ? 矩阵的迹，即矩阵的主对角元素之和。通过对式（ 8）求极小值，便可以得到频率响应函数矩阵的右分式矩阵模型各系数的估计值，即矩阵的估计值。式（ 8）中的加权非线性最小二乘误差函数被定义：NLSWof U of , o D1of ,WofH of ,H o ff ,H of（9）上式中 Wof是一个加权函数。一般地，为了提高估计的质量，我们采用WoH of（ 10）fvar H of式中： var ? 方差，可用相关函数求取。也可使用公式1（ 11）Wo fvar H of来做加权函数的。这两种加权函数都考虑了测量频响函数数据

6、的好坏：测得频响的方差越小，对目标函数的贡献越大。非线性误差函数可以经过一个近似的处理为一个线性的问题。实际上，通过对oNLSf ,右乘 Df ,，则可以得到一个关于参数为线性的方程，此加权LS线性最小二乘（ LS）方程误差of ,为LSNLSf ,Df ,of ,o%1（ 12）WofU of ,oD,H ofD f ,fN%WofzrfBorzrfArH or0这样式（ 12）关于参数为线性，将所有频率点装配成一列， f1,2,L , N f，它可用矩阵形式来表示LS1 ,oLSMXo Yooo（ 13）oJoLS,oNf其中：3欢迎下载。o 和精品文档Wo1zo1 , z11 ,L ,

7、zN1X oMC N fN 1（ 14）WoN fzoN f , z1N f ,L , zNN fWo1zo1 , z11 ,L , zN1%1H oYoMC N f Ni N 1（ 15）WoN fzoN f , z1N f ,L , zNN f%N fH o式中， Kronecker 积。（3）缩减标准方程加权线性最小二乘估计表达式为N oLSHLSl LStrooo1（ 16）NoTTRoSootroSoTToo 1式中：RoRe X oH X oR N1N 1SoHYoN 1Ni N1Re X oRToHRNiN1Ni N1Re Yo Yo同时，目标函数（ 16）等价于l LStrTR

8、e J H J式中， J 是 Jacobian 矩阵，被如下定义X10L0Y10X 2LM Y2N NN 1 NNJOMCofoiM MM00LX NoYNo（ 17）（ 18）为使 l LS值最小，将 l LS对系数矩阵求导，并令其为零lolLS2 Ro oSo0o1,2, L , No（ 19）NoLS2SoToTo0（ 20）o14欢迎下载。精品文档由式（ 19）得到 oRo1So，把它代入式（ 20）得NoSoT Ro 1So（ 21）2ToM0o 1NoSoT Ro1SoRNiN 1 Ni N 1其中， M 2To。o1由式（ 19）和（ 20）得到标准方程，经过整理，此标准方程的表

9、达式为R10L0S110R2L0S22MMOMMMH22Re JJ0（ ）00LRNoSNo22NoNoS1TS2TLSNToToo 1式（21）即为“缩减”标准方程，其中矩阵M 维数为 NiN1NiN1 ，比标准方程式 （22）中的 Re J H J 的维数 NoNi N 1NoNiN1要小的多。（ 4）求解缩减标准方程通过求解“缩减”标准方程，便可得到分母系数矩阵。根据线性方程组的求解理论，先对系数矩阵施加一个约束。假如，设定系数矩阵中的一个系数矩阵块等于正则常数矩阵 （例如设系数矩阵的最后一个矩阵块N1 I Ni ）,在这种前提下，缩减标准方程变为A X B（ 23）其中A M 1: N

10、i N ,1: Ni NB M 1: Ni N Ni N 1: NiN 1系数矩阵的最小二乘估计为?LSXXA 1 B（ 24）I Ni?1?一旦求得了oRo So就可得到所有的分子系数这种方法LS，那么通过LS,o考虑了标准方程的结构特性，比直接求解方程（22）要快得多。确定了分母系数5欢迎下载。精品文档矩阵后，通过求解的伴随矩阵的特征值和特征向量，这样就可以得到了系统的极点和相应的模态参与因子。方程如下0IL0000L00MMOMMV V .（ 25）00L0IA0TA1TLANT2ANT1上式中， V ,C N o N No N ，矩阵 V 的最后 N i 行就是模态参与因子；对角阵的角

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

12、特征矩阵中的对角元ie iTs ，由式（ 26）， i用 Rei i Imi 描述，则Re i i Im iei i %i TseiTsei %iTseiTscos%i sin%（ 27）iTsiTsReii ImieiTs（ 28）arctanImi%（ 29）ReiTsi所以i1lni（ 30）Ts6欢迎下载。精品文档%i1arctanImTsRe由此可求得频率i 和阻尼比ii（ 31）i22i%iii（ 32）ii在求得的频率i 和阻尼比i 包含有结构的固有频率r 和模态阻尼比r ，因此，必须对所有求得的i 和 i 进行有效的分析和选取，以确定系统真实的固有频率和阻尼比。建立稳态图就是一

13、种行之有效的方法。（ 6）建立稳态图在模态分析中，稳态图是帮助实验者分离结构物理极点和数学极点的一个有力工具，如图 1 所示。通过逐渐增大多项式的阶次N ，且进行相应的重复性分析计算可以建立起稳态图。16.7e-3ovoov vso vvovvvs ovvvso vvssvs32ooooo ossvos so ovovvvvvsvvo31o ooovo oovv vvo vvvvsvsvsv30o ooosoovo vvovvvvvvsvvs29ooo oosovov ooovvvvvvsvo28oo oovoovvsovvovvvvvv27oooovovov ooo vvoovvv26o o

14、oovvoov ovvvv vsvvso 25oooosvooosovsvvsvsv24eooovvooovvovvvsvvvs23N) doo ooovsvo voovvsvssv22u/tiglooo ovoov ooovsvvs o21( pmooooo vvoo sooosvsvso20Aooo vsosvovvosvs19oooosvoso vovosvs18oooovvosoooovss17oo ovoosoooovss16oosovosooov ss15oovo osvooo sv14oovvvvooo sv13oovovoooss12oovovovosv11ovovoosv1079.0e-61.00Linear2.05e+3Hz图 1模态分析的稳态图对于不同的 N 值，由式（ 25）得到相应的