天天看點

【雷達通信】非相幹多視處理(CSA)(Matlab代碼實作)

 👨‍🎓個人首頁:​​研學社的部落格​​ 

💥💥💞💞歡迎來到本部落格❤️❤️💥💥

🏆部落客優勢:🌞🌞🌞部落格内容盡量做到思維缜密,邏輯清晰,為了友善讀者。

⛳️座右銘:行百裡者,半于九十。

📋📋📋本文目錄如下:🎁🎁🎁

目錄

​​💥1 概述​​

​​📚2 運作結果​​

​​🎉3 參考文獻​​

​​🌈4 Matlab代碼實作​​

💥1 概述

“在SAR系統中,多個獨立的視可以由飛行載體以不同的方位角通過觀察點時天。

一視由天線沿方位向第一個前向四分之一波束部分産生,下一視則來自下一個四分之一波束,

以此類推。然後,由于來自波束各部分的信号到達雷達接收機是重疊在一起的,所

以在時域或者空域上無法對資料進行分離。然而,具有高方位時間帶寬積的一個實用SAR系

統是将時間和頻率兩者綁定在一起的,在多普勒域内包含了各視的所有資訊。也就是說,具

有較高多普勒頻率的資料一定是由方位向波束前緣觸及到的地形點産生的,而當同一地點出

在方位波束後緣四分之一時,産生了多普勒頻段低四分之一部分。”[1]

📚2 運作結果

【雷達通信】非相幹多視處理(CSA)(Matlab代碼實作)
【雷達通信】非相幹多視處理(CSA)(Matlab代碼實作)
【雷達通信】非相幹多視處理(CSA)(Matlab代碼實作)
【雷達通信】非相幹多視處理(CSA)(Matlab代碼實作)

部分代碼:

Kr = -Kr;                       % 将調頻率Kr改成負值

BW_range = 30.111e+06;          % 脈沖寬度

Vr = 7062;                      % 有效雷達速率

Ka = 1733;                      % 方位調頻率

fnc = -6900;                    % 多普勒中心頻率

Fa = PRF;                       % 方位向采樣率

lamda = c/f0;                   % 波長

T_start = 6.5959e-03;           % 資料窗開始時間

Nr = round(Tr*Fr);              % 線性調頻信号采樣點數

Nrg = Nrg_cells;                % 距離線采樣點數

if b == 1 || b == 2

    Naz = Nrg_lines_blk;         % 每一個資料塊的距離線數

else

    Naz = Nrg_lines;              % 兩個資料塊,總共的距離線數

end

NFFT_r = Nrg;                   % 距離向FFT長度

NFFT_a = Naz;                   % 方位向FFT長度

R_ref = R0;                     % 參考目标選在場景中心,其最近斜距為 R_ref  

fn_ref = fnc;                   % 參考目标的多普勒中心頻率

%%

%

% --------------------------------------------------------------------

% 對原始資料進行補零

% --------------------------------------------------------------------

if b == 1 || b == 2 

    data = zeros(1*2048,3000);

else

    data = zeros(2*2048,3000);

end

data(1:Naz,1:Nrg) = s_echo;

clear s_echo;

s_echo = data;

clear data;

[Naz,Nrg] = size(s_echo);

NFFT_r = Nrg;                   % 距離向FFT長度

NFFT_a = Naz;                   % 方位向FFT長度

% 作圖顯示

figure;

imagesc(abs(s_echo));

title('補零後的原始資料');       % 補零後的原始回波資料(未處理)的幅度圖像

%}

%%

% --------------------------------------------------------------------

% 距離(方位)向時間,頻率相關定義

% --------------------------------------------------------------------

% 距離

tr = 2*R0/c + ( -Nrg/2 : (Nrg/2-1) )/Fr;                % 距離時間軸

fr = ( -NFFT_r/2 : NFFT_r/2-1 )*( Fr/NFFT_r );          % 距離頻率軸

% 方位

ta = ( -Naz/2: Naz/2-1 )/Fa;                            % 方位時間軸

fa = fnc + fftshift( -NFFT_a/2 : NFFT_a/2-1 )*( Fa/NFFT_a );    % 方位頻率軸

% 生成距離(方位)時間(頻率)矩陣

tr_mtx = ones(Naz,1)*tr;    % 距離時間軸矩陣,大小:Naz*Nrg

ta_mtx = ta.'*ones(1,Nrg);  % 方位時間軸矩陣,大小:Naz*Nrg

fr_mtx = ones(Naz,1)*fr;    % 距離頻率軸矩陣,大小:Naz*Nrg

fa_mtx = fa.'*ones(1,Nrg);  % 方位頻率軸矩陣,大小:Naz*Nrg

🎉3 參考文獻

部分理論來源于網絡,如有侵權請聯系删除。

[1]《合成孔徑雷達成像——算法與實作》 ,(美)卡明等著;洪文等譯;電子工業出版社;

[2]《合成孔徑雷達——系統與信号處理》 ,(美)柯蘭德等著;韓傳钊等譯;電子工業出

​​🌈​​4 Matlab代碼實作

繼續閱讀