gas / ZAG

uses crt,gnk;
const ay:array [1..18] of real=(0,-5,11,11,-6,7,-7,11,-8,0,1,3,9,10,9,3,1,0);
ax:array [1..18] of real=(1,4,8,13,19,25,35,45,71,117,135,140,145,150,155,160,165,170);
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,ypk: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(190,ax,ay,x,y);

{q:=2;}

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

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','');

end.