clc
clear
close all
% 参数设置
TBP = 724; % 时间带宽积
T = 42e-6; % 脉冲持续时间
t_0 = 1e-6; % 脉冲回波时延
Nfft = 2^11; % fft长度
% 参数计算
B = TBP/T; % 信号带宽
K = B/T; % 线性调频频率
alpha_os = 1.07; % 过采样率
F = alpha_os*B; % 采样频率
N = 2*ceil(F*T/2); % 采样点数
dt = T/N; % 采样时间间隔
df = F/N; % 采样频率间隔
% 变量设置
t = -T/2:dt:T/2-dt; % 时间变量
f = -F/2:df:F/2-df; % 频率变量
% 信号表达
st = exp(1j*pi*K*t.^2); %chirp信号时域表达式
Sf = fft(st); %chirp信号频域表达式
%窗函数
window = kaiser(N,2.5)'; %时域窗函数
Window = fftshift(window); %频域窗函数
%信号变换
ht = st; %复制信号
ht_window = window.*ht; %加窗
Hf_zero = conj(fft(ht_window,Nfft)); %补零,傅里叶变换,复共轭
%绘图
subplot(211),plot(abs(Hf_zero)),axis([-50 2100 0 35])
subplot(212),plot(1:990,-unwrap(angle(Hf_zero(1:990)))),hold on
plot(1060:2048,-unwrap(angle(Hf_zero(1060:2048)))),axis([-50 2100,-2000 600])
1.
plot(1:990,-unwrap(angle(Hf_zero(1:990)))),hold on
plot(1060:2048,-unwrap(angle(Hf_zero(1060:2048)))),axis([-50 2100,-2000 600])
因为在样本点990-1060之间的匹配滤波能量接近于0,相当于过采样带来的复制信号中的缺失频率