雷达模糊函数matlab代码

clear all;

f0=10e9;
B=300e6;
fs=f0*10;
Tt=1e-7;
tt=linspace（0TtTt*fs）;
k=B/Tt;
u=1*（exp（1j*pi*（2*（f0-B/2）*tt+k*tt.^2）））;%postive chirp

M=length（u）;

F=8;
K=100;
df=F/K/M;
T=1;
N=100;
uamp=abs（u）;
phas=angle（u）;
t=0:M-1;
tscale1=[0 0:M-1 M-1];
dphas=[NaN mod（diff（phas）2*pi）]/2/pi;

% plot the signal parameters
figure（1） clf hold off 
subplot（311）
plot（tscale1[0 abs（uamp） 0]‘linewidth‘1.5）
ylabel（‘ Amplitude ‘）
axis（[-inf inf 0 1.2*max（abs（uamp））]）

subplot（312）
plot（t phas‘linewidth‘1.5）
axis（[-inf inf -inf inf]）
ylabel（‘ Phase [rad] ‘）

subplot（313）
plot（tdphas*ceil（max（t））‘linewidth‘1.5）
axis（[-inf inf -inf inf]）
xlabel（‘ \itt / t_b ‘）
ylabel（‘ \itf * Mt_b ‘）

dtau=ceil（T*M）/N;
tau=round（[0:1:N]*dtau）;
f=[0:1:K]*df; 
f=[-fliplr（f） f];
mat1=spdiags（u‘0M+ceil（T*M）M）;
u_padded=[zeros（1ceil（T*M））uzeros（1ceil（T*M））];
cidx=[1:M+ceil（T*M）];
ridx=round（tau）‘;
index = cidx（ones（N+11）:） + ridx（:ones（1M+ceil（T*M）））;
% mat2 = sparse（u_padded（index））; 
mat2 = （u_padded（index））; 
uu_pos=mat2*mat1;
e=exp（-j*2*pi*f‘*t）;
a_pos=