Question

现在几天我一直想弄清楚什么是错的但是不能，因为我可以找到关于使用Mathgl中的Dens（）函数进行错误处理的非常有限的文档。以下是我的问题。我正在尝试从输入信号构建频谱图。我创建了一个M * N数组来保存我的DFT数据。其中M =没有行=用于该信号的窗口总数和N =每个窗口中样本的总数。有一个X矩阵，其大小（M，1）保存时间数据，另一个矩阵/ y大小（N，1）保存标准化的频率数据（在db中），其中最小db为0，max db = 1 。大小（M，N）的DFT数据数组保持的值范围为数据类型float的0-1。以下是我的代码。 vector fftstore = dSTFT（＆amp; signal，signalLength，windowSize，hopSize）; int num_rows = fftstore.size（）/ windowSize;

  float *m_gl = new float[num_rows * windowSize];
   memset( m_gl,0,num_rows*windowSize*sizeof(float));
   float *m_glstart;
   m_glstart=m_gl;
  //t m_gl[num_rows][windowSize]={0};
   float *begin_loc;
 begin_loc=&fftstore[0];
 float *end_loc;
 end_loc=&fftstore[0]+windowSize;
   //t *dest;
  //st=m_gl;
   float *offset;
  for(int kk=0; kk<num_rows; kk++)
    {

      offset=m_glstart+(kk*windowSize+0);
      //oat offset=&m_gl[kk*windowSize+0];
       std::copy(begin_loc, end_loc, offset);
     begin_loc+=windowSize;
      end_loc+=windowSize;

   }
   //r(vector<float>::iterator it = fftstore.begin(); it != fftstore.end(); it++) {
  //out << *it << endl;
   //

  //converting M_gl to double

  double *m_gldob = new double[num_rows * windowSize];
   memset( m_gldob,0,num_rows*windowSize*sizeof(double));
  for (int id=0;id<num_rows;id++){
    for(int ij=0;ij<4096;ij++){
       m_gldob[id * windowSize + ij]=double(m_gl[id * windowSize + ij]);
    }
  }


  //mgl plotting
  double xdat[num_rows][1];
     memset( xdat,0,num_rows*1*sizeof(int));
   for (int xt=1;xt<=num_rows;xt++){
    xdat[xt-1][1]=((100*(2*xt-1))/(2*xt));
   }

  double ydat[windowSize][1];
   memset( ydat,0,windowSize*1*sizeof(int));
   for (int yf=1;yf<=windowSize;yf++){
    ydat[yf-1][1]=(yf*(1/windowSize));
   }

  double xst;
 xst=xdat[0][0];
 double xend=xdat[num_rows-1][1];
  double yst=ydat[0][0];
  double yend=ydat[windowSize-1][1];
    double *px;
 double *py;
  double *pdat;
   px=&xdat[0][0];
   py=&ydat[0][0];
  pdat=m_gldob;
   mglGraph gr;
   mglData mgl_x;
  mglData mgl_y;
  mglData mgl_dat;



    mgl_dat.Link(m_gldob,num_rows*windowSize,1);//c *
   mgl_dat.Rearrange(num_rows,windowSize);
    mgl_x.Link(px,num_rows,1);
  mgl_y.Link(py,windowSize,1);
  mgl_x.Norm();//c
  mgl_y.Norm();//c
   gr.SetQuality(6);
  gr.Alpha(true);

   gr.SetRange('x',xst,xend);
   gr.SetRange('y',yst,yend);
   gr.SetRange('c',0.0, 1.0);



   gr.Box();
  gr.Dens(mgl_x,mgl_y,mgl_dat);
  gr.Colorbar();
  gr.WriteBMP("/home/koyel/test1.bmp");

  delete[] m_gl;

  return 0;

这段代码运行给我一个如下所示的输出： http://i.stack.imgur.com/8uAq6.png

你能告诉我什么出错吗？谢谢！

Mathgl Dens（）图不起作用

0 个答案: