clear all;
%定義參數
B=10e6;       % B ——帶寬
Fc=10e9;     % Fc——載頻
Fs=10*B;     % Fs——采樣率?
T2=10e-6;   % T1——脈沖寬度
Tr=375e-6;   % Tr——PRT
C=3e8;	

%產生線性調頻脈沖信號
mu=B/T2;
NT2=[Fs*T2];               %采樣點數
t=[-NT2/2:NT2/2-1]*T2/NT2; %取一個脈沖
fai=Fc+mu*t.*t/2;          %相位 
singal=exp(2*pi*j*fai);    %LFM信號，幅度為1
figure(1);
subplot(2,1,1)
plot(t*1e6,real(singal));
title('近距離時線性調頻信號');
xlabel('延時 - us');
ylabel('歸一化幅度');
grid on;

%線性調頻信號頻譜
singalfft=fftshift(fft(singal));
Nf=length(singalfft);
f=[-Nf/2:Nf/2-1]*Fs/Nf;
subplot(2,1,2)
plot(f*1e-6,abs(singalfft)/max(abs(singalfft)));
title('近距離時的信號頻譜圖');
xlabel('頻率f/MHz');
ylabel('歸一化幅度');
grid on;

%模糊函數
Taw=100;
i=0;
for t=[-Taw:Taw-1]*T2/Taw
    i=i+1;
    j=0;
    for Fd=[-Taw:Taw-1]*B/Taw;
        j=j+1;
        val1=1-abs(t)/T2;
        val2=pi*T2*val1;
        val3=Fd+mu.*t+eps;
        val4=val2*val3;
        ambg(j,i)=abs(val1*(sin(val4)/(val4))).^2;
    end
end
t=[-Taw:Taw-1]*T2/Taw;
Fd=[-Taw:Taw-1]*B/Taw;
figure(2);
mesh(t*1e6,Fd*1e-6,ambg);
title('近距離時的模糊函數圖');
xlabel('延時-us');
ylabel('多普普勒頻移fd/MHz');
zlabel('歸一化幅度');
grid on;
%繪圖
%模糊函數等高線圖
figure(3);
contour(t*1e6,Fd*1e-6,ambg);
title('近距離時模糊函數等高線圖');
xlabel('延時-us');
ylabel('多普普勒頻移fd/MHz');
grid on;

%自相關函數
figure(4);
subplot(2,1,1);
plot(t*1e6,abs(ambg(:,Taw+1)));
title('近距離時自相關函數');
xlabel('延時-us');
ylabel('歸一化幅度');
grid on;

%多普勒敏感度
subplot(2,1,2)
plot(Fd*1e-6,abs(ambg(Taw+1,:)));
title('近距離時多普勒敏感度');
xlabel('多普勒頻移-MHz');
ylabel('歸一化幅度');
grid on;