CN101662450A - 查表自适应削峰方法 - Google Patents

Abstract

本发明公开一种查表自适应削峰方法，在跳频通信系统削峰过程中及时自适应调节滤波器的系数，准确构建削峰脉冲对输入信号进行有效削峰处理，包括两个步骤：在其预备的步骤中，将系统总带宽差分分为若干带宽区段，为任意一带宽区段构造一对应的滤波器获得其滤波器系数，存储表征所有带宽区段与相应滤波器系数之间的一一对应关系的数据表；而在其削峰的步骤中，当系统带宽改变时，依当前时刻的信号带宽运算出相应的索引序列的对应元素，再依该元素进一步确定相应的滤波器系数。本发明能根据系统载波频点变化，以查找预存的数据表得到与当前信号匹配合适的滤波器系数，而自适应调整用于构造削峰脉冲的滤波器，从而降低系统复杂度，保证削峰效果较优。

Description

查表自适应削峰方法

技术领域

本发明涉及应用在GSM和EDGE之类的跳频通信系统中的一种查表自适应削峰方法。

背景技术

在现行的主流移动通信系统中，GSM和EDGE都是多载波跳频通信系统的典型代表。

多载波系统相比于单载波系统，由于多个载波信号在时域上的叠加，其具有较高的峰均功率比(PAPR)。为了无失真的传输这些高峰均功率比的信号，对一些部件(如功率放大器)的线性要求很高，若信号的变化范围超出器件的线性范围，则会产生明显的带内失真和带外辐射，从而导致误码率的升高。所以信号峰均功率比的大小对多载波系统性能有着直接的影响。因此，有效地降低系统的峰均功率比就显得十分必要，这将直接关系到整个系统的性能和实现成本。

跳频就是按跳频序列随机地改变一个载波信号占有的频道频率，即改变载波的中心频点的技术。在一个频道组内各跳频序列应是正交的，各信道在跳频传输过程中不能碰撞。跳频技术能减小多径瑞利衰落对信号的影响，改善误码性能，收到频率分集的效果。跳频的另一个作用是干扰分集。蜂窝小区中的同频干扰，由于不相关的跳频，分离了来自许多小区的强干扰，有效地减轻了远近效应影响。另外，跳频使空中接口信息传输愈益随机化，提高了移动无线系统安全性。

综上所述，多载波跳频系统具有如下的特点：

1)较高的峰均功率比(PAPR)；

2)各载波的中心频点随着时间按跳频序列而变化；

3)系统带宽和系统中心频点随着各载波中心频点的变化而相应的改变。

因此，在设计多载波跳频通信系统的削峰方法时需要考虑该系统的特点，才能更为有效地降低多载波跳频通信系统的峰均功率比(PAPR)。

现有一种PC削峰算法可有效降低系统的峰均功率比，其主要原理为利用具有与输入信号相同频谱模板(Spectrum Mask)的滤波器冲激响应来生成削峰脉冲，对输入信号中超过幅度门限值的点进行削峰。原理如图1和图2所示。

图1中的上图为时域信号幅度曲线，Ath为幅度门限，A为检测到的峰值幅度；图1中的下图为构造好的削峰脉冲。

先检测到输入信号峰值，将滤波器冲激响应移位，使其峰值与检测到的信号峰值位置对准，设此时具有与输入信号相同频谱模板(Spectrum Mask)的滤波器冲激响应为H(n)，1≤n≤L，n和L为整数，L表示滤波器冲激响应的阶数，利用公式(1)构造出削峰脉冲：

CP_val(n)＝(A-Ath)*e^jθ*H(n)，1≤n≤L    (1)

其中e为自然对数的底数，j为虚部符号， $j=\sqrt{-1},$ θ为信号峰值的相位，e^jθ＝cosθ+j·sinθ。构造好削峰脉冲之后用原始信号减去削峰脉冲，结果可见图2。

由图2可见，超过幅度门限值Ath的峰值已经被削到门限值以下，并且削下的峰值是圆滑的，这是因为滤波器冲激响应的尖峰是圆滑的，所以本算法相比于其它削峰算法具有较好的频谱性能。

PC削峰算法优点在于：

1)可有效地降低输入信号的峰均功率比；

2)由于使用具有与输入信号相同频谱模板(Spectrum Mask)的滤波器冲激响应进行削峰，因此引入的带外噪声很小；

