Commit 0239a660 authored by Riccardo La Placa's avatar Riccardo La Placa
Browse files

Changed det. stat. threshold by adding mask on ni^s combinations considered

parent ee865032

scsearch.m

+13 −6
Original line number Diff line number Diff line
@@ -195,7 +195,7 @@ tic 


% Try s* and check \nu_s range

g_jj=(((pi*Tseg)^2)/3).*[1; (Tseg^2)/60; (Tseg^4)/1344; (Tseg^6)/172800]; %eq. 22 M2015 + calcoli da eq. 21 M2015

mu_s=0.025; %massimo mismatch sulla griglia coerente da scegliere

mu_s=0.05; %massimo mismatch sulla griglia coerente da scegliere

%s_s=uint8(4);

s_s = 4;

while(1)
@@ -247,9 +247,9 @@ Lambda = zeros(length(nibank),length(f_gr),'gpuArray'); 
F1=gpuArray((0:N-1)./(N*dt_psd)).';

f_ind = zeros(1,length(f_gr),'uint32');

for n = 1:length(f_gr)

    [bles,f_ind(n)] = min(abs(F1-f_gr(n)));

    [~,f_ind(n)] = min(abs(F1-f_gr(n)));

end

clear bles

%clear bles

% nibank = combinations(nis{:}).Variables;

% Lambda = zeros(length(nibank),length(f_gr));

toc
@@ -341,9 +341,11 @@ nisearcher = KDTreeSearcher(gather(nibank),'BucketSize',100); 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

parbank = combinations(porb_gr,a_gr,tasc_gr).Variables;

bestpar = zeros(1,5);

%bestpar = zeros(1,5);

toc





nimask = false(length(nibank));

totlam = zeros(length(parbank),length(f_gr));

tic

% for n=1:length(f_gr)
@@ -397,7 +399,9 @@ for m = 1:M 
        curni = -curni.*curpar(3);

        curni(1) = curni(1) + 1;

        [Idx,D] = knnsearch(nisearcher,curni);



        if (~nimask(Idx))

            nimask(Idx) = 1;

        end

        totlam(i,:) = totlam(i,:)+Lambda(Idx,:);



    end
@@ -439,7 +443,8 @@ toc 


%% Calculate detection threshold for a multi-trial false alarm prob. = pfa

pfa = 0.01;

sigmastar = 2*gammaincinv((1-(1-pfa)^(1/(length(parbank)*length(f_gr)))),M,'upper');

%sigmastar = 2*gammaincinv((1-(1-pfa)^(1/(length(parbank)*length(f_gr)))),M,'upper');

sigmastar = 2*gammaincinv((1-(1-pfa)^(1/(sum(nimask)*length(f_gr)))),M,'upper');

[bru,tto] = max(totlam);

[maxtotlam,maxinde] = max(bru);

bestpar = [f_gr(maxinde), parbank(tto(maxinde),:), maxtotlam];
@@ -455,6 +460,8 @@ save([pathfi,'risultelli_',char(datetime('now','Format','dd_MM')),'_',num2str(Ts 
%% 





daterella = ['Date and time now ',char(datetime('now')),' (UTC), maybe'];

disp(daterella); clear daterella

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%% Fai il grafico delle totlam a posteriori %%%%%%%%%%%%%%%%