testfeatkar.m

一种快速的多叉搜索树

clear;
X = load('mfeat-kar');
for i=1:10
Y((i-1)*200+1:i*200,:) = repmat(i,200,1);
end

train_size = 150;
test_size = 200-train_size;
%select training samples
for i=1:10
    XT((i-1)*train_size+1:i*train_size,:) = X((i-1)*200+1:(i-1)*200+train_size,:);
    YT((i-1)*train_size+1:i*train_size,:) = Y((i-1)*200+1:(i-1)*200+train_size,:);
end
%select test samples
for i=1:10
    XTT((i-1)*test_size+1:i*test_size,:) = X((i-1)*200+train_size+1:i*200,:);
    YTT((i-1)*test_size+1:i*test_size,:) = Y((i-1)*200+train_size+1:i*200,:);
end

N = size(XT,1);
nXdim = size(XT,2);
nYdim = size(YT,2);
q = 5;
knn = 1;
freezeW = 50;
x_quantize = 0;
y_quantize = 0;

IX = randperm(N); 

tid = createIHDRT(nXdim,nYdim,q,knn,freezeW,x_quantize,y_quantize);

for i=1:N 
    ret = addPatternIHDRT(tid, XT(IX(i),:), YT(IX(i),:));
end

disp('resubstitute test');
good = 0;
for i=1:N
    [ret, Xout, Yout, dist] = RetrieveIHDRT(tid, XT(i,:));
    if Yout==YT(i,:) 
        good=good+1;
    end
end
ROC = good/N

disp('cross-test');
good = 0;
for i=1:size(XTT,1)
    [ret, Xout, Yout, dist] = RetrieveIHDRT(tid, XTT(i,:));
    if Yout==YTT(i,:) 
        good=good+1;
    end
end
ROC = good/size(XTT,1)

closeIHDRT(tid);