Sparse and low-rank approximations of large symmetric matrices using biharmonic interpolation