我的FRF怎么了?
36 views (last 30 days)
Show older comments
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.
0 Comments
Accepted Answer
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
More Answers (0)
See Also
s manbetx 845
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
开始狩猎!