3)与其它削峰方法相比较，本算法对信号的EVM(误差向量幅度)性能影响小。

由于这一算法采用滤波器脉冲来削峰，便需要在进行削峰处理之前将滤波器系数设定好，这在载波中心频点不变，带宽不变的通信系统(例如WiMAX系统)中实现时比较简单，只需将具有与输入信号频谱模板相对应的滤波器系数直接设定好即可。

但是，对于GSM和EDGE之类的跳频通信系统来说，其载波中心频点与系统带宽是随时间而变化的，在这种情况下滤波器的频谱特性也应随信号频谱特性的改变而改变，因此，滤波器系数不该是唯一确定的，而应是可变的。

这就必然要求业内需要针对多载波跳频通信系统提出合适的自适应调节滤波器系统的解决方案，以便滤波器能构建出较好的削峰脉冲用于完成信号的削峰操作。

进一步地，如果仅仅考虑自适应调节滤波器系数的实时调节，每次频点改变都需重新计算滤波器系数以进行自适应调节，这样会导致系统额外引入大量的复杂度，因而，需结合系统复杂度进行综合考虑。

发明内容

本发明的目的在于针对现有技术的不足，提供一种查表自适应削峰方法，在跳频通信系统削峰过程中，及时高效地确定一自适应调节滤波器的系数，以便准确构建削峰脉冲对输入信号进行有效削峰处理。

为实现该目的，本发明采取如下技术方案：

本发明查表自适应削峰方法，在跳频通信系统削峰过程中及时自适应调节滤波器的系数，准确构建削峰脉冲对输入信号进行有效削峰处理，其特征在于，包括：

预备步骤：

将系统总带宽差分分为若干带宽区段，为任意一带宽区段构造一对应的滤波器获得其滤波器系数，存储表征所有带宽区段与相应滤波器系数之间的一一对应关系的数据表；

削峰步骤：

1)利用当前输入的系统信号的数字序列判断当前时刻的信号带宽相对于前一相邻时刻是否发生改变，若改变则跳转执行步骤3)，否则执行步骤2)；

2)计算当前系统信号与前一相邻时刻的中心频点的频差，利用该频差对滤波器进行移频处理使其中心频点与当前系统信号的中心频点保持一致，跳转执行步骤4)；

3)确定当前时刻的信号带宽所属的带宽区段，依该带宽区段确定其相应的滤波器系数，并将滤波器进行移频处理使其中心频点与当前系统信号的中心频点保持一致，执行步骤4)；

4)构造出削峰脉冲，以滤波器系数的峰值对准系统信号的峰值，用系统信号减去削峰脉冲。

进一步的，在预备步骤中，先将各带宽区段系数化为唯一性对应的索引序列，每个带宽区段对应索引序列的一个元素，再将索引序列与相应滤波器系数之间建立一一对应关系；而在削峰步骤的子步骤3)中，依当前时刻的信号带宽运算出相应的索引序列的对应元素，再依该元素进一步确定相应的滤波器系数。

具体的，在预备步骤中，存在如下具体步骤：

a)设跳频系统总带宽为Bw。设差分间隔为ΔBw，且Bw/ΔBw＝M为整数，从0-Bw以每ΔBw为差分间隔分为[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]共M+1个带宽区段，得出索引序列[0，1，2，3......M]，完成系数化；

b)设系统采样率为Fs，则归一化带宽为：

[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]/Fs，

分别设计归一化带宽为[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]/Fs的数字低通滤波器，获取其系数并与索引序列中各元素建立一一对应关系；

故而，相应地，在削峰步骤的子步骤2)中，存在具体如下步骤：

设当前带宽为B，对其取整round(B/ΔBw)＝N，0≤N≤M，以确定索引序列中的具体元素，在数据表中检索该具体元素并获得相应的滤波器系数。

由于带宽为0的滤波器并不存在，故而，子步骤b)中索引序列的元素为0的带宽区段的滤波器系数，采用带宽为

代替后获得。

特别的，当所述频差为零时，不需进行移频，滤波器系数保持不变，跳转执行步骤4)。

更为具体的，削峰步骤的子步骤2)中应用如下公式利用频差进行移频：

$H^{t_{k+1}}\left(n\right)=H^{t_k}\left(n\right)*e^{j*\mathrm{\Δ\ω}*n},1\≤n\≤L_{t_{k+1}},L_{t_{k+1}}=L_{t_k},$

