%彈頭雷達1成像
%飛行器雷達1成像
clear
clc

Nr?=?21;
Na?=?Nr;
Br=1.5e9;

pi?=?3.141592653589793;
eps0?=?8.854187817?*?1.d-12;
mu0?=?4.0?*?pi?*?1.0e-7;
cj?=?sqrt（-1）;
c?=?1.0?/?sqrt（eps0?*?mu0）;
RangeReso=c/2.0/Br;

%---------------------------------------------
%讀取圖像數據
fid?=?fopen（‘re.dat‘?‘r‘）;
z1?=?fread（fid?‘double‘）;
fclose（fid）;
fid?=?fopen（‘im.dat‘?‘r‘）;
z2?=?fread（fid?‘double‘）;
fclose（fid）;

cz?=?cj?*?z2（1?:?Nr?*?Na）;
clear?z2;
cz?=?cz?+?z1（1?:?Nr?*?Na）;
clear?z1;

cim?=?reshape（cz?Nr?Na）;??????%到這里cim和BP成像程序里輸出的echo1圖像數據一樣
clear?cz;
%---------------------------------------------

%---------------------------------------------
%對圖像數據插值
[xy]?=?meshgrid（1:1:Nr）;
[x1y1]?=?meshgrid（1:0.1:Nr）;
cim?=?interp2（x?y?cim?x1?y1）;
%這個函數的功能是利用（x?y）對應網格點的已知數據cim?對（x1y1）對應的更密網格上的未知數據進行插值，使圖像更加平滑；
im?=?abs（cim）;
im?=?db（?im?/?（max（max（im））+1.0e-8）?）;
%---------------------------------------------

%---------------------------------------------
%計算坐標軸
temp=（Nr-1）*5+1;???%插值后總點數
x?=?linspace（1?temp?temp）;
x=int16（x）;
x=x-int16（temp）/2;
xaxis=RangeReso/5.0*double（x）;??????%坐標軸
%---------------------------------------------

%---------------------------------------------
%畫圖
figure（1）;
clim?=?[-50];?????????????????????????%clim指定圖像中數值范圍與當前色圖的映射，clim中指定數據會映射到整個色圖范圍上；
figure（1）
imagesc（xaxisxaxisim?clim）;

??axis（[-2?2?-2?2]）
colormap（jet）?????????????%色圖設置為jet型，顏色變化為藍-黃-紅，
shading?interp??????????????????????????%通過插值使顏色漸變
xlabel（‘距離向?（m）‘）;?
ylabel（‘方位向?（m）‘）;
title（‘OGSBL重構回波成像‘）
%?set（gca‘FontName‘‘Times?New?Roman‘‘FontSize‘12）;
%---------------------------------------------

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1851??2018-07-06?11:41??BP\draw_mis_new.m

?????文件?????143157??2018-05-21?10:58??BP\f15_fbp_imag.txt

?????文件?????143157??2018-05-21?10:58??BP\f15_fbp_real.txt

?????文件????????576??2013-03-01?11:15??BP\fighterSC.mat

?????文件???????3250??2018-07-02?22:49??BP\Figure4_23and24.m

?????文件???????4264??2018-07-05?16:34??BP\Figure4_26thru4_28.m

?????文件????????746??2018-06-19?17:09??BP\from_gu_shixiong\f15\cx1.m

?????文件???????5226??2018-06-27?21:04??BP\from_gu_shixiong\f15\DirIntA2.m

?????文件?????143157??2018-05-21?10:58??BP\from_gu_shixiong\f15\f15_fbp_imag.txt

?????文件?????143157??2018-05-21?10:58??BP\from_gu_shixiong\f15\f15_fbp_real.txt

?????文件??????53227??2018-06-19?16:35??BP\from_gu_shixiong\f15\f15_fft_imag.txt

?????文件??????53227??2018-06-19?16:35??BP\from_gu_shixiong\f15\f15_fft_real.txt

?????文件??????53227??2017-11-22?17:33??BP\from_gu_shixiong\f15\f15_vv_imag.txt

?????文件??????53227??2017-11-22?17:33??BP\from_gu_shixiong\f15\f15_vv_real.txt

?????文件???????4089??2018-07-06?15:33??BP\from_gu_shixiong\f15\f15_xiaodegree.m

?????文件????????519??2018-06-19?16:52??BP\from_gu_shixiong\f15\fft_isar.m

?????文件?????170299??2018-07-06?12:19??BP\from_gu_shixiong\f15.zip

?????文件???????3528??2018-07-06?11:38??BP\im.dat

?????文件?????524834??2018-07-06?10:55??BP\imag.txt

?????文件????????193??2013-03-01?11:15??BP\matplot2.m

?????文件???????2462??2018-07-06?12:15??BP\qiu_lunwen.m

?????文件???????3528??2018-07-06?11:38??BP\re.dat

?????文件?????524718??2018-07-06?10:55??BP\real.txt

?????目錄??????????0??2018-06-19?17:17??BP\from_gu_shixiong\f15

?????目錄??????????0??2018-07-06?12:21??BP\from_gu_shixiong

?????目錄??????????0??2019-01-22?15:35??BP

-----------?---------??----------?-----??----

??????????????2035619????????????????????26