progam1.m

总共有三个源代码： progam1:波达方向估计DOA：三个接收器对波达方向的联合估计问题 progam2:利用电子旋转求DOA progam3:波达方向矩阵法对水听器的DOA进行估计 超好

%
clear
rot=2;
for i=1:3
    rot=rot-1;
    theta=0:180;
    theta=theta*pi/180;
    r=(sin(theta)).^0.25;
    phi=theta+rot*pi/4;
    X(:,i)=phi';Y(:,i)=r';
end
figure
polar(X,Y);grid;
A=ones(271,3);A=A*0;
A(1:181,1)=Y(1:181,1);
A(46:226,2)=Y(1:181,2);
A(91:271,3)=Y(1:181,3);
for i=1:271
    s=3*i-2;e=s+2;
    x(s:e)=A(i,:);
end
test=A(70,:);test=test+0.01*randn(size(test));
for i=1:271
    ek=A(i,:)-test;
    Jk(i)=log(ek*ek');
end
figure
plot(Jk);grid
text(70,-11.5,'\leftarrow一个明显的最小值')