sitafun.m

我自己在学习空时编码一书时编写的第一章的仿真程序

function Sita=sitafun(x,lamda,alfa)
for m=1:length(x)
    for k=1:length(x)
        da=alfa./2./10000;
        a=-alfa./2:da:alfa./2;
        d_mk=abs(x(m)-x(k));
        f=exp(-j*2*pi*d_mk./lamda.*sin(a))./alfa.*da;
        Sita(m,k)=real(sum(f));%int eq(1-103)
    end
end