% Doppler shift
        path = A(j)*exp(1i*2*pi*t*Doppler_shift(j)).*WhiteNoise';

        Sig = path.*SigOut;
        SignalOutChan = SignalOutChan + Sig;

%% ***** Инициализация начальных параметров сигнала ****
       Fs=20; %<--- Частота дискретизации;
       f=5;    %<--- Частота сигнала;
       t=0:1/Fs:5; %<--- Время;
%% **************************************************
%%
%% *********** Формирование сигнала ******************
       y=sin(2*pi*t*f); %<--- Исходный сигнал;
%% **************************************************
%%
%%
%% ******* Перевод в частотную область *************
       FF=fft(y);
%% **************************************************
%%
%%
%% ************ Простой сдвиг спектра на 10 шагов ********
       FF2=circshift(FF,[0 10]);
%% **************************************************
%%
%%
%% ******** Перевод во временную область ***************
       y2=ifft(FF2);
%% **************************************************
%%
%%
%% ********* Построение графиков **********************
%   1. График исходного сигнала:    
%       1.1. Построение графика:
                subplot(1,2,1);
                plot(t,y);
%       1.2. Настройка графика:
                title('Исходный сигнал');
%   2. График принятого сигнала:    
%       2.1. Построение графика:
                subplot(1,2,2);
                plot(t,y2);
%       1.2. Настройка графика:
                title('Принятый сигнал после сдвига спектра');
%% ***************************************************

       if(Delay(j)>0)
            k_integ = fix(Delay(j)/ts);
            SigOut = [zeros(1,k_integ - 1),SigIn(1:end-k_integ+1)]*A(j);
        else
            SigOut = SigIn*A(j);
        end

coeffs = c./(real(c).*real(c)+imag(c).*imag(c));

coeffs = c./(sqrt(sum(c.*c)));

        WhiteNoise =  filter(coeffs,1,WhiteNoise);

path = A(j)*exp(1i*2*pi*t*Doppler_shift(j)).*WhiteNoise';

clear;

Fs = 500;
Doppler_spread = 1;
L = 128;

win = hamming(2*L);
c = zeros(1,2*L);
coeff = -4.0*(pi*Doppler_spread/Fs)^2;

for  i=1:L
                tmp = win(L+i) * exp(coeff*(i+0.5)^2);
                c(L+i) = tmp;
                c(L-(i-1)) = tmp;    
end

coeffs = c./norm(c);

plot(coeffs);