active_rgn.m

来自「software to find the threshold of an LDP」· M 代码 · 共 10 行

function [A_r,map_r,K_r] = active_rgn(fp,A,delta,lpad,K_rmax)

foo = sum(cumsum(fp) < (sum(fp)-1e-6));
K_r = foo;

map_r = 1:K_r;

A_r = A(1:(foo+2*lpad+1),map_r);