干涉仪测向

%%该程序仿真利用干涉技术测向的方法,测试不同分数比的影响 clear all; close all; clc; result=zeros(1,100); for Loop=1:500 for loop=1:100 d1=0.58;%% d2=(loop+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*(randn(1)+j*randn(1)); sig2=exp(j*phi2)+0.0001*(randn(1)+j*randn(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+(answer-sita0).^2; Loop end result=(result/Loop).^0.5; figure;plot((10:10:100),result(10:10:100),':*') xlabel('不同分数比'); ylabel('测量误差');