gas / KOSYAK1

uses crt,gnk;
const ay:array [1..10] of real=(10,-5,11,11,-6,7,-7,11,-11,-6);
ax:array [1..10] of real=(0,4,8,13,19,25,35,45,71,117);
T:word=60;
N:word=200;
M:real=0;
v1:real=70;
s:real=0.4;
dy:real=0.2;
var x,y,yk,ykp:array [1..200] of real;i:integer;
ku,nuk,c:word;q,v2,y1,y2:real;

PROCEDURE PV(i,N:word;YK,YPK,dY:real;var C:word;var V2:real);
{CONST
dY=1e-5;}
BEGIN
if i=1 then C:=0;
if abs(YPK-YK)<=dY then inc(C);
if i=N then V2:=C/i
END;


begin
clrscr;
gklf(60,ax,ay,x,y);
{sag(0,v1,m,s,min(y),max(y),ku,q);}
q:=2*s;ku:=round((max(y)-min(y))/q);
for i:=1 to 200 do begin
kvn(ku,q,min(y),y[i],nuk,yk[i]);
yk[i]:=yk[i]+gauss(M,s);
end;
ku:=round((max(yk)-min(yk))/q);
for i:=1 to 200 do begin
kvn(ku,q,min(yk),yk[i],nuk,ykp[i]);
pv(i,N,yk[i],ykp[i],dy,c,v2);
end;

gra(1,0,x,y,x,ykp,'1','2','3');
end.




sag(1,v1,m,s,min(y),max(y),ku,q);

{Єў ­в㥬 бЁЈ­ «}
for i:=1 to 200 do begin
kvn(ku,q,min(y),y[i],nuk,yk[i]);
end;

{ўў®¤Ё¬ Ї®¬Ґег}
for i:=1 to 200 do begin
ykp[i]:=yk[i]+gauss(M,s);
end;

{ku:=round((max(yk)-min(yk))/q);}

for i:=1 to 200 do begin
kvn(ku,q,min(yk),yk[i],nuk,ypk[i]);
pv(i,N,yk[i],ypk[i],dy,c,v2);
end;

gra(1,1,x,y,x,yk,'source signal','discrete signal','');
gra(1,1,x,yk,x,ykp,'discrete signal','gauss noise signal','');
gra(1,1,x,ykp,x,ypk,'gauss noise signal','discrete-noise-discrete','');
gra(1,1,x,y,x,ypk,'source signal','discrete-noise-discrete','');