一种基于小波多分辨率分析的动载荷识别方法技术

本发明专利技术公开了载荷识别领域的一种基于小波多分辨率分析的动载荷识别方法，用以解决目前载荷识别领域研究中存在的问题。该方法包括：首先求解识别参数；其次，基于时域卷积模型，利用小波基函数对载荷进行小波重构，得到小波响应函数；再次，对响应和小波响应函数进行小波变换，得到小波域下的系统响应和小波域下的小波响应函数；最后，计算权系数，反求载荷，完成识别。本发明专利技术可辨识冲击、突变等非平稳载荷，识别精度高；方法对多输入多输出系统中同时包含快慢变化的多路载荷间的相互干扰不敏感，可以区别并辨识多路载荷；提出了多个识别参数的定量/定性确定方法，可用于基于小波分析的动载荷识别方法的参数确定。

【技术实现步骤摘要】
一种基于小波多分辨率分析的动载荷识别方法
本专利技术涉及载荷识别领域，特别涉及一种基于小波多分辨率分析的动载荷识别方法。
技术介绍
工程结构上的动载荷是结构的动力学设计、减振隔振设计、强度校核、动力学优化等的重要原始参数，然而由于技术或经济条件的限制，很多情况下结构所受的动态载荷很难直接进行测量。载荷识别就是在这样的背景下产生的，其可以在许多无法直接量测载荷的情况下通过响应等信息反求载荷。传统载荷识别方法按信号所处域可分为频域识别和时域识别。载荷识别是一类非常基础的问题，研究中还存在一定问题，如频域下的载荷识别要求测量数据的样本具有一定的长度，一般只适用于稳态或平稳随机载荷的识别，对冲击、突变等瞬态载荷的识别不理想，当频率接近共振区时会出现数值不稳定问题；时域下的载荷识别包含递推，对边值条件和初始条件敏感，初始的误差会随着递推推导不断传递造成误差累积，有些方法需要知道一种类型以上的响应，给使用带来困难；小波分析是一种新颖的时频分析方法，其通过伸缩平移运算对信号逐步进行多尺度细化，最终达到高频处时间细分，低频处频率细分，能自动适应时频信号分析的要求，从而可聚焦到信号的任意细节。因此，一些研究将小波分析用于载荷识别来获得更理想的识别结果。基于小波分析的载荷识别方法是在小波的时频域下识别载荷，利用小波如尺度特性、平移特性、紧支性、移位正交性、多分辨率分析等优良性质，可识别包括冲击载荷在内的任意确定性载荷，但由于发展时间短，方法的理论基础并不完善，部分推导的表达式针对特定问题，局限性大，一些参数没有定性或定量的选择依据，一些方法没有利用二维分析特性。
技术实现思路
本专利技术的目的在于，提出一种基于小波多分辨率分析的动载荷识别方法，用以解决目前载荷识别研究中存在的问题。为实现上述目的，本专利技术提出的技术方案是，一种基于小波多分辨率分析的动载荷识别方法，其特征是所述方法包括下列步骤：步骤1：求解识别参数，积分上下限J和I、分解层级j；步骤2：基于时域卷积模型，利用小波基函数对载荷进行小波重构，得到小波响应函数其中，t表示时间，h(t)是单位脉冲响应函数，φ(2jt-k)是小波基函数，k为平移因子，j为分解层级(尺度因子)；步骤3：对响应和小波响应函数进行小波变换，得到小波域下的系统响应和小波域下的小波响应函数其中y(t)是系统的响应；步骤4：反求载荷在小波域下的权系数P＝G-1u，其中u和G分别为小波域下的系统响应u(i)和小波域下的小波响应函数gik的矩阵形式；步骤5：由权系数计算载荷其中P(k)是权系数矩阵P的每一个元素，并完成识别。本专利技术实现了一种基于小波多分辨率分析的动载荷识别方法，其有益效果如下：1给出一种基于小波多分辨率分析的动载荷识别方法，方法可识别非平稳载荷，对冲击、突变等快变瞬态载荷识别能力优秀；2给出的方法对多输入多输出系统中同时包含快慢变化的多路载荷间的相互干扰不敏感，可以区别并辨识多路载荷；3提出了多个识别参数的定量/定性确定方法，可用于基于小波分析的动载荷识别方法的参数确定。附图说明图1本专利技术方法的流程图。图2是实验示意图。图3是悬臂梁节点布置示意图。图4是原载荷和识别载荷的时域图。具体实施方式下面结合附图，对优选实施例作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本专利技术的范围及其应用。本专利技术解决问题的思路是基于小波和其多分辨率分析性质完成动载荷识别：首先求解识别参数；其次，计算小波响应函数；再次，计算小波域下的系统响应和小波域下的小波响应函数；最后，计算权系数，反求载荷，完成识别。图1是本专利技术完成的一种基于小波多分辨率分析的动载荷识别方法的流程图，该方法具体步骤如下：步骤1：量测频响函数(FrequencyResponseFunctions，FRF)。实验示意图如图2所示，首先使用力锤激励法量测悬臂梁系统的频响函数信息，悬臂梁节点布置如图3所示，节点2、3布置传感器，节点1、4不放置激振器，使用力锤敲击节点1，敲击5次取平均值，得到FRF21和FRF31，同理可得FRF24和FRF34，其中数字如21表示第1个节点激励下第2个节点的频响函数。步骤2：计算单位脉冲响应函数(ImpulseResponseFunctions，IRF)。采集的FRF为单边谱。FRF和IRF是傅里叶对子，即IRF经过傅里叶变换后可以得到关于频域对称的FRF。所以首先将FRF扩展成共轭对称形式，然后再由逆傅里叶变换得到IRF。步骤3：施加激励和采集响应。利用振动加速度传感器采集外激励作用下的响应信号。图3中从右至左的4个节点分别设置：节点1施加冲击激励、节点4施加频率10Hz幅值为15N的正弦激励，节点2和3采集响应数据，采样频率4096Hz。共进行3次重复实验，选择最优一组数据作为下节的分析数据，选择依据为冲击特征明显并且毛刺少。步骤4：求解识别参数，积分上下限J和I，分解层级j。步骤4中参数求解的推导说明依托于步骤5至步骤7，因此先介绍步骤5至步骤7。步骤5：计算小波响应函数。小波响应函数是小波基函数通过相同系统得到的输出，综合反映系统特性和小波函数。计算过程为：首先生成db小波基函数(0尺度)，使用数据压缩和重采样生成j尺度下的小波函数，使用平移生成和权系数个数(J-I+1)相等的小波函数族，进而和IRF做卷积可得小波响应函数。其具体步骤包括：步骤51：对于时域卷积模型y(t)＝h(t)*f(t)(1)其中，y(t)是系统的响应，h(t)是单位脉冲响应函数，f(t)是输入载荷。对f(t)使用小波重构，可以得到P(k)＝α(k)2j/2(3)其中，α(k)是尺度系数，P(k)是权系数，I和J是积分下限和上限(两者差值是权系数的个数)，j是分解层级，k是平移因子。这里使用的小波是正交Daubechies(db)小波。步骤52：将式(2)代入式(1)，可以得到其中，是单位冲击响应函数和小波基函数在分解层级j上的卷积，为小波响应函数。步骤6：计算小波域下的系统响应和小波响应函数。对响应和小波响应函数同时进行小波变换，可得小波域下的系统响应和小波域下的小波响应函数，如式(6)所示。其具体步骤为：步骤61：对式(4)两侧同时进行小波变换，步骤62：进行如下替代其中u(i)是小波域下的系统响应，gik是小波域下的小波响应函数。式(6)可化为矩阵形式为u(J-I+1)×1＝G(J-I+1)×(J-I+1)P(J-I+1)×1(10)步骤7：计算载荷，完成识别。由小波域下的系统响应和小波响应函数反求载荷的权系数，如(11)所示，再由权系数计算载荷，完成识别，如式(12)所示。反求载荷在小波域下的权系数的具体方法为，对式(10)进行逆变换，求得载荷在小波域下的权系数P＝G-1u(11)所述由权系数计算载荷，完成识别的具体方法为，将式(11)反求得到权系数代入式(2)，可计算载荷，完成识别。步骤4中求解识别参数，即求解从式(1)到式(11)的推导中的积分上下限J和I、分解层级j，其具体步骤为：步骤41：I和J的差值决定了权系数的个数，对于J，需要保证最后一个基函数的起始时间不大于响应信号的结束时间。对于I，需要保证第一个基函数的结束时间不小于响应信号的起始时间。在这种情况下，基函数才能覆盖响应信号的全部时间，本文档来自技高网...

