doa_phase.m

干涉仪测向具有精度高、速度快的特点

 %%该程序仿真利用干涉技术测向的方法
 clear all;
 close all;
 clc;
 result=zeros(1,100);
 
 for Loop=1:500
 for loop=1:100
 d1=0.58;%%
 d2=(3+0.0015)*d1;%%
 k=floor(d2/d1);%%基线的整数倍
 deltak=d2/d1-k;%%基线的小数部分
 f=1e10;%%信号频率
 C=3e8;%光速
 lambda=C/f;%信号波长
 sita0=30;%%入射角
 sita=((sita0)/180)*pi;
 phi1=2*pi*d1*sin(sita)/lambda;
 phi2=2*pi*d2*sin(sita)/lambda;
 delta_phi10=mod(phi1,2*pi);
 delta_phi20=mod(phi2,2*pi);
sig1=exp(j*(phi1+0.0001*loop*rand(1)));%%给出相位误差
sig2=exp(j*(phi2+0.0001*loop*rand(1)));
delta_phi1=mod(angle(sig1),2*pi);
delta_phi2=mod(angle(sig2),2*pi);
c=floor(delta_phi1*(d2/d1)/(2*pi));
% m=n*k+c;
phi=delta_phi2+2*pi*c-delta_phi1*(d2/d1);
if abs(phi)<=1e-5
    phi=0;
end
 if phi>=0
     n=floor((phi/(2*pi*deltak)+0.5));
     m=n*k+c;
 else
     n=floor((phi+2*pi)/(2*pi*deltak)+0.5);
     m=n*k+c+1;
 end
  temp1(loop)=n*deltak;
 answer1(loop)=abs(asin((2*pi*n+delta_phi1)*lambda/(2*pi*d1))/(2*pi)*360);%利用phi1求解
 answer2(loop)=abs(asin((2*pi*m+delta_phi2)*lambda/(2*pi*d2))/(2*pi)*360);%利用phi2求解
 delta_phi=1/2*(delta_phi2+(2*pi*n+delta_phi1)*d2/d1-2*pi*m);%最小二乘解
 answer(loop)=abs(asin((2*pi*m+delta_phi)*lambda/(2*pi*d2))/(2*pi)*360);
end
result=result+abs(answer-sita0).^2;
Loop
end
result=(result/Loop).^0.5;
figure;plot(0.0001*(10:10:100),result(10:10:100),':*')
xlabel('相位误差');
ylabel('测量误差');