Add power spectrum matlab code