Hidrológiai Közlöny 1974 (54. évfolyam)

8. szám - Szőllősi-Nagy András: Optimális előrejelző-függvény meghatározása a Wiener-féle extrapoláció elmélet alkalmazásával

Szöllősi-Nagy A.: Optimális előrejelző-függvény Hidrológiai Közlöny 1974. 8. sz. 33Í) ( START ) \ INPUT / \jjjnj 4. ábra. A számitógép program folyamatábrája Fig. 4. Block diagram of the computer programme az x(t) és h(t) függvények diszkrét konvolució]Sit hajtja végre [27], Az eljárás formális paraméterei: a'[Ü : p~\ a konvolválandó függvény diszkrét értékeiből alkotott vektor, A[0 : M] a súlyfüggvény vektora, ye[ 0 \ p+u] eredmény vektor, p az x diszkrét értékeinek száma — 1, u a súlyfüggvény ordinátáinak száma. Az eljárásban szereplő aktuális tömbök deklarációjában valamennyi index alsó korlátja nulla kell legyen. A 2. táblázat a számítás teljes programját foglalja 2. táblázat Optimális előrejelző-függvény számításának programja ALGOL nyelven RAZDAN-3 gépi reprezentációban Programme in ALGOL for Computing the optimál forecastiiid function Table 2. begin integer 11, m, i; real fix, fiy; ujra: input (n, in); if n =0 and in=

0 then goto vege else begin array x, y [l:n], txx, fyy, fxy [0:n

—

1], rxx, ryy, rxy [0:n

—

1], hopt [0:in]; boolean korr, pred; input (array x, array y); cov (x, x, fxx, n); cov (y, y, fyy, n); cov (x

, y, 1'xy, n); input (korr); il korr then begin fix: =fxx [0]; fix: =fyy [0]; lor i:>=0 step 1 until n

—

1 do begin rxx[i]:=fxx[i]/fix; ryyfij: —fyy[l]/fiy; rxy[i]: =fxy[i]/sqrt(fix • fiy) end: lines 5; spaces 25; text korrelaciofuggvenyek; lines 5; text k rxx ryy rxy; lines 3; for i: —0 step 1 until entier ri/2 do begin oulput(i:ll,rxx[il:7:4,ryy[i],rxy[il); line end end; lines 20; spaces 12; text optimális előrejelző függvény; lines 5; text t liopt; lines 3; tor i: =—(J step 1 until in do begin hopt[i]:=it fxy[i]«=0 then 0 else fxy[i]/fix; output (i:ll,hopt[i]:10:5); line end; Input (pred); if not pred then goto nine else begin integer p; input (p); begin array csap[0:p],ye[0:p + m]; input (arra? csap); diko (csap, liopt, ye, p, m); lines 20; spaces 10; text ido; spaces 10; text elorejelzett vízhozamok; lines 5; for i:=0 step 1 until p + m do begin output (i: 10,ye[i]: 10:5); line end end end end: nine: goto ujra; vege; end WHFZ;