%近程模式下的單脈沖信號產生，頻譜圖,模糊函數圖,模糊度等高線圖，自相關函數圖和多普勒敏感度：
clear all;close all;
%定義參數
B=10e6;       % B ——帶寬
Fc=10e9;      % Fc——載頻
Fs=10*B;      % Fs——采樣率
T2=10e-6;     % T2——脈沖寬度
Tr=375e-6;    % Tr——帶寬
C=2.997e8;

%產生線性調頻脈沖信號
mu=B/T2;
NT2=[Fs*T2];
t=[-NT2/2:NT2/2-1]*T2/NT2; 
f=Fc+mu*t.*t/2;
sig=exp(2*pi*j*f);
figure(1);
plot(t*1e6,real(sig));
title('近程模式下線性調頻脈沖信號');
xlabel('延時（一個脈寬時間內）/us');
ylabel('歸一化幅度');
grid on;
%線性調頻脈沖信號頻譜
sigfft=fftshift(fft(sig));
Nf=length(sigfft);
f=[-Nf/2:Nf/2-1]*Fs/Nf;
figure(2);
plot(f*1e-6,abs(sigfft)/max(abs(sigfft)));
title('近程模式下的信號頻譜圖');
xlabel('頻率/MHz');
ylabel('歸一化幅度');
grid on;

%模糊函數
Naw=100;
Fdmax=B;
i=0;
for t=[-Naw:Naw-1]*T2/Naw
    i=i+1;
    j=0;
    for Fd=[-Naw:Naw-1]*Fdmax/Naw;
        j=j+1;
        vala=(T2-abs(t))/T2;
        valb=T2-abs(t);
        valc=pi*(Fd+mu.*t);
        vald=valb*valc+eps;
        ambg(j,i)=abs(vala*(sin(vald)/(vald))).^2;
    end
end
Td=[-Naw:Naw-1]*T2/Naw;
Fd=[-Naw:Naw-1]*Fdmax/Naw;
figure(3);
mesh(Td*1e6,Fd*1e-6,ambg);
title('近程模式下的模糊函數圖');
xlabel('延時/us');
ylabel('多普普勒頻移/MHz');
zlabel('歸一化幅度');
grid on;

%模糊度等高線圖
figure(4);
contour(Td*1e6,Fd*1e-3,ambg);
title('近程模式下的模糊度等高線圖');
xlabel('延時/us');
ylabel('多普普勒頻移/KHz');
grid on;

%模塊功能：自相關函數
figure(5);
plot(Td*1e6,abs(ambg(:,Naw+1)));
title('近程模式下的自相關函數');
xlabel('延時/us');
ylabel('歸一化幅度');
grid on;

%多普勒敏感度
figure(6);
plot(Fd*1e-3,abs(ambg(Naw+1,:)));
title('近程模式下的多普勒敏感度');
xlabel('多普勒頻移/KHz');
ylabel('歸一化幅度');
grid on;