直接序列扩频信号的keystone变换捕获

对直扩扩频信号用两种方法：1，sinc函数采样,2，Chirp_ZT变换，分别进行Keystone变换。随后基于FFT实现扩频信号的捕获。

% keystone变换，根据理想的频域距离模型进行计算验证。未产生信号收发流程。
clear all;
close all;

N = 106; % 快时间采样点数
M = 256;  % 慢时间采样点数
L = N;
K = M;

B = 100e6; % 信号带宽
c = 3e8; % 光速
fc = 1e9; % 中心频率

eta = B/（fc*L）;

S_lm = zeros（LM）; % 原始距离向频谱
S_lk = zeros（LM）; % 中间变量
S_lmp = zeros（LM）; % keystone变换后的距离向频谱

l = （-L/2:L/2-1）‘;
m = （-M/2:M/2-1）‘;
% m = （0:M-1）‘;

T = 25e-6; % 脉冲重复周期


vi = [1000]; % 目标速度
R0 = [20]; % 初始位置
fdci = 2*fc*vi/c;
bi = 2*R0/c * B/L;

% 产生目标回波距离向频谱
for ll=-L/2:L/2-1
   for mm=-M/2:M/2-1
       S_lm（ll+1+L/2mm+1+M/2） = sum（exp（-1i*2*pi.*bi*（ll）） .* exp（1i*2*pi*（1+eta*（ll））.*fdci*（mm）*T））+  randn;
       S_lm（ll+1+L/2mm+1+M/2） = sum（exp（-1i*2*pi.*bi*（ll）） .* exp（1i*2*pi*（1+eta*（ll））.*fdci*（mm）*T））;
   end
end



% 求 S（lk‘）
for ll=-L/2:L/2-1
    for kp=-M/2:M/2-1
        S_lk（ll+1+L/2kp+1+M/2） = S_lm（ll+1+L/2:）*exp（-1i*2*pi*（1+eta*ll）/M*（kp）*m）; % 
    end
end

% 求 S（lm‘）
% S_lmp = fft（S_lk[]2）;

for ll=-L/2:L/2-1
    for mp=-M/2:M/2-1
        kp = m;
        S_lmp（ll+1+L/2mp+1+M/2） = S_lk（ll+1+L/2:）*exp（1i*2*pi/M*kp*（mp））;
    end
end

varef = M;
% IFFT求距离像
s_lm = ifft（S_lm）;  % fftshift
s_lmp = ifft（S_lmp）;
fig_position = [1050450310];    
set（gcf‘position‘fig_position）;
fig_fontsize = 12;  % 字号

figure;
contour（abs（s_lm））  % 原始信号
hold on;
figure;
contour（abs（s_lmp）） % keystone变换后的信号

set（gcf‘Color‘‘none‘）  % 背景色
set（gca‘LineWidth‘0.5）; % 坐标轴线宽 该宽度影响gird线宽
set（gca‘FontSize‘fig_fontsize‘FontName‘‘Times New Roman‘）;
% h=legend（‘\itK\rm = 1‘‘\itK\rm = 2‘‘\itK\rm = 4‘‘\itK\rm = 8‘‘Location‘‘Best‘‘Orientation‘‘vertical‘）; % 方向：vertical 或 horizontal
% set（h ‘Box‘ ‘off‘）    % 或者 legend（gca‘boxoff‘）  或者  legend（‘boxoff‘）
% set（h ‘Color‘ ‘none‘） % 图注背景透明
xlabel（‘\fontname{Times New Roman}\fontsize{12}雷达发射脉冲数 \itM‘）;
ylabel（‘\fontname{Times New Roman}\fontsize{12}距离分辨单元‘）;
set（gca‘MinorGridLinestyle‘‘--‘）;
set（gca‘GridLinestyle‘‘--‘）;  % 网格线型
% set（gca‘XGrid‘‘off‘）; % x轴网格
% set（gca‘YGrid‘‘off‘）; % y轴网格
set（gca‘XMinorGrid‘‘off‘）; % x轴子网格
set（gca‘YMinorGrid‘‘off‘）; % y轴子网格
set（gca‘XMinorTick‘‘on‘）; % x轴子网格标识
set（gca‘YMinorTick‘‘on‘）; % y轴子网格标识


s_lmp_FFT = fft（s_lmp[]2）;
figure;
contour（abs（s_lmp_FFT））
[vaktvidkt]=max（max（abs（s_lmp_FFT/M）））;

s_lm_FFT = fft（s_lm[]2）;
figure;
contour（abs（s_lm_FFT））
[vafftvidfft]=max（max（abs（s_lm_FFT）））;
figure; 
imagesc（abs（S_lk））
figure; 
imagesc（angle（S_lmp））

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        632  2018-07-11 07:33  gold_gener.m

     文件       3437  2018-09-17 15:03  keystoneSinc.m

     文件       3040  2018-07-27 20:21  Kstone10MCXFultSpectrum.m

     文件       2673  2018-08-06 16:02  ch4_keystone_contour.m

     文件        363  2018-07-26 15:02  gold_gen.m

----------- ---------  ---------- -----  ----

                10145                    5