《照虎画猫成狗系列》之三，线性信号的近似

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

感觉有点问题： 图上表示的是近似信号的频谱，但实际的代码依然是原来的信号。

先不管，反正自己目前没能力写出来近似信号的表达，留着以后吧,按大神的路走。

本节主要是 POST近似与原信号的比较，

也就是说通过分析说明：POST 近似在一定条件下完全可以代替原信号

 

% 这是一个线性调频脉冲的相位与频率
% 因为线性调频是难以直接精确推导，所以用驻定相位原理进行简化。
% 线性调频信号的频域性质有更清晰的认识


clear;clc;close all;

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

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;   % 时间轴的另一个求法,
tt = linspace(-0.5,0.5,N);

%按理说此处应该是近似信号
st = exp(1i*pi*K*t.^2);
Sf = fftshift(fft(fftshift(st))); % FFT

figure(1) ;
subplot(221);
plot(tt, real(Sf))
axis([-0.1 0.1 -inf inf]); %只看[-0.1 0.1]部分
title('(a)频谱实部'),ylabel('幅度');

%pi*K*t.^2
subplot(2,2,2),plot(tt,abs(Sf));
title('(b)频谱幅度'),ylabel('幅度');


subplot(2,2,3),plot(tt,imag(Sf));axis([-0.1 0.1 -inf inf]);
title('(c)频谱虚部'),xlabel('频率（归一化后）'),ylabel('幅度');

subplot(2,2,4),plot(tt,unwrap(angle(Sf)));
title('(d)频谱相位'),xlabel('频率（归一化后）'),ylabel('相位（弧度）');

重要：

1.  与原信号相比，有相似的线性结构。不同的是PI/4处的相变与斜率符号的改变

2. 包络近似一致。 什么什么源自1.25倍的采样率。（看不懂，估计是翻译的事）

3. 频域相位与时域相位相同，基本是二次的。这表时f=Kt 就是可以用时间直接算出频率。

4. 当TBP大于100时，POST是非常精确的近似。