Требуется помощь со звуком в MATLAB - PullRequest
Новая
       49

Требуется помощь со звуком в MATLAB

1 голос
/ 04 августа 2010

Я пытаюсь написать файл .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 и получаюследующая ошибка:

???Попытка доступа к х (2);индекс вне границ, потому что цифра (х) = 1.Ошибка в ==> мирмах при энергии 12 (кадр) = энергия (кадр) + x (i) ^ 2;

Любая помощь будет принята с благодарностью.

Ответы [ 2 ]

2 голосов
/ 04 августа 2010

Вы должны использовать f вместо x, поскольку f - это фактический сигнал, загруженный из вашего файла .wav. Переменная x, вероятно, была просто другим скаляром в вашем рабочем пространстве, поэтому вы получили сообщение об ошибке.

Есть несколько других исправлений / улучшений, которые следует внести в ваш код. Во-первых, как Пол Р указал , вам нужно исправить способ вычисления numFrames. Во-вторых, energy следует инициализировать как вектор нулей. В-третьих, вы можете уменьшить внутренний цикл for до однострочной векторизованной операции.

Вот как я бы переписал ваш код ( РЕДАКТИРОВАТЬ: основываясь на комментариях, я обновил код для сохранения нескольких дополнительных переменных, вычисленных в цикле ): 

[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
0 голосов
/ 04 августа 2010

Не должна ли эта строка: 

numFrames=(numSamples/1);

быть что-то вроде: 

numFrames=(numSamples / frameWidth);

или, возможно: 

numFrames=((numSamples + frameWidth - 1) / frameWidth);

...