其中，k为整数，t_k与t_k+1分别表示相邻的前后两个时刻，Δω为前后两个中心频点的频差，

表示移频前后的滤波器冲激响应数字序列，n为整数，表示滤波器冲激响应数字序列的序号，L表示滤波器阶数。

而所述预备步骤中，构造某个具体滤波器以获取其系数的具体过程如下：

b1、以如下公式确定逼近理想频响函数H_d(e^jω)：

$H_d\left(e^{\mathrm{j\&omega;}}\right)=\left\{\begin{array}{cc}{e}^{-\mathrm{j\&omega;a}},& 0\&le;\left|\mathrm{\&omega;}\right|\&le;{\mathrm{\&omega;}}_c\\ 0,& {\mathrm{\&omega;}}_c<\left|\mathrm{\&omega;}\right|\end{array}\right.,$

其中，ω为归一化的数字频率，ω_c为系统信号的数字带宽，也即滤波器频谱特性的归一化带宽；

b2、L为低通滤波器的阶数，设α＝(L-1)/2，以如下公式求出逼近理想低通冲激响应H_d(n)：

$H_d\left(n\right)=\frac{1}{2\mathrm{\&pi;}}{\&Integral;}_{-\mathrm{\&pi;}}^{\mathrm{\&pi;}}{H}_d\left(e^{\mathrm{j\&omega;}}\right)e^{\mathrm{j\&omega;n}}\mathrm{d\&omega;}$

$=\frac{1}{2\mathrm{\&pi;}}{\&Integral;}_{-{\mathrm{\&omega;}}_c}^{{\mathrm{\&omega;}}_c}{e}^{-\mathrm{j\&omega;\&alpha;}}{e}^{\mathrm{j\&omega;n}}\mathrm{d\&omega;}=\frac{\sin \left[{\mathrm{\&omega;}}_c(n-\mathrm{\&alpha;})\right]}{\mathrm{\&pi;}(n-\mathrm{\&alpha;})},0\&le;n\&le;L-1,$

n为整数，表示低通滤波器冲激响应数字序列的序号，d为重新构建的低通滤波器的标识符；

b3、确定一窗函数类型，估算窗函数长度即为滤波器阶数L，根据窗函数表达式计算窗函数系数win(n)，0≤n≤L-1；

b4、通过下式计算当前时刻的滤波器冲激响应即滤波器系数

$H_d^{t_{k+1}}\left(n\right)=H_d\left(n\right)*\mathrm{win}\left(n\right)=\frac{\sin \left[{\mathrm{\&omega;}}_c(n-\mathrm{\&alpha;})\right]}{\mathrm{\&pi;}(n-\mathrm{\&alpha;})}*\mathrm{win}\left(n\right),0\&le;n\&le;L-1,$

t_k+1为当前时刻的标识符；

完成上述步骤后，形成满足当前系统信号的数字带宽和频谱模板要求的低通滤波器。

所述窗函数可以为矩形窗、三角窗、汉宁窗、海明窗、布莱克曼窗及凯塞窗中任意一种，但不局限于所述的几种，所有公知的窗函数都可以为本发明所采用。

所述削峰步骤的子步骤3)中获得了与当前时刻的信号带宽相对于的低通滤波器系数之后，采用如下公式进行移频处理：

$H^{t_{k+1}}\left(n\right)={H_d}^{t_{k+1}}\left(n\right)*e^{j*{\mathrm{\&omega;}}_{\mathrm{center}}\left(t_{k+1}\right)*n},1\&le;n\&le;L_{t_{k+1}},$

其中，ω_center(t_k+1)为当前时刻系统信号中心频点，

为移频前后的滤波器冲激响应数字序列。

此外，所述步骤4)中，采用如下公式构造削峰脉冲：

$\mathrm{CP}\_\mathrm{val}\left(n\right)=(A-\mathrm{Ath})*e^{\mathrm{j\&theta;}}*H^{t_{k+1}}\left(n\right),1\&le;n\&le;L_{t_{k+1}}$

其中，A为检测到的峰值幅度，Ath为峰值的幅度门限，e为自然对数的底数，j为虚部符号， $j=\sqrt{-1},$ θ为信号峰值的相位，e^jθ＝cosθ+j·sinθ，n为整数， $1\&le;n\&le;L_{t_{k+1}}.$

与现有技术相比，本发明具有如下优点：

