function [yNoisy,ps] = gaussianSpectrogram(t,fs)
% Generate noisy data and the power spectrum for a Gaussian in the time
% domain.

% Set center and width of the Gaussian, in the time domain.
baselineMu = 3;
baselineSigma = 1;

% Set how much those parameters can randomly vary.
muVariation = 0.4;
sigmaVariation = 1;

% Gaussian time series.
mu = baselineMu + muVariation*randn;
sigma = baselineSigma + sigmaVariation*randn;

y = exp(-(t-mu).^2./sigma^2);

y = addBackground(y,t);
yNoisy = addNoise(y);

ps = computePowerSpectrum(yNoisy,fs);
end