diophatine.cpp

广义预测算法MATLAB转化为VC++的程序,可用于三容水箱的数据的采集.

  #include "matlib.h"
  #pragma hdrstop
  
  #include "diophatine.h"
  
  
  Mm diophatine(Mm a, Mm b, Mm N1, i_o_t, Mm& e__o, Mm& F__o, Mm& g__o, Mm& H__o) {
    begin_scope
    a.setname("a"); b.setname("b"); N1.setname("N1"); 
    dMm(e); dMm(F); dMm(g); dMm(H); dMm(na); dMm(nb); dMm(abar); dMm(f); dMm(h); dMm(j_); dMm(j_diophatine_v0); dMm( \
      i_); dMm(i_diophatine_v1); dMm(ff); dMm(i_diophatine_v2); 
    
    #line 1 "c:/matcom45/diophatine.m"
    call_stack_begin;
    #line 1 "c:/matcom45/diophatine.m"
    // nargin, nargout entry code
    double old_nargin=nargin_val; if (!nargin_set) nargin_val=3.0;
    nargin_set=0;
    double old_nargout=nargout_val; if (!nargout_set) nargout_val=4.0;
    nargout_set=0;
    
    // translated code
    
    #line 2 "c:/matcom45/diophatine.m"
_   na = length(a)-1.0;
    #line 3 "c:/matcom45/diophatine.m"
_   nb = length(b)-1.0;
    #line 4 "c:/matcom45/diophatine.m"
_   a = ctranspose(a);
    #line 4 "c:/matcom45/diophatine.m"
_   b = ctranspose(b);
    #line 5 "c:/matcom45/diophatine.m"
_   abar = conv((BR(1.0),-1.0),a);
    #line 6 "c:/matcom45/diophatine.m"
_   f = -abar(colon(2.0,1.0,abar.size()));
    #line 7 "c:/matcom45/diophatine.m"
_   F = f;
    #line 8 "c:/matcom45/diophatine.m"
_   e(1.0) = 1.0;
    #line 9 "c:/matcom45/diophatine.m"
_   g(1.0) = e(1.0)*b(1.0);
    #line 10 "c:/matcom45/diophatine.m"
_   h = b(colon(2.0,1.0,b.size()));
    #line 11 "c:/matcom45/diophatine.m"
_   H = h;
    #line 12 "c:/matcom45/diophatine.m"
_   j_diophatine_v0 = colon(2.0,1.0,N1); int j_diophatine_i0;
    for (j_diophatine_i0=0;j_diophatine_i0<j_diophatine_v0.cols();j_diophatine_i0++) {
      forelem(j_,j_diophatine_v0,j_diophatine_i0);
      #line 13 "c:/matcom45/diophatine.m"
_     e(j_) = f(1.0);
      #line 14 "c:/matcom45/diophatine.m"
_     i_diophatine_v1 = colon(1.0,1.0,na); int i_diophatine_i1;
      for (i_diophatine_i1=0;i_diophatine_i1<i_diophatine_v1.cols();i_diophatine_i1++) {
        forelem(i_,i_diophatine_v1,i_diophatine_i1);
        #line 15 "c:/matcom45/diophatine.m"
_       ff(i_) = f(i_+1.0)-abar(i_+1.0)*f(1.0);
        #line 16 "c:/matcom45/diophatine.m"
      }
      #line 17 "c:/matcom45/diophatine.m"
_     ff(na+1.0) = -abar(abar.size())*f(1.0);
      #line 18 "c:/matcom45/diophatine.m"
_     f = ff;
      #line 19 "c:/matcom45/diophatine.m"
_     F = (BR(F),semi,
      f);
      #line 20 "c:/matcom45/diophatine.m"
      //end
      
      #line 22 "c:/matcom45/diophatine.m"
      //for j=2:N1
      #line 23 "c:/matcom45/diophatine.m"
_     g(j_) = e(j_)*b(1.0)+h(1.0);
      #line 24 "c:/matcom45/diophatine.m"
_     i_diophatine_v2 = colon(2.0,1.0,nb); int i_diophatine_i2;
      for (i_diophatine_i2=0;i_diophatine_i2<i_diophatine_v2.cols();i_diophatine_i2++) {
        forelem(i_,i_diophatine_v2,i_diophatine_i2);
        #line 25 "c:/matcom45/diophatine.m"
_       h(i_-1.0) = e(j_)*b(i_)+h(i_);
        #line 26 "c:/matcom45/diophatine.m"
      }
      #line 27 "c:/matcom45/diophatine.m"
_     h(nb) = e(j_)*b(b.size());
      #line 28 "c:/matcom45/diophatine.m"
_     H = (BR(H),semi,
      h);
      #line 29 "c:/matcom45/diophatine.m"
    }
    call_stack_end;
    
    // nargin, nargout exit code
    nargin_val=old_nargin; nargout_val=old_nargout;
    
    // function exit code
    a.setname(NULL); b.setname(NULL); N1.setname(NULL); 
    e__o=e; F__o=F; g__o=g; H__o=H; 
    return x_M;
    end_scope
  }
  
  
  Mm diophatine(Mm a, Mm b, Mm N1) {
    begin_scope
    double old_nargin=nargin_val; nargin_val=3.0; nargin_set=1;
    double old_nargout=nargout_val; nargout_val=1.0; nargout_set=1;
    dMm(e__o); dMm(F__o); dMm(g__o); dMm(H__o); 
    diophatine(a, b, N1, i_o, e__o, F__o, g__o, H__o);
    nargout_val=old_nargout;
    nargin_val=old_nargin;
    return(e__o);
    end_scope
  }
  
  Mm diophatine(Mm a, Mm b, Mm N1, i_o_t, Mm& e__o, Mm& F__o) {
    begin_scope
    double old_nargin=nargin_val; nargin_val=3.0; nargin_set=1;
    double old_nargout=nargout_val; nargout_val=2.0; nargout_set=1;
    dMm(g__o); dMm(H__o); 
    diophatine(a, b, N1, i_o, e__o, F__o, g__o, H__o);
    nargout_val=old_nargout;
    nargin_val=old_nargin;
    return x_M;
    end_scope
  }
  
  Mm diophatine(Mm a, Mm b, Mm N1, i_o_t, Mm& e__o, Mm& F__o, Mm& g__o) {
    begin_scope
    double old_nargin=nargin_val; nargin_val=3.0; nargin_set=1;
    double old_nargout=nargout_val; nargout_val=3.0; nargout_set=1;
    dMm(H__o); 
    diophatine(a, b, N1, i_o, e__o, F__o, g__o, H__o);
    nargout_val=old_nargout;
    nargin_val=old_nargin;
    return x_M;
    end_scope
  }