-
Notifications
You must be signed in to change notification settings - Fork 7
/
Copy pathnormalize_eigenvalues.m
92 lines (80 loc) · 2.08 KB
/
normalize_eigenvalues.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
function [eigen_normalized,normalized_eigen_matrix, data_dimention]=normalize_eigenvalues(data,D)
dim=D;
eigen_normalized=data;
[Lamdas_gap2,data_dimention]=estimate_token(eigen_normalized,dim);
N=size(data,2);
s=zeros(1,N);
normalized_eigen_matrix=zeros(dim,size(data,2));
for i=1:dim
%normalized_eigen_matrix(i,:)=(Lamdas_gap2(i,:));
normalized_eigen_matrix(i,:)=(data(dim+1+i,:));
end
%normalized_eigen_matrix(1,:)=(Lamdas_gap2(1,:))+(Lamdas_gap2(2,:));
for i=1:N
s(i)=sqrt(sum(normalized_eigen_matrix(:,i)));
end
for i=1:N
normalized_eigen_matrix(:,i)=normalized_eigen_matrix(:,i)./s(i);
end
for i=1:N
eigen_normalized(dim+2:(2*dim)+1,i)=normalized_eigen_matrix(:,i);
end
end
% dim=D;
%
% if dim==2
% eigen_normalized=data;
% [Lamdas_gap2]=estimate_token(eigen_normalized,dim);
% %l12=(Lamdas_gap2(1,:));
% l1=(data(dim+2,:));
% l2=(data(dim+3,:));
% %l1 = l12 + l2;
% %s = sqrt(sum([l1', l2'].^2,2))';
% s = sqrt((l1+l2).^2);
% l2 = l2 ./ s;
% l1=l1./s;
% eigen_normalized(dim+2,:)=l1;
% eigen_normalized(dim+3,:)=l2;
%
%
% elseif dim==3
% eigen_normalized=data;
% [Lamdas_gap2,data_dimention]=estimate_token(eigen_normalized,dim);
% l12=(Lamdas_gap2(1,:));
% l2=(Lamdas_gap2(2,:));
% l3=(Lamdas_gap2(3,:));
% l1 = l12 + l2;
% s = sqrt(sum([l1', l2', l3'].^2,2))';
% l3 = l3 ./ s;
% l2 = l2 ./ s;
% l1=l1/s;
% eigen_normalized(dim+2,:)=l1;
% eigen_normalized(dim+3,:)=l2;
% eigen_normalized(dim+4,:)=l3;
% elseif dim>3
% eigen_normalized=data;
% [Lamdas_gap2]=estimate_token(eigen_normalized,dim);
% N=size(data,2);
% s=zeros(1,N);
% normalized_eigen_matrix=zeros(dim,size(data,2));
%
% for i=1:dim
% %normalized_eigen_matrix(i,:)=(Lamdas_gap2(i,:));
% normalized_eigen_matrix(i,:)=(data(dim+1+i,:));
%
% end
% %normalized_eigen_matrix(1,:)=(Lamdas_gap2(1,:))+(Lamdas_gap2(2,:));
%
%
% for i=1:N
% s(i)=sqrt(sum(normalized_eigen_matrix(:,i)));
% end
% for i=1:N
% normalized_eigen_matrix(:,i)=normalized_eigen_matrix(:,i)./s(i);
% end
% for i=1:N
% eigen_normalized(dim+2:(2*dim)+1,i)=normalized_eigen_matrix(:,i);
% end
% end
%
%