【脉冲压缩】线性调频+步进频和相位编码信号脉冲压缩仿真【含Matlab源码 2753期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【脉冲压缩】基于matlab线性调频+步进频和相位编码信号脉冲压缩仿真【含Matlab源码 2753期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab信号处理（初级版）

备注：
点击上面蓝色字体付费专栏Matlab信号处理（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab信号处理（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、雷达脉冲压缩简介

1 脉冲压缩的定义
脉冲压缩即pulse compression，它是指发射宽编码脉冲并对回波进行处理以获得窄脉冲，因此脉冲压缩雷达既保持了窄脉冲的高距离分辨力，又能获得宽脉冲的强检测能力。

2 脉冲压缩的主要手段
目前的脉冲压缩的手段主要有线性调频、非线性调频与相位编码等。
1)线性调频
是最简单的脉冲压缩信号，容易产生，而且其压缩脉冲形状和信噪比对多普勒频移不敏感，因而得到了广泛的应用，但是，在利用多普勒频率测量目标方位和距离的情况下很少使用；
2)非线性调频
非线性调频具有几个明显的优点，不需要对时间和频率加权，但是系统复杂。为了达到所需的旁瓣电平，需要对每个幅度频谱分别进行调频设计，因而在实际中很少应用；
3)相位编码
相位编码波形不同于调频波形，它将宽脉冲分为许多短的子脉冲。这些子脉冲宽度相等，其相位通过编码后被发射。根据所选编码的类型，包括巴克码、伪随机序列编码以及多项制编码等。

3 脉冲压缩的产生背景
随着飞行技术的飞速发展，对雷达的作用距离、分辨能力、测量精度和单值性等性能指标提出越来越高的要求。测距精度和距离分辨力对信号形式的要求是一致的，主要取决于信号的频率结构，为了提高测距精度和距离分辨力，要求信号具有大的带宽。而测速精度和速度分辨力则取决于信号的时域结构，为了提高测速精度和速度分辨力，要求信号具有大的时宽。除此之外，为提高雷达系统的发现能力，要求信号具有大的能量。由此可见，为了提高雷达系统的发现能力、测量精度和分辨能力，要求雷达信号具有大的时宽、带宽、能量乘积。但是，在系统的发射和馈电设备峰值功率受限制的情况下，大的信号能量只能靠加大信号的时宽来得到。测距精度和距离分辨力同测速精度和速度分辨力以及作用距离之间存在着不可调和的矛盾。于是在匹配滤波器理论指导下，人们提出了脉冲压缩的概念。
由于发射机效率的限制，雷达真正采用的脉压信号是由调频和相位编码产生的，其中以线性调频和二相编码信号的研究与应用最为广泛。

4 线性调频信号（LFM）
脉冲压缩雷达最常见的调制信号是线性调频（Linear Frequency Modulation）信号,接收时采用匹配滤波器（Matched Filter）压缩脉冲。
LFM信号的数学表达式，式中为载波频率，为矩形信号。是调频斜率，于是，信号的瞬时频率为。
图 典型的LFM信号（a）up-LFM(K>0)（b）down-LFM(K<0)
将(1)式中的up-LFM(信号重写 式中，是信号s(t)的复包络。由傅立叶变换性质，S(t)与s(t)具有相同的幅频特性，只是中心频率不同而以，因此，Matlab仿真时，只需考虑S(t)。由Matlab程序产生信号，并作出其时域波形和幅频特性。
由此可以得到关于线性调频脉冲信号频率特性的两个重要结论:
(1)在满足大时宽带宽积的条件下, 线性调频脉冲信号的振幅谱接近矩形函数, 频谱宽度近似等于信号的调频变化范围B, 与时宽无关。
(2)在满足大时宽带宽积的条件下, 线性调频脉冲信的相位谱具有平方律特性。
以上两点是设计匹配滤波器，进行脉冲压缩处理的主要依据。
注：匹配滤波器是在白噪声背景中检测信号的最佳线性滤波器，其输出信噪比在某个时刻可以达到最大，它是对线性调频信号进行脉冲压缩的主要手段。

5 对脉冲压缩信号的仿真
结合以上分析，用Matlab仿真雷达发射信号，回波信号，和压缩后的信号的复包络特性，其载频不予考虑（实际中需加调制和正交解调环节），仿真信号与系统模型如图3。
仿真程序模拟产生理想点目标的回波，并采用频域相关方法（以便利用FFT）实现脉冲压缩，得到仿真图。
对比图可以看出，采用脉冲压缩后，在接收机的接收端信噪比有了明显的改善，进而可以在满足分辨率的基础上，提高雷达的作用距离。

⛄三、部分源代码

% function LFM_disp
close all;clear; clc
c = 31e8; %%光速 m/s
Bw = 51e6; %%雷达波形 信号带宽
Tau = 501e-6; %%雷达波形 脉冲宽度
f0 = 101e9; %%载波频率
deta_t = 1/(Bw*16); %%保证采样率

time = 0:deta_t:Tau;
N_t = length(time);
time = (0:1:(N_t -1))deta_t ;%%重新调整一下时间采样序列
range = ctime/2; %%距离采样序列
range_cell = mean(diff(range));
f0 = -Bw/2;gamma = Bw/Tau

%%原始信号
x = exp(i2pi*(f0time+1/2gammatime.^2));
x_fft = ifftshift( fft(x)) ;f = linspace(-1/21/deta_t, 1/21/deta_t,N_t);
x_fft_amp = 20log10(abs(x_fft));x_fft_amp_max = max(x_fft_amp ); x_fft_amp = x_fft_amp - x_fft_amp_max ;
figure; subplot(2,1,1);plot(range , real(x),‘k’);hold on; plot(range, imag(x),‘r-.’);xlabel(‘range / m’)
legend(‘I’, ‘Q’);xlabel(‘range /m’);title(‘before range compress’);
subplot(2,1,2);plot(f/1e6,x_fft_amp );xlabel(‘freq / MHz’);ylim([-15,5]);title(‘bandwidth disp’)

%%脉冲压缩 (FFT 实现)
P = nextpow2(2*N_t);
x_append = [x, zeros(1,2^P - N_t)]; fft_x_append = fft(x_append);
x_match_append = [conj(fliplr(x)),zeros(1,2^P-N_t)]; fft_x_match_append = fft(x_match_append );
fft_x_pc = fft_x_append .*fft_x_match_append;
x_pc = ifft(fft_x_pc);
range_pc = (0:1:2^P-1) *range_cell;
figure;plot(range_pc,abs(x_pc));xlabel(‘range /m’);title(‘after range compress’)

⛄四、运行结果

⛄五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社，2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社，2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社，2018.
[4]李树锋.基于完全互补序列的MIMO雷达与5G MIMO通信[M].清华大学出版社.2021
[5]何友,关键.雷达目标检测与恒虚警处理（第二版）[M].清华大学出版社.2011

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除