【技术实现步骤摘要】
本专利技术涉及通信卫星
，具体涉及一种近似精确重构余弦调制滤波器组的设计方法与装置。
技术介绍
余弦调制滤波器组广泛应用于卫星通信领域，灵活实现信号的复用和解复用功能，降低星载信号转发的复杂度。因此，余弦调制滤波器组具有重要的工程应用价值，广大研究者对其设计方法产生了极大兴趣。根据不同的应用，余弦调制滤波器组有精确重构余弦调制滤波器组和近似精确重构余弦调制滤波器组之分。其中对于近似精确重构余弦调制滤波器组的设计，在应用精度允许的范围内放松滤波器组的精确重构性，从而使得设计复杂度相对较低。通常情况下，近似精确重构余弦调制滤波器组的设计被表示为一个非凸的二次约束二次规划问题，因此其设计过程即为求解非凸的二次约束二次规划。迄今为止，这类规划问题非常难于求解，大都求助于一些启发式算法。这类算法的时效性不强，甚至有些情况下找不到问题的全局最优解，因此，近年来，有研究者把待求问题松弛成为一个半定规划问题，然后利用已有算法高效求解该半定规划，从而得到近似精确重构余弦调制滤波器组。但是，这种方法设计得到的滤波器组，其重构误差非常大，可能无法满足应用需求，更为致命的是，这类算法中没有考虑算法的收敛性和数值稳定性，因此，寻找并提出数值稳定且收敛的方法来设计近似精确重构余弦调制滤波器组，对于未来卫星通信系统中星载信号处理的设计而言，具有非常重要的理论指导意义和工程实用价值。
技术实现思路
由于现有的方法设计得...
【技术保护点】
一种近似精确重构余弦调制滤波器组的设计方法，其特征在于，包括：根据余弦调制滤波器组的定义和余弦调制滤波器组的输出与输入关系，计算得到余弦调制滤波器组的近似精确重构条件；根据所述近似精确重构条件，建立近似精确重构余弦调制滤波器组设计的数学模型；根据所述数学模型，计算得到一组原型滤波器系数；根据所述余弦调制滤波器组的定义和所述一组原型滤波器系数，计算得到近似精确重构余弦调制滤波器组。
【技术特征摘要】
1.一种近似精确重构余弦调制滤波器组的设计方法，其特征在于，包括：
根据余弦调制滤波器组的定义和余弦调制滤波器组的输出与输入关系，计算得到余弦
调制滤波器组的近似精确重构条件；
根据所述近似精确重构条件，建立近似精确重构余弦调制滤波器组设计的数学模型；
根据所述数学模型，计算得到一组原型滤波器系数；
根据所述余弦调制滤波器组的定义和所述一组原型滤波器系数，计算得到近似精确重
构余弦调制滤波器组。
2.根据权利要求1所述的方法，其特征在于，所述根据余弦调制滤波器组的定义和余弦
调制滤波器组的输出与输入关系，计算得到余弦调制滤波器组的近似精确重构条件，包括：
所述近似精确重构条件为：
其中，M为余弦调制滤波器组的通道数，k为大于等于0的整数，表示虚数单位，
ω表示角频率，e表示指数函数的底数，H(ej(ω-kπ/M))为H(z)的幅频响应H(ejω)频移了kπ/M，
而H(z)为余弦调制滤波器组的原型滤波器的传递函数，|·|表示取模操作，表示
对所有属于区间[0,π]上的频率ω而言。
3.根据权利要求2所述的方法，其特征在于，所述根据所述近似精确重构条件，建立近
似精确重构余弦调制滤波器组设计的数学模型，包括：
所述数学模型为：
其中：“minimize”和“subjectto”为数学规划术语，分别表示“最小化”和“约束条件”
的意思，为目标函数，表示以最小化原型滤波器的阻带能量为设计目标，
为约束条件，表示在逼近设计目标过程中，必须保持近似精确
重构条件，b＝[2h(L-1),…,2h(1),2h(0)]T，[·]T表示共轭转置操作，h(n)为原型滤波器，
其长度为N，L＝N/2，n为整数，矩阵矩阵矩阵Dn的第(i,
j)个元素为i＝0,1,…,L-1，j＝0,1,…,L-1，δ(·)为
狄拉克函数，表示向下取整操作，RL表示长度为L的实数集。
4.根据权利要求3所述的方法，其特征在于，所述根据所述数学模型，计算得到一组原
型滤波器系数，包括：
根据式(3)、(4)、(5)和(6)计算W(k)、g(k)、和并根据式(7)确定尺度因子dj(k)(j
＝1,2,…,m)，根据式(8)和(9)进行尺度变换，计算A(k)和c(k)；
根据式(10)对A(k)进行QR分解，得到Y(k)和Z(k)；
根据式(12)、(13)和(14)得到y(k)、s(k)和λ(k+1)；
如果k≥kmax或||s(k)||∞≤ε1或||c(k)||∞≤ε2，则得到原型滤波器系数b＝b(k+1)；否则，
令k＝k+1，重新计算式(3)-(14)；
g(k)＝2Ub(k)(4)
d(k)＝diag(d1,d2,…,dm)(8)
Y(k)＝Q1R-T，Z(k)＝Q2(11)
Z(k)TW(k)Z(k)y(k)＝-Z(k)T(g(k)-W(k)Y(k)c(k))(12)
s(k)＝-Y(k)c(k)+Z(k)y(k)(13)
λ(k+1)＝-Y(k)T(W(k)s(k)+g(k))(14)
其中，“diag”为对角矩阵标识符，λ(k)为第k次迭代Lagrange乘子；向量b的初始值b(0)、
式(2)的Lagrange乘子的初始估计值λ(0)、最大迭代次数kmax、决策变量改变量的下限ε1和约
束违反度的容差ε2均为预设值，k＝0，b(k+1)＝b(k)+s(k)，||·||∞表示向量的无穷范数，kmax表
示最大迭代次数，ε1表示决策变量改变量的下限，ε2表示约束违反度的容差。
5.根据权利要求4所述的方法，其特征在于，所述根据所述余弦调制滤波器组的定义和
所述一组原型滤波器系数，计算得到近似精确重构余弦调制滤波器组，包括：
根据式(1)、(2)中滤波器的对偶特性h(n)＝h(N-1-n)和所述一组原型滤波器系数b＝
[2h(L-1),…,2h(1),2h(0)]T，计算得到原型滤波器系数h(n)(n＝1,2,…,N-1)组成的向量
[h(0),h(1),…,h(L-1),h(L-1),…,h(1),h(0)]＝0.5[bL-1,bL-2,…,b1,b0,bT]，从而得到式
(15)和(16)的近似精确重构余弦...
【专利技术属性】
技术研发人员：阳志明，郭兴波，夏冬玉，陈国杰，刘涛，吕学刚，董丽娟，李巍，
申请(专利权)人：中国人民解放军空军装备研究院侦察情报装备研究所，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。