problem in for loop for plotting

Question

Mary 2017-4-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333397-problem-in-for-loop-for-plotting

评论： Mary 2017-4-3

hi, I have an EEG data and its matrix dimention is: 1068*4 now I want to undrestand which column relate to which brain wave, so I use Fourier transform and I write this code:

if true
load 'sig_x.dat'
fs = 500; % Hz
[N, M] = size(sig_x);
n = N/2;
f = (0:n-1)/n * fs/2;
for j = 1:M;
  subplot(M, 1, j)
  ff1 = fft(sig_x(:, j));
  eegfft = abs(ff1);
  figure(3)
  plot(f, eegfft(1:n))
  ylabel(['Wave ' num2str(j)])
end

but I don't get the spectrum of the first column. is there anything wrong whit this code ? or is there something I missed?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Joseph Cheng 2017-4-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/333397-problem-in-for-loop-for-plotting#answer_261564

编辑：Joseph Cheng 2017-4-3

do you mean you're missing the first spectrum because you put the figure after the first instance of subplot? you should move figure(3) before the for loop. with this you create the subplot in whatever is the active figure/axes for the first iteration, jump to figure(3) to then plot the first iteration fft. however with that you overwrite that first iteration with the subsequent subplots.

if you're commenting that there is no frequency response maybe taking the log10(eeft(1:n)) would help you see the smaller frequency responses.