我的FRF怎么了?

36 views (last 30 days)
Manuel Lozano
Manuel Lozano 上23 May 2020
编辑: Alberto Mora 上23 May 2020
Hello,
From a time domain response and excitation signal, I'm trying to compute and plot an FRF function. I've tried using both FFT and the built-in function modalfrf to compute the FRF. However the resulting FRFs look slightly different and differ greatly in magnitude. Notice how the modalfrf curve has a somewhat "downhill" look (orange one) while the one I computed using my code (blue) looks more "flat". Does anyone know what might be the issue here?
%% My FRF code using FFT
tempx=fft(a5);
tempf=fft(F);
x1 =(tempx);
f1=(tempf);
l =长度(x1)/2;
X=x1(1:L);
Fu=f1(1:L);
frf = abs(x./fu);
freq=linspace(0,Nq,length(FRF));
figure
loglog(freq,FRF);
网格
Xlabel('Frecuencia(Hz)')
ylabel('FRF (dB)')
xlim([3 70])
ylim([yinf ysup])
%% Using modalfrf function
[FRF,f]=modalfrf(F,a5,Fs,hann(50000),0.1*50000);
xlim([3 70])
FRF=abs(FRF);
loglog(f,FRF,'Linewidth', 1)
xlim([3 70])
ylim([yinf ysup])
Xlabel("Frecuencia (Hz)")
ylabel("FRF (g/N)")
Thank you in advance.

Accepted Answer

Alberto Mora
Alberto Mora 上23 May 2020
编辑:Alberto Mora 上23 May 2020
The average slope of the two lines suggest that in one case you are integrating the accelartion into velocity (orangeFRF) so you are plotting somthing like (m/s)/N (this aspect could be confimed also by the 3 order of magnitude).
While in the second case (blueFRF) it seems that you are plotting the "accelerance" in g/N.
如果您想检查此方面,只需计算以下内容:
blueFRF / ( j*omega )%整合加速度(bluefrf)[(m/s^2)]一次,以获取[(m/s)/n]中测量速度(OrangeFRF)
Theoretically you should get the orangeFRF. You can also compute the viceversa, so orangeFRF * ( j*omega )
I check my theory at 10Hz by hands, and the result confirm what I suspect:
2e-2/(2*pi*10)^2 = 5.0661E-006%由于我们正在研究规模,因此我们可以考虑j。如果您也对阶段感兴趣,请考虑j。
因此,您可能必须注意将AcceleranceFRF的整合到VelocityFRF中。
Check more carefully this aspect and let me know.
Morover,当您计算FRF时,不要忘记将窗户应用于原始信号。我注意到您正在使用Hann窗口,但是如果您选择了正确的窗口,请检查最新状态。通常,对于带有撞击锤的FRF,最好在两个信号上应用指数窗口。如果您使用来自恒定源的随机噪声(例如电动力振荡器),则最好使用Hann窗口。
Last point, don't forget to check the coherence of your FRF that gives you info about the quality of your measurements. The command " modalfrf ”也给你这个输出。看医生了info.
If this answer help you to solve you problem, please consider to accept the answer as a sign of gratitude.
Best regards
1 Comment
Manuel Lozano
Manuel Lozano 上23 May 2020
解决了问题!非常感谢您的时间和明确的解释!关于窗口,数据是从振动筛测试中绘制的,这就是为什么我使用汉宁窗口的原因。我也将检查连贯函数。
此致。

Sign in to comment.

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

开始狩猎!

Translated by