1、以数字化处理的思维进行信号处理，适于在诸如FPGA、DSP之类的数字芯片中编程实现；

2、通过比较前后两个相邻时刻系统的数字带宽和中心频率的变化，实时做出响应，在带宽不变中心频点改变的情况下即时进行频率搬移操作，而在带宽改变时则通过重新选取滤波器系数使之适应当前系统信号，换言之，保证滤波器的系数可实现自适应调整，以保证系统信号的连续性处理，从而使滤波器可准确地构造出削峰脉冲，保证对系统信号进行削峰后的优良性能；

3、滤波器系数的选取，并非采用实时计算的方式来实现，而是先将各种不同频谱特性的数字低通滤波器系数预存在系统中，在需要时才通过索引的方式进行检索确定相应项，这种查表方式相对于实时计算的方式而言，能大大降低系统的计算量，从而降低系统复杂度，节约系统开销，自然也提高了系统的响应率；

4、由于所预存的数据表中，各滤波器系数的频谱归一化带宽之间的差值是既定的，如ΔBw/Fs，因此与实际信号相比，归一化带宽误差最大为ΔBw/2Fs，这样便为用户提供了灵活性，可以适当调整差分量ΔBw的大小来调整误差范围。

附图说明

图1为传统的削峰方法的信号峰值检测与削峰脉冲的波形示意图；

图2为图1的削峰方法在削峰前后的信号功率曲线示意图；

图3为多载波通信系统的频谱示意图。

具体实施方式

如下结合附图和实施例进行说明：

本发明的查表自适应削峰方法适合以程序的形式集成于FPGA(FieldProgrammable Gate Array，现场可编程门阵列)、DSP(Digital Signal Processor，数字信号处理器)之类的数字芯片中，这种芯片大量应用于移动通信系统中，因而具有广泛的代表性，在本实施例中被用以辅助阐述，但是，本领域内普通技术人员所公知的其它具有可编程能力的数字芯片，以及将来被开发出具有同等能力的新型数字芯片不应被排除在本发明的保护范围之外。

对于多载波跳频通信系统而言，其各个载波的中心频点随着时间按跳频序列而变化，系统带宽和系统中心频点随着各载波中心频点的变化而相应的改变。因此，用于生成削峰脉冲的滤波器冲激响应应当具有自适应调节的功能。本发明查表自适应削峰方法正是基于此而提出的。

为了后续表述的方便，首先介绍多载波跳频通信系统所涉及的系统参数。本文以N载波跳频通信系统为例阐述本发明的内容，N为大于1的整数。在各载波频率没有发生变化的时间段内，一个N载波跳频通信系统即为一个普通的N载波通信系统，其频谱示意图如图3所示。图3的纵坐标为功率谱密度P_ω，横坐标为归一化的数字频率ω：

$\mathrm{\&omega;}=\frac{f}{f_s}\&CenterDot;2\mathrm{\&pi;}$

其中，f表示模拟频率，f_s表示系统的采样率。图3中每载波的载波数字带宽都为B，从左至右各个载波的数字中心频点依次为ω₁，ω₂，…ω_N。图3所示意的多载波通信系统的系统数字带宽为B_W，B_W＝ω_N-ω₁+B，系统的数字中心频点为ω_center，ω_center＝(ω₁+ω_N)/2。需要特别指出的是当B＜＜B_W时，B_W≈ω_N-ω₁。

对于N载波跳频通信系统，上述的各载波的中心频点随着时间按跳频序列而变化，系统数字带宽和系统数字中心频点随着各载波中心频点的变化而相应的改变，即都为时间t的函数：ω₁(t)，ω₂(t)，…ω_N(t)，B_W(t)＝ω_N(t)-ω₁(t)+B，ω_center(t)＝[ω₁(t)+ω_N(t)]/2。

假设t_k，t_k+1为载波频点跳变的相邻两个时刻(k为整数，仅起标识时刻作用，如k＝0，下同)，其对应的系统参数分别为：

1)t_k时刻：

载波频点：ω₁(t_k)，ω₂(t_k)，…ω_N(t_k)，

系统数字带宽：B_W(t_k)＝ω_N(t_k)-ω₁(t_k)+B，

系统数字中心频点：ω_center(t_k)＝[ω₁(t_k)+ω_N(t_k)]/2。

2)t_k+1时刻：

载波频点：ω₁(t_k+1)，ω₂(t_k+1)，…ω_N(t_k+1)，

