clc
clear
close all
%% 距离向参数
R_eta_c = 850e+3; % 景中心斜距
%% 方位向参数
V_r = 7100; % 等效速度
T_a = 0.64; % 合成孔径时间
K_a = 2095; % 方位向调频率
theta_r_c = -0.3; % 斜视角
%% 参数计算
Delta_f_dop = abs(K_a*T_a); % 多普勒带宽
t_eta_c = -R_eta_c*sind(theta_r_c)/V_r; % 波束中心穿越时刻
%% 参数设置
alpha_a_s = 1; % 方位向过采样率
F_a = alpha_a_s*Delta_f_dop; % 方位向采样率PRF
N_a = 2*ceil(F_a*T_a/2); % 方位向采样点数
dt = T_a/N_a; % 采样时间间隔
df = F_a/N_a; % 采样频率间隔
%% 变量设置
t_eta = (-T_a/2:dt:T_a/2-dt)+t_eta_c; % 时间轴
%% 信号表达
RCM_1 = -V_r*sind(theta_r_c)*(t_eta-t_eta_c+T_a/2); % 线性分量
RCM_2 = (1/2)*(V_r^2*cosd(theta_r_c)^2/R_eta_c)*(t_eta-t_eta_c).^2;
%% 绘图
plot(RCM_1,t_eta),hold on
plot(RCM_1+RCM_2,t_eta),set(gca,'ydir','reverse')
axis([-5 35 0 1.2])
