function?value_ISLR=ISLR（line）
%?Copyright（c）2004Beijing?Institute?of?TechnologyRadar?Research?Lab
%?All?rights?reserved.
%?
%?Filename:	ISLR.m
%?Description:	the?ISLR?of?range?line?or?azimuth?line?after?compression
%
%?Author:	zss
%?Date:	2005.4.17
%?
%?Modification?History
%?--------------------
%
%?測試程序
%?c?=?3e8;
%?wave?=?0.05;??????????????
%?Tp?=?20e-6;????????
%?B?=?150e6;?????????
%?fs?=?200e6;?
%?Kr?=?B?/?Tp;
%?t?=?-Tp/2?:?1/fs?:?（Tp/2?-?1/fs）;
%?L?=?length（t）;
%?chirp?=?exp（1i?*?pi?*?Kr?*?t.^2）;
%?YY?=fftshift（fft（chirp））;
%?figureimagesc（abs（YY））;
%?chirp_ref?=?conj（fftshift（fft（chirp）））;
%?line?=?fftshift（ifft（fftshift（（YY.*chirp_ref））））;
%?figureimagesc（abs（line））;
%?figureplot（abs（line））;
%?求某一條距離線或方位線的峰值旁瓣比，信號頻譜一般為偶數點，一般為8倍插值
interp?=?8;
fft_line?=?fft（line）;
%?figure（1）
%?plot（abs（fft_line））;
[rowcol]?=?size（line）;
fft_line_interp?=?zeros（1interp?*?col）;
%?信號的頻譜在兩邊
fft_line_interp（1?:?col/2）?=?fft_line（1?:?col/2）;
fft_line_interp（interp?*?col?-?col/2?+?1?:?interp?*?col）?=?fft_line（col/2?+?1?:?col）;
%?信號的頻譜在中央
%?fft_line_interp（col?*?7/2?+?1?:?col?*?9/2）?=?fft_line（1?:?col）;
result_interp?=?interp?*?ifft（fft_line_interp）;?%乘以插值系數是為了插值后頻譜的幅度不變
temp??=?abs（result_interp）;
[value_mainlobeposition_mainlobe]?=?max（temp）;
db_result?=?20?*?log10（temp?/?max（temp））;
%?figure
%?plot（db_result）;
right_flag?=?1;
while?right_flag
????for?k?=?1?:?interp?*?col
????????if?temp（position_mainlobe?+?k）?<=?value_mainlobe?*?sqrt（2）/2
????????????right_posi?=?k;
????????????right_flag?=?0;
????????????break;
????????end
????end
end
right_width?=?right_posi?-1?+?（temp（position_mainlobe?+?right_posi?-?1）?-?...
??????????value_mainlobe?*?sqrt（2）/2）?/?（temp（position_mainlobe?+?right_posi?-?1）?-?...
??????????temp（position_mainlobe?+?right_posi））;

left_flag?=?1;
while?left_flag
????for?k?=?1?:?interp?*?col
????????if?temp（position_mainlobe?-?k）?<=?value_mainlobe?*?sqrt（2）/2
????????????left_posi?=?k;
????????????left_flag?=?0;
????????????break;
????????end
????end
end?
left_width?=?left_posi?-?1?+?（temp（position_mainlobe?-?left_posi?+?1）?-?...
??????????value_mainlobe?*?sqrt（2）/2）?/?（temp（position_mainlobe?-?left_posi?+?1）?-?...
??????????temp（position_mainlobe?-?left_posi））;?????
%?width_mainlobe?=?left_width?+?right_width;
left_width=2.5*left_width;
right_width=2.5*right_width;
energy_total?=?sum（temp.^2）;
energy_mainlobe?=?sum（temp（（position_mainlobe?-?ceil（left_width））?:?（position_mainlobe?+?ceil（right_width）））.^2）;
value_ISLR?=?-10?*?log10（energy_mainlobe?/?（energy_total?-?energy_mainlobe））;
%?return;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2703??2018-11-11?17:31??新建文件夾\ISLR.m
?????文件?????????867??2018-11-14?13:14??新建文件夾\NLFM.m
?????文件????????2184??2018-11-11?17:27??新建文件夾\PSLR.m
?????目錄???????????0??2018-11-24?16:59??新建文件夾\