-
Notifications
You must be signed in to change notification settings - Fork 0
/
digRec.m
57 lines (46 loc) · 1.2 KB
/
digRec.m
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
%function dig = digRec(test, recdata)
%calculates the closest digit to test
%by solving min(||Ui' * b - Si * yi||)
%where yi = (Ui * Si)\b and b = format(test)
%
%input:
% test ... test image
% recdata ... path to folder containing matrices U and S
%
%output:
% dig ... recognised digit
function dig=digRec(test, recdata)
%format test image
b = format(test);
%assume the digit is 0
dig = 0;
%set U and S path
upath = [recdata "/U0.mat"];
spath = [recdata "/S0.mat"];
%import lernset
load(upath);
load(spath);
%solve for y1, U1*S1*y1 = b
y = (U * S)\b;
%get ||U1'*b - S1*y1|| and set as minimum value
min = norm(U' * b - S * y);
for i=1:9
%set U and S path
upath = [recdata "/U" num2str(i) ".mat"];
spath = [recdata "/S" num2str(i) ".mat"];
%import lernset
load(upath);
load(spath);
%solve for yi, Ui*Si*yi = b
y = (U * S)\b;
%get ||U1'*b - S1*y1|| and set as compare value
tmp = norm(U' * b - S * y);
%if distance between b and Ai*xi smaller than min
if tmp < min
%set as min
min = tmp;
%set closest digit to digit represented by b
dig = i;
endif
endfor
endfunction