常规、LFM、相位编码、频率编码等雷达调制信号代码（MATLAB）

介绍

有关傅里叶变化介绍的优秀博文：点这里真的超级有惊喜
雷达信号调制种类：

常规脉冲信号

参考博文：
博文1
博文2

优秀博文：主要参考
matlab代码：
加高斯白噪声：awgn(x,SNR)

clc
clear all
close all
fs = 300;%采样率
T = 1/fs;
N = 600;
t = (0:N-1)*T; % 时间序列

A_danpin = 10;
f_danpin = 80;
y_danpin = A_danpin*cos(2*pi*f_danpin*t)
%y_danpin_noise = y_danpin + randn(size(t));     %添加噪声
y_danpin_noise = awgn(y_danpin ,2);%添加高斯白噪声
figure(1);
subplot(2,1,1)
%plot(t,y_danpin_noise)
plot(t,y_danpin_noise);
xlabel('t/s')
ylabel('幅度 / V')

n = 2^nextpow2(N);%采样点数，采样点数越大，分辨的频率越精确，N>=L，超出的部分信号补为0
Y = fft(y_danpin_noise,n)/n*2;%除以N乘以2才是真实幅值，N越大，幅值精度越高
f = fs/n*(0:1:n-1); %频率
P = abs(Y);%幅值
subplot(2,1,2)
plot(f(1:n/2),P(1:n/2));%函数fft返回值的数据结构具有对称性,因此我们只取前一半
xlabel('频率(Hz)')
ylabel('幅值 ')
title('幅值频谱')
grid on;

LFM

优秀博文：点这里有惊喜
主要参考：点这里有惊喜

clc
clear all
close all
T = 1;%采样时间
f = 1e5;%采样率
t = 0:1/f:(T-1/f); % 采样点
n = length(t); % 采样点数
% 线性调频信号 s(t)=a(t)cos[2πf0 t+πkt^2]，a(t)是包络，f0是调频频率
A_lfm = 10;
f_lfm = 1000;
k_lfm = 800000;
y_lfm = A_lfm*cos(2*pi*f_lfm*t+pi*k_lfm*t.^2);
figure(3)
plot(t(1:500),y_lfm(1:500));
xlabel('t/s');
title('线性调频信号波形');
axis([0 0.005 -A_lfm A_lfm]);

注：要使得时域图扩大范围，只需要改变采样点数，比如 t (1:500)【采500个点】可以改成t（1:1000）【采1000个点】。

LFM频谱：

clc
clear all
close all
B = 2e8;
T = 1e-6;%采样时间
fs = 4*B;%采样率
K = B/T;
N =  round( T / (1/fs) );%采样点数
t = linspace(