tu1-20.m

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

clear;
DOA=[pi/180 pi/36 pi/18 pi/9 pi/6];
%DOA=pi/6;
dlabuda=0:0.1:50;
for s=1:length(DOA)
    p=1./(DOA(s));
    %p=1./DOA;
    for k=1:length(dlabuda)
        xiangguan(k)=quadl(@beijifun,-DOA(s)/2,DOA(s)/2,[],[],dlabuda,k,p);
    end
    plot(dlabuda,real(xiangguan));
    axis([0,50,-0.4,1.0]);
    hold on;
    grid on;
    drawnow;
end
%xlabel('dlabuda');
%ylabel('相关系数');
%legend('DOA=1度','DOA=5度','DOA=10度','DOA=20度','DOA=30度');