tu1-21-1.m

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

alfa_r=[1,5,10,20,30];
alfa_r=alfa_r./180.*pi;
da=0.001;
for n=1:length(alfa_r)
    a=-pi:da:pi;%for fig 1.21
    %a=-alfa_r(n)./2 : da: alfa_r(n)./2;%for fig 1.20
    k=1./(2.*sqrt(3));
    sigma=alfa_r(n).*k;
    p_a=1./sqrt(2*pi*sigma.^2).*exp(-a.^2./(2.*sigma.^2));%for fig 1.21
    %p_a=1./alfa_r(n);%for fig 1.20
    dmklamda=0:0.1:50;% that is dmk/lamda
    sitamk(length(dmklamda))=0;%init
    for s=1:length(dmklamda)
        sitamk(s)=sum(exp(-j.*2.*pi.*dmklamda(s).*sin(a)).*p_a.*da);%integral
    end
    plot(dmklamda,real(sitamk));
    axis([0 50 -0.2 1.2]);
    hold on;
    grid on;
    drawnow;
end