Plotting of FRF curves

Question

Shrinivas Gombi 2011-6-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/10362-plotting-of-frf-curves

I want to plot FRF vs Frequency graph.I have a matrix 'a' of acceleration & a force matrix 'f' both of 30000X1 size, collected in 3secs with a time step of 0.0001 sec.I am using the following steps to plot FRF. a=[ ]; f=[ ]; ffa=fft(a); fff = fft(f); frf, H=ffa./fff; aH =abs(H); I want to plot FRF vs Frequency graph using log(aH) on Y-axis. Pl tell me how to specify x-axis?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Rick Rosson 2011-7-11

Please format your code:

http://www.mathworks.com/matlabcentral/about/answers/

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

bym 2011-6-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/10362-plotting-of-frf-curves#answer_14283

see my answer here http://www.mathworks.com/matlabcentral/answers/10306-frf-plotting

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Shrinivas Gombi 2011-6-28

Dear Sir,

I had only asked this question in that link also.Thank u very much for ur answer.But I have tried many ways of setting x- axis but not getting correctly.I request u to kindly provide me the exact answer.

Prof.S L Gombi

请先登录，再进行评论。

Answer 2

bym 2011-6-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/10362-plotting-of-frf-curves#answer_14418

在 MATLAB Online 中打开

Frequency resolution is sampling frequency divided by FFT points. So for a one sided spectrum, it is (directly from the documentation)

NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);

you can then plot as

semilogy(f,aH)

or

loglog(f,aH)

3 个评论
显示 1更早的评论隐藏 1更早的评论

bym 2011-6-29

L = 30000; y = aH; NFFT = 32768; Fs = 10000

Shrinivas Gombi 2011-7-5

Dear SIr,

Kindly check my steps:

a = [30000X1 readings];

ffa = fft(a);

F =[30000X1 readings];

fff = fft(F);

H = ffa./fff;

aH = abs(H);

L = 30000;

NFFT = 32768;

Fs = 10000;

y = aH;

Y = fft(y,NFFT)/L;

f = Fs/2*linspace(0,1,NFFT/2+1);

semilogy(f,aH)

After the Execution I am getting the following comment

??? Error using ==> semilogy

Vectors must be the same lengths.

Kindly let me know where is the mistake and how to correct it.As I required to submit my PhD thesis at the earliest, pl help me in solving this problem.Waiting for your reply

Prof> Shrinivas Gombi

请先登录，再进行评论。