clear all;  
c=3*10^8;
f=10e9;
v=c/f;
d=v/2;
L=8;
l=0:1:L-1; 
angle=-3*pi/5:0.01:3*pi/5;
theta1=-5*pi/180;
a1=exp(-j*2*pi*l*sin(theta1)*d/v);
b1=exp(-j*2*pi*l.'*sin(angle)*d/v);
A1=conj(a1)*b1;
m=max(abs(A1));
log_A1=20*log10(abs(A1)/m);
figure(1);
plot(angle*180/pi,log_A1,'b');
grid on;
title('兩個對稱的波束');
xlabel('角度/度');
ylabel('增益/dB');
hold on;                                                           %再畫一個
theta2=5*pi/180;
a2=exp(-j*2*pi*l*sin(theta2)*d/v);
b2=(exp(-j*2*pi*l.'*sin(angle)*d/v));
A2=conj(a2)*b2;
m=max(abs(A2));
log_A2=20*log10(abs(A2)/m);
plot(angle*180/pi,log_A2,'r');
legend('-5度','5度');
%極坐標表示
figure(2);
polar(angle,abs(A1),'--b');                          %極坐標下的表示不需要使用對數(shù)的計算！
hold on;
polar(angle,abs(A2),'--r');   
title('在極坐標下的差波束天線方向圖');
%直角坐標
figure(3);
plot(angle*180/pi,abs(A2)-abs(A1),'r');
grid on;
title('在直角坐標下的差波束天線方向圖');
%差波束目標跟蹤角度
theta=-7.2*pi/180;
a1=exp(-j*2*pi*l*sin(theta1)*d/v);
b1=exp(-j*2*pi*l.'*sin(theta)*d/v);
A1=abs(conj(a1)*b1);
a2=exp(-j*2*pi*l*sin(theta2)*d/v);
b2=(exp(-j*2*pi*l.'*sin(theta)*d/v));
A2=abs(conj(a2)*b2);
delta=A2-A1;
while(abs(delta)>=0.2)
    if(delta>0)
        theta1=theta1+0.1*pi/180;
        theta2=theta2+0.1*pi/180;
        a1=exp(-j*2*pi*l*sin(theta1)*d/v);
        b1=exp(-j*2*pi*l.'*sin(theta)*d/v);
        A1=abs(conj(a1)*b1);
        a2=exp(-j*2*pi*l*sin(theta2)*d/v);
        b2=(exp(-j*2*pi*l.'*sin(theta)*d/v));
        A2=abs(conj(a2)*b2);
        delta=A2-A1;
    else if(delta<0)
        theta1=theta1-0.1*pi/180;
        theta2=theta2-0.1*pi/180;
        a1=exp(-j*2*pi*l*sin(theta1)*d/v);
        b1=exp(-j*2*pi*l.'*sin(theta)*d/v);
        A1=abs(conj(a1)*b1);
        a2=exp(-j*2*pi*l*sin(theta2)*d/v);
        b2=(exp(-j*2*pi*l.'*sin(theta)*d/v));
        A2=abs(conj(a2)*b2);
        delta=A2-A1;
        end
    end
end
%重新畫出波束調(diào)整后的差波束方向圖
%極坐標表示
a1=exp(-j*2*pi*l*sin(theta1)*d/v);
b1=exp(-j*2*pi*l.'*sin(angle)*d/v);
A1=abs(conj(a1)*b1);
a2=exp(-j*2*pi*l*sin(theta2)*d/v);
b2=(exp(-j*2*pi*l.'*sin(angle)*d/v));
A2=abs(conj(a2)*b2);
figure(4);
polar(angle,abs(A1),'--b');                          %極坐標下的表示不需要使用對數(shù)的計算！
hold on;
polar(angle,abs(A2),'--r');   
title('在極坐標下的差波束天線方向圖');
%直角坐標
figure(5);
plot(angle*180/pi,abs(A2)-abs(A1),'r');
grid on;
title('在直角坐標下的差波束天線方向圖(跟蹤完成)');