errorapart.cpp

对于转台控制系统中的精测误差分析与补偿算法

_   a(2.0,4.0) = cos(alpha)*transpose(cos(alpha*2.0));
    #line 66 "d:/matcom45/samples/errorapart.m"
_   a(2.0,5.0) = cos(alpha)*transpose(sin(alpha*2.0));
    #line 67 "d:/matcom45/samples/errorapart.m"
_   a(3.0,1.0) = a(1.0,3.0);
    #line 68 "d:/matcom45/samples/errorapart.m"
_   a(3.0,2.0) = a(2.0,3.0);
    #line 69 "d:/matcom45/samples/errorapart.m"
_   a(3.0,3.0) = sin(alpha)*transpose(sin(alpha));
    #line 70 "d:/matcom45/samples/errorapart.m"
_   a(3.0,4.0) = sin(alpha)*transpose(cos(alpha*2.0));
    #line 71 "d:/matcom45/samples/errorapart.m"
_   a(3.0,5.0) = sin(alpha)*transpose(sin(alpha*2.0));
    #line 72 "d:/matcom45/samples/errorapart.m"
_   a(4.0,1.0) = a(1.0,4.0);
    #line 73 "d:/matcom45/samples/errorapart.m"
_   a(4.0,2.0) = a(2.0,4.0);
    #line 74 "d:/matcom45/samples/errorapart.m"
_   a(4.0,3.0) = a(3.0,4.0);
    #line 75 "d:/matcom45/samples/errorapart.m"
_   a(4.0,4.0) = cos(alpha*2.0)*transpose(cos(alpha*2.0));
    #line 76 "d:/matcom45/samples/errorapart.m"
_   a(4.0,5.0) = cos(alpha*2.0)*transpose(sin(alpha*2.0));
    #line 77 "d:/matcom45/samples/errorapart.m"
_   a(5.0,1.0) = a(1.0,5.0);
    #line 78 "d:/matcom45/samples/errorapart.m"
_   a(5.0,2.0) = a(2.0,5.0);
    #line 79 "d:/matcom45/samples/errorapart.m"
_   a(5.0,3.0) = a(3.0,5.0);
    #line 80 "d:/matcom45/samples/errorapart.m"
_   a(5.0,4.0) = a(4.0,5.0);
    #line 81 "d:/matcom45/samples/errorapart.m"
_   a(5.0,5.0) = sin(alpha*2.0)*transpose(sin(alpha*2.0));
    #line 82 "d:/matcom45/samples/errorapart.m"
_   f(1.0) = transpose(ax)*ones(size(ax));
    #line 83 "d:/matcom45/samples/errorapart.m"
_   f(2.0) = transpose(ax)*transpose(cos(alpha));
    #line 84 "d:/matcom45/samples/errorapart.m"
_   f(3.0) = transpose(ax)*transpose(sin(alpha));
    #line 85 "d:/matcom45/samples/errorapart.m"
_   f(4.0) = transpose(ax)*transpose(cos(alpha*2.0));
    #line 86 "d:/matcom45/samples/errorapart.m"
_   f(5.0) = transpose(ax)*transpose(sin(alpha*2.0));
    #line 87 "d:/matcom45/samples/errorapart.m"
_   errlong = inv(a)*transpose(f);
    
    #line 89 "d:/matcom45/samples/errorapart.m"
_   alpha = transpose(pitch*2.0*3.14159265359);
    #line 90 "d:/matcom45/samples/errorapart.m"
_   a(1.0,1.0) = length(nn);
    #line 91 "d:/matcom45/samples/errorapart.m"
_   a(1.0,2.0) = cos(alpha)*transpose(ones(size(alpha)));
    #line 92 "d:/matcom45/samples/errorapart.m"
_   a(1.0,3.0) = sin(alpha)*transpose(ones(size(alpha)));
    #line 93 "d:/matcom45/samples/errorapart.m"
_   a(1.0,4.0) = cos(2.0*alpha)*transpose(ones(size(alpha)));
    #line 94 "d:/matcom45/samples/errorapart.m"
_   a(1.0,5.0) = sin(2.0*alpha)*transpose(ones(size(alpha)));
    #line 95 "d:/matcom45/samples/errorapart.m"
