%
%本函數在5.5節中，用于仿真雷達系統，測試LFM波形的徑向距離分辨率。
%
function LFM_radar(T,B,Rmin,Rmax,R,RCS)
if nargin==0
    T=10e-6;                                       %發射脈寬
    B=30e6;                                        %調頻帶寬
    Rmin=10000;Rmax=15000;                         %距離范圍
    R=[10500,10510,11000,11008,12000,12005,13000,13002];  %目標位置向量，相對于雷達
    RCS=[1 1 1 1 1 1 1 1];                           %目標相對有效反射面
end
%參數設置
R=R-Rmin;                                         %目標位置向量，相對距離門下限Rmin
C=3e8;                                               %光速
K=B/T;                                                %調頻斜率
Rrec=Rmax-Rmin;                              %接收門寬/m
Trec=2*Rrec/C;                                  %接收門寬/s
Fs=5*B;Ts=1/Fs;                                %計算機仿真的采樣頻率和周期
N0=ceil(T/Ts);                                    %發射脈寬對應的采樣點數
N=ceil(Trec/Ts);                           %接收窗對應的采樣點數
%產生高采樣率chirp信號
Rate=10;
t1=linspace(0,T,Rate*N0);
Chirp=exp(j*pi*K*t1.^2);
%產生回波信號                            
M=length(R);
Srt(1:N)=0;               %初始化接收信號
for i=1:1:M
    temp(1:N)=0;
    tau=2*R(i)/C;
    K0=ceil(tau/Ts)                %第i個目標回波的相對起始位置
    M0=fix((Ts*K0-tau)*Rate/Ts)
    temp(K0:K0+N0-1)=RCS(i)*Chirp(M0+1:Rate:M0+Rate*(N0-1)+1);
    Srt=Srt+temp;                              %將所有的目標回波疊加
end
%作出相關圖形
figure(1)
subplot(211)
t=linspace(0,Trec,N);
plot(t*10^6,real(Srt));  
xlabel('\fontsize{9}時間(us)');title('\fontsize{9}脈沖壓縮前的回波');
subplot(212)
t1=linspace(-T/2,T/2,N0);
h=exp(-j*pi*K*t1.^2);                               %匹配濾波器單位沖激響應
Sot=conv(Srt,h);                                    %匹配濾波
Z=abs(Sot);Z=Z/max(Z);
Z=20*log10(Z+1e-6);
dist=linspace(0,Rrec,N);
plot(dist+Rmin,Z(N0/2:N+N0/2-1));
axis([Rmin,Rmax,-60,0]);
xlabel('\fontsize{9}目標相對雷達的距離(m)');title('\fontsize{9}脈沖壓縮后的回波');
%作出一張圖，只有脈壓后的波形，用于放大觀察
figure(2)
dist2=linspace(0,Rrec,N);
plot(dist2+Rmin,Z(N0/2:N+N0/2-1));
axis([Rmin,Rmax,-60,0]);
xlabel('\fontsize{9}目標相對雷達的距離(m)');title('\fontsize{9}脈沖壓縮后的回波');