系统数字带宽：B_W(t_k+1)＝ω_N(t_k+1)-ω₁(t_k+1)+B，

系统数字中心频点：ω_center(t_k+1)＝[ω₁(t_k+1)+ω_N(t_k+1)]/2。

系统初始状态可以看成是t_k为负无穷，t_k+1为0时刻的特殊情况。

假设在t_k时刻，根据系统数字带宽B_W(t_k)、系统数字中心频点ω_center(t_k)以及系统频谱模板要求设定的滤波器冲激响应为

$1\&le;n\&le;L_{t_k},$ n和

为整数，

表示此时滤波器的阶数。

从t_k时刻到t_k+1时刻，系统的载波频点发生了变化，此时的滤波器冲激响应

也需要做相应的改变， $1\&le;n\&le;L_{t_{k+1}},$ n和为整数，

表示此时滤波器的阶数。

首先，在预备的阶段，需先建立一表征不同频谱特性的数字低通滤波器系数与带宽区段之间关系的数据表存在于存储单元之中，供后续实时削峰的步骤使用。预备阶段的处理步骤如下：

a)设跳频系统总带宽为Bw。设差分间隔为ΔBw，且Bw/ΔBw＝M为整数，从0-Bw以每ΔBw为差分间隔分为[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]共M+1个带宽区段，得出索引序列[0，1，2，3......M]，完成系数化；

b)设系统采样率为Fs，则归一化带宽为：

[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]/Fs，

分别设计归一化带宽为[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]/Fs的滤波器，获取其系数并与索引序列中各元素[0，1，2，3......M]建立一一对应关系；

c)以数据表的形式表达该一一对应关系，并将该数据表存储在芯片的存储单元(或其它存储设备)中，以供后用。

注意在此一预备步骤中，索引序列[0，1，2，3......M]中各元素已与差分间隔的多个带宽区段相对应，只需求出系统当前信号的带宽与差分间隔ΔBw的比值，即可确定相应的索引序列的元素，从而可获得对应的滤波器系数。

需要指出的是，由于归一化带宽为0的滤波器在实际中是不存在的，因此，在上述子步骤b)中，需采用

代替其中的0值，然后进行参与滤波器的构建，最终得出滤波器系数。

在预备的步骤b)中进行的滤波器构建过程，可以采用多种公知的方式，本实施例中所应用的方式的具体过程如下：

b1、首先由如下公式确定逼近理想低通频响函数H_d(e^jω)：

$H_d\left(e^{\mathrm{j\&omega;}}\right)=\left\{\begin{array}{cc}{e}^{-\mathrm{j\&omega;a}},& 0\&le;\left|\mathrm{\&omega;}\right|\&le;{\mathrm{\&omega;}}_c\\ 0,& {\mathrm{\&omega;}}_c<\left|\mathrm{\&omega;}\right|\end{array}\right.,$

其中ω_c为信号的数字带宽，即滤波器频谱特性的归一化带宽，d为重新构建的低通滤波器的标识符；

b2、因滤波器阶数为L，设α＝(L-1)/2，接着由如下公式求出逼近理想低通冲激响应H_d(n)：

$H_d\left(n\right)=\frac{1}{2\mathrm{\&pi;}}{\&Integral;}_{-\mathrm{\&pi;}}^{\mathrm{\&pi;}}{H}_d\left(e^{\mathrm{j\&omega;}}\right)e^{\mathrm{j\&omega;n}}\mathrm{d\&omega;}$

$=\frac{1}{2\mathrm{\&pi;}}{\&Integral;}_{-{\mathrm{\&omega;}}_c}^{{\mathrm{\&omega;}}_c}{e}^{-\mathrm{j\&omega;\&alpha;}}{e}^{\mathrm{j\&omega;n}}\mathrm{d\&omega;}=\frac{\sin \left[{\mathrm{\&omega;}}_c(n-\mathrm{\&alpha;})\right]}{\mathrm{\&pi;}(n-\mathrm{\&alpha;})},0\&le;n\&le;L-1;$

b3、然后选择窗函数类型，几种常用的窗函数有：矩形窗、三角窗、汉宁(Hanning)窗、海明(Hamming)窗、布莱克曼(Blackman)窗、凯塞(Kaiser)窗，这些窗函数的表达式为本领域所公知，故不行赘述；

b4、窗函数长度即为滤波器阶数L，根据窗函数表达式得到窗函数系数win(n)，0≤n≤L-1，通过如下公式计算数字低通滤波器冲激响应(即滤波器系数)：

$H_d^{t_{k+1}}\left(n\right)=H_d\left(n\right)*\mathrm{win}\left(n\right)=\frac{\sin \left[{\mathrm{\&omega;}}_c(n-\mathrm{\&alpha;})\right]}{\mathrm{\&pi;}(n-\mathrm{\&alpha;})}*\mathrm{win}\left(n\right),0\&le;n\&le;L-1.$

此外，必须重申，构造低通滤波器的方法也不受上述所列限制。本领域内普通技术人员在通读本发明后，应当可举一反三，通过多种等效替换的方式，而本发明为简洁起见，对各知公知的等效替换步骤不行赘述。

在通过上述预备阶段预存数据表之后，便可启动系数进行实时削峰，其具体过程如下：

1)首先判断系统数字带宽是否改变，即B_W(t_k+1)是否等于B_W(t_k)，如果相等则执行步骤2)，否则跳转执行步骤3)；

