comp_iwfac.cc

来自「Matlab时频分析工具箱,希望能对大家有所帮助啊」· CC 代码 · 共 25 行

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

DEFUN_DLD (comp_iwfac, args, ,
  "Computes inverse window factorization.\n\
  Usage: c=comp_iwfac(gf,L,a,M);\n\
  Yeah.")
{

  ComplexMatrix gf = args(0).complex_matrix_value();
  const int L = args(1).int_value();
  const int a = args(2).int_value();
  const int M = args(3).int_value();

  const int R=gf.rows()*gf.columns()/L;

  ComplexMatrix g(L,R);
  
  iwfac((ltfat_complex*)gf.data(),L,R,a,M,(ltfat_complex*)g.data());

  return octave_value (g);
}