%%?小前斜角CS成像算法
clc
clear?all
close?all

%%?基本參數(shù)
C?=?299792458?;?%?電磁波傳播速度，m/s

%?幾何關(guān)系
Vel?=?100?;????%?速度，m/s
Rref?=?5e3?;????%?中心距離，m
BeamWidAzi?=?2?;????%?方位波束寬度，度
Theta?=?10?;????%?前斜角（波束指向與速度法向的夾角，deg）

%?雷達(dá)參數(shù)?
Lambda?=?0.03?;???%?波長(zhǎng)，m
Prf?=?400?;?%?重頻，Hz
Br?=?300e6?;????%?帶寬，Hz
Fs?=?360e6?;????%?采樣率，Hz
Tp?=?2e-6?;?%?脈寬，s
SignOfChirpSlope?=?1?;??%?調(diào)頻斜率符號(hào)，1表示正調(diào)頻，-1表示負(fù)調(diào)頻

%?點(diǎn)目標(biāo)坐標(biāo)
TarX?=?[0?0?20]?+?Rref*cosd（Theta）?;
TarY?=?[0?20?0]?+?Rref*sind（Theta）?;
%?TarX?=?[0?]?+?Rref*cosd（Theta）?;
%?TarY?=?[0?]?+?Rref*sind（Theta）?;
TarPosi?=?[TarX（:）?TarY（:）]?;
TarNum?=?length（TarX）?;

%?錄取數(shù)據(jù)參數(shù)
Na?=?1024?;????%?方位積累點(diǎn)數(shù)
Nr?=?1024?;?%?距離向點(diǎn)數(shù)
Rmin?=?4.95e3?;????%?起始采樣距離，m
tr?=?2*Rmin/C+（0:Nr-1）/Fs?;??%?距離時(shí)間

%?航跡坐標(biāo)?雷達(dá)沿Y正方向飛行，天線沿X方向照射
RadarX?=?ones（Na1）*0?;
RadarY?=?（-Na/2:Na/2-1）‘/Prf*Vel?;?
RadarPosi?=?[RadarX?RadarY]?;


%%?仿真生成回波
disp（‘正在生成回波。。。‘）?;
Echo?=?zeros（NaNr）?;?
for?kPrt?=?1?:?Na
????for?kTar?=?1?:?TarNum
????????TarR?=?sqrt（sum（?（RadarPosi（kPrt:）-TarPosi（kTar:））.^2?））?;?%?點(diǎn)目標(biāo)的斜距，m
????????TarTau?=?2*TarR/C?;?????%?點(diǎn)目標(biāo)的時(shí)延，s
????????AziAng?=?acosd（dot（TarPosi（kTar:）-RadarPosi（kPrt:）[cosd（Theta）?sind（Theta）]）./TarR）?;??%?點(diǎn)目標(biāo)視線與波束中心的夾角，deg
????????if?abs（AziAng）<=BeamWidAzi/2??%?點(diǎn)目標(biāo)在方位波束內(nèi)
????????????trTar?=?tr-TarTau-Tp/2?;????????????????????%?每個(gè)點(diǎn)目標(biāo)的距離時(shí)間
????????????Echo（kPrt:）?=?Echo（kPrt:）?+?exp（-1j*4*pi*TarR/Lambda）?...
????????????????*rectpuls（trTarTp）?.*?exp（1j*pi*SignOfChirpSlope*Br/Tp*（trTar）.^2）?;
????????end
????end
end
figureimagesc（Rmin+（0:Nr-1）*C/Fs/2（-Na/2:Na/2-1）*Vel/Prfreal（Echo））?;
xlabel（‘距離（km）‘）ylabel（‘方位位置（m）‘）title（‘回波實(shí)部‘）?;
figureimagesc（Rmin+（0:Nr-1）*C/Fs/2（-Na/2:Na/2-1）*Vel/Prfimag（Echo））?;
xlabel（‘距離（km）‘）ylabel（‘方位位置（m）‘）title（‘回波虛部‘）?;

%%?方位FFT
disp（‘正在方位FFT。。。‘）?;
ta?=?（-Na/2:Na/2-1）‘/Prf?;??%?方位時(shí)間，s
fa?=?（（-Na/2:Na/2-1）‘/Na*Prf）?;????%?方位頻率，Hz
Fdc?=?2*Vel*sind（Theta）/Lambda?;????%?多普勒中心，Hz
AziFft?=?zeros（NaNr）?;
for?kR?=?1?:?Nr
????AziFft（:kR）?=?fftshift（?fft（?Echo（:kR）?...
????????.*exp（-1j*2*pi*Fdc*ta）?）?）?;
end
figureimagesc（Rmin+（0:Nr-1）*C/Fs/2fa+Fdcreal（AziFft））?;
xlabel（‘距離（km）‘）ylabel（‘方位頻率（Hz）‘）title（‘方位FFT（實(shí)部）‘）?;

%%?距離CS
disp（‘正在距離CS。。。‘）?;
KsRef?=?（SignOfChirpSlope*Br/Tp）./（1-（SignOfChirpSlope*Br/Tp）*Rref*cosd（Theta）?...
????*2*Lambda/C^2*（Lambda*（fa+Fdc）/2/Vel）.^2./（sqrt（1-（Lambda*（fa+Fdc）/2/Vel）.^2））.^3）?;
Cs?=?1./sqrt（1-（Lambda*（fa+Fdc）/2/Vel）.^2）-1?;
RanCS?=?zeros（NaNr）?;
for?kFa?=?1?:?Na
????RanCS（kFa:）?=?AziFft（kFa:）?.*?exp（1j*pi*KsRef（kFa）*Cs（kFa）?...
????????*（tr-Tp/2-2/C*Rref*cosd（Theta）*（1+Cs（kFa）））.^2）?;
end
figureimagesc（Rmin+（0:Nr-1）*C/Fs/2fa+Fdcreal（RanCS））?;
xlabel（‘距離（km）‘）ylabel（‘方位頻率（Hz）‘）title（‘距離CS（實(shí)部）‘）?;

%%?距離FFT
disp（‘正在距離FFT。。。‘）?;
fr?=?（-Nr/2:Nr/2-1）/Nr*Fs?;??%?距離頻率軸
RanFft?=?fftshift（fft（RanCS[]2）2）?;
figureimagesc（fr/1e6fa+Fdcreal（RanFft））?;
xlabel（‘距離頻率（MHz）‘）ylabel（‘方位頻率（Hz）‘）title（‘距離FFT（實(shí)部）‘）?;

%%?距離壓縮、徙動(dòng)校正、二次距離壓縮
disp（‘正在距離壓縮、徙動(dòng)校正、二次距離壓

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????4681??2017-07-09?15:56??Main_CS.m

?????文件???????4604??2017-07-13?22:18??Main_OmegaK.m

?????文件???????3937??2017-06-25?23:44??Main_RD.m

?????文件?????490391??2017-07-09?16:01??合成孔徑雷達(dá)CS算法點(diǎn)目標(biāo)仿真實(shí)驗(yàn).docx

?????文件?????380731??2017-07-13?22:19??合成孔徑雷達(dá)OmegaK算法點(diǎn)目標(biāo)仿真實(shí)驗(yàn).docx

?????文件?????203702??2017-06-25?23:49??合成孔徑雷達(dá)RD算法點(diǎn)目標(biāo)仿真實(shí)驗(yàn).docx

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

??????????????1088046????????????????????6