2)系统数字带宽既未改变，即行计算两个时刻系统数字中心频点的频差Δω＝ω_center(t_k+1)-ω_center(t_k)。

若Δω≠0，则系统数字中心频点已发生改变，利用如下公式进行移频，将滤波器冲激响应的中心频点移至与当前系统信号的中心频点保持一致：

$H^{t_{k+1}}\left(n\right)=H^{t_k}\left(n\right)*e^{j*\mathrm{\&Delta;\&omega;}*n},1\&le;n\&le;L_{t_{k+1}},L_{t_{k+1}}=L_{t_k},$

反之，若Δω＝0，则系统数字带宽及系统数字中心频点均未发生改变，仍采用原来的滤波器冲激响应，即：

$H^{t_{k+1}}\left(n\right)=H^{t_k}\left(n\right),1\&le;n\&le;L_{t_{k+1}},L_{t_{k+1}}=L_{t_k};$

上述两种情况被处理后，完成移频工作，便可跳转至步骤4)进行具体的削峰操作；

3)在系统带宽已发生变化时，便需获取满足系统数字带宽B_W(t_k+1)和系统频谱模板要求的低通滤波器冲激响应(滤波器系数)，也即已预存在数据表中的 ${H_d}^{t_{k+1}}\left(n\right),1\&le;n\&le;L_{t_{k+1}}.$ 具体如下：

如果改变后的系统带宽为B，利用取整函数进行计算：

round(B/ΔBw)＝N，0≤N≤M，如果N＞0，则查表在系统中选择归一化带宽为NΔBw/Fs的滤波器系数，如果N＝0，则查表在系统中选择归一化带宽为

的滤波器系数，这样，利用数据表的一一对应关系，便可通过N检索到相应的数字低通滤波器系数

由于低通滤波器的数字中心频点为零频，而当前时刻系统频谱的中心频点不一定为零频，如果直接用

进行削峰，会对信号性能带来很大影响，故需要将其搬频至此时的系统信号的系统数字中心频点ω_center(t_k+1)，形成更适合本发明使用的带通滤波器，执行如下移频公式：

$H^{t_{k+1}}\left(n\right)={H_d}^{t_{k+1}}\left(n\right)*e^{j*{\mathrm{\&omega;}}_{\mathrm{center}}\left(t_{k+1}\right)*n},1\&le;n\&le;L_{t_{k+1}}.$

这样，便保证重新构造的滤波器的系统数字中心频点与当前系统信号的数字中心频点保持一致，用于后续步骤4)的削峰操作。

4)本实施例因基于传统的削峰方法进行优化，在步骤2)或3)中得到满足当前系统带宽和频谱要求的滤波器冲激响应

之后，构造削峰脉冲和削峰的方式与传统削峰方法中的方式是相同的，即检测输入信号峰值，将滤波器冲激响应移位，使其峰值与检测到的信号峰值位置对准，然后构造削峰脉冲进行削峰。

因而，其主要原理仍然是利用与输入信号的系统数字带宽和频谱模板(Spectrum Mask)相同的滤波器的冲激响应(也即滤波器系数)来生成削峰脉冲，对原系统信号中超过预设的幅度门限值的点进行削峰。故可参阅图1和图2所揭示的公知的构造削峰脉冲的方法，构造削峰脉冲的方程式同于背景技术所引用的公式：

