cv04-katastrofa

poc = 1;
for x1 = 2:0.2:3.8
  for y1 = 0:0.2:1.8
    z = (x1.*x1+x1.*y1)./(x1-y1);
    VstupSiete(poc,1)=x1;
    VstupSiete(poc,2)=y1;
    VstupSiete(poc,3)=z;
    poc=poc+1;
   end;
end;

if (max(abs(VstupSiete(:,3))))>1
  N_VstupSiete = VstupSiete / max(abs(VstupSiete))
else
  N_VstupSiete = VstupSiete;
end;

PocetVstup=2;
Pocet1Skryta=11;
PocetSkryta=7;

% nahodna inicializacia vah
Vahy1 = rand(pocetVstup+1, pocet1Skryta);
Vahy1(3,:) = ones(1);
Vahy2 = rand(pocet1Skryta+1, pocet2Skryta);
Vahy2(pocet1Skryta+1,:) = ones(1);
Vahy3 = rand(pocet2Skryta+1, 1);
Vahy3(pocet2Skryta+1,:) = ones(1);

Potencial1Skrytej=ones(1,pocet1Skryta)
Potencial2Skrytej=ones(1,pocet2Skryta)
PotencialVystupnej=onew(1,1)

% inicializacia VSTUPY
Vystupy0 = ones(1,3);
Vystupy0(3)=-1;
Vystupy1 = ones(1,pocet1Skryta=1);
Vystupy1(pocet1Skryta+1)=-1;
Vystupy2 = ones(1,pocet2Skryta+1);
Vystupy2(pocet2Skryta+1)=-1;
Vystupy3 = ones(1,1);

VysledkySiete=ones(100,1)

for a=1:50
  for r=1:100

    Vystupy0(1)=N_VstupSiete(r,1);
    Vystupy0(2)=N_VstupSiete(r,2);

    for i=1:Pocet1Skryta
      Potencial1Skrytej(1,i)=Vahy1(:,i)*Vystupy0(1,:)
      Vystupy1(1,i)=tanh(Potencial1Skrytej(1,i)):
    end;

    for i=1:Pocet2Skryta
      Potencial2Skrytej(1,i)=Vahy2(:,i)*Vystupy1(1,:)
      Vystupy2(1,i)=tanh(Potencial2Skrytej(1,i)):
    end;

    PotencialVystupnej(1,1)=Vahy3(:,i)*Vystupy2(1,:)
    Vystupy3(1,1)=tanh(PotencialVystupnej(1,1)):

    ChybaNaVystupe = ones(1,1);
    ChybaNaVystupe = (1-tanh(PotencialVystupnej(1,1))*tanh(PotencialVystupnej(1,1)))*(N_VstupSiete(r,3)-Vystupy3);

    Chyba2Vrstva = ones(1,Pocet2Skryta);
    for i=1:Pocet2Skryta
      Chyba2Vrstva(1,i) = (1-tanh(Potencial2Skrytej(1,i))*tanh(Potencial2Skrytej(1,i)))*(Vahy3(1,i)*ChybaNaVystupe(1,1));
    end;

    pomocnySucin = ones(1, pocet1Skryta)
    Chyba1Vrstva = ones(1,pocet1Skryta);
    for i=1:pocet1Skryta
      for j=1:pocet2Skryta
        pomocnySucin(1,i)=Vahy2(i,j)*Chyba2Vrstva(1,j)
      end;

      Chyba1Vrstva = (1,i) = (1-tanh(Potencial2Skrytej(1,i))*tanh(Potencial2Skrytej(1,i)))*(Vahy3(1,i)*pomocnySucin(1,1));
    end;

    % zmena vah
    UciaciPomer = 0.1;
    for i=1:(pocetVstup+1)
      for j=1:pocet1Skryta
        pomocnySucin(i,j) = UciaciPomer*Chyba1Vrstva(1,j)*Vystupy0(1:j));
        noveVahy(i,j)=Vahy1(i,j)+pomocnySucin);
      end;
    end;

    celkovaChyba = N_VstupSiete(r,3)*N_VstupSiete(r,3) - Vystupy3(1,1))*0.5;

  end;

  chybaSiete = chybaSiete + celkovaChyba;

end;