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

%%?系統參數
tao?=?1e-6/37.5;??????%?時寬，對應帶寬37.5MHz
PRT?=?2e-6;?????????%?脈沖重復周期?2us
detf?=?30e6;??????????%?步進頻率間隔?30MHz
c?=?3e8;??????????????%?光速?
N?=?50;???????????????%?步進頻率跳變次數?
B?=?N?*?detf;?????????%?雷達合成帶寬?1.5GHz
Rmax?=?PRT?*?c?/?2;???%?最大不模糊距離

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

nprt?=?ceil（（PRT-tao）*fs）;?????????%?1個prt內采樣的點數?ceil?向上取整
ntao?=?ceil（tao*fs）;???????????????%?1個脈寬內采樣的點數

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

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

%%?坐標軸參數設置
range_axis?=?0?:?Ru?:?（ts_num*nprt-1）*Ru;
range_axis?=?range_axis?-?3;

%%?生成回波信號
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?根據數據定?
????end
????%%%%%%%%%??高斯型的回波?%%%%%%%%%%
????x1（n+1:）?=?[zeros（1ntarget）exp（-1j*2*pi*（f0+n*detf）*DELAY）...
?????????????????*Func_Gauss（ntao1）zeros（1nprt-ntao-ntarget）];?
end????%%%?結束N個脈沖，產生回波

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

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