一种卫星多波束天线指向标校信号簇设计方法与流程

1.本发明涉及一种基于粒子群算法的卫星多波束天线指向标校信号簇设计方法，适用于卫星多波束天线指向的标校，属于卫星通信技术领域。


背景技术：

2.多波束天线技术已在国内外绝大多数卫星通信系统中得到应用，以提高系统容量。但由于受到太阳和月亮的牵引、地球引力场不均匀、太阳辐射压力等因素的影响，也会造成波束指向的偏离误差较大。为了保证系统的使用性能，需要在卫星采用三轴姿态动态偏置方法基础上，使用星地一体化波束标校进一步提高波束的指向精度，满足系统使用要求。
3.卫星多波束天线指向的标校可采用下行标校体制。下行波束标校通常由信关站或卫星发射标校信号簇(包括东、西、南、北四个波束)，地面设置波束标校接收站，接收站对不同波束的标校信号能量进行处理，得到卫星天线指向偏差的测量值，并通过控制卫星的姿态或卫星多波束天线的波束成形参数进行指向调整。


技术实现要素：

4.本发明针对针对卫星多波束天线指向偏差的测量与校准问题，提出了一种卫星多波束天线指向标校信号簇设计方法。该方法通过粒子群算法，选择适用于卫星多波束天线指向高精度测量的标校信号簇。本发明设计的标校信号簇具有高功率隔离度、低测量误差等特点，并易于工程实现。
5.本发明所采取的技术方案为：
6.一种卫星多波束天线指向标校信号簇设计方法，包括如下步骤：
7.(1)根据卫星多波束天线标校波束立体方向图、hadamard矩阵以及卫星信道特点，给出地面站接收到的标校信号的表达式；
8.(2)基于卫星多波束天线指向校准的归一化能量差值测向原理，结合卫星多波束天线指向高精度测量与工程实现需求，设计目标函数；
9.(3)根据步骤(2)设计的目标函数，通过粒子群算法进行标校信号簇的设计。
10.进一步的，步骤(1)的具体方式为：
11.(1-1)卫星采用连续并行的方式发射标校信号簇，并由地面站进行接收；根据卫星多波束天线标校波束的立体方向图，得到东、西、南、北四个波束的传输增益分别为ke(θ)、kw(θ)、ks(θ)、kn(θ)，其中，θ为方向图中测量范围内某个点与四个波束等功率交叠点的夹角；
12.(1-2)生成一个hadamard矩阵：
[0013][0014]
式中，m＝2n，n≥5，为kronecker算子，hm的每一个行向量为一组walsh码，码长为m；
[0015]
(1-3)将地面站接收到的标校信号表示为：
[0016]
r(t)＝ke(θ)w1(t)*h(t-τ1) kw(θ)w2(t)*h(t-τ2) ks(θ)w3(t)*h(t-τ3) kn(θ)w4(t)*h(t-τ4) n(t)
[0017]
式中，w1(t)、w2(t)、w3(t)、w4(t)分别为东、西、南、北四个波束使用的walsh码，码型从hadamard矩阵中进行选择，但不包括hadamard矩阵的第一行；h(t)为卫星信道的传递函数；τ1、τ2、τ3、τ4为东、西、南、北四个波束的传输延时；n(t)为高斯白噪声。
[0018]
进一步的，步骤(2)的具体方式为：
[0019]
(2-1)对地面站接收到的标校信号进行离散采样，得到rk(θ，τ1，τ2，τ3，τ4)，其中，θ为在测量角度范围内随机选取的20组不同的值，τ1、τ2、τ3、τ4中任意两个的差值为0到100ns间的随机数；将离散的walsh码表示为w
ik
，构建用于实现接收设备信号捕获的目标函数：
[0020][0021]
其中，c(rk(θ)，w
ik
，k)为rk(θ)与w
ik
的互相关运算，-m≤k≤m；标校信号簇设计的目标是使f1最小化；
[0022]
(2-2)构建用于实现接收设备信号同步的目标函数：
[0023][0024]
f2越小，越易于接收机设备的信号同步；
[0025]
(2-3)针对卫星多波束天线指向的高精度测量中，标校信号簇码型间的高功率隔离度要求，构建目标函数：
[0026][0027]
其中，cov(w
ik
，w
jk
)为协方差矩阵，k＝1、2、3.......m；f3越小，标校信号簇的码型间的功率隔离度越大；
[0028]
(2-4)为了避免传输信道中滤波器对信号功率的影响，标校信号的功率要求集中在滤波器的通带内，对此，构建目标函数：
[0029][0030][0031]
其中，[pi，ni]＝max
i＝1、2、3、4
|fft(w
ik
)|，k＝1、2、3.......m；pi为最大点的功率值，ni为最大点的位置，li为东、西、南、北四个波束使用的码型在hadamard矩阵中的行数；f4与f5越小，传输信道中滤波器对信号功率的影响越小；
[0032]
(2-5)为了实现卫星多波束天线指向标校信号簇设计的综合最优，构建目标函数：
[0033]
f＝＝af1 bf2 cf3 df4 ef5
[0034]
其中，a、b、c、d、e为权值。
[0035]
进一步的，步骤(3)中，通过粒子群算法搜索目标函数f的最小值，进行标校信号簇的设计，具体方式如下：
[0036]
(3-1)设置粒子群参数：粒子数目j＝50，空间维度d＝4，迭代次数n＝100，位置矢量初值为p
j，d
(1)，速度矢量初值为v
j，d
(1)；设置每个个体的最佳适应度初值为1000，种族历史最佳适应度为正无穷；
[0037]
(3-2)根据目标函数计算所有粒子当前的适应度值：
[0038]
f(n)＝＝f＝af1 bf2 cf3 df4 ef5
[0039]
其中，n为迭代次数，权值为：
[0040]
a＝1、b＝1、c＝1、d＝1、e＝100；
[0041]
(3-3)评估粒子个体的适应度，更新粒子群历史最佳可行解s
jd
与整个群体的历史全局最佳可行解gd，j＝1，2，.......50，d＝1，2，3，4；
[0042]
(3-4)更新粒子群的速度及其矢量位置：
[0043]vj，d
(n 1)＝wv
j，d
(n) c1r1(s
jd-p
jd
) c2r2(g
d-p
jd
)
[0044]
p
j，d
(n 1)＝p
j，d
(n 1) v
j，d
(n 1)
[0045]
其中，惯性权重w为0.9，学习因子c1与c2为0.5，r1与r2为介于[0，1]区间的均匀随机数；
[0046]
(3-5)判断是否满足迭代终止条件，满足条件则算法结束，不满足条件则回到步骤(3-2)继续迭代。
[0047]
本发明与背景技术相比具有如下优点：
[0048]
1、本发明根据卫星多波束天线标校波束的立体方向图构建地面站接收到的标校信号表达式。
[0049]
2、本发明所设计的标校信号簇具有高功率隔离度、低测量误差的特点，适用于卫星多波束天线的高精度测量。
[0050]
3、本发明所设计的标校信号簇采用连续并行发射的方式，标校站同时接收四路标校信号并进行能量差值计算，可以消除天线指向、通道增益波动等带来的影响。
附图说明
[0051]
图1是本发明实施例方法的流程图。
具体实施方式
[0052]
下面结合附图对本发明做详细说明。
[0053]
如图1所示，一种卫星多波束天线指向标校信号簇设计方法，包括如下步骤：
[0054]
(1)根据卫星多波束天线标校波束立体方向图、hadamard矩阵以及卫星信道特点，给出地面站接收到的标校信号的表达式。
[0055]
步骤(1)的具体方式如下：
[0056]
(1-1)卫星多波束天线标校波束包括东、西、南、北四个波束。卫星采用连续并行的方式发射标校信号簇，由地面站进行接收。根据卫星多波束天线标校波束的立体方向图，可以得到东、西、南、北四个波束的传输增益分别为ke(θ)、kw(θ)、ks(θ)、kn(θ)，θ为方向图中测量范围内某个点与四个波束等功率交叠点的夹角。
[0057]
(1-2)生成一个hadamard矩阵
[0058][0059]
式中，m＝2n，为kronecker算子。hm的每一个行向量为一组walsh码，码长为m。
[0060]
(1-3)地面站接收到的标校信号可表示为：
[0061]
r(t)＝ke(θ)w1(t)*h(t-τ1) kw(θ)w2(t)*h(t-τ2) ks(θ)w3(t)*h(t-τ3) kn(θ)w4(t)*h(t-τ4) n(t)
[0062]
式中，w1(t)、w2(t)、w3(t)、w4(t)分别为东、西、南、北四个波束使用的walsh码，码型从hadamard矩阵中进行选择(不包括第一行)。h(t)为卫星信道的传递函数；τ1、τ2、τ3、τ4为东、西、南、北四个波束的传输延时；n(t)为高斯白噪声。
[0063]
(2)基于卫星多波束天线指向校准的归一化能量差值测向原理，结合高精度测量与工程可实现需求设计目标函数。
[0064]
步骤(2)的具体方式如下：
[0065]
(2-1)对地面站接收到的标校信号进行离散采样，得到rk(θ，τ1，τ2，τ3，τ4)，θ为在测量角度范围内随机选取的20组不同的值，τ1、τ2、τ3、τ4中任意两个间的差值为0到100ns间的随机数；离散的walsh码可以表示为w
ik
。为了实现接收设备的信号捕获，构建目标函数：
[0066][0067]
c(
rk
(θ)，w
ik
，k)为rk(θ)与w
ik
的互相关运算，-m≤k≤m。标校信号簇设计的目标是使f1最小化。
[0068]
(2-2)为了实现接收设备信号同步，构建目标函数：
[0069][0070]
f2越小，越易于接收机设备的信号同步。
[0071]
(2-3)卫星多波束天线指向的高精度测量中，标校信号簇的码型间要求具有高功率隔离度，因此构建目标函数：
[0072][0073]
其中，cov(w
ik
，w
jk
)为协方差矩阵，k＝1、2、3.......m。f3越小，标校信号簇的码型间的功率隔离度越大。
[0074]
(2-4)为了避免传输信道中滤波器对信号功率的影响，标校信号的功率要求集中在滤波器的通带内，可以构建目标函数：
[0075][0076]
其中，[pi，ni]＝max
i＝1、2、3、4
|fft(w
ik
)|，k＝1、2、3.......m。pi为最大点的功率值，ni为最大点的位置，li为东、西、南、北四个波束使用的码型在hadamard矩阵中的行数。f4与f5越小，传输信道中滤波器对信号功率的影响越小。
[0077]
(2-5)为了实现卫星多波束天线指向标校信号簇设计的综合最优，目标函数可以表示为：
[0078]
f＝af1 bf2 cf3 df4 ef5
[0079]
其中，a、b、c、d、e为权值。可通过粒子群算法搜索目标函数f的最小值，进行标校信号簇的设计。
[0080]
(3)根据上面设计的目标函数，通过粒子群算法进行标校信号簇的设计。
[0081]
步骤(3)的具体方式如下
[0082]
(3-1)设置粒子群参数：粒子数目j＝50；空间维度d＝4；迭代次数n＝100；位置矢量初值为p
j，d
(1)；速度矢量初值为v
j，d
(1)；设置每个个体的最佳适应度初值为1000，种族历史最佳适应度为正无穷。
[0083]
(3-2)根据目标函数计算所有粒子当前的适应度值。
[0084]
f(n)＝f＝af1 bf2 cf3 df4 ef5
[0085]
其中，n为迭代次数，权值为：
[0086]
a＝1、b＝1、c＝1、d＝1、e＝100。
[0087]
(3-3)评估粒子个体的适应度，更新粒子群历史最佳可行解s
jd
与整个群体的历史全局最佳可行解gd。j＝1，2，.......50，d＝1，2，3，4。
[0088]
(3-4)更新粒子群的速度及其矢量位置：
[0089]vj，d
(n 1)＝wv
j，d
(n) c1r1(s
jd-p
jd
) c2r2(g
d-p
jd
)
[0090]
p
j，d
(n 1)＝p
j，d
(n 1) v
j，d
(n 1)
[0091]
其中，惯性权重w为0.9；学习因子c1与c2为0.5；r1与r2为介于[0，1]区间的均匀随机数。
[0092]
(3-5)判断是否满足迭代终止条件，满足条件则算法结束，不满足条件则重复步骤(3-2)。
[0093]
总之，本发明方法基于卫星多波束天线指向标校的归一化能量差值测向原理，通过使用卫星多波束天线立体方向图与hadamard矩阵建立分析数据源；结合卫星多波束天线指向高精度测量与工程可实现需求设计了目标函数；通过粒子群算法选择适用于卫星多波束天线指向高精度测量的标校信号簇。本发明设计的信号簇具有高功率隔离度、低测量误差等特点，并易于工程实现。