clc
clear?all
fs=40e9;
fc=7.5e9;
B=3e9;
N=1024;
M=512;
Z=512;

c=3e08;
deta_d=0.01;
tc=gauspuls（‘cutoff‘fcB/fc[]-40）;
t=-tc:1/fs:tc;
ps=gauspuls（tfcB/fc）*1000;
s=zeros（MN）;
ss=zeros（MN）;
zb=[100600
????200600
????300600];
d=0.3;
mb=3;
for?i=1:1:mb
????m0=zb（i1）;
????n0=zb（i2）;
for?m=1:1:M
????r=sqrt（（（m0-m）*deta_d）^2+（（n0-1）/fs*c/2）^2）-d;
????n1=floor（r*2/c*fs+1）;
????n2=floor（sqrt（（r+d）^2-d^2）*2/c*fs+1）;
????s2=zeros（1N+62）;
????for?n=n1:1:n2
????if?n>N&atan（abs（（m0-m）*deta_d/（（n0-1）/fs*c/2）））>10/180*pi;
????????ss（m:）=zeros（1N）;
????else
????????
????????s2（n+16:n+46