clear
# 1 класс: (1,1),(2,4),(7,18),(12,15)
# 2 класс: (-3,5),(-4,9),(-5,12),(-15,2)
# 3 класс: (7,-3),(14,-18),(3,-22),(7,-13)
# 4 класс: (-1,-1),(-4,-9),(-3,-22),(-12,-15)
points = [1,2,7,12, -3,-4,-5,-15, 7,14,3,7, -1,-4,-3,-12;
1,4,18,15, 5,9,12,2, -3,-18,-22,-13, -1,-9,-22,-15];
Y = [1,0,0,0;
1,0,0,0;
1,0,0,0;
1,0,0,0;
0,1,0,0;
0,1,0,0;
0,1,0,0;
0,1,0,0;
0,0,1,0;
0,0,1,0;
0,0,1,0;
0,0,1,0;
0,0,0,1;
0,0,0,1;
0,0,0,1;
0,0,0,1;
];
[ht, wd] = size(points);
for i = 1:ht + 1
for j = 1:wd
if (i == 3)
H1(i, j) = 1;
else
H1(i, j) = points(i, j);
end
end
end
#disp(H1);
s1 = 3;
s2 = 4;
for i = 1:s2
for j = 1:s1
W(i,j) = 1;
end
end
#disp(W);
Z2 = W * H1;
#disp(Z2);
N = wd;
function HI = softmax(Z2)
[N, s2] = size(Z2);
for j = 1:N
sum = 0;
for k = 1:s2
sum = sum + exp(Z2(j, k));
end
for i = 1:s2
HI(j, i) = exp(Z2(j, i)) / sum;
end
end
endfunction;
H2 = softmax(Z2);
for i = 1:N
sum = 0;
for j = 1:s2
sum = sum + H2(j, i);
end
disp(sum);
end
#disp(H2);