Преобразование Фурье нескольких строк

Question

Я хочу принять БПФ из 10 измерений. Это означает, что у меня есть 10 строк в моих данных, каждая строка имеет размерность [1 * 2000]. В приведенном ниже коде я делаю неправильный путь. Может ли кто-нибудь сказать мне, какую ошибку я делаю?

load('fb2010');         % loading the data 
x = fb2010(3:1:15,:); 
y_filt = filter(b,a,x);       % filtering the received signal 
%%%%%%% Fourier transform 
nfft = length(y_filt); 
for i = 1:10 
res(i,:) = fft(y_filt(i,:),nfft)/ nfft;   %%%% Taking first fft and normalizing it 
end 
res2 = res(:,1:nfft/2+1);     %%%% taking single sided spectrum 
f = fs/2*linspace(0,1,nfft/2+1);   % choosing correct frequency axes 
figure, plot(f,abs(res2)); 

Answer 1

Функция filter по умолчанию передает передаточную функцию по первому измерению. Таким образом, вы фильтруете свои столбцы вместо своих строк. Чтобы получить правильный результат, используйте следующую строку:

filter(b,a,x,[],2); 

Затем вы можете опустить петлю с помощью fft с третьим аргументом, чтобы определить размер он работает вместе. Это будет следующая строка:

res = fft(y_filt,nfft,2); 

Answer 2

doc fft:

Y = fft(X,n,dim) , где тусклым является измерение. Полагаю, вы хотите dim = 2. Это избавится от цикла и должно работать. Я просто использовал его сам.