FisherAnalysis/libSDR_SCSF.m at master · ShinonomeLaboratory/FisherAnalysis · 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
function [w] = libSDR_SCSF(label_y,data_x,beta)
%去看libSDR函数的说明
size_y=size(label_y);
size_x=size(data_x);
if size_y(1)==1
    size_y=size_y';
    data_x=data_x';
end
size_y=size(label_y);
size_x=size(data_x);

label_1=label_y(1);
lst_c1=find(label_y==label_1);
lst_c2=find(label_y~=label_1);

if length(lst_c1)<=2 || length(lst_c2)<=2
    disp('某一类的样本数量过少...');
end

m1=sum(data_x(lst_c1,:))./length(lst_c1);
m2=sum(data_x(lst_c2,:))./length(lst_c2);

C1=data_x(lst_c1,:);
C2=data_x(lst_c2,:);

D1=C1-repmat(m1,length(lst_c1),1);
D2=C2-repmat(m2,length(lst_c2),1);

Sw=D1'*D1+D2'*D2;
for i=1:100
    if rank(Sw)<size_x(2)
        Sw=Sw+eye(size_x(2)).*(beta/100);
    else
        break;
    end
end
disp(['辅助矩阵系数：',num2str(beta*(i-1))])
w=Sw\(m2-m1)';
w=w./sqrt(sum(w.^2));
w=w(:);

end