$\mathrm{CP}\_\mathrm{val}\left(n\right)=(A-\mathrm{Ath})*e^{\mathrm{j\&theta;}}*H^{t_{k+1}}\left(n\right),1\&le;n\&le;L_{t_{k+1}},$

其中，A为检测到的峰值幅度，Ath为峰值的幅度门限，e为自然对数的底数，j为虚部符号， $j=\sqrt{-1},$ θ为信号峰值的相位，e^jθ＝cosθ+j·sinθ，n为整数， $1\&le;n\&le;L_{t_{k+1}},$ tk+1仅起标识作用，用于表示系统的当前时刻。

综上所述，本发明查表自适应削峰方法，能根据系统载波频点变化，以查找预存的数据表得到与当前信号匹配合适的滤波器系数，而自适应调整用于构造削峰脉冲的滤波器系数，从而保证削峰效果优于传统技术。

Claims (10)

1、一种查表自适应削峰方法，在跳频通信系统削峰过程中及时自适应调节滤波器的系数，准确构建削峰脉冲对输入信号进行有效削峰处理，其特征在于，包括：

预备步骤：

将系统总带宽差分分为若干带宽区段，为任意一带宽区段构造一对应的滤波器获得其滤波器系数，存储表征所有带宽区段与相应滤波器系数之间的一一对应关系的数据表；

削峰步骤：

1)利用当前输入的系统信号的数字序列判断当前时刻的信号带宽相对于前一相邻时刻是否发生改变，若改变则跳转执行步骤3)，否则执行步骤2)；

2)计算当前系统信号与前一相邻时刻的中心频点的频差，利用该频差对滤波器进行移频处理使其中心频点与当前系统信号的中心频点保持一致，跳转执行步骤4)；

3)确定当前时刻的信号带宽所属的带宽区段，依该带宽区段确定其相应的滤波器系数，并将滤波器进行移频处理使其中心频点与当前系统信号的中心频点保持一致，执行步骤4)；

4)构造出削峰脉冲，以滤波器系数的峰值对准系统信号的峰值，用系统信号减去削峰脉冲。

2、根据权利要求1所述的查表自适应削峰方法，其特征在于：在预备步骤中，先将各带宽区段系数化为唯一性对应的索引序列，每个带宽区段对应索引序列的一个元素，再将索引序列与相应滤波器系数之间建立一一对应关系；而在削峰步骤的子步骤3)中，依当前时刻的信号带宽运算出相应的索引序列的对应元素，再依该元素进一步确定相应的滤波器系数。

3、根据权利要求2所述的查表自适应削峰方法，其特征在于，在预备步骤中，存在如下具体步骤：

a)设跳频系统总带宽为Bw。设差分间隔为ΔBw，且Bw/ΔBw＝M为整数，从0-Bw以每ΔBw为差分间隔分为[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]共M+1个带宽区段，得出索引序列[0，1，2，3......M]，完成系数化；

b)设系统采样率为Fs，则归一化带宽为：

[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]/Fs，

分别设计归一化带宽为[0，ΔBw，2ΔBw，3ΔBw，...，MΔBw]/Fs的数字低通滤波器，获取其系数并与索引序列中各元素建立一一对应关系；

相应地，在削峰步骤的子步骤2)中，存在具体如下步骤：

设当前带宽为B，对其取整round(B/ΔBw)＝N，0≤N≤M，以确定索引序列中的具体元素，在数据表中检索该具体元素并获得相应的滤波器系数。

4、根据权利要求3所述的查表自适应削峰方法，其特征在于：子步骤b)中索引序列的元素为0的带宽区段的滤波器系数，采用带宽为

代替后获得。

5、根据权利要求1所述的查表自适应削峰方法，其特征在于：步骤2)中，当所述频差为零时，不需进行移频，滤波器系数保持不变，跳转执行步骤4)。

6、根据权利要求1至5中任意一项所述的查表自适应削峰方法，其特征在于，步骤2)中应用如下公式利用频差进行移频：

$H^{t_{k+1}}\left(n\right)=H^{t_k}\left(n\right)*e^{j*\mathrm{\&Delta;\&omega;}*n},1\&le;n\&le;L_{t_{k+1}},L_{t_{k+1}}=L_{t_k},$

其中，k为整数，t_k与t_k+1分别表示相邻的前后两个时刻，Δω为前后两个中心频点的频差，

