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;                 %fumaichong kuandu
delta_f=0.5             %pinlv bujin liang
Tr=5;                 %maichong 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.*rect

Z=fudu.*rect;
surf(X,Y,Z)