最小二乘复频域法

最小二乘复频域法（PolyMax）SX1201069虞刚PolyMax模态识别方法，属于多自由度时域识别法，也称作多参考点最小二乘复频域法（Polyreferenceleastsquarescomplexfrequencydomainmethod），是最小二乘复频域法（LSCF）的多输入形式，是一种对极点和模态参预因子进行整体估计的多自由度法，一般首先通过实验建立稳态图，以判定真实的模态频率、阻尼和参预因子；建立可以线性化的直交矩阵分式模型，然后基于正则方程缩减最小二乘问题，得到压缩正则方程，于是模态参数可以通过求解最小二乘问题得到。该方法集合了多参考点法和LSCF方法的优点，可以得出非常清晰的稳态图，并且密集空间可以被分离出来，尤其在模态较密集的系统（动力总成系统），或者FRF数据受到严重噪声污染的情况下仍可以建立清晰的稳态图，识别出高度密集的模态，对每一个模态的频率、阻尼和振型都有很好的识别精度，是国际最新发展并流行的基于传递函数的模态分析方法。其基本思想如下：（1）建立频率响应函数模型多参考点最小二乘复频域识别技术（PRLSCF或PolyMAX）要以频响函数矩阵作为识别的初始数据，其数学模型采用右矩阵分式模型来描述。在频域中，系统输出o（o=1,2...,N0，其中N0为输出点数）和全部输入的关系可用右矩阵分式模型（RMFD）来描述，右矩阵分式模型的表达式为HS）=US）D（n>i（1）式中：H（o）eCm’一理论频响函数的第o行，N.是输入点数，即激励数；U（©）£CixN一分子多项式行向量；D（©）eCnxn.一分母多项式矩阵。且U。（J和D。（J可以表示成如下形式：U侦）=艺（Z侦）.B）（o=1,2...,N）（2）r=0D侦）=艺（Z侦）.A）（3）r=0

式中：N—多项式阶次其中分母系数矩阵气GRNxN和分子系数行向量BoGRg是待估计的参数。所有这些系数合并为一个矩阵。（4）其中P二o其中P二or、B_o0BD(.4o1\GR3+1)xn，:，BIoNJr、AA°P=<1>GRn(n+1^xN.O:，AN（5）式（2）和式（3）中出现的多项式基函数Z（©），一般地，有以下两种选择:r侦）』i.对于连续时域模型，可取为(6)式中：卜壬电—缩放因子,用来提高方程的数值状况。侦）』ii.对于离散时域模型，可取为ZCo）=e（-e-rs）（7）r式中：T一采样周期。通常采用离散时域模型。（2）参数的线性化通过试验测量出的频率响应函数矩阵虱JGC叩N用H（JgClxNj表示实测频响矩阵的第O行O=1,2,...,No,f=1,2,...,Nf）那么关于参数矩阵0的非线性最小二乘（NLS）目标函数可表示为£（0）=我/七trtnlsC,0》•（NLSCd,0））（8）nlsofofO=1f=1通过试验测量出的频率响应函数矩阵虱JGC叩N式中：・H一矩阵的复共扼转置；tr（・）一矩阵的迹，即矩阵的主对角元素之和。

通过对式（8）求极小值，便可以得到频率响应函数矩阵的右分式矩阵模型各系数的估计值，即0矩阵的估计值。G,0)G,0)=wG舄G,0)一H，》=W5G,p).*G"HG》OfOfOfOfOfOf-Of（9）上式中W（o）是一个加权函数。一般地，为了提高估计的质量，O上式中W（o）是一个加权函数。一般地，为了提高估计的质量，Of我们采用W(of)=HT

varHof（10）式中：var{•}一方差，可用相关函数求取。也可使用公式（11）(o)=（11）f来做加权函数的。这两种加权函数都考虑了测量频响函数数据的好坏：测得频响的方差越小，对目标函数的贡献越大。非线性误差函数可以经过一个近似的处理为一个线性的问题。实际上，通过对£NLS0）右乘D、a），则可以得到一个关于参数为线性的方程，此加权线性最小二乘（LS）方程误差£LS（of,0）为£ls(o,0)=£nls(o,0).D(o,a)=W(o)(U(o,p).D-i(o,a)-H(o)).D(o,0)(12)ofofofoff=WCo)丈(zCo).B—zCo)H(o).A)ofrforrforr=0这样式（12）关于参数为线性，将所有频率点装配成一列，f=1,2,...,七，它可用矩阵形式来表示£LS(0)=」O£LS(o,0)

