Question

我想在这个问题前加上这样一个事实，即在编码方面我很缺乏经验，所以这个问题的解决方案可能比我一直在尝试的容易得多。我有一个定义为vector<float> phase;的矢量'phas'，它有7987200个元素，我想把这个矢量重新排列成60个元素的133120个矢量（称为line2定义为vector<long double> line2;）。然后应将60的每个向量一个接一个地放置在定义为vector< vector<long double> > RFlines2;和RFlines2.resize(7987200);的向量'RFlines2'的向量中。我想用'phas'的元素填充60个元素向量中的每一个，用128分隔。例如，60个元素的第一个向量将用phas [0]，phas [128]，phas [256]，...填充。 .phas [7680]。然后用ps [1]，phas [129]，phas [257]，...... phas [7681]等填充60的第二个载体。我目前的代码如下：

    for(int x = 0; x<133120; x++){
        if((x == 128 || x == 7680+128 || x == (7680*a)+128)){
            x = 7680*a;
            a = a + 1;
        }
        int j = x;
        for(int i = 0; i<60;i++){
            line2.pushback(i);
            line2[i] = phas[j];
            j = j + 128;
        }
        cout<<"This is x: "<<x<<endl;
        RFlines2[x] = line2;
        line2.clear();
        }

然而，在外循环的128次迭代之后（已经创建了128个60个向量并且已经使用了76个来自phas的元素），我需要x值跳转到7680以避免从已经存在的phas中放置元素因为当x = 128时，60的下一个矢量的第一个元素将是相位[128]，它已被用作第一个矢量60的第二个元素。然后在另一个128 x迭代之后使用60的下一个矢量，我需要x值跳到15,360等等。上面的代码是我最近的尝试，但是当我尝试在RFlines2中的每个60矢量上执行fftw时，如下所示：

    int c = 0;
   for(int x = 0; x < 133120; x++){
    //cout<<x<<endl;

   fftw_plan p2;

   inter = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * W);
   outter = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * W);

    /* cast elements in line to type fftw_complex */
    for (int i = 0; i <60; i++) {
    //cout<<i<<endl;
    //inter[i][0] = phas[i];
    //inter[x][0] = zlines[x]; 
    inter[i][0] = RFlines2[x][i];
    inter[i][1] = 0;
  }

    p2 = fftw_plan_dft_1d(60, inter, outter, FFTW_FORWARD, FFTW_ESTIMATE);
    fftw_execute(p2);
    //inter[x][0].clear();
    for(int u = 0; u<60;u++){
        if(u == 0){
        cout<<' '<<outter[0][0]<<' '<<c++<<endl;
        }
    }
    fftw_free(inter);
    fftw_free(outter);
    fftw_destroy_plan((p2));
   }

程序在显示外部[0] [0] 128次后崩溃。任何想法如何解决这一问题？另外，如果我说的任何内容都没有意义，请告诉我，我会尽力澄清。提前谢谢！

-Mike

Answer 1

我不知道为什么你的代码崩溃，因为我在这里看不到整个代码。但我将建议一种分散数据和管理载体的方法。

（但有一个重要的警告：你不应该使用向量（至少不是向量的向量）来完成这项任务;你最好使用一维向量并自己管理二维索引。但这是一个表现的东西，并且不影响正确性。）

我建议您填写RFLines2 :(我还没有尝试过此代码，因此可能无效。）

// first, build the memory for RFLines2...
vector<vector<long double>> RFLines2 (133120, vector<long double>(60));

// assuming a "phase" vector...
for (unsigned i = 0; i < 7987200; ++i)
{
    unsigned const row = (i / (128 * 60)) * 128 + (i % (128 * 60)) % 128;
    unsigned const col = (i % (128 * 60)) / 128;
    RFLines[row][col] = phase[i];
}

这样你不需要line2中间人。

其余代码“ 应该 ”有效。（顺便说一句，我根本不理解for 上的内部u循环。你想在那里做什么？）

用c ++重组向量

1 个答案: