PatternRecognitionSamples/assignmentII_2.m at master · SevdanurGENC/PatternRecognitionSamples · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
clear all;
close all;
clc;

W1 = [ 8,4; 6,4; 7,5; 7,3; ]
W2 = [ 10,-2; 4,-2; 7,-1; 7,-3;]

%% A Sikki %%

nokta_w1 = [ 8,6,7,7; 4,4,5,3; ];
figure;
plot(nokta_w1(1,:),nokta_w1(2,:),'*r');
nokta_w2 = [ 10,4,7,7; -2,-2,-1,-3;]
hold on;
plot(nokta_w2(1,:),nokta_w2(2,:),'*b');
axis([-15 15 -15 15 -15 15]);
grid on;

%% B Sikki %%

ortalamaW1 = mean(W1)
ortalamaW2 = mean(W2)

kovaryansW1 = cov(W1)
kovaryansW2 = cov(W2)

%% C Sikki %%

E1 = [1/2 , 0; 0, 1/2];
E2 = [4.5, 0; 0 0.5];

Wb1 = - (1/2) * E1';
Wb2 = - (1/2) * E2';

wk1 = E1' * ortalamaW1(1,1);
wk2 = E2' * ortalamaW1(1,1);

wk10 = [- (1/2) * ortalamaW1(1,1)' * E1' * ortalamaW1(1,1)] - [ (1/2) * (log(det(E1)) + log(exp(1))) ];
wk20 = [- (1/2) * ortalamaW2(1,1)' * E2' * ortalamaW2(1,1)] - [ (1/2) * (log(det(E2)) + log(exp(1))) ];

for i=1:length(W1)
    G11 = W1(i,1)' * Wb1 * W1(i,1) + wk1' * W1(i,1) + wk10;
    G12 = W1(i,2)' * Wb1 * W1(i,2) + wk1' * W1(i,2) + wk10;
end

for i=1:length(W1)
    G21 = W2(i,1)' * Wb2 * W2(i,1) + wk2' * W2(i,1) + wk20;
    G22 = W2(i,2)' * Wb2 * W2(i,2) + wk2' * W2(i,2) + wk20;
end

%% D Sikki %%


%% E Sikki %%

prototype = [9,1]

for i=1:length(W1)
sonuc(i,:) = [(prototype(1,1) - W1(i,1))^2  + (prototype(1,2) - W1(i,2))^2];
end

NNMatrixW1 = sonuc(:,:)
vektorAraligi1 = min(NNMatrixW1)

for i=1:length(W2)
sonuc(i,:) = [(prototype(1,1) - W2(i,1))^2  + (prototype(1,2) - W2(i,2))^2];
end

NNMatrixW2 = sonuc(:,:)
vektorAraligi2 = min(NNMatrixW2)

prototype = [0,1]      %% 2.satir'a ait.

for i=1:length(W1)
sonuc(i,:) = [(prototype(1,1) - W1(i,1))^2  + (prototype(1,2) - W1(i,2))^2];
end

NNMatrixW1 = sonuc(:,:)
vektorAraligi1 = min(NNMatrixW1)

for i=1:length(W2)
sonuc(i,:) = [(prototype(1,1) - W2(i,1))^2  + (prototype(1,2) - W2(i,2))^2];
end

NNMatrixW2 = sonuc(:,:)
vektorAraligi2 = min(NNMatrixW2)

hold on;
prototypeNokta1 = [9; 1];
plot(prototypeNokta1(1,:),prototypeNokta1(2,:),'Or');

hold on;
prototypeNokta2 = [0; 1];
plot(prototypeNokta2(1,:),prototypeNokta2(2,:),'Ob');