Курсовые / curstau / CLM

               {      ”Ё«мва  ЉЂ‹ЊЂЌЂ   }

           {_X(k+1/k+1)=_X(k+1/k)+Kr*[y(k+1)-_y(k+1/k)] }


  unit CLM;
    interface
        uses ALG,GB,CANON;
        procedure RICC (n:integer;l:real;A,B,Po,Q:matr;var P,R:matr);
        procedure CALM (n:integer;A,C:matr;var Kr:matr);

  implementation


    procedure EST (n,k:integer;A,B,C:matr;Yo:vect;u,y:real;var _X:vect);
             var i,j                      :integer;
                 F,Ft,FI,P1,Po            :matr;
                 Xkk,Xk1k,Xk1k1,Yk1k,Kr   :vect;

    begin
      OBSF(n,k,A,C,F);
      for i:=1 to k do for j:=1 to n do Ft[j,i]:=F[i,j];
      MUL(n,k,n,Ft,F,P1);
      INV(n,P1,Po];
      MUL(n,n,k,Po,Ft,FI);
      MULV(n,k,FI,Yo,Xkk);
      MULV(n,n,A,Xkk,Xk1k);
      MULV(1,n,C,Xk1k,Yk1k);
      CALM(n,A,C,Kr);
      for i:=1 to n do
      _X[i]:=Xk1k[i]+B[i,1]*u+Kr[i,1]*(y-Yk1k[1]);

    end;
    {-----------------------------------------------------}

  end.