_   a(2.0,1.0) = a(1.0,2.0);
    #line 96 "d:/matcom45/samples/errorapart.m"
_   a(2.0,2.0) = cos(alpha)*transpose(cos(alpha));
    #line 97 "d:/matcom45/samples/errorapart.m"
_   a(2.0,3.0) = cos(alpha)*transpose(sin(alpha));
    #line 98 "d:/matcom45/samples/errorapart.m"
_   a(2.0,4.0) = cos(alpha)*transpose(cos(alpha*2.0));
    #line 99 "d:/matcom45/samples/errorapart.m"
_   a(2.0,5.0) = cos(alpha)*transpose(sin(alpha*2.0));
    #line 100 "d:/matcom45/samples/errorapart.m"
_   a(3.0,1.0) = a(1.0,3.0);
    #line 101 "d:/matcom45/samples/errorapart.m"
_   a(3.0,2.0) = a(2.0,3.0);
    #line 102 "d:/matcom45/samples/errorapart.m"
_   a(3.0,3.0) = sin(alpha)*transpose(sin(alpha));
    #line 103 "d:/matcom45/samples/errorapart.m"
_   a(3.0,4.0) = sin(alpha)*transpose(cos(alpha*2.0));
    #line 104 "d:/matcom45/samples/errorapart.m"
_   a(3.0,5.0) = sin(alpha)*transpose(sin(alpha*2.0));
    #line 105 "d:/matcom45/samples/errorapart.m"
_   a(4.0,1.0) = a(1.0,4.0);
    #line 106 "d:/matcom45/samples/errorapart.m"
_   a(4.0,2.0) = a(2.0,4.0);
    #line 107 "d:/matcom45/samples/errorapart.m"
_   a(4.0,3.0) = a(3.0,4.0);
    #line 108 "d:/matcom45/samples/errorapart.m"
_   a(4.0,4.0) = cos(alpha*2.0)*transpose(cos(alpha*2.0));
    #line 109 "d:/matcom45/samples/errorapart.m"
_   a(4.0,5.0) = cos(alpha*2.0)*transpose(sin(alpha*2.0));
    #line 110 "d:/matcom45/samples/errorapart.m"
_   a(5.0,1.0) = a(1.0,5.0);
    #line 111 "d:/matcom45/samples/errorapart.m"
_   a(5.0,2.0) = a(2.0,5.0);
    #line 112 "d:/matcom45/samples/errorapart.m"
_   a(5.0,3.0) = a(3.0,5.0);
    #line 113 "d:/matcom45/samples/errorapart.m"
_   a(5.0,4.0) = a(4.0,5.0);
    #line 114 "d:/matcom45/samples/errorapart.m"
_   a(5.0,5.0) = sin(alpha*2.0)*transpose(sin(alpha*2.0));
    #line 115 "d:/matcom45/samples/errorapart.m"
_   f(1.0) = transpose(ax)*ones(size(ax));
    #line 116 "d:/matcom45/samples/errorapart.m"
_   f(2.0) = transpose(ax)*transpose(cos(alpha));
    #line 117 "d:/matcom45/samples/errorapart.m"
_   f(3.0) = transpose(ax)*transpose(sin(alpha));
    #line 118 "d:/matcom45/samples/errorapart.m"
_   f(4.0) = transpose(ax)*transpose(cos(alpha*2.0));
    #line 119 "d:/matcom45/samples/errorapart.m"
_   f(5.0) = transpose(ax)*transpose(sin(alpha*2.0));
    #line 120 "d:/matcom45/samples/errorapart.m"
