function[amf]=af_sfw(Tp,Tr,N,delta_f)

delta_tau=Tp/10; % 时间轴间隔

taumin=-(Tp+Tr); % 时间轴最小值

delta_fd=delta_f/100; % 频率轴间隔

B=N*delta_f;

tau=taumin:delta_tau:-taumin; % 用于绘图的时间轴

fd=-B/2:delta_fd:B/2; % 用于绘图的频率轴

eps=0.0000001;

nf=length(fd);

nt=length(tau);

amf=zeros(nf,nt);

....

....

amf=abs(amf)/max(max(abs(amf)));

[m,n]=find(amf==1.0);

figure(1);mesh(tau*1e6,fd*1e-6,amf);

xlabel('Delay - seconds');

ylabel('Doppler - Hz');

zlabel('Ambiguity function');

title('步进频率脉冲信号的模糊函数');

figure(2);contour(tau*1e6,fd*1e-6,amf,1,'b');

xlabel('Delay - seconds');

ylabel('Doppler - Hz');

title('步进频率脉冲信号的模糊度图');

figure(3);plot(tau*1e6,amf(m,:));grid on;

xlabel('Delay - seconds');

ylabel('Ambiguity function');

title('步进频率脉冲信号的距离模糊函数主瓣');

figure(4);plot(fd*1e-6,amf(:,n));grid on;

xlabel('Doppler - Hz');

ylabel('Ambiguity function');

title('步进频率脉冲信号的多普勒模糊函数主瓣');