% Математическая модель по определению fдоп по фазовому спектру 
 
format long 
clear all;
% Исходные данные:
j=sqrt(-1);    % Мнимая единица
F_prom=63;     % Промежуточная частота (МГц)
fs=28;         % Частота дискретизации (МГц)
dt=1/fs;       % Шаг дискретизации (мкс)
F_rep=0.2;     % Частота повторения импульса (МГц)
T_rep=1/F_rep; % Период повторения импульса (мкс)
tau_imp=1;     % Длительность импульса (мкс)
N_imp=input('Количество импульсов в пачке: ');                 % Количество импульсов в пачке
F_Doppler=input('Доплеровский сдвиг частоты (кГц): ')*10^(-3); % Доплеровский сдвиг частоты (кГц)
tau_delay=input('Задержка сигнала (мкс): ');                   % Задержка сигнала (мкс)                                                 
q_dB=input('q_дБ: ');                                          % Отношение сигнал-шум в дБ
K=input('Число реализаций: ');                                 % Число реализаций
     
q=10^(q_dB/20);
phase=2*pi*rand;                  % Случайная начальная фаза сигнала
n_delay=round(tau_delay*fs);      % Количество отсчетов в задержке строба
N_sample_1=T_rep*fs;              % Количество отсчетов в одном периоде повторения
N_sample_imp=tau_imp*fs;          % Количество отсчетов в одном импульсе
N_sample_signal=N_sample_1*N_imp; % Количество отсчетов в сигнале
 
h=waitbar(0,'Please wait...');
for k=1:K
waitbar(k/K)
figure(1)
% Формирование когерентной пачки радиоимпульсов
t=(0:(N_sample_signal+n_delay)-1)*dt;           % Интервал времени сигнала
x=cos(2*pi*(F_prom+F_Doppler)*t+phase);         % Формирование гармонического сигнала
U_1=zeros(1,N_sample_1); U_1(1:N_sample_imp)=1; % Формирование одного импульса
U_imp=ones(1,N_imp);                            % Массив содержащий, столько единиц, сколько импульсов в пачке
U=kron(U_imp,U_1);                              % Последовательность импульсов
U(end+n_delay)=0;
U=x.*U;                                         % Когерентная пачка радиоимпульсов
U=circshift(U.',n_delay).';
plot(t,U)
xlabel('t (мкс)')
ylabel('Амлитуда')
title ('Когерентная пачка радиоимпульсов')
grid on
 
figure(2)
% Формирование аддитивной смеси сигнала и шума
sigma_signal=std(U);                   % СКО сигнала
subplot(211)
sigma_noise=sigma_signal/q;            % СКО шума
noise=sigma_noise*randn(1,length(t));  % Формирование шума
plot(t,noise)                          % График шума
xlabel('t (мкс)')
title ('Шум')
grid on
subplot(212)
signal_with_noise=U+noise;             % Аддитивная смесь сигнала и шума
plot(t,signal_with_noise)              % График аддитивной смеси сигнала и шума
xlabel('t (мкс)')
ylabel('Амлитуда')
title ('Аддитивная смесь сигнала и шума')
grid on
 
FileName1=fopen('signal.txt','w');              % Открытие текстового файла statistics. Флаг "a" добавляет значения в текстовый файл
fprintf(FileName1,'%f\n',signal_with_noise);    % Запись данных в текстовый файл statistics
fclose(FileName1);

% Определение доплеровского сдвига частоты по фазовому спектру
clear all;
[FileName,PathName] = uigetfile('*.txt','Выбор файла'); % Выбор txt-файла
if FileName==0 % Обработка нажатия клавиши Cancel в окне выбора файла
    return
end
 
figure(1)
% Аддитивная смесь сигнала и шума
real_signal=textread([PathName,FileName]);