comp_wfac.cc

Matlab时频分析工具箱,希望能对大家有所帮助啊

#include <octave/oct.h>
#include "config.h"
#include "../src/wfac.c"
#include "../src/gcd.c"

DEFUN_DLD (comp_wfac, args, ,
  "Computes window factorization.\n\
  Usage: c=comp_wfac(g,a,M);\n\
  Yeah.")
{

  ComplexMatrix g = args(0).complex_matrix_value();
  const int a = args(1).int_value();
  const int M = args(2).int_value();
  const int L = g.rows();
  const int R = g.columns();

  const int b = L/M;

  int h_a, h_m;

  const int c=gcd(a, M,&h_a, &h_m);
  const int p=a/c;
  const int q=M/c;
  const int d=b/p;

  ComplexMatrix gf(p*q*R,c*d);
  
  wfac((ltfat_complex*)g.data(),L,R,a,M,(ltfat_complex*)gf.data());

  return octave_value (gf);
}