电子设计工程
Electronic Design Engineering
第26卷Vol.26第1期No.12018年1月Jan.2018
收稿日期：2016-11-10
稿件编号：201611080
基金项目：国家自然科学基金（61401278）
作者简介：刘豪（1992—），男，湖北钟祥人，硕士研究生。研究方向：数字信号处理、阵列天线综合。
阵列天线由多个天线阵元组成，利用波的干涉
原理使阵列在某一特定方向的增益增强。而在某些特定场合，如卫星与地面通信或者雷达天线，由于距离较远杂波干扰大，所以对低副瓣的要求极高。如果对阵列天线的激励不作任何加权处理，其第一副瓣电平理论值大约为-13.4dB ，达不到相控阵雷达对副瓣电平的要求。切比雪夫加权阵列就是一种能兼顾
主瓣波束宽度和旁瓣高度的数字波束形成的最优阵列。它的特点是：在给定的旁瓣高度下能提供最窄的主瓣波束，在给定主瓣波束宽度下能提供最低的旁瓣高度。传统的切比雪夫加权阵列仅局限于一维线阵，而很多应用中使用的都是二维平面阵列。未解决这一问题，本文提出了两种均匀面阵的切比雪夫权值计算方法，即二维切比雪夫多项式近似法和二维窗函数法，将一维切比雪夫阵列推广至二维，并仿真验证其可行性和有效性。
1切比雪夫线阵波束形成
切比雪夫波束形成主要依托于切比雪夫多项式的性质。切比雪夫多项式是微分方程
任豪事件全过程(1-x 2
)d 2T n (x )d x
2
-x d T n (x )
d x +n 2T n (x )=0的解，其n 阶多项式可表示为：当||x >1时，T n (x )=ch(n ∙ch -1x )；当
||x ≤1，
T n (x )=cos(n ∙cos -1x )。由于当x 的取值范围在[-1，1]时，T n (x )的值均在[-1，1]之间震荡，且
T n (x )=0的解的个数为n ，
均位于[-1，1]之间，并分布在对称于x =0的位置，当x >1时，T n (x )递增。切比雪夫多项式满足如下的递推关系式：
基于切比雪夫加权的面阵波束形成方法
刘豪1，2，余金培1，2
，梁广2
（1.中科院上海微系统与信息技术研究所上海200050；2.上海微小卫星工程中心上海201203）摘要：针对传统的切比雪夫加权阵列仅局限于一维线阵，而很多应用如卫星天线、雷达天线等使用的都是二维平面阵列，本文首先推导了一维切比雪夫阵列权值的计算公式，并仿真验证其对旁瓣抑制的效果。接着重点提出了两种均匀面阵的切比雪夫权值计算方法，即二维切比雪夫多项式近似法和二维窗函数法，仿真验证了该方法的正确性和有效性。关键词：切比雪夫；阵列天线；旁瓣抑制；波束形成中图分类号：TN828.5
文献标识码：A
文章编号：1674-6236（2018）01-0140-04
Beamforming for planar arrays using Chebyshev weighting method
LIU Hao 1，2，YU Jin⁃pei 1，2
，LIANG Guang 2
（1.Shanghai Institute of Microsystem and Information Technology ，Chinese Academy of Sciences ，Shanghai 200050，China ；2.Shanghai Engineering Center for Microsatellites ，Shanghai 201203，China ）Abstract:Traditional Chebyshev weighting is just applicable to one-dimensional arrays ，but in many application scenarios ，such as Satellite antenna or radar antenna ，which always use planar two-dimensional arrays.The paper firstly derive calculating formula for one-dimensional Chebyshev weighting ，and simulate its good effect on Sidelobe Suppression.Then give two methods of calculating Chebyshev weighting for planar two-dimensional arrays ，namely two-dimensional Chebyshev polynomial
approximation and two-dimensional window function method ，and simulation results show the effectiveness of the proposed methods.
Key words:Chebyshev ；antenna array ；sidelobe suppression ；beamforming -
-140
{
T 0(x )=1,  T 1(x )=x ,n =0,1
T n (x )=2xT n -1(x )-T n -2(x ),n ≥2
（1）
根据（1）式可以很方便的推导出各阶切比雪夫多项式的具体表达式。
阵元间距为d 的N 元均匀线阵，假定N 为偶数，则其方向图可表示为：
f ()θ=2∑n =1
N 2I n cos[(n -0.5)∙β]
（2）
一千个伤心的理由>我不是沉默的羔羊
式中，β=2πd ∙sin θ/λ，λ为信号波长。由切比雪夫多项式的性质可知，若N 元阵列方向图与N -1阶切比雪夫多项式曲线一致，则该阵列具有最佳方向性，所有旁瓣高度相同且可控，由此可设计切比雪夫加权天线阵列。假设要求阵列方向图的主旁瓣比为S 。则有T N ()x 0=S ，可解得x 0=ch(1n
ch -1S )。为了
使切比雪夫多项式与阵列天线方向图函数对应，令x =x 0∙cos(β/2)，
则结合f ()θ与T N ()x 可解出阵列系数I n 。具体公式如下：
切比雪夫阵列零点值：
θn =2cos -1
ìíîü
ýþ1x 0cos éëêùûú()2n -1π2∙N -2,
n =1,2,⋯,N 2
-1
（3）根据f ()θn =0，给定第N /2个阵元的幅值为I N/2，
则第1到N /2-1个阵元的幅值可由式（3）求得[1]：
éëêêêêêêù
ûú
úúúúúI 1I 2⋮I N 2-1=-I N 2éë
êêêêê
ê
cos(θ1/2)cos(3∙θ1/2)⋯cos(θ2/2)cos(3∙θ2/2)⋯⋮cos(θN 2-1/2)cos(3∙θN 2-1/2)
我们都要好好的歌词ùû
úú
ú
úúúcos[(N -3)∙θ1/2]cos[(N -3)∙θ2/2]⋮⋯cos[(N -3)∙θN 2-1/2]-1
éëêêêêêùû
úúúúú
cos [](N 2-0.5)∙θ1cos [](N 2-0.5)∙θ2⋮cos éëêùûú(N 2-0.5)∙θN 2-1（4）
2切比雪夫面阵波束形成
2.1
二维切比雪夫多项式近似法
假设均匀面阵x 方向和y 方向阵元间距分别为
d x 、d y ，θ,φ分别代表方位角和俯仰角，I m ,n 代表阵
元的初始复激励。则阵列方向图为：
f (u ,v )=∑m =0M -1∑n =0
N -1I m ,n ∙
exp[-j ∙k c (md x u +nd y v )]
（5）
式中：u =sin θcos φ，v =sin θsin φ。根据工程实际需要，选择x 方向、y 方向等间距的4×4均匀面阵，则由于均匀平面阵列的对称性，可将方向图化简为：
f u ,v =4∙(I 1,1cos 3u 2cos 3v 2+I 1,2cos u 2cos 3v 2
+
I 2,1cos 3u 2cos v 2+I 2,2cos u 2cos v 2
)（6）
考虑用两个切比雪夫多项式相乘来模拟切比雪
夫阵列，则可构造二维切比雪夫多项式，如下：
T m ,n ()x ,y =T m (x )∙T n (y )
（7）
取m =3，n =3，主旁瓣比S ，三阶切比雪夫多项式
为T 3(x )=4x 3-3x ，则T 3,3()x 0,y 0=(4x 03-3x 0)(4y 03-
3y 0)=S 由于阵列x 方向和y 方向对称，故x 0=y 0，令
x =x 0∙cos u 2,y =y 0cos v 2
，带入上式并与（5）式联合可
求得系数I 1,1,I 1,2,I 2,1,I 2,2，再由对称性可求得各个阵
元系数。
由于该方法将两个一维函数相乘模拟二维函数，所以二维切比雪夫多项式近似法仅在
θ=0∘和φ=0∘时达到最优，在波束指向角为其它任意
角度的时候，用该方法求得阵列系数的方向图旁瓣高度虽有一定的优化，但并未达到理想的抑制旁瓣
的效果。并且当阵元数很大时，公式的推导和计算量很大，实现起来较为麻烦。下面引入二维窗函数法实现一维到二维的转换。
2.2二维窗函数法
根据一维N 阶切比雪夫多项式可轻松推导出N
阶切比雪夫系数，该系数是包含N 个元素的一维向量。对于规则的平面阵列，可考虑使用一维窗函数来构造二维窗函数。通常的构造方法有两种：一种方法是将两个一维窗函数直接相乘，得到二维窗函数，如下式所示：
w R (m ,n )=w 1(m )∙w 2(n )T
（8）
式中w 1(m )和w 2(n )是一维窗序列（函数），
w R (m ,n )是二维目标窗序列
（函数）。显然w R (m ,n )具有矩形支持区域。此方法由于是把两个一维矢量相乘得到参数矩阵，所以根据二维窗的性质，该方法不具有最优的效果。
另一种方法是通过以一维窗函数w (x )为轴旋转
w (x )得到圆形轴对称的二维窗函数，如式：
w R (m ,n )=w (m 2+n 2)（9）
根据其傅里叶变换分析可知，如果一维的窗函数具有优良的性质，则其构造的二维窗函数亦具有
刘豪，等基于切比雪夫加权的面阵波束形成方法
-
-141
《电子设计工程》2018年第1期
较好的性质。也可将切比雪夫窗在频域旋转，然后进行傅里叶反变换得到时域窗函数，该方法公式推导较为麻烦，本文不展开详细推导。常用的二维海明窗、凯森窗等都是通过其一维窗的旋转变换获得的。文中采用时域旋转法来构造二维切比雪夫窗函数，可利用MATLAB 对时域旋转做一些改进，计算二维切比雪夫系数矩阵。限于篇幅，本文不对具体计算二维切比雪夫系数的过程加以展开，仅给出通过两种窗函数法计算的阵列系数。
3仿真及性能分析
张朝阳老婆3.1
切比雪夫线阵仿真
选取8个全向天线组成阵元间距为d 的均匀线阵，阵元间距与波长比d /λ=0.5，旁瓣增益为-40dB ，则根据前面推导的公式计算各个阵元切比雪夫系数为：
I 1:I 2:I 3:I 4:I 5:I 6:I 7:I 8=
1:2.86:5.20:6.84:6.84:5.20:2.86:1
波束指向角分别为θ=0°及θ=20°时，仿真结果如
图1和图2
所示。
图1波束指向角为
0°
图2波束指向角为20°观察图1和图2可知，切比雪夫阵列对波束的旁瓣有很好的抑制效果，并且副瓣电平是可控的，即可通过设置目标副瓣电平来计算对应的阵元权值系数。通过对比图1和图2可知，切比雪夫加权可与相
位控制同时进行，在调整波束指向角的同时也可以压低旁瓣。但在抑制旁瓣的同时，也会展宽主瓣波束，所以工程应用中需结合实际情况来对二者取折中。
3.2切比雪夫面阵波束形成仿真
结合工程中实际使用的4×4均匀面阵，x 方向
和y 方向阵元间距d x =d y =0.0192m ，阵列天线的中心频点为8.36GHz ，旁瓣高度为-20dB 。分别应用二维切比雪夫多项式近似法以及二维窗函数法计算阵列各个阵元的系数如表1和表2所示。
表1二维切比雪夫多项式近似法阵元系数阵元序号阵元激励
1,4,13,161.0000
2,3,14,151.6667
5,8,9,121.6667
6,7,10,112.7778
表2
二维切比雪夫窗函数法阵元系数阵元序号阵元激励
1,4,13,16
1.0000
2,3,14,151.9914
5,8,9,121.9914
6,7,10,113.3046
MATLAB 仿真结果如图3所示。
图3上面三幅图是俯仰角θ及方位角φ联合的三维方向图，下面三幅图是俯仰角θ=0∘的一个剖面，
可以较直观的反映旁瓣抑制度。左图是4×4均匀面阵的各全向阵元的幅度激励均为1，即采用等幅激励
情况下的方向图；中图是采用二维切比雪夫多项式近似法计算得到的阵列系数，仿真得到的方向图；右图采用的是二维窗函数法中的旋转法计算得到的阵列系数，据此仿真的方向图。通过对比图3左、中、右三幅图可以发现，采用二维切比雪夫多项式近似法的旁瓣高度为-18dB 左右，采用切比雪夫窗函数法的旁瓣高度略小于-20dB ，效果稍好于二维切比雪夫多项式法。而不采用切比雪夫加权时的旁瓣高度大约是-12dB ，明显劣于切比雪夫阵列。所以以上两种切比雪夫波束形成方法都对均匀平面阵列的旁瓣有较好的抑制效果。但是，和一维切比雪夫阵列仿真的方向图一样，主瓣波束均有不同程度的展宽。而且由于阵元数较少，所以对旁瓣的抑制程度不如线阵明显。
同时可以发现，二维切比雪夫阵列的旁瓣的可控性明显劣于一维切比雪夫阵列，主要原因是因为一维向二维转换的过程中，并不能维持切比雪夫多项
式的[-1，1]区间函数值以及零点分布的性质。所以关于二维切比雪夫阵列还有一定的解析优化空间。
张一山的父母
4结束语
文中首先分析了切比雪夫多项式的性质，并将
-
-142
切比雪夫多项式与一维线性阵列天线的方向图函数结合，推导了切比雪夫阵列系数的公式，完成了切比雪夫加权阵列的仿真。重点提出了二维切比雪夫多项式近似法以及二维切比雪夫窗函数法，通过此两种方法将切比雪夫多项式从一维推广至二维，实现了切比雪夫面阵波束形成，通过MATLAB 仿真验证了二维切比雪夫波束形成的可行性，并且较好的旁瓣抑制效果。随着阵列天线的进一步发展，阵列的多元化已成为必然的趋势，针对特殊阵列或者不规则阵列的旁瓣抑制等还有待研究改进。
参考文献：
[1]张小飞，陈华伟，仇小锋，等.阵列信号处理及MATLAB 实现[M].北京：电子工业出版社，2015.
[2]张友益，徐朝阳.相位加权低副瓣阵列天线综合技术研究[J].船海电子对抗，2014，37（1）:71-77.[3]管自新.星载多波束接收天线自适应宽带校正方法[J]，遥测遥控，2012，33（2）:27-31.
[4]浦锦毅.存在功率控制误差的多波束低轨卫星通
信系统容量分析[J].遥测遥控，2013，34（4）:58-64.
[5]薛正辉.阵列天线分析与综合[M].北京：北京航空航天大学出版社，2011.
[6]何保荣，李晓歌.一种阵列天线赋形波束综合优
化方法[J].华中师范大学学报，2015，49（3）:368-372.
[7]孙冬暖.阵列天线的方向图综合研究与分析[D].
西安：西安电子科技大学，2014.
[8]夏欢.稳健波束形成及方向图综合研究[D].南京：南京大学，2015.
[9]王雨阳.数字波束形成的工程实现方法[J].雷达与对抗，2011，31（1）:39-42.
[10]赵星惟.星载多波束接收天线阵列校准技术研究[D].上海：中国科学院上海微系统与信息技术研究所，2011.
[11]杜勃伦.基于窗函数法的二维FIR 滤波器仿真[J].中国科技博览，2015，6（5）:388.[12]蒋迅，王淑红.切比雪夫和切比雪夫多顼式的故事[J].科学，2016，68（4）：54-58.[13]李亮.阵列天线方向图综合方法研究[D].哈尔滨：哈尔滨工程大学，2013.[14]化小会，张秋生.切比雪夫多项式在代数信号处理中的应用[J].西南师范大学学报，2016，4（15）：
108-114.
[15]AHammami R ，Ghayoula A ，Gharsallah.Antenna
array synthesis with Chebyshev-Genetic Algorithm method[J].Communications ，Computingand Control Applications （CCCA ）.2011：6031445.
[16]赵冬艳，熊童满，陈韶华.一种简化的平面阵波束形成方法[J].声学技术，2013（S1）:79-80.[17]马晓峰.一种星载阵列天线的高效波束赋形算法
[J].电波科学学报，2016，31（3）：
479-485.图3二维切比雪夫阵列仿真方向图
刘豪，等基于切比雪夫加权的面阵波束形成方法
-
-143