tu1-21.m

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

clear;
DOA=[pi/180 pi/36 pi/18 pi/9 pi/6];
dlamda=0:0.1:50;
k1=1./(2.*sqrt(3));
for s=1:length(DOA)
    sigma=DOA(s).*k1;
    for n=1:length(dlamda)
        xiangguan(n)=quadl(@beijifun1,-pi,pi,[],[],dlamda,sigma,n);
    end
    plot(dlamda,real(xiangguan));
    axis([0,50,-0.2,1.2]);
    hold on;
    grid on;
    drawnow;
end
xlabel('d/λ');
ylabel('θ 相关系数');
legend('σ=1°k','σ=5°k','σ=10°k','σ=20°k','σ=30°k');