Systems and methods manage electrical loads in a grid by applying Robust principal component analysis (R-PCA) to decompose annual load profiles into low-rank components and sparse components; extracting one or more predetermined features; constructing a similarity graph; selecting submodular cluster centers through the constructed similarity graph; determining a cluster assignment based on selected centers; and applying the clustering assignment for load analysis.