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（211）
plot（t*1e6real（singal））;
title（‘近距離時線性調頻信號‘）;
xlabel（‘延時?-?us‘）;
ylabel（‘歸一化幅度‘）;
grid?on;

%線性調頻信號頻譜
singalfft=fftshift（fft（singal））;
Nf=length（singalfft）;
f=[-Nf/2:Nf/2-1]*Fs/Nf;
subplot（212）
plot（f*1e-6abs（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;