-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathprocessOvO.m
109 lines (100 loc) · 2.67 KB
/
processOvO.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
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
102
103
104
105
106
107
108
function processOvO()
%PROCESSOVM Summary of this function goes here
% Detailed explanation goes here
close all; clear all; clc;
load('OvO.mat');
dA = [0:9];
iA = [1:5]; %1 is done in 10s, so this is 15x multiplier
kA = [0.000001 0.000003 0.00001 0.00003 0.0001 0.0003 0.001 0.003 0.01 0.03];
sA = [100 300 1000]; %100 is done in 10s, so this is x+3x+10x = 14x multiplier
pO = zeros(4,size(dA,2)*size(iA,2)*size(kA,2)*size(sA,2));
p = zeros(4,size(dA,2)*size(iA,2)*size(kA,2)*size(sA,2));
STr1 = STr;
STe1 = STe;
HTr1 = HTr;
HTe1 = HTe;
i = 1;
for Q=1:size(iA,2)
for W=1:size(dA,2)
for E=1:size(kA,2)
for R=1:size(sA,2)
p(1:1,i) = dA(W);
p(2:2,i) = iA(Q);
p(3:3,i) = kA(E);
p(4:4,i) = sA(R);
i = i+1;
end
end
end
end
i = 1;
for Q=1:size(dA,2)
for W=1:size(iA,2)
for E=1:size(kA,2)
for R=1:size(sA,2)
pO(1:1,i) = dA(Q);
pO(2:2,i) = iA(W);
pO(3:3,i) = kA(E);
pO(4:4,i) = sA(R);
i = i+1;
end
end
end
end
for I=1:size(p,2)
d = find(p(1,:) == pO(1,I));
i = find(p(2,:) == pO(2,I));
k = find(p(3,:) == pO(3,I));
s = find(p(4,:) == pO(4,I));
nI = intersect(intersect(intersect(d,i),k),s);
STr1(:,I) = STr(:,nI);
STe1(:,I) = STe(:,nI);
HTr1(:,I) = HTr(:,nI);
HTe1(:,I) = HTe(:,nI);
end
% processDataTable(reorg(STr1),p);
processDataTable(reorg(STe1),p,'SVM');
% processDataTable(reorg(HTr1),p);
processDataTable(reorg(HTe1),p,'Hebbian');
% max(STe(5,:))
% max(STe(8,:))
% find(STe(5,:) == max(STe(5,:)))
% find(STe(5,:) == max(STe(5,751:900)))
% find(STe(8,:) == max(STe(8,751:900)))
% help reshape
% S = reshape(STe,150,80)
% pR = reshape(p,150,40)
% help reshape
% pR = reshape(p.',150,[])
% pR = reshape(p.',150,[])'
% S = reshape(STe.',150,[])'
% max(STe(41:48,:))
% max(mean(STe(41:48,:)))
% mean(S(41:48,:))
% max(mean(S(41:48,:)))
% find(mean(S(41:48,:) == max(mean(S(41:48,:)))))
% OvM.mat
% find(mean(S(41:48,:)) == max(mean(S(41:48,:))))
% find(mean(S(71:78,:)) == max(mean(S(71:78,:))))
% max(mean(S(41:48,:)))
% max(mean(S(41:48,140:150)))
% find(mean(S(41:48,:) == max(mean(S(41:48,:)))))
% mean(S(41:48,140:150))
% mean(S(41:47,140:150))
% max(mean(S(41:47,140:150)))
% max(mean(S(41:48,:)))
% max(mean(S(71:77,:)))
% find(S == max(mean(S(71:77,:))))
% find(S(71:77,:) == max(mean(S(71:77,:))))
% find(mean(S(71:77,:) )== max(mean(S(71:77,:))))
% find(pR(31,:)==300)
% t = find(pR(31,:)==300)
% s(t)
% S(t)
% S(t,:)
% S(:,t)
% St = s(:,t)
% St = S(:,t)
% max(mean(St(41:48,:)))
% find(mean(St(41:48,:) )== max(mean(St(41:48,:))))
end