Fit gaussian process to time series data

I have time series data for a signal for a particular time period t. The signal profile is different in each data representing the uncertainty. I want to fit a Gaussian process to this data so that I can sample a time series from this Gaussian process. How can I do that in Matlab?