《照虎画猫成狗系列》之 五，回波信号（三）有频移压缩

《合成孔径雷达成像——算法与实现》之【8】仿真图3.8_destiny0321的博客-CSDN博客

本节目地，在无噪信号中加入频移，看一下情况

% 这是一个线性调频脉冲的相位与频率

clear;clc;close all;

T = 7.24e-6;    %持续时间为 7.24 us 
B = 5.8e6;      % 带宽 5.8MHz

K = B/T;        % 调频率？？ 这是什么意思？
ratio = 5;      % 过采样倍数： 比如你是 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;   % 时间轴的另一个求法,

%st = exp(1i*pi*K*t.^2)+0.75*randn(1,N);      % 生成带有高斯噪声的信号
st = exp(1i*pi*K*(t-1e-6).^2);     % 生成信号-频率移动了一点（相当于时间移动了一点）

ht = conj(fliplr(st));      % 匹配滤波器-
%ht = exp(-1i*pi*K*t.^2);     % 匹配滤波器-？？？为什么这个不能用
out = fftshift(ifft(fft(st).*fft(ht)));  % 压缩

%下面计算幅值
Z = abs(out);
Z = Z/max(Z);
Z = 20*log10(Z);   %这是扩展，这时暂时去掉


figure(1) ;
subplot(221);
plot(t*1e6, real(st))
title('信号实部')

subplot(2,2,2),plot(t*1e6,Z);axis([-1 1 -30 0]);
title('(c)压缩后的信号'),ylabel('幅度(dB)');

subplot(2,2,3),plot(t*1e6,real(out) );
title('(b)压缩后的信号'),xlabel('相对于t_{0}时间(\mus)'),ylabel('幅度');

subplot(2,2,4),plot(t*1e6,abs(angle(out)));axis([-1 1 -5 5]);
title('(d)压缩后信号的相位(经扩展)'),xlabel('相对于t_{0}时间(\mus)'),ylabel('相位(弧度)');

subplot(2,2,4),plot(t*1e6,angle(out));axis([-1 1 -5 5]);
title('(d)压缩后信号的相位(经扩展)'),xlabel('相对于t_{0}时间(\mus)'),ylabel('相位(弧度)');

重要：

为什么0频落后，压缩信号的脉冲峰值，却提前了？

我改了时间偏移  st = exp(1i*pi*K*(t+1e-6).^2)， 发现还是提前了....没想通！！