I am using the plomb function to estimate the PSD of an unevenly spaced velocity time series. The results look reasonable and are comparable to a pwelsh PSD; however, I would like to improve the results using a multitaper method and reduce some of the variability. Any recommendations on how to combine multitaper with the Matlab plomb function? Or any other suggestions for improving the plomb results in general? Thanks.