LFM信号模糊函数Matlab代码

%% 窄带模糊函数
% 任意信号NBAF
%%
clear all;clc;
close all
C=1500;
%%
FS=5e3;
WL=2;  
A=1;
B = 50;          % 带宽
T = 0.5;         % 信号脉宽
f0 = 1000;       % 信号的初始频率
K = B/T;         % LFM信号的调制系数
t = -T/2:1/FS:T/2 ; 
x = A*cos（2*pi*（f0）*t+pi*K*t.^2）;
data1=[zeros（1（WL/2-T/2）*FS） x zeros（1（WL/2-T/2）*FS）];
% data1=awgn（data1-5）;
% figure（1）;plot（（1:length（data1））/FSdata1）;title（‘data1‘）;hold on

t=-WL/2:1/FS:WL/2; 
fd=-20:0.2:20; 
for i=1:length（fd） 
    E=exp（-j*2*pi*fd（i）.*t）;
    Y（i:）=data1.*E;
end

p=length（data1）*2;
for j=1:length（fd） 
    tmp=fft（xcorr（Y（j:）data1）p）;
    % envelope
    h = [1; 2*ones（fix（（p-1）/2）1）; ones（1-rem（p2）1）; zeros（fix（（p-1）/2）1）];
    M（j:） = tmp（:）.*h;
    M（j:） = abs（ifft（M（j:）p））‘;             
    [Mmax（j）]=max（M（j:））/p*2;
end
M=M/max（max（M））;
t=-WL:1/FS:WL; 
figure;imagesc（tfdM）; colorbar;title（‘信号NBAF‘）xlabel（‘时间 s‘）ylabel（‘频移 Hz‘）;
[ttffdd]=meshgrid（tfd）;
% figure;mesh（ttffddM（1:length（fd）1:length（t）））;t