c++ - 将 Vector 拆分为 block - 奇怪的结果

Question

我得到了一个 MATLAB 函数，它接受一个一维 vector 和两个大小，然后该函数将数据分成 block ，最后将它们存储在一个二维 vector 中。我一直在编写此函数的 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 .. 13 14 15 等..

附言

iniMatrix 只是一个浮点 vector 的类型定义..

另一个注意事项，变量:

n

maxblockstart

lastblockstart

numblocks

在 Matlab 程序和 C++ 中都具有相同的值，所以我认为这与 for 循环有关..

有人有什么建议吗？

Answer 1

让我看看我是否正确理解了您想要的算法。 说 n = 10 N = 3 M = 2

这应该产生 0、2、4、6 对吗？ (因为 8 不适合，因为 8+M 超出范围 0 <= x < n)。

所以，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;
}

试一试...

请注意，将结果与 MATLAB 进行比较可能会造成混淆，因为 C++ 索引是从零开始的。因此尝试以下操作；

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);