T=10e-3; % Длительность сигнала, 10 мс
S0=5; % Амплитуда сигнала, 5 В
f0=1e3; % Центральная частота, 1 кГц
td=1/(20*f0); % Шаг дискретизации

% Lt=round(T/td); % Длина сигнала
% t=zeros(1,Lt); % инициализация вектора времени
% for i=1:Lt % формирование вектора времени в цикле
%     t(i)=(i-1)*td;
% end

t=(0:td:2*T-td); % формирование вектора времени, короткая реализация цикла
Lt=length(t); % Длина сигнала
S=zeros(1,Lt);
S(1:Lt/2)=S0; % Видеоимпульс с прямоугольной огибающей
s=S.*cos(2*pi*f0*t); % Радиоимпульс с прямоугольной огибающей, 
% .* - поэлементное перемножение двух векторов

Sf=fft(S); % Спектр видеоимпульса
sf=fft(s); % Спектр радио импульса

sigma=20; % СКО шума
mo=0; % МО шума
n=sigma*randn(1,Lt)+mo; % шум

u=s+n; % смесь сигнала и шума
u_cor=xcorr(u,s); % корреляционная обработка смеси
s_cor=xcorr(s,s); % обработка сигнала
n_cor=xcorr(n,s); % обработка шума
Lcor=length(u_cor);

% построение графиков
figure(1),
subplot(3,1,1), plot(t,s), grid on
subplot(3,1,2), plot(t,[s;u]), grid on
subplot(3,1,3), plot(1:Lcor,[u_cor; s_cor; n_cor]), grid on

figure(2),
subplot(2,2,1), plot(abs(Sf)), grid on
subplot(2,2,2), plot(fftshift(abs(Sf))), grid on
subplot(2,2,3), plot(abs(sf)), grid on
subplot(2,2,4), plot(fftshift(abs(sf))), grid on