C'est tout simple :


[ligne,colonne]=find(matrice > (mean(matrice(:) + 3*std(matrice(:)) ) ) )