0

1Dベクトルと2つのサイズを受け取るMATLAB関数が与えられました。この関数は、データをブロックに分割し、最後に2Dベクトル内に格納します。この関数のC++バージョンを作成していますが、C ++関数の(誤った)結果がMATLAB関数の(正しい)結果と一致しません。

Matlab関数:

function f = block(v, N, M)

% This function separates the vector
% into blocks.  Each block has size N.
% and consecutive blocks differ in
% their starting positions by M
%
% Typically
%   N = 30 msec (600 samples)
%   M = 10 msec (200 samples)

  n = length(v);
  maxblockstart = n - N + 1;
  lastblockstart = maxblockstart - mod(maxblockstart-1 , M);

  % Remove the semicolon to see the number of blocks
  % numblocks = (lastblockstart-1)/M + 1
  numblocks = (lastblockstart-1)/M + 1;

  %f = zeros(numblocks,N);

  for i = 1:numblocks
  for j = 1:N
   f(i,j) = ((i-1)*M+j);
  end 
 end

この例では、((i-1)* M + j)の結果を出力しているだけで、MatLabでは次の結果が得られます(例)。

1 201 401 601 .. 1001 1201 .. 1401 .. 1601 .. 1801 ..

そして、これが私のC++関数です。

vector<iniMatrix> Audio::subBlocks(vector<float>& theData, int N, int M)
{
    // This method splits the vector into blocks
    // Each block has size N.
    // and consecutive blocks differ 
    
    int n = theData.size();
    int maxblockstart = n - N+1;
    int lastblockstart = maxblockstart - mod(maxblockstart-1, M);
    
    int numblocks = (lastblockstart-1)/M + 1;
    
    vector<float> subBlock;
    vector<iniMatrix> block;
    
    for(int i=1; (i < numblocks); i++)
    {
        for(int j=1; (j < N); j++)
        {
            cout << ((i-1)*M+j);
        }
        
    }

    return block;
}

私がこれから得た結果:

1 2 3 4 .. 7 8 9 .. 131415など。

PS

iniMatrixは、floatのベクトルの単なるtypdefです。

別の注意、変数:

n

maxblockstart

lastblockstart

numblocks

MatlabプログラムとC++ではすべて同じ値なので、forループと関係があると思います。

誰か提案がありますか?

4

2 に答える 2

1

ご希望のアルゴリズムを正しく理解できたかどうかを確認させてください。n = 10 N = 3 M = 2 とします。

これは 0、2、4、6 になるはずですよね?(8+M が範囲 0 <= x < n の外にあるため、8 は収まらないため)。

したがって、maxblockstart は 7 = 10 - 3、つまり n - N にする必要があります。

lastblockstart は 6 = 7 - 7 % 2 である必要があります。つまり、maxblockstart - maxblockstart % M

numblocks は 4 = 6/2 + 1、つまり lastblockstart/M + 1 である必要があります

次のようにコードを変更すると、正しい結果が得られるようです (これは紙の上でのみ機能し、コンパイルや実行は試していません.....):

vector<iniMatrix> Audio::subBlocks(vector<float>& theData, int N, int M)
{
  int n = theData.size();
  int maxblockstart = n - N;
  int lastblockstart = maxblockstart - (maxblockstart % M);

  int numblocks = (lastblockstart)/M + 1;

  vector<iniMatrix> block;    

  for(int i=0; (i < numblocks); i++)
  {
     vector<float> subBlock;
      for(int j=0; (j < N); j++)
      {
          subBlock.push_back(theData[i*M+j]); //cout << (i*M+j);
      }  
      block.push_back(subBlock);
  }

  return block;
}

試してみる...

C++ のインデックスはゼロから始まるため、結果を MATLAB と比較すると混乱する可能性があることに注意してください。したがって、次のことを試してください。

1)行cout << (i*M+j);を次のように変更しますcout << theData[i*M+j];

2) 次のテストを試してください。

vector<float> test;
for(int i=0; i<=10000; i++)
  test.push_back(i);

Audio::subBlocks(test, 1023, 200);
于 2012-09-18T20:28:13.163 に答える
1

よくわかりませんが、出力の 1D ベクトルの値をどこに書きますか?

多分あなたはこれを変更する必要があります:

  for i = 1:numblocks
  for j = 1:N
   f(i,j) = v((i-1)*M+j);
  end 
 end
于 2012-09-18T20:27:04.110 に答える