本次实验的目的是仿真发射信号、回波信号、中频信号以及对中频信号做range_fft
重要公式:
雷达发射信号:
S
T
(
t
)
=
A
T
e
x
p
(
j
2
π
(
(
f
c
−
B
W
2
)
t
+
B
W
2
T
t
2
)
)
S_T(t)=A_Texp(j2\pi ((f_c-\frac{BW}{2})t+\frac{BW}{2T}t^2))
ST(t)=ATexp(j2π((fc−2BW)t+2TBWt2))
雷达接收信号:
S
R
(
t
)
=
A
R
e
x
p
(
j
2
π
(
(
f
c
−
B
W
2
)
(
t
−
τ
)
+
B
W
2
T
(
t
−
τ
)
2
)
)
S_R(t)=A_Rexp(j2\pi ((f_c-\frac{BW}{2})(t-\tau)+\frac{BW}{2T}(t-\tau)^2))
SR(t)=ARexp(j2π((fc−2BW)(t−τ)+2TBW(t−τ)2))
雷达接收信号延迟:
τ
=
2
∗
d
c
\tau = \frac{2*d}{c}
τ=c2∗d
中频信号:
I
F
=
S
T
(
t
)
∗
S
T
∗
(
t
)
IF = S_T(t)*S_T^*(t)
IF=ST(t)∗ST∗(t)
仿真程序:
%% 超参数
c = 3e8; %光速
fc = 76.5e9; %发射信号载频 中心频率
bw = 500e6; %发射信号带宽
Tr = 10e-6; %扫频时间 也就是周期
N = 256; %采样点
Fs = 25.6e6; %采样率
M = 256; %chirp的数目
k = bw/Tr; %chirp斜率
index = 1:1:N;
%% 发射信号
t = 0:1/Fs:Tr-1/Fs; %时间向量 确定256个点在一个Tr中的每个时刻
t = t - Tr/2; %将fc作为中心频率
St = 10*exp((1i*2*pi)*(fc*t+k/2*t.^2)); %发射信号
%% 回波信号
distance = 50; %目标距离雷达50m的距离
t_d = 2 * distance / c; %目标距离雷达的延迟
A = 0.8; %回波信号衰减的比例值
Sr = A*10*exp((1i*2*pi)*(fc*(t-t_d)+k/2*(t-t_d).^2)); %回波信号
%% 求回波信号的共轭
Sr_conj = conj(Sr); %求回波信号的共轭
%% 中频信号
IF = St .* Sr_conj; %求中频信号
%% 对中频信号做FFT
range_win = hamming(N);
range_win = range_win';
temp = IF .* range_win;
IF_fft = fft(temp);
freq_bin = (index - 1) * Fs / N;
range_bin = freq_bin * c / 2 / k;
%% 画图
figure;
subplot(2,2,1);plot(index,abs(St));title('发射信号St');
subplot(2,2,2);plot(index,abs(Sr));title('回波信号Sr');
subplot(2,2,3);plot(index,abs(IF));title('中频信号IF');
subplot(2,2,4);plot(range_bin,abs(IF_fft));title('rangefft');
实验结果: