一种仅有响应输出情况下获取系统特征函数和信号特征值的方法，属动态测试领域中参数识别的方法。该方法是采用各响应点的互谱密度函数代替频响函数进行时频滤波和频域参数识别的方法，包括步骤（１）对不同测量输出点信号进行互谱密度函数进行分析的计算；（２）根据互谱计算结果在时频域进行非正交小波分析计算；（３）反傅立叶变换求时频分析系数；（４）通过加矩形窗进行时频滤波；（５）求滤波后的输出信号的互谱作为识别用的系统函数；（６）进行曲线拟合求系统参数；本方法能提高系统参数的识别精度，能准确识别模态参数，简单方便，适用于在运行状态下的大型复杂机械设备、高层建筑、桥梁等大型土木工程设施的动态分析、性能验证和故障诊断。
【技术实现步骤摘要】
本专利技术涉及的一种基于互谱函数和时频滤波在仅有响应输出情况下获取信号特征的方法，属于动态测试ip参数识别领域。 背录技术仅有输出信号的特征提取问题是运行情况下机械系统、动态测试和故障诊断 技术的关键技术。利用环境激励作为试验的振动激励源进行摸态识别最大优点是(1)省去了研制以及安装专用激励装置的费用和时间；(2)可同时激励所有 操纵面，同时激出对称和反对称模态，省去霈要多次进行的试验；(3)大大减少 试验次数和周期。运行状态下的大型复杂机械设备、高层建筑、桥梁等大型土木工程设施的动 态分析、设计性能验证和故障诊断中不便使用激励设备，与其他激励(如扫频和 脉冲激励)相比环境激励强度较弱、信噪比低、使得系统参数识别的分散度大大 增加，系统的特征提取十分困难受噪声影响很大。需要采用新的信号处理和参数 识别理论以获得系统特征提取的精度等。综上所述，研究环境激励下仅有输出信 号的信号特征提取方法有着迫切性需求和重大的应用价值。
技术实现思路
本方法的目的在于提供一种新的时频滤波技术和频域识别技术相结合的特 征参数识别算法，针对环境激励下仅有输出信号强度较弱、信噪比低、使得系统 参数识别的分散度大大增加的情况，提髙系统参数的识别精度。专利技术其特征在于: 解决由于环境激励情况下输入力不可测，无法得到频响函数(FRF),因而无法釆 用传统EMA频域识别方法。首先考虑到互谱密度函数与频响函数有相似的表达 式。提出互谱密度函数-时频空间域分解方法。通过随机响应数据估计功率互谱 密度函数(CSD)，接着，由CSD代替频响函数进行基于非正交小波的时频滤波， 然后采用正交多项式方法估算出结构模态频率和阻尼。这种方法简单方便，结果 证明不仅可以准确识别模态参数，还可以实现自动识别，因而特别适用于在线监、 由于环境激励情况下输入力不可测，无法得到进行参数识别的频响函数 (FRF),因而无法采用传统EMA频域识别方法。在包括桥梁在内的环境激励的复 杂结构运行模态分析应用中，存在噪声模态的干扰，以及相应的结构模态区分和 精确识别，James等人于1995年曾提出NexT方法，即在白噪声环境激励下结 构两点之间响应的互相关函数和脉冲激励下的响应函数有相似的表达式。后来 发展了基于脉冲响应函数的Morlet小波参数识别算法，但是所有时域模态识别方法都对测试噪声、信号处理泄漏误差、模态截断、非线性等因素十分敏感。为 了克服时域识别的缺点，希望采用频域参数的识别方法来克服噪声的影响，因此 提出互谱功率密度函数-时频空间域分解方法。具体步骤是 (1).对不同测量点输出信号进行互功率谱密度分析的计算根据信号处理理论，设任意两个随机输出信号分别为义(/7), /(/7)， r=Ut，r为采样周期，w为样本点数，是自然l^t为采样间隔。其离散傅立叶变换J(O)、/(W)，则其互功率谱密度定义为s一)+[竭.n)] (i)建立在傅立叶变化基础上的互功率谱密度估计，其基本过程为首先计算样 本信号序列义()，y()的离散傅立叶变换，然后取的变换结果与的变 换结果的共轭相乘，并除以样本序列的个数W,作为真实功率谱的一个估计，可 表示为，4l (2) 式中W ( )表示F(w)的共轭，Ww)分别为)，/(/7)的离散傅立叶变换，即= gx(K乂2卿，；k(必)=^J:K)e力2咖，^o, i, 2，……，N-1釆用统计平均方法减少随机误差，将长度为n的数据分为若干段，分别求出每一段的互功率谱，然后加以平均，同时采用选取的每段数据有部分重叠。通 过对分段数目、分段大小及数据重叠数目进行权衡，增加平均次数达到最佳优化。功率谱密度的单位是被测信号物理单位的平方除以频率(Hz)。(2).根据步骤(1)中的互功率谱密度计算结果在时频域进行非正交小波 分析计算；对于任意的函数i a)ef(及)的连续小波函数定义为 i i一i/ f 一 ft (3)在这里cp^^cK^)是小波的基函数。6是时间变量。fl是尺度参数。该， O变换限于Heisenberg不确定性原理所限定的分辨率以内。Heisenberg不确定性 原理为Ar」S1/2。这里4 、」t表示频窗和时窗的各一半宽度，单位分别为弧 度/秒和秒。基于随机激励的下的输出信号为工程背景。选择Morlet小波基函数p(0 = ;e_、—w。采用无变化尺度参数的Morlet小波(=常数)作为基函数，来进行时频分析。采样信号巻积小波滤波器得到小波系数。在构造频域滤波器时，采用形式为Morlet小波，其中/分别为力,/2…/变化的n个小波基函数。/!-厶为分析频带，它们之间的间隔为A。在频域，尺度系数a=l,中心频率为/的Morlet小波滤波器的相应小波系数是<formula>formula see original document page 7</formula>这里①y(co)s①(co，/)是Morlet小波的傅立叶变换。等式(4)是复数的小波基函数，在频域滤波会带来相位失真。信号重构相应也会相位失真。为了避 免这个问题，仅用Morlet小波的实部构造滤波器，使得重构的信号与原始信号 同相位。Morlet小波实部的傅立叶变换为-①々)二e 2 +e 2 (5)(3).反傅立叶变换求时频分析系数； 对应的时域系数<formula>formula see original document page 7</formula>、(w)、 OV(Q))分别为信号互谱密度函数和小波实部的傅立叶变换。实的小波滤波器组对于每一个中心频率都提供了一个有限脉冲响应(FIR)线性相位滤 波器。使得原始信号和重构信号的相位一致。(6)式用来求得时域信号，以便进 行时频空间滤波。设个小波基函数用来进行小波的重构，其中n为自然数。等式(5)用来构 造FIR滤波器组。对于给定的频率^，小波系数通过滤波器组乘原始信号的互谱密度函数可得到<formula>formula see original document page 8</formula>(4) .通过加矩形窗进行时频滤波；对步骤(3)中的公式(6)中的小波系数加矩形窗滤波，通过加矩形窗滤波，将矩形窗外部的噪声信号滤掉，这样得到的是良好的系统特征函数，互功率谱密 度函数。(5) .求滤波后的输出信号的互功率谱密度作为识别用的系统函数；从步骤(3)中的等式(7)用广义逆求得滤波后的互谱<formula>formula see original document page 8</formula>在每一个频率0。处求得&y(叫)。该方法用来滤除输入信号和响应信号的失 真和干扰，提取想要的信号特征。(6).进行曲线拟合求系统参数；在时频空间滤波去除了噪声和非平稳的影响后，得到较好的系统功率谱密度 函数，该函数代替平稳信号输出情况下随机激励的频响函数，这样就能采用传统 的频域识别算法进行系统特征参数的识别。利用对噪声不敏感的有理分式正交多 项式识别方法(RFOP)对重构互谱密度函数进行识别，求得模态频率和阻尼比。 其基本做法是将结构的频率响应函数展开为有理分式形式，将总方差定义为目本文档来自技高网...
【技术保护点】
一种仅有响应输出情况下获取系统特征函数和信号特征值的方法，其特征在于，采用互功率谱密度函数和时频空间域分解方法，具体步骤是：（１）．对不同测量点输出信号进行互功率谱密度分析的计算根据信号处理理论，设任意两个随机输出信号分别为 ｘ（ｎ），ｙ（ｎ），Ｔ＝Ｎ.△ｔ，Ｔ为采样周期，Ｎ为样本点数，是自然数，△ｔ为采样间隔，其离散傅立叶变换Ｘ（ω）、Ｙ（ω），则其互功率谱密度定义为Ｓ↓［ｘｙ］（ω）＝２／ＴＥ［Ｘ（ω）.Ｙ↑［＊］（ω）］（１）建立在傅立 叶变化基础上的互功率谱密度估计，其基本过程为：首先计算样本信号序列ｘ（ｎ），ｙ（ｎ）的离散傅立叶变换，然后取ｘ（ｎ）的变换结果与ｙ（ｎ）的变换结果的共轭相乘，并除以样本序列的个数Ｎ，作为真实功率谱的一个估计，表示为Ｓ↓［ｘｙ］（ω） ＝２／ＮＴ＊Ｘ↓［ｉ］（ω）.Ｙ↓［ｉ］↑［＊］（ω）（２）式中Ｙ＊（ω）表示Ｙ（ω）的共轭，Ｘ（ω）、Ｙ（ω）分别为ｘ（ｎ），ｙ（ｎ）的离散傅立叶变换，即Ｘ（ω）＝＊ｘ（ｎ）ｅ↑［－ｊ２πｎω／Ｎ］，Ｙ（ω）＝＊ｙ（ｎ ）ｅ↑［－ｊ２πｎω／Ｎ］，ｎ＝０，１，２，……，Ｎ－１采用统计平均方法减少随机误差，将长度为Ｎ的数据分为若干段，分别求出每一段的互功率谱，然后加以平均；同时采用选取的每段数据有部分重叠，通过对分段数目、分段大小及数据重叠数目进行权 衡，增加平均次数达到最佳优化，功率谱密度的单位是被测信号物理单位的平方除以频率；（２）．根据步骤（１）的互功率谱密度计算结果在时频域进行非正交小波分析计算；对于任意的函数Ｒ（ｔ）∈Ｌ↑［２］（Ｒ）的连续小波函数定义为： Ｗ（ａ，ｂ）＝〈Ｒφ↓［ａ，ｂ］〉＝｜ａ｜↑［－１／２］∫↓［Ｒ］Ｒ（ｔ）φ↑［＊］（＊＊＊）ｄｔ（３）在这里φ↓［ａ，ｂ］（ｔ）＝φ（＊＊＊）是小波的基函数，ｂ是时间变量，ａ是尺度参数，该变换限于Ｈｅｉｓｅｎｂｅ ｒｇ不确定性原理所限定的分辨率以内，Ｈｅｉｓｅｎｂｅｒｇ不确定性原理为：△↓［ω］.△↓［ｔ］≥１／２，这里△↓［ω］、△↓［ｔ］表示频窗和时窗的各一半宽度，单位分别为弧度／秒和秒，基于随机激励的下的输出信号为工程背景，选择Ｍｏｒｌ ｅｔ小波基函数φ（ｔ）＝１／＊ｅ↑［－ｔ↑［２］／２］ｅ↑［－ｉｆｔ］，采用无变化尺度参数ａ即ａ＝常数的Ｍｏｒｌｅｔ小波作为基函数...
【技术特征摘要】
1.一种仅有响应输出情况下获取系统特征函数和信号特征值的方法，其特征在于，采用互功率谱密度函数和时频空间域分解方法，具体步骤是(1).对不同测量点输出信号进行互功率谱密度分析的计算根据信...
【专利技术属性】
技术研发人员：岳林，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：84[中国|南京]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。