function [probAintersectB probAnotB probBnotA sim]= hardCode2() % hardcoded example for simulating purposes % A ~ zipf(3,1), B~ 1-zipf(3,1) % similarity = counting (alpha = 1, beta = gamma = 0.5) f1 = zeros(64,3); f2 = zeros(64,3); sim = zeros(64,1); A= []; for i=1:64 bin = dec2bin(i,6); Features1 = bin(1:3); Features2 = bin(4:6); f1 = bin2dec(Features1); f2 = bin2dec(Features2); for j=1:3 if(Features1(j) == '1') features1(j) =1; elseif(Features1(j) == '0') features1(j)=0; end if(Features2(j) == '1') features2(j) =1; elseif(Features2(j) == '0') features2(j)=0; end end %disp('features1'); %disp(features1); %disp('features2'); %disp(features2); [AintersectB(i,:) AnotB(i,:) BnotA(i,:)] = featurescomp(features1, features2); %probAintersectB(i) = calcTrueProbIntersect(AintersectB(i,:), 'zipf', 0, ... % [3 1]); %[probAnotB(i)] = calcTrueProbExclusive('A', AnotB(i,:), 'zipf', 0, [3 1]); %[probBnotA(i)] = calcTrueProbExclusive('B', BnotA(i,:), 'zipf', 0,[3 1]); %disp('probAintersectB'); %disp(probAintersectB); %disp('probAnotB'); %disp(probAnotB); %disp('probBnotA'); %disp(probBnotA); %disp('sim'); %Sim = simLin(probAintersectB(i), probAnotB(i), probBnotA(i)); Sim = simCounting(AintersectB(i,:), AnotB(i,:), BnotA(i,:), 1, 0.5, 0.5); %disp(simLin(probAintersectB, probAnotB, probBnotA)); sim(i) = Sim; A = [A; f1 f2 Sim AintersectB(i,:) AnotB(i,:) BnotA(i,:)]; end success = xlswrite('zipf-counting-3', A); i = 1:1:64; %plot(i,probAintersectB,'go', i, probAnotB, 'ro', i, probBnotA, 'bo', i, sim, 'ko'); %subplot(4,1,1), plot(i,probAintersectB,'o'); %subplot(4,1,2), plot(i,probAnotB,'o'); %subplot(4,1,3), plot(i,probBnotA,'o'); %subplot(4,1,4), plot(sim,'o');