本发明涉及电磁波技术，特别涉及赋形波束合成技术。
背景技术：
针对阵列赋形波束(shapedbeampattern,sbp)合成问题，现有方法主要分为两大类：进化类算法和凸优化算法。这两类方法通过迭代搜索所有可能的阵列权系数解空间，得到期望赋形波束对应的阵列权系数，能同时适用于均匀阵列和非均匀阵列。与此同时，其缺点也是明显的，即当阵列个数增加时，以上方法的计算量也随之激增，使得这些方法的计算量大大增加，严重影响阵列天线的实时扫描性能。基于均匀分布的阵列，到目前为止，还没有解析方法对赋形波束的阵列权系数进行直接优化，以得到赋形波束，如平顶方向图、余割平方方向图等。
在对赋形波束，如平顶方向图和余割平方方向图，进行合成时，现有方法主要通过迭代方法搜索阵列权系数的解空间来获得最优的阵列权系数，这些方法通常能够获得相应的最优解。但是在处理大型阵列的赋形波束合成问题时，这些方法的计算量大，计算时间长，不利于阵列天线波束的实时扫描。
技术实现要素：
本发明所要解决的技术问题是，提供一种通过解析方法得到赋形波束对应的阵列权系数从而完成波束合成的方法。
本发明为解决上述技术问题所采用的技术方案是，一种基于笔形波束(pencilbeampattern,pbp)的阵列赋形波束解析合成方法，设计的阵列权系数的计算只与阵列单元个数和角度采样相关，可以对阵列权系数进行直接描述，即通过解析方法来计算阵列权系数以完成波束合成，包括步骤：
1)初始化步骤：对n个阵列单元进行等距布置，各阵列单元的位置为dn，n＝1…n，对阵列天线波束方向图的整个空间进行网格划分得到各波束方向θ，设置赋形波束的主瓣范围θml，设置期望赋形波束fdesired,sbp；
2)通过现有的解析方法，如chebychev切比雪夫方法、taylor泰勒方法等计算得到确定笔形波束对应的阵列权系数wpbp；
3)根据笔形波束对应的阵列权系数wpbp构造笔形波束fpbp(θ,θc)；
4)利用笔形波束构建主瓣θml内各采样点
5)计算赋形波束对应的加权系数
6)得到合成的赋形波束fsbp(θ)与阵列中各个阵元的阵列权系数ωn。
本发明的有益效果是，大大加快了阵列赋形波束的合成速度，适应于解决大型阵列波束赋形问题。
附图说明
图1不同中心角度平顶方向图赋形结果。
具体实施方式
以线阵列天线为例，平面阵列天线的理论过程以此类推。假定天线具有均匀布特性的n个阵元，各阵列单元的位置为dn，n＝1…n。当n为偶数时，dn＝(n-n-0.5)×0.5λ，当n为奇数时，dn＝(n-n-1)×0.5λ。阵列天线在接收信号时，阵列天线的合成电场强度可以描述为：
其中ωn、
对上式进行向量化处理，得到：
其中ae(θ,θc)＝a(θ,θc)⊙e(θ)，⊙表示点乘，阵因子向量a(θ,θc)＝[a1(θ,θc)…an(θ,θc)]h，远场电场强度向量e(θ)＝[e1(θ)…en(θ)]h，阵列权系数
阵列天线的赋形波束可以表示为：
其中，
对均匀阵列而言，对给定副瓣要求的笔形波束，wpbp可以由chebychev或taylor等解析方法直接给出，并且不同方向，即不同笔形波束的指向
fsbp(θ)＝fpbp(θ)wsbp(4)
其中，阵列波束向量
假设期望的赋形波束可以描述为：
此时，合成赋形波束可以写成：
ωn,pbp表示在生成笔形波束时，第n个阵元对应的权系数。
这样，n个阵列的第n个阵元的阵列权系数为：
κ为空间波数。
实施例的具体流程如下：
步骤1:初始化
a)对n个单元进行半波长等距布置，当n为奇数时，单元位置为：
dn＝(n-n-1)×0.5λ；当n为偶数时，dn＝(n-n-0.5)×0.5λ；
b)对线性波束方向图的整个空间，即[0°,180°]进行角度网格划分，波束方向步长设置为△θ＝1°，则整个空间可划分为l＝181个不同的角度；
c)设置赋形波束的主瓣和副瓣角度范围分别为：θml和θsl，令
d)设置期望赋形波束fdesired,sbp；
e)设置期望赋形波束的副瓣电平sll，如sll＝-20db。
步骤2：利用chebychev或taylor方法计算笔形阵因子
步骤3：利用wpbp和公式(2)构造笔形波束fpbp(θ,θc)：
步骤4：构建主瓣θml内各采样点
步骤5：利用公式(5)计算赋形波束对应的加权系数
步骤6：根据公式(7)确定合成赋形波束fsbp(θ)，并根据公式(6)返回总体阵列权系数ωn。
实验验证
采用如下41阵元均匀分布的线阵列天线对本发明所设计的方法进行实验验证：
表141阵元阵列天线的位置信息
本发明设计的方法的结果比较：
例1:对平顶方向图扫描特性的仿真:主瓣宽度设为40o，主瓣中心分别为θc＝0°,20°,40°,60°，副瓣电平设置为：-20db。
仿真结果如图1所示，本发明设计的方法能够有效合成具有不同中心角度θc的平顶方向图，从而证明了本方案设计的方法的有效性。同时，由于本方案设计的方法可以直观地表示出阵列的权系数，因此，该方法是一类解析方法，可以不通过迭代的方式而直接求得。
1.一种基于笔形波束的阵列赋形波束解析合成方法，其特征在于，包括以下步骤：
1)初始化步骤：对n个阵列单元进行等距布置，各阵列单元的位置为dn，n＝1…n，对阵列天线波束方向图的整个空间进行角度网格划分得到各波束方向θ，设置赋形波束的主瓣范围θml，设置期望赋形波束fdesired,sbp；
2)确定笔形波束对应的阵列权系数wpbp，
3)构造笔形波束fpbp(θ,θc)：
其中，θc为主瓣波束的中心指向，h表示共轭转置，ae(θ,θc)＝a(θ,θc)⊙e(θ)，⊙表示点乘，阵列向量a(θ,θc)＝[a1(θ,θc)…an(θ,θc)]h，an(θ,θc)为第n个阵因子，远场电场强度向量e(θ)＝[e1(θ)…en(θ)]h，en(θ)为第n个阵元的远场电场强度；
4)利用笔形波束构建主瓣θml内各采样点
5)计算赋形波束对应的加权系数
6)得到合成的赋形波束fsbp(θ)与阵列中各个阵元的阵列权系数ωn为：
其中，κ为空间波数。
2.如权利要求1所述方法，其特征在于，
