1

オーディオ トラックからエネルギー機能を抽出する .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 であるため、インデックスが範囲外です。エラー ==> 12 の myrms energy(frame)=energy(frame)+x(i)^2;

どんな助けでも大歓迎です。

4

2 に答える 2

2

fの代わりにを使用する必要があります。実際の信号は .wav ファイルからロードされるためですx。f変数xはおそらくワークスペース内の他のスカラーであったため、見たエラーが発生しました。

コードに加える必要があるその他の修正/改善がいくつかあります。まず、Paul R が指摘したように、 の計算方法を修正する必要がありますnumFrames。次に、energyゼロのベクトルとして初期化する必要があります。3 番目に、内側の for ループを 1 行のベクトル化操作に減らすことができます。

コードを書き直す方法は次のとおりです（編集：コメントに基づいて、ループで計算されたいくつかの余分な変数を保存するためにコードを更新しました）：

[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
于 2010-08-04T14:16:22.950 に答える
0

この行はすべきではありません：

numFrames=(numSamples/1);

次のようになります。

numFrames=(numSamples / frameWidth);

またはおそらく：

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

?

于 2010-08-04T14:16:26.890 に答える