MATLAB 后向投影成像算法

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