计算余弦距离

load('lfw_feature_jing.mat');
fid = fopen('lfw_all.txt','wt');
a=[];
label=[];
for i=1:2:12000
    feature1=feats(i,:);
    feature2=feats(i+1,:);
    x=[feature1;feature2];
    distance_cosine=1-pdist(x,'cosine');
    a=[a distance_cosine];
    [cos,index]=sort(a);
    new_label=[cos,index]';
    fprintf(fid,'%g\n',distance_cosine);
end

fclose(fid);

load('mydata.mat');
acc1=[];
for i=1:6000
    b=a(1,i);
    m=0;
    n=0;
    for j=1:3000
        if b>=a(1,j)
            m=m+1;
        end
    end
    for j=3001:6000
        if b<=a(1,j)
            n=n+1;
        end
    end
    sum=m+n;
    b;
    acc=sum/6000;
    acc1=[acc1 acc];
    [t,p]=sort(acc1);
end