progam1.m

来自「总共有三个源代码： progam1:波达方向估计DOA：三个接收器对波达方向的」· M 代码 · 共 29 行

%
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一个明显的最小值')