資源簡介
利用步進頻波形應(yīng)用于穿墻雷達,實現(xiàn)目標定位研究算法。此算法為當前熱門重點算法,為后續(xù)壓縮感知算法可作為鋪墊。
代碼片段和文件信息
%計算并畫出SFW的合成HRR像
%該函數(shù)使用長度為步長兩倍的逆傅里葉變化IFFT
function?[hl]?=?hrr_profile(nscatscat_rangescat_rcsndeltafprfvrnotewinid)
%nscat組成目標散射體的個數(shù)
%scat_range包含散射體距離的向量
%scat_rcs?包含散射體RCS的向量
%n步進數(shù)
%deltaf步進頻率
%prf?SFW的PRF
%v?目標速度
%rnote?像起始距離
%winid?>?0為漢明窗??<0不加窗
%hl?距離像:是指目標散射中心在雷達視線上的投影,反映目標各散射中心間的徑向位置關(guān)系,體現(xiàn)目標的部分結(jié)構(gòu)特征,是雷達感知、識別目標的重要途徑。
c?=?3e8;?%光速
num_pulses?=?n;?%步進數(shù)
SNR_dB?=?40;??
nfft?=?256;?%IFFT點數(shù)
carrier_freq?=?9.5e9;?%載波頻率
freq_step?=?deltaf;?%步進頻率
V?=?v;?%雷達速度
PRI?=?1./prf;?%脈沖重復間隔
if?(nfft>2*num_pulses)
????num_pulses?=?nfft/2;
end
%初始化矩陣
Inphase?=?zeros((2*num_pulses)1);
Quadrature?=?zeros((2*num_pulses)1);?
Inphase_tgt?=?zeros(num_pulses1);
Quadrature_tgt?=?zeros(num_pulses1);
IQ_freq_domain?=?zeros((2*num_pulses)1);
Weighted_I_freq_domain?=?zeros((num_pulses)1);
Weighted_Q_freq_domain?=?zeros((num_pulses)1);
Weighted_IQ_freq_domain?=?zeros((2*num_pulses)1);
Weighted_IQ_time_domain?=?zeros((2*num_pulses)1);
abs_Weighted_IQ_time_domain?=?zeros((2*num_pulses)1);
dB_abs_Weighted_IQ_time_domain?=?zeros((2*num_pulses)1);
taur?=?2.*rnote/c;?%起始時間延遲
for?jscat?=?1:nscat
????ii?=?0;
????for?i?=?1:num_pulses
????????ii?=?ii?+?1;
????????rec_freq?=?((i-1)*freq_step);?%fi
????????%實數(shù)分量????利用歐拉公式
????????Inphase_tgt(ii)?=?Inphase_tgt(ii)?+?sqrt(scat_rcs(jscat))*cos(-2*pi*rec_freq*(2.*scat_range(jscat)/c-2*(V/c)*((i-1)*PRI?+?taur/2?+?2*scat_range(jscat)/c)));
????????%虛部分量
????????Quadrature_tgt(ii)?=?Quadrature_tgt(ii)?+?sqrt(scat_rcs(jscat))*sin(-2*pi*rec_freq*(2.*scat_range(jscat)/c-2*(V/c)*((i-1)*PRI?+?taur/2?+?2*scat_range(jscat)/c)));
????end
end
if?(winid?>=?0)?%加漢明窗
????window(1:num_pulses)?=?hamming(num_pulses);
else
????window(1:num_pulses)?=?1;
end
Inphase?=?Inphase_tgt;
Quadrature?=?Quadrature_tgt;
Weighted_I_freq_domain(1:num_pulses)?=?Inphase(1:num_pulses).*window‘;
Weighted_Q_freq_domain(1:num_pulses)?=?Quadrature(1:num_pulses).*window‘;
Weighted_IQ_freq_domain(1:num_pulses)?=?Weighted_I_freq_domain?+?Weighted_Q_freq_domain*j;???%虛部和實部的合并
Weighted_IQ_freq_domain(num_pulses:2*num_pulses)?=?0.+0.i;
Weighted_IQ_time_domain?=?(ifft(Weighted_IQ_freq_domain));?%換算時域上
abs_Weighted_IQ_time_domain?=?abs(Weighted_IQ_time_domain);
dB_abs_Weighted_IQ_time_domain?=?20*log10(abs_Weighted_IQ_time_domain)?+?SNR_dB;
%計算未定義的范圍窗口大小
Ru?=?c/2/deltaf;
hl?=?dB_abs_Weighted_IQ_time_domain;
numb?=?2*num_pulses;
delx_meter?=?Ru/numb;
xmeter?=?0:delx_meter:Ru-delx_meter;
plot(xmeterdB_abs_Weighted_IQ_time_domain‘k‘);
xlabel(‘相對距離/m‘);
ylabel(‘一維距離像‘);
grid
????????
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????2831??2018-03-19?15:03??hrr_profile.m
?????文件???????1619??2018-04-23?20:32??Wall_echo.m
-----------?---------??----------?-----??----
?????????????????4450????????????????????2
評論
共有 條評論