clc
clear
close all参数设置
距离向参数设置
R_eta_c = 20e3; % 景中心斜距
Tr = 2.5e-6; % 发射脉冲时宽
Kr = 20e12; % 距离向调频率
alpha_os_r = 1.2; % 距离过采样率
Nrg = 320; % 距离线采样数
距离向参数计算
Bw = abs(Kr)*Tr; % 距离信号带宽
Fr = alpha_os_r*Bw; % 距离向采样率
Nr = round(Fr*Tr); % 距离采样点数（脉冲序列长度）
方位向参数设置
c = 3e8; % 光速
Vr = 150; % 等效雷达速度
Vs = Vr; % 卫星平台速度
Vg = Vr; % 波束扫描速度
f0 = 5.3e9; % 雷达工作频率
Delta_f_dop = 80; % 多普勒带宽
alpha_os_a = 1.25; % 方位过采样率
Naz = 256; % 距离线数
theta_r_c = 3.5; % 波束斜视角
方位向参数计算
lambda = c/f0; % 雷达工作波长
eta_c = -R_eta_c*sind(theta_r_c)/Vr;% 波束中心偏移时间
f_eta_c = 2*Vr*sind(theta_r_c)/lambda;% 多普勒中心频率
La = 0.886*2*Vs*cosd(theta_r_c)/Delta_f_dop;% 实际天线长度
Fa = alpha_os_a*Delta_f_dop;% 方位向采样率
Ta = 0.886*lambda*R_eta_c/(La*Vg*cosd(theta_r_c));% 目标照射时间
R0 = R_eta_c*cosd(theta_r_c);% 最短斜距
Ka = 2*Vr^2*cosd(theta_r_c)^3/(lambda*R0);% 方位向调频率
theta_bw = 0.886*lambda/La; % 方位向3dB波束宽度
theta_syn = Vs/Vg*theta_bw; % 合成角
Ls = R_eta_c*theta_syn; % 合成孔径
其他参数计算
rho_r = c/2/Bw; % 距离向分辨率
rho_a = La/2; % 方位向分辨率
Trg = Nrg/Fr; % 发射脉冲宽度
Taz = Naz/Fa; % 目标照射时间
d_t_tau = 1/Fr; % 距离向采样时间间隔
d_t_eta = 1/Fa; % 方位向采样时间间隔
d_f_tau = Fr/Nrg; % 距离向采样频率间隔
d_f_eta = Fa/Naz; % 方位向采样频率间隔目标设置
设置目标点距离景中心的距离
A_r = -50;A_a = -50;
B_r = -50;B_a = +50;
C_r = +50;C_a = B_a+(C_r-B_r)*tand(theta_r_c);
坐标
A_x = R0+A_r;A_y = A_a;
B_x = R0+B_r;B_y = B_a;
C_x = R0+C_r;C_y = C_a;
N_position = [A_x,A_y;B_x,B_y;C_x,C_y];
波束中心穿越时刻
N_target = 3;
Target_eta_c = zeros(1,N_target);
for i = 1:N_targetDelta_Y = N_position(i,2)-N_position(i,1)*tand(theta_r_c);Target_eta_c(i) = Delta_Y/Vs;
end
绝对零多普勒时刻
Target_eta_0 = zeros(1,N_target);
for i = 1:N_targetTarget_eta_0(i) = N_position(i,2)/Vs;
end变量设置
时间变量：以景中心绝对零多普勒时刻作为方位向零点
t_tau = (-Trg/2:d_t_tau:Trg/2-d_t_tau)+2*R_eta_c/c; % 距离时间变量
t_eta = (-Taz/2:d_t_eta:Taz/2-d_t_eta)+eta_c; % 方位时间变量
r_tau = (t_tau*c/2)*cosd(theta_r_c); % 距离变量
% r_tau = (t_tau*c/2); % 距离变量
频率变量
f_tau = -Fr/2:d_f_tau:Fr/2-d_f_tau; % 距离频率变量
f_tau = f_tau-round((f_tau-0)/Fr)*Fr; % 将频率折叠入(-Fr/2,Fr/2)，距离可观测频率变量
f_eta = -Fa/2:d_f_eta:Fa/2-d_f_eta; % 方位频率变量
f_eta = f_eta-round((f_eta-f_eta_c)/Fa)*Fa; % 将频率折叠入f_eta_c附近(-Fa/2,Fa/2)范围，方位可观测频率变量
坐标设置
[t_tauX,t_etaY] = meshgrid(t_tau,t_eta); % 距离时间X轴，方位时间Y轴
[f_tauX,f_etaY] = meshgrid(f_tau,f_eta); % 距离频域X轴，方位频域Y轴信号设置，原始回波生成
tic % 计时，与toc搭配使用
wait_title = waitbar(0,'开始生成回波数据 ...');
pause(1);
st_tt = zeros(Naz,Nrg);
for i = 1:N_targetR_eta = sqrt(N_position(i,1)^2+Vs^2*(t_etaY-Target_eta_0(i)).^2);% 瞬时斜距，还有近似公式可以尝试A0 = [1,1,1,1]*exp(+1j*0); % 后向散射系数wr = (abs(t_tauX-2*R_eta/c)<=Tr/2); % 距离向包络wa = sinc(0.886*atan(Vs*(t_etaY-Target_eta_c(i))/N_position(i,1))/theta_bw).^2;% 方位向包络，用波束穿越时刻
% wa = sinc(0.886*(atan(Vs*(t_etaY-Target_eta_0(i))/N_position(i,1))+theta_r_c)/theta_bw).^2;st_tt_target = A0(i)*wr.*wa.*exp(-1j*4*pi*f0*R_eta/c)....*exp(1j*pi*Kr*(t_tauX-2*R_eta/c).^2);st_tt = st_tt+st_tt_target;pause(0.001);time = toc;Display_Data = num2str(roundn(i/N_target*100,-1));Display_Str = ['Computation Progress',Display_Data,'%',' --- ',...'Using Time: ',num2str(time)];waitbar(i/N_target,wait_title,Display_Str); % 三参数：进度，句柄，展示的话
end
pause(1);
close(wait_title);
toc一次距离压缩
方式三：根据脉冲频谱特性直接在频域生成频域匹配滤波器
window = kaiser(Nrg,2.5)'; % 时域窗
Window = fftshift(window); % 频域窗
% figure,plot(window)
% figure,plot(Window)
Hrf = (abs(f_tauX)<=Bw/2).*Window.*exp(1j*pi*f_tauX.^2/Kr);
Sf_ft = fft(st_tt,Nrg,2);
Srf_tf = Sf_ft.*Hrf;
srt_tt = ifft(Srf_tf,Nrg,2);方位向FFT
Saf_tf = fft(srt_tt,Naz,1);绘图
H = figure();
set(H,'position',[50,50,600,300]);
subplot(121)
imagesc(real(Saf_tf)),set(gca,'YDir','normal'),xlim([0,150])
xlabel('距离时间(采样点)'),ylabel('方位时间(采样点)'),title('(a)实部')
subplot(122)
imagesc(abs(Saf_tf)),set(gca,'YDir','normal'),xlim([0,150])
xlabel('距离时间(采样点)'),ylabel('方位时间(采样点)'),title('(b)幅度')
