?%%該程序仿真利用干涉技術測向的方法測試不同分數比的影響
?clear?all;
?close?all;
?clc;
?result=zeros（1100）;
?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（phi12*pi）;
delta_phi20=mod（phi22*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（‘測量誤差‘）;
?
?
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1459??2007-06-11?21:18??work\DOA_k.m

?????文件???????1432??2007-06-12?16:46??work\DOA_phase.m

?????目錄??????????0??2007-06-18?22:27??work

-----------?---------??----------?-----??----

?????????????????3109????????????????????4