close all;clear all;clc;
x=-5:0.1:5;          %shijian
y=-5:0.1:5;          %pinlv
[X,Y]=meshgrid(x,y);
%%%%%%%canshu sheji%%%%%%%%%%%%%
N=10;                %zimaichong geshu
T=5;                 %zimaichong kuandu
delta_f=0.5             %pinlv bujin liang
Tr=10;                 %maijian kuandu
%%%%%%%kai_rect%%%%%%%%%%%%%%%%%
%a=exp(j*pi*Y.*(T-X));
%b=pi*Y.*(T-abs(X))+eps;
%c=((sin(b).*(T-abs(X)))./(T*b));
%rect=a.*c;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
xiangwei=exp(j*pi*((N-1-p)*(Y-p*delta_f)*Tr-(N-1+p)*delta_f*X));
d=pi*((Y-p*delta_f)*Tr-delta_f*X);
f=sin((N-abs(p))*d);
fudu=f./d;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Z=0;
for p=-(N-1):(N-1)
    Z=Z+xiangwei.*fudu.*kai_rect(X-p*Tr,Y-p*delta_f,T)
end

surf(X,Y,Z)