2

基本上,我试图将一个巨大的 1D 向量拆分为通过函数传递的给定大小的块。该函数应该返回一个二维向量,然后我可以访问不同的块。我找到了一个合适的算法来做到这一点,但是,它在 Matlab 中,我不明白如何将元素放在 2D 向量中。

MatLab 代码:

function f = block(v, N, M)

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

numblocks = (lastblockstart-1)/M + 1;

f = zeros(numblocks,N);

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

这是我在 C++ 中的尝试(如果不好,请见谅):

vector<iniMatrix> subBlocks(vector<int>& 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 - (maxblockstart-1 % M);

int numblocks = (lastblockstart-1)/M + 1;

vector<int> subBlock;
vector<iniMatrix> block;

for(unsigned i=0; (i < numblocks); i++)
{
    for(unsigned j=0; (j < N); j++)
    {
        subBlock.push_back(theData[(i-1*M+j)]);

        block.push_back(subBlock);
    }
}

return block;
}

这段代码可以编译,但是,当我试图输出块的大小时,我得到:分段错误:11 ..有什么想法吗?

通过函数传递的数据是:N = 600 M = 200

我希望有人可以帮助我,谢谢:)

4

1 回答 1

5

在 c 和 c++ 中,数组索引从零开始。如果您使用数组长度函数,请记住这一点。因此,您应该将i-1乘数替换为i并从零开始计数。循环条件:

for (unsigned j=1; (j < N); j++)

将开始1和结束于N-1-- 总共的N-1时间。但,

for (unsigned j=0; (j < N); j++)

将开始0和结束于N-1-- 总共N时间。

在 Matlab 中:for 循环从第一个边界开始,然后在第二个边界结束

如果你觉得必须从 index 开始1

for (unsigned j=1; (j < N+1); j++)

将进行N迭代,同时仍从1. 但请记住,当您在 C/C++ 中声明数组时,第一个元素的索引为零。

于 2012-08-16T05:00:19.463 回答