ЛР 3 вариант 14 Код Matlab LR3

var = 14;

f_0 = 180 + var; % несущая частота
delta_f = 25 + var; %полоса сигнала
f_s = 800 + var; %частота дискретизации
N_i = 400 + var; % дляительность сигнала в отсчётах
f_dp = 0; % частота Доплера
phi_0 = 2*pi/var; % начальная фаза
N = 1024;

var_2 = mod(14, 4); % остаток от деления = 2, => ЛЧМ pass сигнал по центру

T = 1/fs;
ti = N_i*T;

sg = @(t) cos(2*pi*(f_0*t+1/2*delta_f/ti*t.^2+phi_0/(pi*2))); % сигнал
signal = [sg(0:T:ti/2), zeros(1,N-N_i+1), sg(-ti/2:T:-T)];
Re = real(signal); % действительный спектр
Im = imag(signal); % мнимый спектр

imp_resp = @(t) cos(2*pi*(f_0*t-1/2*delta_f/ti*t.^2)); % импульсная характеристика
impulse_resp = [imp_resp(0:T:ti/2), zeros(1,N-N_i+1), imp_resp(-ti/2:T:-T)];
Re_imp = real(impulse_resp); % действительный спектр импульсной характеристики
Im_imp = imag(impulse_resp); % мнимый спектр импульсной характеристики

dpf = abs(fft(signal)); % амплитудный спектр сигнала
dpf_imp = abs(fft(impulse_resp)); % амплитудный спектр импульсной характеристики

phase_spec = angle(fft(signal)); % фазовый спектр сигнала
phase_spec_imp = angle(fft(impulse_resp)); % фазовый спектр импульсной характеристики

my_plot_LR3 ('Действительный спектр сигнала', Re, 'Мнимый спектр сигнала', Im, 'Действительный спектр импульсной характеристики', Re_imp, 'Мнимый спектр импульсной характеристики', Im_imp, 'Амплитудный спектр сигнала', dpf, 'Фазовый спектр сигнала', phase_spec, 'Амплитудный спектр импульсной характеристики', dpf_imp, 'Фазовый спектр импульсной характеристики', phase_spec_imp)

% спектры произведения
multiply = fft(signal).*fft(impulse_resp); % произведения БПФ сигнала и БПФ импульсной характеристики
Re_mult = real(multiply); % действительный спектр произведения
Im_mult = imag(multiply); % мнимый спектр произведения
dpf_mult = abs(multiply); % амплитудный спектр произведения
phase_mult = angle(multiply); % фазовый спектр произведения

% сжатый сигнал
OBPF = ifft(multiply); % ОБПФ произведения БПФ сигнала и БПФ импульсной характеристики
Re_OBPF = real(OBPF); % действительный спектр сжатого сигнала
Im_OBPF = imag(OBPF); % мнимый спектр сжатого сигнала
dpf_OBPF = abs(OBPF); % модуль сжатого сигнала

my_plot_LR3 ('Действительный спектр сжатого сигнала', Re_OBPF, 'Мнимый спектр сжатого сигнала', Im_OBPF, 'Модуль сжатого сигнала', dpf_OBPF, 'No picture', 0, 'Действительный спектр произведения', Re_mult, 'Мнимый спектр произведения', Im_mult, 'Амплитудный спектр произведения', dpf_mult, 'Фазовый спектр произведения', phase_mult)