-
Notifications
You must be signed in to change notification settings - Fork 0
/
compute_Euclidean_bhattacharyya.m
35 lines (24 loc) · 1.24 KB
/
compute_Euclidean_bhattacharyya.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
function [ index_Euclidean_bhattacharyya,extend_distance_bhattacharyya ] = compute_Euclidean_bhattacharyya(extend_Euclidean_helix,extend_Euclidean_stick,helix_number,stick_number )
% for i=1:helix_number
% sort_Euclidean_helix(i,:)=sort(extend_Euclidean_helix(i,:));
% sort_Euclidean_stick(i,:)=sort(extend_Euclidean_stick(i,:));
%
% end
for i=2:helix_number
p= extend_Euclidean_helix(i,2:helix_number);
for j=2:stick_number
q= extend_Euclidean_stick(j,2:helix_number);
extend_distance_bhattacharyya(i,j)=bhattacharyya(p,q); % compute bhattacharyya distance
end
end
% insert index of row and column for computing correct index after removing
extend_distance_bhattacharyya(2:helix_number,1)=extend_Euclidean_helix(2:helix_number,1);
extend_distance_bhattacharyya(1,2:helix_number)=extend_Euclidean_stick(1,2:helix_number);
for i=2:helix_number
d=extend_distance_bhattacharyya(i,2:helix_number);
[minimum,index_bat] = min(d(:));
index_Euclidean_bhattacharyya(i-1).num_helix=extend_distance_bhattacharyya(i,1);
index_Euclidean_bhattacharyya(i-1).num_stick=extend_distance_bhattacharyya(1,index_bat+1);
index_Euclidean_bhattacharyya(i-1).min=minimum;
end
end