Исходные данные:
Амплитуда 1 В
Количество импульсов 2
Задержка 0 мкс
ско 0 (шума нет)
Доплеровская частота 25
Подскажите как мне тут быть?
Код: Выделить всё
% Исходные данные
U_m=input('Амплитуда сигнала (В): '); % Амплитуда сигнала
F_prom=63; % Промежуточная частота (МГц)
T_rep=10; % Период повторения импульсов (мкс)
F_rep=1/T_rep; % Частота повторения (МГц)
tau_imp=1; % Длительность импульса (мкс)
N_imp=input('Количество импульсов в пачке: '); % Количество импульсов в пачке
T_pack=N_imp*T_rep; % Длительность пачки радиоимпульсов (мкс)
F_Doppler=input('Доплеровский сдвиг частоты (кГц): ')*10^(-3); % Доплеровский сдвиг частоты (МГц)
fs=28; % Частота дискретизации (МГц)
dt=1/fs; % Шаг дискретизации (мкс)
tay_delay=input('Задержка отраженного сигнала (мкс): '); % Задержка отраженного сигнала (мкс)
N_sample_1=T_rep*fs; % Количество отсчетов в одном периоде
N_sample_imp=tau_imp*fs; % Количество отсчетов в импульсе
sigma_noise=input('СКО шума: ');
figure(1)
% Отр. от цели сигнал. Сигнал представляет собой ког. пачку импульсов.
t=0:0.001:T_pack; % Интервал времени пачки импульсов
d=[0.5:10:T_pack]; % Массив значений центров импульсов
y=U_m*pulstran(t-tay_delay,d,'rectpuls',1); % Формирование последовательности прямоугольных импульсов
b=cos(2*pi*(F_prom+F_Doppler)*t); % Формирование гармонического сигнала т.е. ВЧ заполнение в импульсе
x=y.*b; % Получение радиоимпульса путем умножения огибающей сигнала и гармонического сигнала
plot(t,x),grid % График принятого отраженного сигнала
ylim([-U_m-0.5,U_m+0.5]) % Интервал значений по оси y
figure(2)
% Аналоговый полосовой фильтр f0=63 МГц, полоса пропускания 14 МГц.
Center_frequency=63; % Центральная частота аналогового полосового фильтра
[b,a]=butter(5,[2*pi*(Center_frequency-7) 2*pi*(Center_frequency+7)],'s'); % Синтез коэффициентов фильтра
f=0:0.1:90; % Массив значений частоты (МГц)
h=freqs(b,a,2*pi*f); % Комплексный коэффициент передачи
plot(f,abs(h)) % График АЧХ полосового фильтра
xlim([30 90]) % Интервал значений по оси x
H=tf(b,a); % Передаточная функция соответствующая своим коэффициентам числителя и знаменателя
% Формирование шума. Нормальный случайный процесс.
T=0:dt:T_pack-dt; % Массив дискретного времени (мкс)
noise=sigma_noise*randn(1,length(T)); % Формирование шума
noise_f=lsim(H,noise,T); % Шум на выходе аналогового фильтра
noise_exit=noise_f'; % Изменение вектора столбца на вектор строку
figure(3)
% Оцифрованный (отраженный) сигнал. Сигнал на выходе АЦП
y_digitize=U_m*pulstran(T-tay_delay,d,'rectpuls',1);
b_digitize=cos(2*pi*(F_prom+F_Doppler)*T);
x_digitize=y_digitize.*b_digitize+noise_exit; % Сигнал-шум на выходе АЦП
plot(T,x_digitize),grid % График сигнал-шум на выходе АЦП
xlim([0 1])
figure(4)
% Вычисление спектра от дискретного (отраженного) сигнала на выходе АЦП.
NFFT=N_sample_1*N_imp; % Количество дискретных отсчетов сигнала
df=fs/NFFT; % Шаг по частоте в спектре (МГц)
f=(0:NFFT-1)*df; % Диапозон частот в спектре (МГц)
Spectr=abs(fft(x_digitize,NFFT)); % Вычисление спектра
plot(f,Spectr),grid % График дискретного спектра
xlim([0,fs]) % Интервал значений по оси x
figure(5)
% Цифровой полосовой фильтр f0=21 МГц, полоса пропускания 10 МГц.
Center_frequency=21; % Центральная частота цифрового полосового фильтра
[b1,a1]=butter(6,[2*pi*(Center_frequency-5) 2*pi*(Center_frequency+5)],'s'); % Синтез коэффициентов фильтра
f=0:0.1:40; % Массив значений частоты (МГц)
h1=freqs(b1,a1,2*pi*f); % Комплексный коэффициент передачи цифрового фильтра
plot(f,abs(h1)) % График АЧХ полосового фильтра
xlim([0 40]) % Интервал значений по оси x
H1=tf(b1,a1); % Передаточная функция соответствующая своим коэффициентам числителя и знаменателя
signal_noise=lsim(H1,x_digitize,T); % Сигнал-шум на выходе цифрового фильтра
signal_noise1=signal_noise'; % Изменение вектора столбца на вектор строку
figure(6)
plot(T,signal_noise1),grid % График сигнал-шум на выходе АЦП
figure(7)
% Вычисление спектра от дискретного (отраженного) сигнала на выходе цифрового фильтра.
f1=(0:NFFT-1)*df; % Диапозон частот в спектре (МГц)
Spectr1=abs(fft(signal_noise1,NFFT)); % Вычисление спектра
plot(f1,Spectr1),grid % График дискретного спектра
xlim([0,fs]) % Интервал значений по оси x