io.cpp

vc代用matcom编程

  #include "matlib.h"
  #pragma hdrstop
  
  #include "io.h"
  
  
  Mm io(Mm x, Mm imf) {
    begin_scope
    x.setname("x"); imf.setname("imf"); 
    dMm(ort); dMm(lx); dMm(n); dMm(s); dMm(i_); dMm(i_io_v0); dMm(j_); dMm(j_io_v1); 
    
    #line 1 "d:/matcom45/samples/emd/io.m"
    call_stack_begin;
    #line 1 "d:/matcom45/samples/emd/io.m"
    // nargin, nargout entry code
    double old_nargin=nargin_val; if (!nargin_set) nargin_val=2.0;
    nargin_set=0;
    double old_nargout=nargout_val; if (!nargout_set) nargout_val=1.0;
    nargout_set=0;
    
    // translated code
    
    
    #line 3 "d:/matcom45/samples/emd/io.m"
    // ort = IO(x,imf) computes the index of orthogonality
    #line 4 "d:/matcom45/samples/emd/io.m"
    //
    #line 5 "d:/matcom45/samples/emd/io.m"
    // inputs : - x    : analyzed signal
    #line 6 "d:/matcom45/samples/emd/io.m"
    //          - imf  : empirical mode decomposition
    
    #line 8 "d:/matcom45/samples/emd/io.m"
_   lx = size(imf,2.0);
    #line 9 "d:/matcom45/samples/emd/io.m"
_   n = size(imf,1.0);
    
    #line 11 "d:/matcom45/samples/emd/io.m"
_   s = 0.0;
    
    #line 13 "d:/matcom45/samples/emd/io.m"
_   i_io_v0 = colon(1.0,1.0,n); int i_io_i0;
    for (i_io_i0=0;i_io_i0<i_io_v0.cols();i_io_i0++) {
      forelem(i_,i_io_v0,i_io_i0);
      #line 14 "d:/matcom45/samples/emd/io.m"
_     j_io_v1 = colon(1.0,1.0,n); int j_io_i1;
      for (j_io_i1=0;j_io_i1<j_io_v1.cols();j_io_i1++) {
        forelem(j_,j_io_v1,j_io_i1);
        #line 15 "d:/matcom45/samples/emd/io.m"
_       if (istrue(i_!=j_)) {
          #line 16 "d:/matcom45/samples/emd/io.m"
_         s = s+abs(sum(times(imf(i_,c_p),imf(j_,c_p)))/sum(power(x,2.0)));
          #line 17 "d:/matcom45/samples/emd/io.m"
        }
        #line 18 "d:/matcom45/samples/emd/io.m"
      }
      #line 19 "d:/matcom45/samples/emd/io.m"
    }
    
    #line 21 "d:/matcom45/samples/emd/io.m"
_   ort = 0.5*s;
    
    call_stack_end;
    
    // nargin, nargout exit code
    nargin_val=old_nargin; nargout_val=old_nargout;
    
    // function exit code
    x.setname(NULL); imf.setname(NULL); 
    return ort;
    end_scope
  }