Real Gaussian Windowed Chirp (Time Domain)

fs = 8000;
x = chirp([0:1/fs:0.1],1000,1,2000);
M = length(x);
n=(-(M-1)/2:(M-1)/2)';
w = exp(-n.*n./(2*sigma.*sigma));
xw = w(:) .* x(:);