表示移频前后的滤波器冲激响应数字序列，n为整数，表示滤波器冲激响应数字序列的序号，L表示滤波器阶数。

7、根据权利要求1至5中任意一项所述的查表自适应削峰方法，其特征在于，所述预备步骤中，构造某个具体滤波器以获取其系数的具体过程如下：

b1、以如下公式确定逼近理想频响函数

$H_d\left(e^{\mathrm{j\&omega;}}\right)=\left\{\begin{array}{cc}{e}^{-\mathrm{j\&omega;a}},& 0\&le;\left|\mathrm{\&omega;}\right|\&le;{\mathrm{\&omega;}}_c\\ 0,& {\mathrm{\&omega;}}_c<\left|\mathrm{\&omega;}\right|\end{array}\right.,$

其中，ω为归一化的数字频率，ω_c为系统信号的数字带宽，也即滤波器频谱特性的归一化带宽；

b2、L为低通滤波器的阶数，设α＝(L-1)/2，以如下公式求出逼近理想低通冲激响应H_d(n)：

$H_d\left(n\right)=\frac{1}{2\mathrm{\&pi;}}{\&Integral;}_{-\mathrm{\&pi;}}^{\mathrm{\&pi;}}{H}_d\left(e^{\mathrm{j\&omega;}}\right)e^{\mathrm{j\&omega;n}}\mathrm{d\&omega;}$

$=\frac{1}{2\mathrm{\&pi;}}{\&Integral;}_{-{\mathrm{\&omega;}}_c}^{{\mathrm{\&omega;}}_c}{e}^{-\mathrm{j\&omega;\&alpha;}}{e}^{\mathrm{j\&omega;n}}\mathrm{d\&omega;}=\frac{\sin \left[{\mathrm{\&omega;}}_c(n-\mathrm{\&alpha;})\right]}{\mathrm{\&pi;}(n-\mathrm{\&alpha;})},0\&le;n\&le;L-1,$

n为整数，表示低通滤波器冲激响应数字序列的序号，d为重新构建的低通滤波器的标识符；

b3、确定一窗函数类型，估算窗函数长度即为滤波器阶数L，根据窗函数表达式计算窗函数系数win(n)，0≤n ≤L-1；

b4、通过下式计算当前时刻的滤波器冲激响应即滤波器系数

$H_d^{t_{k+1}}\left(n\right)=H_d\left(n\right)*\mathrm{win}\left(n\right)=\frac{\sin \left[{\mathrm{\&omega;}}_c(n-\mathrm{\&alpha;})\right]}{\mathrm{\&pi;}(n-\mathrm{\&alpha;})}*\mathrm{win}\left(n\right),0\&le;n\&le;L-1,$

t_k+1为当前时刻的标识符；

完成上述步骤后，形成满足当前系统信号的数字带宽和频谱模板要求的低通滤波器。

8、根据权利要求7所述的查表自适应削峰方法，其特征在于：所述窗函数为矩形窗、三角窗、汉宁窗、海明窗、布莱克曼窗及凯塞窗中任意一种。

9、根据权利要求1至5中任意一项所述的查表自适应削峰方法，其特征在于，所述削峰步骤的子步骤3)中获得了与当前时刻的信号带宽相对于的低通滤波器系数之后，采用如下公式进行移频处理：

$H^{t_{k+1}}\left(n\right)={H_d}^{t_{k+1}}\left(n\right)*e^{j*{\mathrm{\&omega;}}_{\mathrm{center}}\left(t_{k+1}\right)*n},1\&le;n\&le;L_{t_{k+1}},$

其中，ω_center(t_k+1)为当前时刻系统信号中心频点，

为移频前后的滤波器冲激响应数字序列。

10、根据权利要求1至5中任意一项所述的查表自适应削峰方法，其特征在于，所述步骤4)中，采用如下公式构造削峰脉冲：

$\mathrm{CP}\_\mathrm{val}\left(n\right)=(A-\mathrm{Ath})*e^{\mathrm{j\&theta;}}*H^{t_{k+1}}\left(n\right),1\&le;n\&le;L_{t_{k+1}}$

其中，A为检测到的峰值幅度，Ath为峰值的幅度门限，e为自然对数的底数，j为虚部符号，θ为信号峰值的相位，e^jθ＝cosθ+j·sinθ，n为整数，

Images