2

我需要编写一个四列矩阵（“g_Grid.r”、“g_Grid.t”、“g_Grid.b”、“g_Grid.ph”）

通常，我使用文件流写入文件：

ofstream fout; 
fout.open("GRID.dat");
 for (int m=0;m<N_PH;++m)
    {   
    for (int k=0;k<N_B;++k)
        {
          for (int j=0;j<N_T;++j)
            {   
            for (int i=0;i<N_R;++i)
                {
                  fout << setprecision(32) << g_Grid.r[i]<<" "<<g_Grid.t[j]<<" "<<g_Grid.b[k]<<" "<<g_Grid.ph[m]<< endl;

                }
            }
        }
    }
 fout.close();

它工作正常，但现在我正在处理非常大（长）的矩阵，编写格式化输出（“.txt”）需要很长时间。因为我真的不需要格式化文件，所以它更喜欢把它写成二进制文件。

问题：给定四个向量/数组（r，t，b，ph），如何将表写入二进制文件？（我后来用matlab读取这个文件，也需要把它读成二进制表）

编辑

通常，我使用以下简单代码编写二进制文件：

ofstream myFile (fileName, ios::out | ios::binary);
double val;
for (int m=0;m<N_PH;++m)
    {   
    for (int k=0;k<N_B;++k)
        {
          for (int j=0;j<N_T;++j)
            {   
            for (int i=0;i<N_R;++i)
                {
                  val = g_N.Amp[m][k][j][i];
                  myFile.write(reinterpret_cast<const char*>(&val), sizeof(val));
                }
            }
        }
    }
myFile.close();

但是这个只会产生“一列”，而稍后将由 matlab 读取。

4

1 回答 1

5

要将二进制数据写入文件，首先需要以二进制模式打开文件，需要使用fstream::write. 它接受一个char *参数，并从中写入多个字节 -char *如果数据不是数组char或指向的指针，则该参数需要强制转换char。

为了避免多次调用写入，最好将数据填充到一个结构中。

所以，以下应该做你需要的：

ofstream fout; 
fout.open("GRID.dat", ios::binary);
 for (int m=0;m<N_PH;++m)
    {   
    for (int k=0;k<N_B;++k)
        {
          for (int j=0;j<N_T;++j)
            {   
            for (int i=0;i<N_R;++i)
                {
                  struct X
                  {
                      double a, b, c, d;
                  } x;

                  x.a = g_Grid.r[i];
                  x.b = g_Grid.t[j];
                  x.c = g_Grid.b[k];
                  x.d = g_Grid.ph[m];
                  fout.write(reinterpret_cast<char *>(&x), sizeof(x));
                }
            }
        }
    }
 fout.close();
于 2013-06-12T08:55:46.470 回答