%%?同距離選大距離像拼接
clc;
close?all;
clear?all

%%?系統(tǒng)參數(shù)
tao?=?1e-6/37.5;??????%?時(shí)寬，對(duì)應(yīng)帶寬37.5MHz
PRT?=?2e-6;?????????%?脈沖重復(fù)周期?2us
detf?=?30e6;??????????%?步進(jìn)頻率間隔?30MHz
c?=?3e8;??????????????%?光速?
N?=?50;???????????????%?步進(jìn)頻率跳變次數(shù)?
B?=?N?*?detf;?????????%?雷達(dá)合成帶寬?1.5GHz
Rmax?=?PRT?*?c?/?2;???%?最大不模糊距離

f0?=?9.025e9;?????????%?載頻起始頻率9.025GHz
fs?=?50e6;????????????%?回波采樣頻率50MHz???????

nprt?=?ceil（（PRT-tao）*fs）;?????????%?1個(gè)prt內(nèi)采樣的點(diǎn)數(shù)?ceil?向上取整
ntao?=?ceil（tao*fs）;???????????????%?1個(gè)脈寬內(nèi)采樣的點(diǎn)數(shù)

detR?=?c/（2*B）;???????%?細(xì)化后的距離分辨率
Rfs?=?c/2/fs;?????????%?采樣距離分辨率
Rtao?=?c*tao/2;???????%?單脈沖分辨距離
Ri?=?c/2/detf;????????%?單點(diǎn)不模糊距離
Ru?=?c/2/detf/N;
ts_num?=?Rfs/Ru;??????%每個(gè)采樣點(diǎn)代表的有效數(shù)據(jù)段，60%的有效，重疊20%，冗余20%
ttao_num?=?round（Rtao/Ru）;???%每個(gè)脈寬代表的有效數(shù)據(jù)段，80%的有效，冗余20%
Overlapping=ttao_num-ts_num;??%?每個(gè)脈沖重疊數(shù)據(jù)段

snr?=?2;??????%回波信號(hào)疊加信噪比大小
R0?=?30;??????%目標(biāo)所在距離
V0?=0;????????%假設(shè)為靜止目標(biāo)回波
ntao?=?3;?????%回波一般都有展寬，從而ntao設(shè)為3，而不是2

%%?坐標(biāo)軸參數(shù)設(shè)置
range_axis?=?0?:?Ru?:?（ts_num*nprt-1）*Ru;
range_axis?=?range_axis?-?3;

%%?生成回波信號(hào)
x1?=?zeros（Nnprt）;
for?n?=?0:N-1
????t?=?n*PRT;
????R（n+1）?=?R0-V0*t;
????DELAY?=?2*R（n+1）/c;
????temp?=?abs（DELAY*fs）;
????if?（?abs（temp-round（temp））????????ntarget?=?round（temp）;
????else
????????ntarget?=?ceil（temp）-1;???%?ceil/floor?根據(jù)數(shù)據(jù)定?
????end
????%%%%%%%%%??高斯型的回波?%%%%%%%%%%
????x1（n+1:）?=?[zeros（1ntarget）exp（-1j*2*pi*（f0+n*detf）*DELAY）...
?????????????????*Func_Gauss（ntao1）zeros（1nprt-ntao-ntarget）];?
end????%%%?結(jié)束N個(gè)脈沖，產(chǎn)生回波

%?%?echo?=?x1;
echo?=?awgn（x1snr‘measured‘‘linear‘）;????%回波信號(hào)中加入噪聲

%%?一維成像
echo_ifft?=?ifft（echoN）;??????%IFFT處理
figure;
mesh（abs（echo_ifft））;
xlabel（‘采樣點(diǎn)‘）;ylabel（‘距離分辨單元‘）;
zlabel（‘幅度‘）;
title（‘各頻點(diǎn)IFFT相參脈壓結(jié)果‘）