Ошибка при попытке открыть аудиофайл с помощью FFT в matlab

Question

Я пытаюсь использовать алгоритм FFT для построения графика простого аудиофайла с величиной против частоты, он работает, но я пытаюсь открыть окно сигнала с целью устранения утечки частоты на боковых лепестках основных частот. Он просто дает мне ошибку

здесь код без кадрирования (который работает):

[x,fs]=audioread('E:\voice_hum.wav'); 

x1=x; 

plot_points = 100000; 
h=fft(x1,plot_points); 
freq = h(1:50001); %dump data after half FS 
mag = abs(freq); 
freq = 0.5 * fs *(0:50000)/50000; %frequency on x-axis 
plot(freq,mag) 

вот код, когда я попытался реализации окно Хэмминга:

[x,fs]=audioread('E:\voice_hum.wav'); 

x1 = x(1:100); 

w = hamming(100); 
plot_points = 100000; 

x1 = x1.*w'; %window the signal 
h=fft(x1,plot_points); 
freq = h(1:50001); %dump data after half FS 

mag = abs(freq); 
freq = 0.5 * fs *(0:50000)/50000; 

plot(freq,mag) 

его дает мне эту ошибку:

Error using .* Matrix dimensions must agree.

Error in lab8b (line 12) x1 = x1.*w'; %window the signal

Я немного смущен, кто-то знает, делать?

Answer 1

Вы используете умножение точек между w (вектор-столбец) и x1 (вектор-столбец, если у вас есть только один аудиоканал).

Для умножения двух векторов с умножением точек (то есть a.*b) они должны иметь одинаковую размерность, в то время как вы используете a.*b', который умножает число на столбец и строку.

Таким образом, вы можете использовать x1'*w (matrix multiplication) или x1.*w с тем же результатом.

Answer 2

Используйте вместо этого x1 = x1.*w;. Всегда помните, чтобы проверить размер вашего массива перед операцией.