Курсовые / 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.