MATLAB中的音频需要帮助

时间:2010-08-04 14:10:40

标签: matlab audio waveform wave energy

我正在尝试编写一个.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); index超出界限因为numel(x)= 1。   ==>中的错误12岁的巨人       能量(帧)=能量(帧)+ X(I)^ 2;

非常感谢任何帮助。

2 个答案:

答案 0 :(得分:2)

您应该使用f而不是x,因为f是从.wav文件加载的实际信号。变量x可能只是您工作区中的其他标量,这就是您收到错误的原因。

应该对您的代码进行一些其他更正/改进。首先,作为Paul R pointed out,您需要更正计算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

答案 1 :(得分:0)

不应该这一行:

numFrames=(numSamples/1);

类似于:

numFrames=(numSamples / frameWidth);

或可能:

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