《合成孔徑雷達成像——演算法與實現》之【9】模擬圖3.10
阿新 • • 發佈:2019-01-29
加窗是為了增加目標響應函式的峰值旁邊比,可分為時域窗和頻域窗。
% SAR_Figure_3_10
% 2016.10.12
close all;clear all;clc
T = 7.24e-6; % 訊號持續時間
B = 14.2e6; % 訊號頻寬
K = B/T; % 調頻率
ratio = 1.25; % 過取樣率
Fs = ratio*B; % 取樣頻率
dt = 1/Fs; % 取樣間隔
N = ceil (T/dt); % 取樣點數
t = ((0:N-1)-N/2)/N*T; % 時間軸
f = ((0:N-1)-N/2)/N*Fs; % 頻率軸
st = exp(1i*pi*K*t.^2); % 生成訊號
Sf = fft(st); % FFT
Hf = exp(1i*pi*f.^2/K); % 頻域匹配濾波器
Out = Sf.*Hf; % 頻域匹配濾波
window = kaiser(N,2.5)'; % 時域窗
Window = fftshift(window); % 頻域窗
st_window = st.*window; % 時域加窗
Out_Window = Out.*Window; % 頻域匹配濾波後加窗
figure,set(gcf,' Color','w');
subplot(2,2,1),plot(t*1e6,window);axis tight
title('時域窗函式'),ylabel('幅度');
subplot(2,2,2),plot(Window);axis tight
title('頻域窗函式');
subplot(2,2,3),plot(t*1e6,real(st_window));axis tight
title('訊號實部'),xlabel('時間(\mus)'),ylabel('幅度');
subplot(2,2,4),plot(real(Out_Window));axis tight
title('頻譜實部' ),xlabel('頻率(單元)');