Я пытаюсь записать.m файл для извлечения энергетических функций из аудиотрека, но я, кажется, испытываю затруднения при его реализации:
% Formula for calculating RMS
[f, fs, nb] = wavread('Three.wav');
frameWidth=441; %10ms
numSamples=length(x);
numFrames=(numSamples/1);
energy(frame)=0;
for frame=1:numFrames,
startSample=(frame-1)*frameWidth+1;
endSample=startSample+frameWidth-1;
% Calculate frame energy
for i=startSample:endSample
energy(frame)=energy(frame)+x(i)^2;
end
end
Я петляю в MATLAB и получаю следующую ошибку:
??? Предпринятый для доступа к x (2); индексируйте за пределы потому что numel (x) =1. Ошибка в ==> myrms в 12 энергиях (кадр) =energy (кадр) +x (i) ^2;
Любая справка очень ценилась бы.
Вы должны использовать f
вместо x
, поскольку f
- это фактический сигнал, загруженный из вашего .wav файла. Переменная x
, вероятно, была просто каким-то другим скаляром в вашей рабочей области, поэтому вы и получили ошибку, которую увидели.
Есть еще несколько исправлений/улучшений, которые следует внести в ваш код. Во-первых, как указал Paul R, вам нужно исправить способ вычисления numFrames
. Во-вторых, energy
должен быть инициализирован как вектор нулей. В-третьих, вы можете сократить внутренний цикл for до однострочной векторной операции.
Вот как я бы переписал ваш код (EDIT: Основываясь на комментариях, я обновил код, чтобы сохранить несколько дополнительных переменных, вычисляемых в цикле):
[y, fs, nb] = wavread('Three.wav'); %# Load the signal into variable y
frameWidth = 441; %# 10 msec
numSamples = length(y); %# Number of samples in y
numFrames = floor(numSamples/frameWidth); %# Number of full frames in y
energy = zeros(1,numFrames); %# Initialize energy
startSample = zeros(1,numFrames); %# Initialize start indices of frame
endSample = zeros(1,numFrames); %# Initialize end indices of frame
for frame = 1:numFrames %# Loop over frames
startSample(frame) = (frame-1)*frameWidth+1; %# Starting index of frame
endSample(frame) = frame*frameWidth; %# Ending index of frame
frameIndex = startSample(frame):endSample(frame); %# Indices of frame samples
energy(frame) = sum(y(frameIndex).^2); %# Calculate frame energy
end
Не должна ли эта строка:
numFrames=(numSamples/1);
быть чем-то вроде:
numFrames=(numSamples / frameWidth);
или, возможно:
numFrames=((numSamples + frameWidth - 1) / frameWidth);
?