o1£LSONf,0)[=&£LS(0)=」O£LS(o,0)

o1£LSONf,0)[=&y心=j[poOOlaIOla（13）其中:W0)「Z(",z(",•••,zo1Lo111NWoeCNfx(N+1)（14）-Wo1](气),...,zn(气)]®H0)eCNWoeCNfx(N+1)（14）-Wo1](气),...,zn(气)]®H0)eCNfxn(n+1)(15)WoNfNNf式中，®—Kronecker积。（3）缩减标准方程加权线性最小二乘估计表达式为/(0)=以GL(0)》•gLS(0))LSooo=1=&o=1(VcJ「RS■rp「\1PTat•00」0>111-o\二StTaJ（16）式中:R=Re(XhX)eR(n+»(n+1)同时，目标函数（16）等价于I（0）=trGt•Re（JhJ）.0）LS（17）式中，J是Jacobian矩阵，被如下定义yly2:.o・・・••••••0X2:.xlo一一一JL00…Xno"」为使孔（0）值最小，将、G）对系数矩阵po和a求导，并令其为零(0)=2(RP+Sa)=0(o=1,2,...,N)QPLSooooo—£(0)=2仔(StP+Ta)=0河LSLo=oooJ(19)(20)由式（19）得到P广-R-iSo・a，把它代入式（20）得区(T-StR-iS)ooooLo=1-(21)其中，-St其中，-StR-iS)oooo1-o=1-eRn.(n+i)xn(n+i)。由式（19）和（20）得到标准方程，经过整理，此标准方程的表达式为-r-r0••♦0S-11p0R••♦0s122.p•.2-..•5•-r00••♦rsNoNopNNoStSt••♦St1oTa_12Nooo=1=2Re(Jh.J)・。=0(22)式（21）即为“缩减”标准方程，其中矩阵M维数为N（N+1）xN（N+1），比标准方程式（22）中的Re（JHJ）的维数（N+N）（N+1）x（N+N）（N+1）要小的多。式（21）即为“缩减”标准方程，（4）求解缩减标准方程通过求解“缩减”标准方程，便可得到分母系数矩阵a。根据线性方程组的求解理论，先对系数矩阵a施加一个约束。假如，设定系数矩阵^中的一个系数矩阵块等于正则常数矩阵（例如设系数矩阵a的最后一个矩阵块a（N+1）=、）,在这种前提下，缩减标准方程变为TOC\o"1-5"\h\zA-X=B（23）其中A=M（1:NN,1:NN）系数矩阵a的最小二乘估计为\o"CurrentDocument"人「X1/、a=<>X=A-1B（24）LIvNjJ一旦求得了alS那么通过p广-r-1So・a就可得到所有的分子系数6LSo这种方法考虑了标准方程的结构特性，比直接求解方程（22）要快得多。确定了分母系数矩阵a后，通过求解以的伴随矩阵的特征值和特征向量，这样就可以得到了系统的极点和相应的模态参与因子。方程如下

0I...000•0•...0••0••0•0・•...0•I—At—At…—At—At0iN—2-V=KA(25)N—1-1上式中，V,AeCN°Nf，矩阵V的最后N.行就是模态参与因子；对角阵A的角元记录为Ai（i=1,2,・..,NN）由不稳定的数学极点和稳定的物理结构点两部分组成记稳定的物理结构极点为A=e-小，通过对这些物理结构极进行转换，便可得出结构r(26)的固有频率s,和模态阻尼比匚，；关系式如下人,人*=b土is或人0I...000•0•...0••0••0•0・•...0•I—At—At…—At—At0iN—2-V=KA(25)N—1-1(26)（5）计算频率点和阻尼比点根据信号与系统基本理论中对系统稳定性的描述：系统的全部极点落于s域左半平面（不包括虚轴），且满足有界输入有界输出原则，系统是稳定的。复特征矩阵A中的复特征值总是以共轭对的形式出现，同时也包含实数（虚轴上），在求解频率点七和阻尼比点匚，时，对于每个共轭对只取其中一个进行分析，且不考虑实数。复特征矩阵A中的对角元△**，由式（26），A用Re（A）+i•Im（A）描ii述，则Re(A)+i-Im(ARe(A)+i-Im(A)=e-Gi+isRs=e-叮•e—isTs=e-。t「cos(sT)—isin(sT)"i（27）所以|Re(A)+i-Im(△)=e-叩（28）arctan/Im(A)'I【（29）（30）~1f&=—-arctaniTs'Im(A)'IRemIXi/(31)~1f&=—-arctaniTs'Im(A)'IRemIXi/(31)(32)iii(32)c=二i①i在求得的频率气和阻尼比c，包含有结构的固有频率®r和模态阻尼比cr，因此，必须对所有求得的气和C，进行有效的分析和选取，以确定系统真实的固有频率和阻尼比。建立稳态图就是一种行之有效的方法。建立稳态图在模态分析中，稳态图是帮助实验者分离结构物理极点和数学极点的一个有力工具，如图1所示。通过逐渐增大多项式的阶次N，且进