CN105242111B - 一种采用类脉冲激励的频响函数测量方法 - Google Patents
一种采用类脉冲激励的频响函数测量方法 Download PDFInfo
- Publication number
- CN105242111B CN105242111B CN201510594076.6A CN201510594076A CN105242111B CN 105242111 B CN105242111 B CN 105242111B CN 201510594076 A CN201510594076 A CN 201510594076A CN 105242111 B CN105242111 B CN 105242111B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- frequency
- sut
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005316 response function Methods 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000005284 excitation Effects 0.000 title claims abstract description 32
- 238000012360 testing method Methods 0.000 claims abstract description 62
- 238000005259 measurement Methods 0.000 claims abstract description 37
- 230000003595 spectral effect Effects 0.000 claims abstract description 14
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000013459 approach Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 21
- 238000005070 sampling Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims 2
- 230000031068 symbiosis, encompassing mutualism through parasitism Effects 0.000 description 7
- 238000005086 pumping Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000011888 foil Substances 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明涉及一种采用类脉冲激励的频响函数测量方法,属于电系统和机械系统频响函数的测量技术领域;该方法包括:确定频率范围,构造类脉冲激励信号作为被测系统的输入信号,测量被测系统输入信号与被测系统输出信号,采用中心零点窗分析法对输入输出信号进行加窗傅里叶变换得到加窗奈奎斯特频谱,利用奈奎斯特频谱计算得到被测系统的频响函数。该方法可实现频响函数的快速测量,且在测量强共振系统的频响函数时,可获得其他方法难以获得的高精度。
Description
技术领域
本发明属于电系统和机械系统频响函数的测量技术领域,特别涉及预期频带的频响函数的高效测量以及强共振频响函数的精确测量。
背景技术
系统的频响函数是系统的响应(输出)信号频谱函数与激励(输入)信号频谱函数之比,可以反映出在频率范围内系统的动力学/电学特性。频响函数方法是测量机械系统和电系统的特性的最重要的手段之一。在快速傅里叶变换(FFT)方法出现之后,为了提高测量效率,人们广泛采用宽频带激励。对于一个系统来说,在测量频响函数的过程中应关注波形畸变,能量泄露以及信噪比等问题。测量出的频响函数的质量是衡量测量方法好坏的最直接的标准。
目前一般采用宽频带激励的频响函数测量方法,该方法更关注于余弦激励的叠加,因为这种激励比较灵活且具有可行性。对于采用这种激励的频响函数测量方法来说,可以保证每个频率点处的能量足够大,降低信噪比。可以证明,当测量时间有限时,随机相位的余弦信号对畸变和干扰具有较强的抵抗力。虽然整周期的随机相位取样能够消除泄露误差,但利用这种方法很难精确测量强共振处的频响函数。
在实际工程中,通常通过加窗来限制积分,同时也会带来能量泄露。当整周期采样时,如果采用矩形窗,则可以消除能量泄露,如果是非整周期采样,结果就很容易受到能量泄露的影响。类脉冲激励法可以在非整周期截取的情况下得到系统精确的频响函数。
在信号的数字化采样过程中,由于被采的模拟信号中存在共生信号,因此易引起频率混淆。由于共生信号频率都大于奈奎斯特频率,因此必须在采样之前通过模拟的低通滤波方式,滤除被采模拟信号中的共生信号,才能保证在采样过程中,不发生频率混淆现象。因此当发现被采信号中存在共生频率,或者不能确定被采信号中是否存在共生频率时,就要对被采信号实施模拟低通滤波,滤除已知或潜在的共生信号,以确保被采信号中没有共生信号,再进行采样,这样的滤波称为抗混滤波。由此需要的抗混滤波器应是模拟低通滤波器。有多种形式的模拟低通滤波器,巴特沃斯模拟低通滤波器因其具有非常平直的通带特性和较宽的线性相位带宽,经常被用作抗混滤波器。偶阶巴特沃斯模拟低通滤波器可以用较少的元器件构成较高阶次的滤波器,在抗混滤波器中被大量应用。
悬臂梁是一个连续弹性体,具有无限多个自由度,即有无限多个固有频率和主振型。在一般情况下,梁的振动是无限多个主振型的叠加,如果给悬臂梁施加一个大小合适的激振力,其频率正好等于悬臂梁的某阶固有频率,就会产生共振,对于这一阶固有频率的确定的振动形态叫做这一阶的主振型,这时其他各阶振型的影响可以忽略不计。
发明内容
本发明的目的是为克服已有技术的不足之处,提出了一种采用类脉冲激励的频响函数测量方法,该方法可实现频响函数的快速测量,且在测量强共振系统的频响函数时,可获得其他方法难以获得的高精度。
本发明提出的一种采用类脉冲激励测量频响函数的方法,其特点在于,对被测系统(机械系统或电系统)输入一个由在一定频率范围内均匀分布的余弦信号组成的类脉冲激励信号,通过测量输入输出信号的加窗奈奎斯特频谱,可以计算出被测系统在该频段的频响函数。
该方法具体包括以下步骤:
1)确定频率范围:确定所需的被测系统频响特性的频率范围[fl,fh],对于电系统或机械系统,频率范围为0~5000Hz;
2)构造类脉冲激励信号作为被测系统的输入信号:
21)确定类脉冲信号的频率间隔Δf:
其中,T为类脉冲信号的周期,Δf范围为0.005Hz~20Hz,T的范围为0.05s~200s;
22)测量的频响函数频率范围为[fL,fH],2Δf≤fL<fH-Δf,fH<fsmax/nH取fL和fH之间(包含fL和fH)的整数频率点:
[fL,fH]:fL≤f≤fH,且f∈I
其中,fH-fL的值应能被Δf整除,fsmax为测量仪器的最大采样频率,nH为每两个频率间隔点间的采样点数,nH的取值范围为[5,50];
23)用一系列振幅相同相位为零的余弦信号进行叠加构造类脉冲激励信号,作为被测系统输入信号,如式(1):
式(1)中,cm(t)为余弦波,fm为第m个余弦波的频率,fm=(m-1)Δf+fL,M为余弦波的总数,当m=M时,fm=fH,A为输入信号振幅,取0V<A<40V;
3)测量被测系统输入信号与被测系统输出信号:
将式(1)中得到的类脉冲信号作为机械系统或电系统的输入信号,测量被测系统输入信号以及被测系统输出信号,记测量到的被测系统输入信号为x(t),被测系统输出信号为y(t),测量时采样频率fs满足fs=nHfH;
4)采用中心零点窗分析法对输入输出信号进行离散加窗傅里叶变换得到加窗奈奎斯特频谱:
对测量得到的被测系统输入输出信号进行加窗处理,式(2)为窗函数,然后进行傅里叶变换得到被测系统输入输出信号的加窗无理频谱Xw(f)和Yw(f),如式(3)和式(4):
式(2)中β为0.54/0.46或1,零点为窗函数的中点,Tw为窗宽,Tw≤T,x(t)和y(t)分别为输入和输出信号,W0为正则化系数,W0=βTw/(β+1),设零点在窗的中间,则Xw(f)的相位与x(t)中点的相位一致,Yw(f)的相位与y(t)中点的相位一致。λ(λ≥1&λ∈I)为补零倍数;
设x(tk)为测得的x(t)的离散信号,y(tk)为测得的y(t)的离散信号,离散加窗傅里叶变换后得到的输入输出信号的加窗奈奎斯特频谱Xw(fn)和Yw(fn)分别为:
式(5)和式(6)中fn为离散化频率,fn=n/λTw,频率间隔由1/λTw决定,Δt=1/fs,N=λTw/Δt,0≤n≤N-1,n∈I,tk为离散时间,tk=kΔt,0≤k≤N-1,k∈I,FFT[·]为快速傅里叶变换;
5)利用加窗奈奎斯特频谱计算得到的被测系统的频响函数如式(7)所示:
式(7)中,Xw(fn)和Yw(fn)为输入信号和输出信号的加窗奈奎斯特频谱,由式(5)和式(6)得到,窗宽Tw=rwT(0<rw≤1),fn为离散化频率。
本发明的特点及有益效果:
本发明所采用的激励信号是利用余弦波叠加的方式合成的谐波型信号,在制作信号时,在目标频率范围内采用定频率间隔的余弦信号进行叠加,可使激励信号的能量在频率范围内均匀分布。使输入信号在时域上接近脉冲信号。
本发明在对测量到的电系统/机械系统激励信号x(tk)和响应信号y(tk)进行信号处理时,利用傅里叶变换将信号从时域变换到频域,得到加窗奈奎斯特频谱Xw(fn)和Yw(fn)。两者相除即可得到目标频段频响函数。
本发明在计算频响函数时,即使非整周期截取信号,也不会带来泄露误差。
本发明通过对系统输入信号的制作,对系统输入输出信号的测量、计算,实现了电系统/机械系统频响函数的准确测量。
附图说明
图1为本发明采用类脉冲激励测量频响函数的具体方法流程框图。
图2为本发明方法实施例1(抗混滤波器)的实现装置示意图。
图3为利用本发明的类脉冲激励法测得的抗混滤波器频响函数实验结果曲线,其中(a)为抗混滤波器频响函数幅值谱,(b)为抗混滤波器频响函数相位谱。
图4为本发明方法实施例2(悬臂梁振动)的实现装置示意图。
图5为利用本发明的类脉冲激励方法测得的悬臂梁振动频响函数实验结果曲线,其中(a)为悬臂梁振动频响函数幅值谱,(b)为悬臂梁振动频响函数相位谱。
图6为本发明方法实施例2(悬臂梁振动)在共振频率处频响函数的精细测量结果,其中(a)为悬臂梁振动频响函数共振频率处幅值谱,(b)为悬臂梁振动频响函数共振频率处相位谱。
具体实施方式
下面结合附图对本发明的具体实施方式做详细说明。
本发明提出的一种采用类脉冲激励测量频响函数的方法,其特点在于,对被测系统(机械系统或电系统)输入一个由在一定频率范围内均匀分布的余弦信号组成的类脉冲激励信号,通过测量被测系统输入输出信号,可以计算出被测系统在该频段的频响函数。
本发明方法采用具体步骤如图1所示:
1)确定频率范围:
确定被测的电系统/机械系统的频率范围[fl,fh],对于电系统或机械系统,频率范围为0~5000Hz;
2)构造类脉冲激励信号作为被测系统的输入信号:
21)确定类脉冲信号的频率间隔Δf:
其中,T为类脉冲信号的周期,在选择Δf时,如果被测系统的频响函数有强烈的共振,频率间隔要取得小些,以准确测量共振峰值;如果被测系统频响函数比较平滑,频率间隔可以取得大些;
本发明的频率间隔范围为0.005Hz~20Hz,T的范围为0.05s~200s;
22)测量的频响函数频率范围为[fL,fH],2Δf≤fL<fH-Δf,fH<fsmax/nH,取fL和fH之间(包含fL和fH)的整数频率点:
[fL,fH]:fL≤f≤fH,且f∈I
其中,fH-fL的值应能被Δf整除,fsmax为测量仪器的最大采样频率,nH为每两个频率间隔点间的采样点数,为了保证频率为fH的余弦分量的每个周期内有足够数量的测试点,nH的取值范围为[5,50];
23)用一系列振幅相同相位为零的余弦信号进行叠加构造出类脉冲激励信号,作为被测系统输入信号,如式(1):
式(1)中,cm(t)为余弦波,fm为第m个余弦波的频率,fm=(m-1)Δf+fL,M为余弦波的总数,当m=M时,fm=fH,A为输入信号振幅,取0V<A<40V;
3)测量被测系统输入信号与被测系统输出信号:
将式(1)中的类脉冲激励信号作为被测系统的输入信号,测量被测系统输入信号以及被测系统输出信号,记测量到的被测系统输入信号为x(t),被测系统输出信号为y(t),测量时采样频率fs=nHfH。
4)采用中心零点窗分析法对输入输出信号进行离散加窗傅里叶变换得到加窗奈奎斯特频谱:
对测量得到的被测系统输入输出信号进行加窗处理,式(2)为窗函数,然后进行傅里叶变换得到被测系统输入输出信号的加窗无理频谱Xw(f)和Yw(f),如式(3)和式(4):
式(2)中β取0.54/0.46或1,零点为窗函数的中点,Tw为窗宽,Tw≤T,x(t)和y(t)分别为被测系统输入输出信号,W0为正则化系数,W0=βTw/(β+1),设时间零点在窗的中间,则Xw(f)的相位与x(t)中点的相位一致,Yw(f)的相位与y(t)中点的相位一致。λ(λ≥1&λ∈I)为补零倍数(补零对计算结果没有影响,但可以使结果更加精细化)。
实际测量中得到的被测系统输入输出信号均为离散信号。设x(tk)为测得的x(t)的离散信号,y(tk)为测得的y(t)的离散信号,离散加窗傅里叶变换后得到的被测系统输入输出信号加窗奈奎斯特频谱Xw(fn)和Yw(fn)分别为:
式(5)和式(6)中fn为离散化频率,fn=n/λTw,频率间隔由1/λTw决定,Δt=1/fs,N=λTw/Δt,0≤n≤N-1,n∈I,tk为离散时间,tk=kΔt,0≤k≤N-1,k∈I,FFT[·]为快速傅里叶变换;
5)利用加窗奈奎斯特频谱计算被测系统的频响函数如式(7)所示:
式(7)中,Xw(fn)和Yw(fn)为输入信号和输出信号的加窗奈奎斯特频谱,由式(5)和式(6)得到,窗宽Tw=rwT(0<rw≤1),fn为离散化频率。
本方法还可包括:
6)在得到被测系统的频响函数曲线后,如果需要得到某段频率更加精确的测量值,则可以降低Δf值,重复上述步骤1)到步骤5),得到该段频率更加精细的频响函数测量值。
本发明提出的方法实施例1为测量抗混滤波器频响函数,抗混滤波器由三个萨伦凯低通滤波器串联组成,中心频率为500Hz。实施例2为测量悬臂梁的频响函数,所用悬臂梁为一端固定,另一端自由的钢梁。
实施例1采用的测量系统如图2所示,测量系统采用常规设备组成,包括NI数据采集器,电源,抗混滤波器,计算机等。利用计算机制作类脉冲激励信号,该类脉冲激励信号通过NI数据采集器输入到抗混滤波器中,同时用数据采集器采集抗混滤波器的输入信号x(t)输出信号y(t),在计算机中按本发明方法编程进行抗混滤波器输入输出信号的信号处理以及频响函数的计算。
实施例1的具体测量流程如图1所示,具体步骤如下:
1)确定被测系统频率范围:
实施例1中频率范围为[10Hz,1500Hz]
2)设计类脉冲激励信号:
21)确定类脉冲信号的频率间隔Δf:
实施例1中频率间隔为5Hz,信号周期T=1/Δf=0.2s;
22)测量的频响函数频率范围为[10Hz,1500Hz]
23)构造类脉冲激励信号作为被测系统的输入信号:
用一系列振幅相同相位为零的余弦信号进行叠加构造出类脉冲激励信号,作为输入信号,如式(1):
式(1)中,cm(t)为余弦波,fm为第m个余弦波的频率,fm=(m-1)Δf+fL,M为余弦波的总数,当m=M时,fm=fH,A为输入信号振幅,取A=3.4V;
3)测量被测系统输入信号与被测系统输出信号:
将式(1)中的类脉冲激励信号作为机械系统或电系统的输入信号,测量该系统输入信号和输出信号,记测量到的被测系统输入信号为x(t),被测系统输出信号为y(t),测量时采样频率fs=30000Hz
4)采用中心零点窗分析法对输入输出信号进行加窗傅里叶变换得到加窗奈奎斯特频谱:
对测量得到的被测系统输入输出信号进行加窗处理,式(2)为窗函数,然后进行傅里叶变换得到被测系统输入输出信号的加窗无理频谱Xw(f)和Yw(f),如式(3)和式(4):
式(2)中β取0.54/0.46,零点为窗函数的中点,Tw为窗宽,Tw=T,x(t)和y(t)分别为被测系统输入输出信号,W0为正则化系数,W0=βTw/(β+1),补零倍数λ=10。
实际测量中得到的输入输出信号均为离散信号。设x(tk)为测得的x(t)的离散信号,y(tk)为测得的y(t)的离散信号,离散加窗傅里叶变换后得到的输入输出信号的加窗奈奎斯特频谱Xw(fn)和Yw(fn)分别为:
式(5)和式(6)中fn为离散化频率,fn=n/λTw,频率间隔由1/λTw决定,Δt=1/fs,N=λTw/Δt,0≤n≤N-1,n∈I,tk为离散时间,tk=kΔt,0≤k≤N-1,k∈I,FFT[·]为快速傅里叶变换;
5)计算被测系统的频响函数:
采用中心零点窗离散傅里叶变换可以得到输入输出信号的加窗奈奎斯特频谱Xw(fn)和Yw(fn),其中fn为离散化频率,则离散频响函数H(fn)可由式(7)得到:
6)不需要对结果进行局部精细测量。
实施例1的测量结果如图3,图中虚线为抗混滤波器频响函数的测量值,实线为理论值,(a)为抗混滤波器频响函数幅值谱,(b)为抗混滤波器频响函数相位谱。由图3可以看出,利用本发明方法得到的抗混滤波器频响函数测量结果与理论值吻合,说明本发明的采用类脉冲激励的频响函数测量方法十分有效。
实施例2采用的测量系统如图4所示,测量系统采用常规设备组成,包括NI数据采集器,加速度传感器,压电陶瓷应变片,电源,悬臂梁,计算机等。利用计算机制作类脉冲激励信号,通过NI数据采集器输入到压电陶瓷应变片中,应变片在悬臂梁的固定端给悬臂梁提供一个类脉冲激励,加速度传感器装在悬臂梁的悬臂端,用数据采集器采集固定端的输入信号x(t)和悬臂端的输出信号y(t),在计算机中按本发明方法编程进行输入信号和输出信号的信号处理以及系统频响函数的计算。
实施例2的具体测量流程如图1所示,与实施例1有以下几点不同:
1)在实施例1的第1)步中,频率范围为[1Hz,50Hz]。
2)在实施例1的21)步中,频率间隔Δf=0.1Hz。
3)在实施例1的23)步中,A=37.4V。
4)在实施例1的3)步中,采样频率fs=2000Hz。
5)在实施例1的4)步中,补零倍数λ=20,窗宽Tw=0.8T。
6)在实施例1的6)步中,实施例2需对共振峰值处的频响函数曲线进行精细测量,即图5中的31Hz~41Hz处,减小Δf的值,分别取Δf=0.05Hz、0.025Hz、0.02Hz,测量悬臂梁的频响函数。
实施例2的测量结果如图5,可以看出,实施例中的悬臂梁的共振频率约为36Hz,为了精确测量共振频率处的幅值,减小Δf进行重复测量,结果如图6所示,图6中,图(a)为幅值谱,图(b)为相位谱,曲线1为Δf=0.1Hz的频响函数,曲线2为Δf=0.05Hz的频响函数,曲线3为Δf=0.025Hz的频响函数,曲线4为Δf=0.02Hz的频响函数。由图6可以看出,随着Δf的降低,悬臂梁共振峰处的频响函数幅值逐渐增大,当Δf≤0.025Hz,幅值不再增大,说明此时频响函数的幅值是准确的值,因此利用本专利提出的采用类脉冲激励的频响函数的测量方法可以准确测量出悬臂梁的频响函数的精确值。
Claims (2)
1.一种采用类脉冲激励测量频响函数的方法,其特点在于,包括以下步骤:
1)确定频率范围:确定所需的被测系统频响特性的频率范围[fl,fh],对于电系统或机械系统,频率范围为0~5000Hz;
2)构造类脉冲激励信号作为被测系统的输入信号:
21)确定类脉冲信号的频率间隔Δf:
<mrow>
<mi>&Delta;</mi>
<mi>f</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>T</mi>
</mfrac>
</mrow>
其中,T为类脉冲信号的周期,Δf的范围为0.005Hz~20Hz,T的范围为0.05s~200s;
22)测量的频响函数频率范围为[fL,fH],2Δf≤fL<fH-Δf,fH<fsmax/nH,
取fL和fH之间,包含fL和fH的整数频率点:
[fL,fH]:fL≤f≤fH,且f∈I
其中,fH-fL的值能被Δf整除,fsmax为测量仪器的最大采样频率,nH为每两个频率间隔点间的采样点数,nH的取值范围为[5,50];
23)用一系列振幅相同相位为零的余弦信号进行叠加构造出类脉冲激励信号,作为被测系统输入信号,如式(1):
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>A</mi>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mi>m</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中,cm(t)为余弦波,fm为第m个余弦波的频率,fm=(m-1)Δf+fL,M为余弦波的总数,当m=M时,fm=fH,A为输入信号振幅,取0V<A<40V;
3)测量被测系统输入信号与被测系统输出信号:
将式(1)中的类脉冲激励信号作为机械系统或电系统的输入信号,测量被测系统输入信号以及被测系统输出信号,记测量到的被测系统输入信号为x(t),被测系统输出信号为y(t),测量时采样频率fs满足fs=nHfH;
4)采用中心零点窗分析法对输入输出信号进行离散加窗傅里叶变换得到加窗奈奎斯特频谱:
对测量得到的被测系统输入输出信号进行加窗处理,式(2)为窗函数,然后进行傅里叶变换得到被测系统输入输出信号的加窗无理频谱Xw(f)和Yw(f),如式(3)和式(4):
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<mi>&beta;</mi>
<mo>+</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mi>t</mi>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&beta;</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>&le;</mo>
<mi>t</mi>
<mo><</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mn>2</mn>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>t</mi>
<mo>&Element;</mo>
<mi>e</mi>
<mi>l</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>X</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>W</mi>
<mn>0</mn>
</msub>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<msub>
<mi>&lambda;T</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>i</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Y</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>W</mi>
<mn>0</mn>
</msub>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<msub>
<mi>&lambda;T</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>i</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>f</mi>
<mi>t</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)中β取0.54/0.46或1,零点为窗函数的中点,Tw为窗宽,Tw≤T,x(t)和y(t)分别为被测系统输入输出信号,W0为正则化系数,W0=βTw/(β+1),设时间零点在窗的中间,则Xw(f)的相位与x(t)中点的相位一致,Yw(f)的相位与y(t)中点的相位一致;λ,λ≥1&λ∈I,为补零倍数;
设x(tk)为测得的x(t)的离散信号,y(tk)为测得的y(t)的离散信号,离散加窗傅里叶变换后得到的被测系统输入输出信号加窗奈奎斯特频谱Xw(fn)和Yw(fn)分别为:
<mrow>
<msub>
<mi>X</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&pi;T</mi>
<mi>w</mi>
</msub>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
</mrow>
</msup>
<mrow>
<msub>
<mi>W</mi>
<mn>0</mn>
</msub>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
<mo>&lsqb;</mo>
<mi>w</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Y</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&pi;T</mi>
<mi>w</mi>
</msub>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
</mrow>
</msup>
<mrow>
<msub>
<mi>W</mi>
<mn>0</mn>
</msub>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
<mo>&lsqb;</mo>
<mi>w</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式(5)和式(6)中,fn为离散化频率,fn=n/λTw,频率间隔由1/λTw决定,Δt=1/fs,N=λTw/Δt,0≤n≤N-1,n∈I,tk为离散时间,tk=kΔt,0≤k≤N-1,k∈I,FFT[·]为快速傅里叶变换;
5)利用加窗奈奎斯特频谱计算得到被测系统的频响函数如式(7)所示:
<mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Y</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>X</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
式(7)中,Xw(fn)和Yw(fn)为输入信号和输出信号的加窗奈奎斯特频谱,由式(5)和式(6)得到,窗宽Tw=rwT,0<rw≤1,fn为离散化频率。
2.如权利要求1所述方法,其特征在于,该方法还包括:
6)在得到被测系统的频响函数曲线后,降低Δf值,重复上述步骤1)到步骤5),得到该段频率更加精细的频响函数测量值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510594076.6A CN105242111B (zh) | 2015-09-17 | 2015-09-17 | 一种采用类脉冲激励的频响函数测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510594076.6A CN105242111B (zh) | 2015-09-17 | 2015-09-17 | 一种采用类脉冲激励的频响函数测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105242111A CN105242111A (zh) | 2016-01-13 |
CN105242111B true CN105242111B (zh) | 2018-02-27 |
Family
ID=55039826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510594076.6A Active CN105242111B (zh) | 2015-09-17 | 2015-09-17 | 一种采用类脉冲激励的频响函数测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105242111B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108845000A (zh) * | 2018-08-20 | 2018-11-20 | 四川大学 | 一种脉冲场指纹法测量管道缺陷的方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108490438B (zh) * | 2018-03-13 | 2021-11-16 | 深圳大学 | 利用无线射频技术对室内物体一次性成像的方法及系统 |
CN111596611A (zh) * | 2020-06-09 | 2020-08-28 | 上海交通大学 | 一种数控机床动力学特性测试分析系统 |
CN114113780A (zh) * | 2021-12-06 | 2022-03-01 | 国网湖南省电力有限公司 | 一种线性系统频域特性测试方法及系统 |
CN117909801A (zh) * | 2023-11-10 | 2024-04-19 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种工业机器人关节的少标签振动数据的故障诊断方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3808526A (en) * | 1972-06-12 | 1974-04-30 | Solartron Electronic Group | Frequency response analyser |
US4544883A (en) * | 1982-09-21 | 1985-10-01 | Michael Grimsland | Method and device for rapidly measuring frequency response of electronic test objects |
CN1969319A (zh) * | 2004-04-21 | 2007-05-23 | 诺基亚公司 | 信号编码 |
CN101420516A (zh) * | 2008-12-05 | 2009-04-29 | 硅谷数模半导体(北京)有限公司 | 用于影像缩放的自适应方法和装置 |
CN101509942A (zh) * | 2009-02-20 | 2009-08-19 | 华南理工大学 | 一种超声波换能器频率响应特性的检测方法 |
CN101561342A (zh) * | 2009-05-15 | 2009-10-21 | 北京工业大学 | 分时快速稳态正弦扫频激振频响函数测量系统及方法 |
DE102008060385A1 (de) * | 2008-10-01 | 2010-04-08 | Rohde & Schwarz Gmbh & Co. Kg | Zeitoptimierte Ermittlung des Frequenzganges von in einem Messobjekt erzeugten Harmonischen |
EP2314217A1 (en) * | 2009-10-23 | 2011-04-27 | Tallinn University of Technology | Method and device for fast measurement of frequency response with scalable short chirp signals |
CN102099856A (zh) * | 2008-07-17 | 2011-06-15 | 弗劳恩霍夫应用研究促进协会 | 具有可切换旁路的音频编码/解码方案 |
CN102903370A (zh) * | 2011-07-27 | 2013-01-30 | 株式会社东芝 | 多速率系统的高速频率响应识别法和装置 |
CN104065603A (zh) * | 2014-06-14 | 2014-09-24 | 温州大学 | 升余弦发送成型及多径环境下的正交相移键控信号自动识别方法 |
CN104079517A (zh) * | 2014-06-19 | 2014-10-01 | 哈尔滨工程大学 | 一种基于脉冲函数加权处理的连续相位调制实现方法 |
CN104198810A (zh) * | 2014-08-14 | 2014-12-10 | 深圳市爱普泰科电子有限公司 | 系统频率响应的测量方法 |
-
2015
- 2015-09-17 CN CN201510594076.6A patent/CN105242111B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3808526A (en) * | 1972-06-12 | 1974-04-30 | Solartron Electronic Group | Frequency response analyser |
US4544883A (en) * | 1982-09-21 | 1985-10-01 | Michael Grimsland | Method and device for rapidly measuring frequency response of electronic test objects |
CN1969319A (zh) * | 2004-04-21 | 2007-05-23 | 诺基亚公司 | 信号编码 |
CN102099856A (zh) * | 2008-07-17 | 2011-06-15 | 弗劳恩霍夫应用研究促进协会 | 具有可切换旁路的音频编码/解码方案 |
DE102008060385A1 (de) * | 2008-10-01 | 2010-04-08 | Rohde & Schwarz Gmbh & Co. Kg | Zeitoptimierte Ermittlung des Frequenzganges von in einem Messobjekt erzeugten Harmonischen |
CN101420516A (zh) * | 2008-12-05 | 2009-04-29 | 硅谷数模半导体(北京)有限公司 | 用于影像缩放的自适应方法和装置 |
CN101509942A (zh) * | 2009-02-20 | 2009-08-19 | 华南理工大学 | 一种超声波换能器频率响应特性的检测方法 |
CN101561342A (zh) * | 2009-05-15 | 2009-10-21 | 北京工业大学 | 分时快速稳态正弦扫频激振频响函数测量系统及方法 |
EP2314217A1 (en) * | 2009-10-23 | 2011-04-27 | Tallinn University of Technology | Method and device for fast measurement of frequency response with scalable short chirp signals |
CN102903370A (zh) * | 2011-07-27 | 2013-01-30 | 株式会社东芝 | 多速率系统的高速频率响应识别法和装置 |
CN104065603A (zh) * | 2014-06-14 | 2014-09-24 | 温州大学 | 升余弦发送成型及多径环境下的正交相移键控信号自动识别方法 |
CN104079517A (zh) * | 2014-06-19 | 2014-10-01 | 哈尔滨工程大学 | 一种基于脉冲函数加权处理的连续相位调制实现方法 |
CN104198810A (zh) * | 2014-08-14 | 2014-12-10 | 深圳市爱普泰科电子有限公司 | 系统频率响应的测量方法 |
Non-Patent Citations (5)
Title |
---|
Frequency Response Measurement for SAW Devices Based on Fourier Analysis of Impulse Response;KEN-YA HASHIMOTO等;《IEEE TRANSACTIONS ON SONICS AND ULTRASONICS》;19850131;第SU-32卷(第1期);76-77页 * |
Impulse method to calculate the frequency response of the electromagnetic forces on whirling cage rotors;A. Tenhunen等;《IEE Proc.-Electr. Power Appl》;20031007;第150卷(第6期);752-756页 * |
Smoothing the LPM Estimate of the Frequency Response Function via an Impulse Response Truncation Technique;Mikaya L. D. Lumori等;《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》;20140131;第63卷(第1期);214-220页 * |
频响函数测量的脉冲激励方法;钱济国;《振动、测试与诊断》;19940630;第14卷(第2期);32-39页 * |
频率响应函数实时处理程序的设计;周德忠等;《武汉冶金科技大学学报(自然科学版)》;19990330;第22卷(第1期);60-62页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108845000A (zh) * | 2018-08-20 | 2018-11-20 | 四川大学 | 一种脉冲场指纹法测量管道缺陷的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105242111A (zh) | 2016-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105242111B (zh) | 一种采用类脉冲激励的频响函数测量方法 | |
CN103308804B (zh) | 基于快速k-s变换电能质量扰动信号时频参数提取方法 | |
CN101403774B (zh) | 一种基于非同步采样的谐波分析方法 | |
US20090164142A1 (en) | Method and system for use in analyzing vibrations of a variable speed rotating body | |
CN102288807A (zh) | 一种测量电网电压闪变的方法 | |
CN108875706A (zh) | 基于滑动平均与能量归集的海洋结构时频分析方法 | |
CN106018956A (zh) | 一种加窗谱线插值的电力系统频率计算方法 | |
Xu et al. | A low-cost and field portable electromechanical (E/M) impedance analyzer for active structural health monitoring | |
CN114354112B (zh) | 一种叶片多阶耦合振动疲劳分析方法 | |
CN102353839A (zh) | 基于多层前馈神经网络的电力系统谐波分析方法 | |
CN104050147A (zh) | 将时域信号转换成频域信号的方法与系统 | |
CN109446552B (zh) | 多轴相关随机激励下结构疲劳寿命时域计算方法 | |
CN103983849A (zh) | 一种实时高精度的电力谐波分析方法 | |
CN101917163B (zh) | 一种改善非正弦周期信号电液激振控制波形的方法 | |
CN105675126A (zh) | 一种用于检测多频多源复杂稳定声场声压的新方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN109635399A (zh) | 一种振动加速度信号的加窗积分转换方法 | |
CN108181486A (zh) | 加速度信号的处理方法及装置 | |
CN104678170A (zh) | 一种基于谐波分析仪的电力谐波分析方法和谐波分析仪 | |
RU2363005C1 (ru) | Способ спектрального анализа полигармонических сигналов и устройство для его реализации | |
CN108491777A (zh) | 一种基于Morlet母小波的Lamb波空间采样信号拟合方法 | |
Shulin et al. | Application of improved EMD algorithm for the fault diagnosis of reciprocating pump valves with spring failure | |
CN115728513A (zh) | 一种新型冲击加速度测量装置及其测量方法 | |
Chen et al. | A phase difference measurement method based on the extended kalman filter for Coriolis mass flowmeters | |
CN111551785B (zh) | 基于无迹卡尔曼滤波的频率与谐波检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |