《合成孔径雷达成像——算法与实现》之【9】仿真图3.10_destiny0321的博客-CSDN博客
本节就是讨论,回波信号加窗处理的影响。
% 这是一个线性调频脉冲的相位与频率
clear;clc;close all;
T = 7.24e-6; %持续时间为 7.24 us
B = 14.2e6; % 带宽 5.8MHz
K = B/T; % 调频率?? 这是什么意思?
ratio = 1.25; % 过采样倍数: 比如你是 100hz 则用 500hz进行采样
Fs=B*ratio; % 采样频率
dt=1/Fs; % 每次采样的间隔时间
N = ceil(T/dt); % 采样点数
t = ((0:N-1) - N/2)/N * T; % 时间轴
t1 = ((0:N-1) - N/2) * dt; % 时间轴的另一个求法,
f = ((0:N-1)-N/2)/N*Fs; % 频率轴
st = exp(1i*pi*K*t.^2);
st_fft = fft(st);
hf = exp(1i*pi*f.^2/K);
Out = st_fft.*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('频率(单元)');
重点:
kaiser有个参数,可以根据不同应用进行调节!!