【技术保护点】
一种基于小波多分辨率分析的动载荷识别方法，其特征是所述方法包括：步骤1：求解识别参数，积分上下限J和I、分解层级j；步骤2：基于时域卷积模型，利用小波基函数对载荷进行小波重构，得到小波响应函数其中，t表示时间，h(t)是单位脉冲响应函数，φ(2jt‑k)是小波基函数，k为平移因子，j为分解层级(尺度因子)；步骤3：对响应和小波响应函数进行小波变换，得到小波域下的系统响应和小波域下的小波响应函数其中y(t)是系统的响应；步骤4：反求载荷在小波域下的权系数P＝G‑1u，其中u和G分别为小波域下的系统响应u(i)和小波域下的小波响应函数gik的矩阵形式；步骤5：由权系数计算载荷其中P(k)是权系数矩阵P的每一个元素，并完成识别。

【技术特征摘要】
1.一种基于小波多分辨率分析的动载荷识别方法，其特征是所述方法包括：步骤1：求解识别参数，包括积分上J、积分下限I及分解层级j；其具体步骤为：步骤101：I和J的差值决定了权系数的个数，对于J，需要保证最后一个基函数的起始时间不大于响应信号的结束时间；对于I，需要保证第一个基函数的结束时间不小于响应信号的起始时间；在这种情况下，基函数才能覆盖响应信号的全部时间，所以在分解层级j下，有2-jtl≤tend(1)2-j(tf+2N-1)≥tbegin(2)其中，tl是最后一个基函数的起始时间，tf是第一个基函数的起始时间，tbegin是响应信号的起始时间，tend是响应信号的结束时间；求解得：tl≤2jtend(3)tf≥2jtbegin-(2N-1)(4)J＝max(tl)J∈Z(5)I＝min(tf)I∈Z(6)步骤102：依据分解层级、权系数个数和...

【专利技术属性】
技术研发人员：李宗，褚福磊，秦朝烨，阎绍泽，
申请(专利权)人：清华大学，
类型：发明
国别省市：北京;11