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

{ “а ў­Ґ­ЁҐ ђ€ЉЉЂ’€ }


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

implementation

{ аҐиҐ­ЁҐ га ў­Ґ­Ёп ђЁЄЄ вЁ }

procedure RICC (n:integer;l:real;A,B,Q:matr;var P,R:matr);
var i,j :integer;
At,AtP,AtPA,AtPB,AtPBR,
Bt,BtP,BtPA,BtPB,Q1,P1,P2,P3,P4,DP :matr;
lb,ndp :real;
begin
for i:=1 to n do for j:=1 to n do
begin
At[i,j]:=A[j,i];Bt[i,j]:=B [j,i];
Q1[i,j]:=Q[i,j];P [i,j]:=Q1[i,j];
end;
repeat
MUL(1,n,n,Bt,P,BtP);
MUL(1,n,n,BtP,B,BtPB);
MUL(1,n,n,BtP,A,BtPA);
lb:=1/(l+BtPB[1,1]);
for i:=1 to n do for j:=1 to n do R[i,j]:=lb*BtPA[i,j];
MUL(n,n,n,At,P,AtP);
MUL(n,n,n,AtP,A,AtPA);
MUL(n,n,1,AtP,B,AtPB);
MUL(n,1,n,AtPB,R,AtPBR);
for i:=1 to n do for j:=1 to n do P1[i,j]:=-AtPBR[i,j];
ADD(n,n,Q,AtPA,P2);
ADD(n,n,P2,P1,P3);
for i:=1 to n do for j:=1 to n do P4[i,j]:=-P3[i,j];
ADD(n,n,P,P4,DP);
ndp:=0;
for i:=1 to n do for j:=1 to n do ndp:=ndp+sqr(DP[i,j]);
ndp:=sqrt(ndp);
for i:=1 to n do for j:=1 to n do P[i,j]:=P3[i,j];
until ndp<0.00001;
end;
{----------------------------------------}

{ дЁ«мва Љ «¬ ­  }

procedure CALM (n:integer;A,C:matr;var Kr:matr);
var i,j :integer;
At,CA,CAt,B,E,Q,P,R :matr;
begin
MUL(1,n,n,C,A,CA);
for i:=1 to n do for j:=1 to n do
begin
At[i,j]:=A [i,j]; A[i,j]:=At[j,i]; CAt[i,j]:=CA[j,i];
B [i,j]:=CAt[i,j]; E[i,j]:=0; E[i,i]:=1;Q[i,j]:=0;
end;
RICC(n,1,A,B,Q,P,R);
for i:=1 to n do for j:=1 to n do
begin
P[i,j]:=P[j,i]; Kr[i,j]:=R[j,i];
end;
end;
{----------------------------------------}


end.