%%?窄帶模糊函數
%?任意信號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