正负斜率LFM信号的仿真（信号生成，模糊函数矩阵，波形隔离度）

%% 正负斜率LFM信号在同时极化雷达中的应用
%% 初始化
clc
clear
B=2e7;
taup=1e-6;
tbp= B*taup;
u=2*pi*B/taup;
npoints = 5 * B * taup + 1; %可调整时宽带宽积来调整抽样频率，数值越大，抽样频率越高，频谱恢复得越好
delt = linspace（-taup/2. taup/2. npoints）; %时间轴
Ichannal = cos（u .* delt.^2 / 2.）; %信号的实部
Qchannal = sin（u .* delt.^2 / 2.）; %信号的虚部
LFMh=（1/sqrt（taup））*（Ichannal + sqrt（-1） .*Qchannal）;
LFMv=（1/sqrt（taup））*（Ichannal - sqrt（-1） .*Qchannal）;%信号的复解析形式
LFMhFFT = fftshift（fft（LFMh））;
LFMvFFT = fftshift（fft（LFMv））;
sampling_interval = taup / npoints;
freqlimit = 0.5 / sampling_interval;
freq = linspace（-freqlimitfreqlimitnpoints）;
%% 信号的幅度谱
% figure（1）;subplot（311）;
% plot（freqabs（LFMhFFT）/max（abs（LFMhFFT））‘r‘）;
% axis（[-4e7 4e7 0 1]）;
% subplot（312）;
% title（‘spectrum for LFM signal‘）;
% plot（freqabs（LFMvFFT）/max（abs（LFMvFFT））‘b‘）;
% axis（[-4e7 4e7 0 1]）;
% subplot（313）;
% plot（deltIchannal‘k‘）;
% axis（[-taup/2 taup/2 -1 1]）;
% title（‘real part of LFM signal‘）;
% 两路正交信号具有同样的幅度谱
%% 正负斜率LFM信号模糊函数矩阵
% eps=1e-8;
% [delt1freq1]=meshgrid（deltfreq）;
% syms t
% f1=（1/taup）*exp（-sqrt（-1）*2*pi*u*t^2）;
% f2=（1/taup）*exp（sqrt（-1）*2*pi*u*t^2）;
% karma1=zeros（npointsnpoints）;
% karma2=zeros（npointsnpoints）;
% for i=1:length（npoints）
%     for j=1:length（npoints）
%         karma1（ij）=int（f1t[0.5*（-taup+abs（delt（j））-freq（i）/u） 0.5*...
%             （taup-abs（delt（j））-freq（i）/u）]）;  
%         karma2（ij）=int（f2t[0.5*（-taup+abs（delt（j））+freq（i）/u） 0.5*...
%             （taup-abs（delt（j））+freq（i）/u）]）;
%     end
% end
% Ahh=exp（-sqrt（-1）*pi*u*delt1.^2）.*exp（sqrt（-1）*2*pi*（freq1+u*delt1）.*delt1）.*...
%     sin（pi*（freq1+u*delt1）.*（taup-abs（delt1））+eps）./（pi*（freq1+u*delt1）*taup+eps）; %自模糊函数
% Avv=exp（sqrt（-1）*pi*u*delt1.^2）.*exp（sqrt（-1）*2*pi*（freq1-u*del