_   errshort = inv(a)*transpose(f);
    #line 121 "d:/matcom45/samples/errorapart.m"
    //end
    
    call_stack_end;
    
    // nargin, nargout exit code
    nargin_val=old_nargin; nargout_val=old_nargout;
    
    // function exit code
    x1.setname(NULL); x2.setname(NULL); x3.setname(NULL); x4.setname(NULL); x5.setname(NULL); x6.setname(NULL); x7.setname( \
      NULL); x8.setname(NULL); x9.setname(NULL); x10.setname(NULL); x11.setname(NULL); x12.setname(NULL); x13.setname( \
      NULL); x14.setname(NULL); x15.setname(NULL); x16.setname(NULL); x17.setname(NULL); x18.setname(NULL); x19.setname( \
      NULL); x20.setname(NULL); x21.setname(NULL); x22.setname(NULL); x23.setname(NULL); x24.setname(NULL); x25.setname( \
      NULL); x26.setname(NULL); x27.setname(NULL); x28.setname(NULL); x29.setname(NULL); x30.setname(NULL); x31.setname( \
      NULL); x32.setname(NULL); x33.setname(NULL); x34.setname(NULL); x35.setname(NULL); x36.setname(NULL); x37.setname( \
      NULL); x38.setname(NULL); x39.setname(NULL); x40.setname(NULL); x41.setname(NULL); x42.setname(NULL); x43.setname( \
      NULL); x44.setname(NULL); x45.setname(NULL); x46.setname(NULL); x47.setname(NULL); x48.setname(NULL); x49.setname( \
      NULL); x50.setname(NULL); x51.setname(NULL); 
    cx__o=cx; errlong__o=errlong; errshort__o=errshort; 
    return x_M;
    end_scope
  }
  
  
  Mm errorapart(Mm x1, Mm x2, Mm x3, Mm x4, Mm x5, Mm x6, Mm x7, Mm x8, Mm x9, Mm x10, Mm x11, Mm x12, Mm x13, Mm  \
    x14, Mm x15, Mm x16, Mm x17, Mm x18, Mm x19, Mm x20, Mm x21, Mm x22, Mm x23, Mm x24, Mm x25, Mm x26, Mm x27, Mm  \
    x28, Mm x29, Mm x30, Mm x31, Mm x32, Mm x33, Mm x34, Mm x35, Mm x36, Mm x37, Mm x38, Mm x39, Mm x40, Mm x41, Mm  \
    x42, Mm x43, Mm x44, Mm x45, Mm x46, Mm x47, Mm x48, Mm x49, Mm x50, Mm x51) {
    begin_scope
    double old_nargin=nargin_val; nargin_val=51.0; nargin_set=1;
    double old_nargout=nargout_val; nargout_val=1.0; nargout_set=1;
    dMm(cx__o); dMm(errlong__o); dMm(errshort__o); 
    errorapart(x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16, x17, x18, x19, x20, x21, x22, \
       x23, x24, x25, x26, x27, x28, x29, x30, x31, x32, x33, x34, x35, x36, x37, x38, x39, x40, x41, x42, x43, x44, \
       x45, x46, x47, x48, x49, x50, x51, i_o, cx__o, errlong__o, errshort__o);
    nargout_val=old_nargout;
    nargin_val=old_nargin;
    return(cx__o);
    end_scope
  }
  
  Mm errorapart(Mm x1, Mm x2, Mm x3, Mm x4, Mm x5, Mm x6, Mm x7, Mm x8, Mm x9, Mm x10, Mm x11, Mm x12, Mm x13, Mm  \
    x14, Mm x15, Mm x16, Mm x17, Mm x18, Mm x19, Mm x20, Mm x21, Mm x22, Mm x23, Mm x24, Mm x25, Mm x26, Mm x27, Mm  \
    x28, Mm x29, Mm x30, Mm x31, Mm x32, Mm x33, Mm x34, Mm x35, Mm x36, Mm x37, Mm x38, Mm x39, Mm x40, Mm x41, Mm  \
    x42, Mm x43, Mm x44, Mm x45, Mm x46, Mm x47, Mm x48, Mm x49, Mm x50, Mm x51, i_o_t, Mm& cx__o, Mm& errlong__o)  \
    {
    begin_scope
    double old_nargin=nargin_val; nargin_val=51.0; nargin_set=1;
    double old_nargout=nargout_val; nargout_val=2.0; nargout_set=1;
    dMm(errshort__o); 
    errorapart(x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16, x17, x18, x19, x20, x21, x22, \
       x23, x24, x25, x26, x27, x28, x29, x30, x31, x32, x33, x34, x35, x36, x37, x38, x39, x40, x41, x42, x43, x44, \
       x45, x46, x47, x48, x49, x50, x51, i_o, cx__o, errlong__o, errshort__o);
    nargout_val=old_nargout;
    nargin_val=old_nargin;
    return x_M;
    end_scope
  }