MATLAB仿真应用_第9章

1、第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验 9.1 全向天线的波束下倾全向天线的波束下倾 9.2 天线阵的波达方向估计天线阵的波达方向估计9.3 天线阵的波束形成天线阵的波束形成第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验9.1 全向天线的波束下倾全向天线的波束下倾9.1.1 全向天线的波束下倾的基本原理广播用的发射机天线建在天线塔上，希望覆盖广大的地域。在覆盖区的边沿，由于距离远，地球表面弯曲，电磁波场强衰减很快，因此为了有效地覆盖既定的区域，保证区域内场强不低于特定的值，通常采用波束下倾的方法，将轴向排列的半

2、波振子天线，通过调节天线的轴向间距、馈电的相位，使得轴向天线阵列的方向图实现波束下倾。第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-1所示是轴向排列的半波振子天线的结构图。阵元之间的距离是d，垂直轴线与电磁波辐射方向的夹角是，相邻阵元馈入信号的相位差是，相邻阵元发出电磁波到达同一地点的程差是dcos，相邻阵元发出电磁波到达同一地点的相位差参看公式(9-1)。(9-1） 2cos).1 ()1(20deeeEEnjjj第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-1 轴向排列的半波振子天线阵元间的相位差第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验9.1

3、.2 仿真实现根据以上讨论，编写出绘制轴向排列天线阵列的方向图的程序。程序9-1 %均匀轴向排列（8阵元）线状天线阵lam=1;%波长t=0:.01:2*pi;d=.6;%阵元间距a1=-1;%阵元馈入信号的相位差即第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验ps=（d*2*pi*cos（t）/lam）-a1;%阵元发出电磁波到达同一地点的相位差f=abs（1+（exp（j*（ps）+（exp（j*2*（ps）+（exp（j*3*（ps）+.（exp（j*4*（ps）+（exp（j*5*（ps）+（exp（j*6*（ps）+（exp（j*7*（ps）;%8阵元信号的叠加，然后取绝对

4、值T=t+pi/2;polar（T，f)%绘极坐标图程序运行后得出的方向图如图9-2所示。 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-2 轴向阵列波束下倾的仿真图第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验9.2 天线阵的波达方向估计天线阵的波达方向估计 9.2.1 原理全向天线不仅利用率不高，而且对各种信号不加区别地接收，降低了通信质量。定点无线通信采用定向天线，大幅度地改善了通信质量。面对众多移动用户的公众通信网基站和专用移动通信网，采用天线指向即波束可变的天线（智能天线），可以使移动通信的通信质量得到很大的改善。第第9章章 天线及智能天线仿真试验天线及智能

5、天线仿真试验为使天线的波束指向可控，甚至形状可控，采用阵列天线是合适的。在距离通信源足够远的空间里，可以将到达的电磁波视为平面波。对于等距离直线阵天线，由于调制在载波上的基带信号码元宽度与波速的乘积远大于天线阵列的尺寸，因此多个天线阵元上的信号的幅度可视为不变，而它们的载波的相位差则取决于其相互位置、尺寸、波长和到达方向。图9-3所示是智能天线的原理框图。天线接收的无线电信号中有许多成分，其中有我们关心的S信号。天线阵列各个阵元接收的电磁波信号因为阵元排列位置的不同带来相位差。第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验经过特定参数的加权控制器w处理后，进一步改变了各个阵元输出信号

6、的相位和幅度。处理的目标是使得阵元输出的信号和Y中的S成分具有最大输出。用S信号作为基准信号，反馈控制单元的功能就是将输出信号Y与基准信号S的差值（即误差信号），作为调节控制加权控制器w参数的依据。反馈控制的结果是使减小，Y中的S成分加大,也就是说,天线阵列接收方向图指向了S信号的方向。图9-4所示是等距离直线阵智能天线的原理框图，图9-5所示是等距离直线阵的三维图。 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-3 智能天线原理框图第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-4 均匀直线阵原理图第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-5

7、 均匀直线阵三维图第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验以下通过两个简单的例子介绍智能天线中波达方向估计和波束形成的原理。在图9-4中，我们把天线阵元沿x轴排列，从0到M-1。若有一平面波以角（入射线与z轴的夹角）和角（入射线与x轴的夹角）入射到阵列上，第K号阵元上产生的信号为xK，它与0号阵元的相位差是29： cossin2dKAk第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验 式中，与d分别是入射波的波长和阵元的间距，AK亦称阵因子。计入阵因子的影响，第K号阵元的输出是AKxK，即uK。为了使天线阵的输出满足需要，在每个阵元上，用加权因子wK进行控制。这样第K号

