《合成孔径雷达成像——算法与实现》之【5】仿真图3.6_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);
ht = conj(fliplr(st)); % 匹配滤波器
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 -1 1]);
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,angle(out));axis([-1 1 -5 5]);
title('(d)压缩后信号的相位(经扩展)'),xlabel('相对于t_{0}时间(\mus)'),ylabel('相位(弧度)');
重要的东西:
在 ifft 之后,虚部全为 0 ,所以对偶数旁瓣相位是0,偶对旁瓣相位为pi(这句没看懂)
原程序有个ABS,但我感觉应该没有,就去掉了,但跟书上的图不一样了!!
感觉我不对的,去参考大神的代码