8、阵元上输出的信号为wKAKxK，即wK wK 。若到达天线阵的信号是N个，则天线阵的输出是N个信号在M个阵元上的输出的叠加。将问题简化为xy平面的二维问题（sin=1），并用解析式表达如下：X(n)=x1(n),x2(n),，xN(n) (9-2） 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验其中，为第1个信号的入射角。A=A1,A2,AN (9-4）W=W1,W2,WM (9-5）y(t)=WHAXH=WHU (9-6） 其中，U=u0,u1,uM-1H。 求多个信号到达的方向（波达方向）的方法有多种，下面讨论其中的两种方法及其仿真试验的结果。 1第第9章章 天线及智能天线仿真试

9、验天线及智能天线仿真试验9.2.2 求波达方向估计的方法 1. Capon法 Capon法亦称最小方差无畸变响应MVDR（MinimumVarianceDistortionLessResponse）。天线阵列中的阵元数决定了阵列方向图设计中的自由度数。Capon法将阵列中可控的自由度用来形成期望的波束形状，达到对有用信号进行提升和对无用信号进行抑制的目的，并将其优化问题表达为 (9-7） WRWnyEUUHwwmin)(min2第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验其约束条件为WHA( )=1。可以证明上式的解为 代入式（9-6），可以得到相应的功率为 Vn=vN+1vN+2

10、vM(Vn是RUU相应的噪声特征矢量）因为A与Vn的正交性，分母很小，峰值很大,这样可以得出MUSIC法的空间谱8为0)()()(11ARAARWUUHUUp)()(11ARAPUUHp(9-8） 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验其中, 称为噪声子空间的正交投影估计。以下是对七单元线形天线阵在四信号输入情况下的编程仿真设信号1从/4方向入射,信号2从/3方向入射,信号3从/6方向入射,信号4从3/4方向入射。下面是用两种方法求波达方向估计的应用程序。 HnIVnV第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验程序9-2%七单元线形天线阵的波达方向估计d=1;%

11、天线阵元的间距lma=2;%信号中心波长q1=1*pi/4;q2=pi/3;q3=pi/6;q4=3*pi/4;%四个输入信号的方向A1=exp（-2*pi*j*d*0:6*cos（q1）/lma）;%求阵因子A2=exp（-2*pi*j*d*0:6*cos（q2）/lma）;第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验A3=exp（-2*pi*j*d*0:6*cos（q3）/lma）;A4=exp（-2*pi*j*d*0:6*cos（q4）/lma）;A=A1，A2，A3，A4;%得出A矩阵n=1:1900;v1=.015;%四信号的频率v2=.05;v3=.02;v4=.035

12、;d=1.3*cos（v1*n）;1*sin（v2*n）;1*sin（v3*n）;1*sin（v4*n）;%输入信号矢量 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验U=A*d;%总的输入信号U1=（U）;c=cov（U*U1）;%总输入信号的协方差矩阵s，h=eig（c）;%求协方差矩阵的特征矢量及特征值Vn=s（:，5:7）;%取出与零特征值对应的特征矢量ci=inv（c）;%求协方差矩阵的逆矩阵q1b=pi/180:pi/180:pi;forn=1:length（q1b）q1a（n）=q1b（n）;A1a=exp（-2*pi*j*1*0:6*cos（q1a（n）/lma）;P

13、music（n）=（A1a）*A1a*（inv（A1a）*Vn*（Vn）*A1a）;%应用MUSIC法估计输出第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验Pcap（n）=inv（A1a）*ci*（A1a）;%应用Capon法估计输出T（n）=q1a（n）;P1=abs（Pmusic）;P2=abs（Pcap）;endfigure（1)%绘出应用MUSIC法估计的波达方向图polar（T，P1）figure（2)%绘出应用MUSIC法估计的波达方向图T1=T*180/pi;第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验semilogy（T1，P1）;gridfigure（3

14、)%绘出应用Capon法估计的波达方向图polar（T，P2）figure（4)%绘出应用Capon法估计的波达方向图T1=T*180/pi;semilogy（T1，P2）;grid图9-6图9-9所示是程序运行后显示的结果。 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-6 Capon法作出的波达方向估计（极坐标）第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-7 Capon法作出的波达方向估计（直角坐标）第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-8 MUSIC法作出的波达方向估计（极坐标）第第9章章 天线及智能天线仿真试验天线及智能天线仿真

15、试验图9-9 MUSIC法作出的波达方向估计（直角坐标） 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验从对以上四图的分析中可以得到如下结论：两种方法估计得都比较准确。在/6,/4,/3,3/4处有着尖锐的方向图线。MUSIC法的方向图线的幅度更大。第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验9.3 天线阵的波束形成天线阵的波束形成我们以等距离圆阵为例来讨论天线阵的波束形成。图9-10所示是等距离圆线阵的三维图。第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验图9-10 等距离圆线阵的三维图第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验我们把天线阵元顺序

16、定为从OB起顺时针排列为0到M-1。若有一平面波以角入射到阵列上，第K号阵元上产生的信号为xK，它与到达阵元中心的波前的相位差是 ；与r分别是入射波的波长和阵列圆的半径，AK亦称阵因子。为了使天线阵的输出满足需要，在每个阵元加上加权因子wK控制。这样第K号阵元上输出的信号为wKAKxK。若到达天线阵的信号是N个，天线阵的输出是N个信号在M个阵元上的输出的叠加。用解析式表达如下： )cos(2exp()(kKrjA第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验AK=AK1,AK2,AMHX(n)=x0(n),x1(n),xN-1 (n)H (9-11） 式中, 是K阵元以OB为基准顺时

17、针画出的角度。A=A1,A2,AN(9-12）W=w0,w1,wM-1H (9-13）y(t)=WHAX=WHU (9-14）为了求得多个信号到达的方向（波达方向），可以采用上述的Capon、MUSIC两种方法。波束形成可以采用下面的方法：当有多（N）个信号输入时，其中有1个信号是我们关心的，N-1个信号是需要抑制的。方程组（9-15）描述了上述需求的约束条件（四个信号输入中，第一个信号是我们关心的，其余的信号是需要抑制的）。 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验 WHU=1,0,0,0T (9-15）根据信号波达方向的知识U(U=AX) 及约束条件求解方程组（9-15），

18、可以得到 W=w0,w1,wM-1 (9-16）代入式（9-14）可以得到阵列输出的方向特性。下面是四个输入信号八单元圆阵列天线阵的波达方向估计和波束形成。设信号1从2/3方向入射,信号2从/3方向入射,信号3从3/2方向入射,信号4从5.8/6方向入射。下面是用两种算法求波达方向估计及波束形成的应用程序。第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验程序9-3 clear;q1=pi*2/3;q2=pi*1/3;q3=pi*3/2;q4=pi*5.8/6;%四个输入信号的方向E=1;lam=1;%信号波长r=1;%天线阵的半径m1=1;m2=1;m3=1;m4=1;m5=1;m6=

19、1;m7=1;m8=1;%信号幅度a=0:7第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验forn=1:length（a）A1（n）=exp（-j*pi*cos（q1-（pi*a（n）/4）/lam）; A2（n）=exp（-j*pi*cos（q2-（pi*a（n）/4）/lam）; A3（n）=exp（-j*pi*cos（q3-（pi*a（n）/4）/lam）;A4（n）=exp（-j*pi*cos（q4-（pi*a（n）/4）/lam）; 第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验endA=A1，A2，A3，A4;%得出A矩阵n=1:1900;v1=.06;v2=.

20、02;v3=.03;v4=.073;D=1*cos（v1*n）;1*sin（v2*n）;1*sin（v3*n）;1*square（v4*n）;%四个输入信号U=A*D;%总的输入信号第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验U1=（U）;c=cov（U*U1）;%总输入信号的协方差矩阵s，z=eig（c）;%求协方差矩阵的特征矢量及特征值Vn=s（:，5:8）;%取出与零特征值对应的特征矢量ci=inv（c）;bb=1000;Wopte=Abb;%求解线性方程组（求Wopte）q1b=2*pi:-2*pi/180:2*pi/180;b=0:7;第第9章章 天线及智能天线仿真试验天线及智能天线仿真试验fort=1:length（q1b） h（t）=q1b（t