背景技术
在OFDM系统中,通常在频域内将给定的信道分成若干个独立的子信道,在每个子信道上使用一个子载波进行调制,并且各子载波并行传输。由于这一系列子信道信号在时域累加后重叠起来形成多载波信号,这样会在某些时刻出现较大的尖峰脉冲,很容易造成较大的峰平比(PAPR,Peak-to-average-power ratio),即信号的最大峰值功率与平均功率之比较大。而现有的功率放大器无法对频率范围较大的信号进行线性放大,因此存在较高PAPR的信号在通过功率放大器时会产生较大的频谱扩展和严重的带内失真,造成子信道间的相互干扰,从而影响OFDM系统的性能。
为了降低PAPR,通常在通过天线执行发送操作之前,利用非线性峰值抑制技术对待处理的多载波信号进行迭代限幅滤波处理或者峰值抵消处理。其中,迭代限幅滤波处理的核心思想为:预先设置门限值,以迭代的方式对时域OFDM信号中幅值超过该门限值的信号点进行直接消除,并保持消除前该信号点的相位,而将幅值小于或者等于该门限值的信号点直接输出。但是迭代限幅滤波处理会带来一定程度的信号失真、带外辐射以及邻道干扰,并且其复杂度随着迭代次数的增加会迅速升高,因此若从提高频谱性能以及降低复杂度的角度出发,则峰值抵消处理成为一种较佳的选择。
图1示出了现有峰值抵消处理方法的流程图。参见图1,该方法包括:
在步骤101中,计算作为待处理信号的时域OFDM信号中各采样点的限幅比率,并生成限幅序列。
本步骤中,假设时域OFDM信号为s=[s(0),s(1),...,s(n),...,s(JN-1)],其中J为过采样倍数,N为子载波数目,n为[0,JN-1]闭区间内的整数。对于采样点s(n),对应的限幅比率为:
公式1
其中An为采样点s(n)的限幅比率,|s(n)|为s(n)的幅值,Th为预先设置的目标门限。
然后,将所有采样点的限幅比率组成限幅序列A=[A0,A1,...,An,...,AJN-1]T。
在步骤102中,利用限幅序列和预先设置的抵消函数,生成抵消信号。
这里,抵消函数g(n)具有有限的频域带宽,并且根据下述公式2,得到抵消信号k(n):
公式2
在步骤103中,将时域OFDM信号与抵消信号叠加,完成峰值抵消。
至此,结束峰值抵消过程。
假设目标门限Th为1分贝(dB),其对应的幅值为10Th/20=1.12;采用的抵消函数为加矩形截断窗长为3的sinc函数,其函数的时域表达式为g(n)={1,0.64,0,0,0,0,0.64},并且过采样因子J为2。
当时域OFDM信号为s=[-1.30,-1.34,-1.17,-0.52,-0.53,0.41,-1.45,0.55]时,其平均功率为1,峰值功率为2.1025,
按照上述的峰值抵消处理方法,首先按照公式 计算信号s中各个采样点对应的限幅比率,得到限幅序列A=[0.18,0.22,0.05,0,0,0,0.33,0]。然后,将这一限幅序列与前述的抵消函数g(n)进行圆周卷积,得到抵消信号 此后再将得到的抵消信号叠加于信号s上,得到峰值抵消后的时域信号 此时的PAPR值为2.39dB,比原信号s的PAPR值3.23dB有所降低。
在上述的峰值抵消处理方法中,针对每个采样点均计算对应的限幅比率,并按照计算出来的限幅比率和抵消函数确定出抵消信号。但是,由于抵消函数g(n)体现出的是时域多点信号,该函数不仅将对应采样点的峰值抵消至目标门限,同时也给该采样点周围的其他采样点的幅值带来影响,即产生多点时域干扰。该干扰与被干扰的采样点上的信号叠加后,存在幅值增大的可能性,从而导致新的高幅值点的出现,即出现峰值再升(peak regrowth)。例如,上述时域OFDM信号s中的最后一个采样点s(7)=0.55,但是该采样点在峰值抵消后得到
,其幅值因峰值抵消处理而增大。
再者,目标门限越小,峰值抵消过程中产生的多点干扰越多,峰值再升的概率越高。此时就有可能出现再升的峰值大于该采样点原幅值的情况,那么峰值抵消处理操作并没有真正实现PAPR的降低,从而导致PAPR抑制性能较差。
发明内容
本发明提供一种多载波信号处理方法,能够提高PAPR抑制性能。
在本发明的多载波信号处理方法中,包括:
A.将待处理信号中幅值高于预设的比较门限且未被峰值抵消处理的一个采样点作为对象采样点,根据预设的目标门限对所述对象采样点进行单点峰值抵消处理,所述比较门限大于或者等于所述目标门限;
B.利用对所述对象采样点进行单点峰值抵消处理的结果更新所述待处理信号,并在更新后的待处理信号中存在未被峰值抵消的采样点时,返回执行所述步骤A。
较佳地,该方法进一步包括:从所述预设的比较门限中选择当前比较门限;
步骤A所述将待处理信号中幅值高于预设的比较门限且未被峰值抵消处理的一个采样点作为对象采样点,根据预设的目标门限对所述对象采样点进行单点峰值抵消处理为:
将所述待处理信号中未与当前比较门限相比较的一个采样点作为对象采样点,在所述对象采样点的幅值大于所述当前比较门限时,根据所述目标门限和预先设置的抵消函数计算出所述对象采样点的抵消信号;
将所述对象采样点的抵消信号叠加至所述待处理信号,得到对所述对象采样点进行单点峰值抵消处理的结果。
其中,所述将所述待处理信号中未与当前比较门限相比较的一个采样点作为对象采样点为:
按照采样点的排列顺序或者与所述排列顺序相反的顺序,从所述待处理信号的各采样点中选择未与当前比较门限相比较的一个采样点,作为对象采样点。
其中,所述将所述待处理信号中未与当前比较门限相比较的一个采样点作为对象采样点为:
按照随机的方式从所述待处理信号的各采样点中选择未与当前比较门限相比较的一个采样点,作为对象采样点。
其中,所述按照随机的方式从所述待处理信号的各采样点中选择未与当前比较门限相比较的一个采样点为:
预先设置用于表示各采样点被比较的顺序的指针位置查询表,该指针位置查询表以各采样点在所述待处理信号中的序号为元素,各元素的排列顺序为随机设置;读取所述指针位置查询表,并按照该指针位置查询表中的元素排列顺序,从所述待处理信号的各采样点中选择对应的采样点。
其中,所述根据所述目标门限和预先设置的抵消函数计算出所述对象采样点的抵消信号为:
按照
计算出所述对象采样点的限幅比率,其中s(m)为对象采样点值,|s(m)|为对象采样点的幅值,α
m为对象采样点的限幅比率,Th为所述目标门限;
对所述抵消函数进行圆周移位,使得该抵消函数的最高点位于所述对象采样点处;
将所述对象采样点的限幅比率与圆周移位后的抵消函数相乘,得到所述抵消信号。
其中,步骤B所述待处理信号中存在未被峰值抵消的采样点为:
所述待处理信号中存在未与当前比较门限比较的采样点。
较佳地,该方法进一步包括:
当所述待处理信号中不存在未与当前比较门限比较的采样点时,判断是否所有预设的比较门限均已被使用,如果是,则结束本流程;否则,从未被使用的比较门限中选择一个比较门限作为当前比较门限,并返回执行所述将所述待处理信号中未与当前比较门限相比较的一个采样点作为对象采样点。
其中,所述从所述预设的比较门限中选择当前比较门限为:从预设的比较门限中选择数值最大的比较门限,作为当前比较门限;
所述从未被使用的比较门限中选择一个比较门限作为当前比较门限为:从所述未被使用的比较门限中选择数值最大的比较门限,作为当前比较门限。
较佳地,该方法进一步包括:预先设置峰值抵消最大迭代次数;
所述结束本流程之前,进一步包括:判断迭代次数是否达到所述峰值抵消最大迭代次数,如果是,则结束本流程;否则,将迭代次数加1,根据此时待处理信号的平均功率更新所述比较门限和目标门限,从更新后的比较门限中选择当前比较门限,并返回执行所述将所述待处理信号中未与当前比较门限相比较的一个采样点作为对象采样点。
较佳地,该方法进一步包括:
在所述待处理信号中不存在未被峰值抵消的采样点时,对传输至接收端的所述待处理信号的时域信号进行限幅噪声重构处理。
其中,所述进行限幅噪声重构处理为:
对所述传输至接收端的待处理信号的时域信号进行快速傅立叶变换FFT处理、解调、调制以及快速傅立叶逆变换IFFT处理,得到时域的估计信号;
按照步骤A所述单点峰值抵消处理时所采用的参数,对该估计信号进行单点峰值抵消处理;
根据所述估计信号和该估计信号的单点峰值抵消处理结果,重构出峰值抵消噪声信号,并从所述传输至接收端的待处理信号的时域信号中消除重构的峰值抵消噪声信号,得到噪声消除结果。
较佳地,所述得到噪声消除结果之后,进一步包括:
在限幅噪声重构处理迭代次数未达到预先设置的限幅噪声重构处理最大迭代次数时,将所述噪声消除结果作为所述传输至接收端的待处理信号的时域信号,将限幅噪声重构处理迭代次数加1,并返回执行所述对所述传输至接收端的待处理信号的时域信号进行快速傅立叶变换FFT处理。
本发明还提供一种多载波信号处理装置,能够提高PAPR抑制性能。
在本发明的多载波信号处理装置中,包括:峰值抵消模块和存储模块,其中,
所述峰值抵消模块用于从所述存储模块中读取待处理信号,将待处理信号中幅值高于预设的比较门限且未被峰值抵消处理的一个采样点作为对象采样点,根据预设的目标门限对所述对象采样点进行单点峰值抵消处理,并利用单点峰值抵消处理结果更新所述待处理信号,保存在存储模块中;
所述存储模块用于保存预设的比较门限和目标门限以及待处理信号,所述的比较门限大于或者等于目标门限。
较佳地,所述峰值抵消模块包括:控制子模块和计算子模块,其中,
所述控制子模块用于读取所述存储模块所保存预设的比较门限和目标门限,当存在未被使用的比较门限时,选择当前比较门限,并将目标门限和所选择的当前比较门限发送给计算子模块,通知计算子模块启动单点峰值抵消处理,并且接收来自于计算子模块的更新后的待处理信号;
所述计算子模块在控制子模块的通知下,根据目标门限,对待处理信号中幅值高于当前比较门限的采样点逐一进行单点峰值抵消处理和用单点峰抵消处理结果更新待处理信号,在完成所有采样点的处理后,将更新后的待处理信号发送给控制子模块。
较佳地,所述存储模块进一步保存预先设置的峰值抵消最大迭代次数;
所述控制子模块进一步用于从所述存储模块中读取峰值抵消最大迭代次数,在迭代次数未达到该峰值抵消最大迭代次数时,将迭代次数加1,利用来自于所述计算子模块的更新后的待处理信号的平均功率更新所述预设的比较门限和目标门限,将更新结果作为下次迭代的比较门限和目标门限,发送给存储模块,并继续从前述读取存储模块所保存的比较门限和目标门限的操作开始执行。
较佳地,该装置进一步包括:
限幅噪声重构处理模块,用于对传输至接收端的所述待处理信号的时域信号进行限幅噪声重构处理。
较佳地,所述限幅噪声重构处理模块包括:信号准备子模块、峰值抵消子模块和噪声去除子模块,其中,
所述信号准备子模块用于对接收到的所述待处理信号的时域信号进行快速傅立叶变换FFT处理、解调、调制以及快速傅立叶逆变换IFFT处理,得到时域的估计信号;
所述峰值抵消子模块用于按照所述待处理信号在发射端进行峰值抵消处理时所采用的参数,对该估计信号进行单点峰值抵消处理;
所述噪声去除子模块用于根据所述估计信号和该估计信号的单点峰值抵消处理结果,重构出峰值抵消噪声信号,从所述传输至接收端的待处理信号的时域信号中消除重构的峰值抵消噪声信号,得到噪声消除结果。
较佳地,所述噪声去除子模块进一步在限幅噪声重构处理迭代次数未达到预先设置的限幅噪声重构处理最大迭代次数时,将限幅噪声重构处理迭代次数加1,将噪声去除结果作为所述传输至接收端的待处理信号的时域信号发送给信号准备子模块,通知所述信号准备子模块返回执行所述对所述传输至接收端的待处理信号的时域信号进行快速傅立叶变换FFT处理。
由上述方案可见,本发明对高于比较门限的采样点进行单点峰值抵消处理,即每次仅抵消一个点的峰值,并利用峰值抵消处理结果来更新待处理信号,此后再对更新后的待处理信号中的采样点进行单点的峰值抵消操作。这样能够更为准确地反映已处理采样点的峰值抵消操作给全部待处理信号带来的干扰,并且能够及时调整未处理采样点对应的峰值抵消程度,从而有效避免峰值再升情况的出现,提高PAPR抑制性能。
再者,这里预先设置有比较门限和目标门限,其中的比较门限用于确定对象采样点是否需要进行峰值抵消处理,而目标门限用于对所确定的采样点进行单点峰值抵消处理。在比较门限大于或者等于目标门限的情况下,能够实现对带来更多干扰的采样点的优先处理,降低后续过程中峰值抵消操作给待处理信号带来的干扰,从而保证峰值抵消操作的有效性,进一步提高PAPR抑制性能。
具体实施方式
为使本发明的目的、技术方案更加清楚明白,以下参照附图并举实施例,对本发明做进一步的详细说明。
本发明在多载波信号处理过程中,预先设置一个目标门限和至少一个比较门限,利用每一个比较门限对待处理多载波信号中的采样点进行单点峰值抵消。
图2为本发明多载波信号处理方法的示例性流程图。参见图2,该方法包括:
在步骤201中,将待处理信号中幅值高于预设的比较门限且未被峰值抵消处理的一个采样点作为对象采样点,根据预设的目标门限对该对象采样点进行单点峰值抵消处理,所述比较门限大于或者等于所述目标门限;
在步骤202中,利用对所述对象采样点进行单点峰值抵消处理的结果更新所述待处理信号,并在更新后的待处理信号中存在未被峰值抵消的采样点时,返回执行所述步骤201。
图3为本发明多载波信号处理装置的示例性结构图,参见图3,该装置包括:
峰值抵消模块和存储模块,其中峰值抵消模块用于从存储模块中读取待处理信号,将待处理信号中幅值高于预设的比较门限且未被峰值抵消处理的一个采样点作为对象采样点,根据预设的目标门限对该对象采样点进行单点峰值抵消处理,并利用单点峰值抵消处理结果更新所述待处理信号,保存在存储模块中;存储模块用于保存预设的比较门限和目标门限以及待处理信号,所述的比较门限大于或者等于目标门限。
在上述的多载波信号处理方案中,对高于比较门限的采样点进行单点峰值抵消处理,即每次仅抵消一个点的峰值,并利用峰值抵消处理结果来更新待处理信号,此后再对更新后的待处理信号中的采样点进行单点的峰值抵消操作。这样能够更为准确地反映已处理采样点的峰值抵消操作给全部待处理信号带来的干扰,并且能够及时调整未处理采样点对应的峰值抵消程度,从而有效避免峰值再升情况的出现,提高PAPR抑制性能。
再者,这里预先设置有比较门限和目标门限,其中的比较门限用于确定对象采样点是否需要进行峰值抵消处理,而目标门限用于对所确定的采样点进行单点峰值抵消处理。在比较门限大于或者等于目标门限的情况下,能够实现对带来更多干扰的采样点的优先处理,降低后续过程中峰值抵消操作给待处理信号带来的干扰,从而保证峰值抵消操作的有效性,进一步提高PAPR抑制性能。
下面将详细描述本发明的多载波信号处理方案。
实施例1
图4示出了本发明实施例1中多载波信号处理方法的流程图。假设本实施例中目标门限为Th,预先设置了K个比较门限{Th_comp0,Th_comp1,…,Th_compk,…,Th_compK-1},并且待处理信号为时域OFDM多载波信号s=[s(0),s(1),...,s(n),...,s(JN-1)],其中J为过采样倍数,N为子载波数目。参见图4,该方法包括:
在步骤401中,从K个比较门限中选择当前比较门限。
本步骤中在选择当前比较门限时,可以按照降序的顺序进行选择,即先选择较大的比较门限,后选择较小的比较门限;当然,也可以按照其他常用的顺序进行选择。
在步骤402中,将待处理信号中未与当前比较门限比较的一个采样点作为对象采样点。
本步骤中选取对象采样点的方式可以有多种,例如:按照采样点的排列顺序或者与该排列顺序相反的顺序,从待处理信号中选择采样点,或者,按照随机的方式从待处理信号的各采样点中选择采样点。一种随机方式的实例是,预先设置用于表示各采样点被比较的顺序的指针位置查询表INDEX,该查询表中以各采样点在所述待处理信号中的序号为元素,各元素的排列顺序为随机设置,例如:INDEX={2,6,3,0,4,1,7,5},那么按照该查询表对待处理信号中各采样点的处理顺序为[s(2),s(6),s(3),s(0),s(4),s(1),s(7),s(5)]。
在步骤403中,判断对象采样点的幅值是否大于当前比较门限,如果是,则执行步骤404;否则,执行步骤407。
若对象采样点的幅值小于或者等于当前比较门限,表明对象采样点目前在可接受的幅值范围内,因此无需进行单点峰值抵消。
在步骤404中,根据目标门限生成对象采样点的限幅比率。
这里根据下述的公式3计算出对象采样点的限幅比率:
公式3
其中,s(m)为对象采样点值,|s(m)|为对象采样点的幅值,αm为对象采样点的限幅比率,Th为目标门限。
在步骤405中,根据限幅比率和预先设置的抵消函数计算得到对象采样点的抵消信号。
本实施例中,预先设置出具有有限频域带宽并且时域为有限长序列的抵消函数g(n),例如时域截断的理想低通信号。换言之,可以通过带限信号h(n)乘以时域窗函数序列WNs(n)得到抵消函数g(n)。该抵消函数可以表示为:g(n)=h(n)WNs(n),其中(Ns+1)为截断窗长。本实施例中的带限信号h(n)可以是时域为sinc函数的理想低通信号,或者是时域为余弦函数傅立叶变换结果的余弦信号等。这些函数的频域带宽可以与待处理的OFDM信号的带宽不同;此外,时域窗函数序列WNs(n)可以是矩形窗、汉宁(Hanning)窗、海明(Hamming)窗、布拉克曼(blackman)窗、凯撒(Kaiser)窗等。
本步骤中在计算对象采样点的抵消信号时,对抵消函数g(n)进行圆周移位,使得该抵消函数的最高点位于对象采样点处,即沿圆周位移m位,得到g(n-m),其中m为对象采样点在待处理信号中的位置减去1,例如,对象采样点为待处理信号中的第4个采样点,则m为3;此后,将对象采样点的限幅比率与圆周移位后的抵消函数相乘,得到抵消信号αmg(m-n)。
在步骤406中,利用对象采样点的抵消信号更新待处理信号。
本步骤中,将步骤405中得到的抵消信号直接叠加至选出对象采样点的待处理信号中,即步骤402中的待处理信号,作为对象采样点的单点峰值抵消处理结果。具体而言,此时的待处理信号为:c(m)+αmg(m-n)。
以上的步骤404至406即为本实施例中根据目标门限和抵消函数计算抵消信号从而实现单点峰值抵消的过程。
在步骤407中,判断更新后的待处理信号中是否存在未与当前比较门限比较的采样点,如果是,则返回执行步骤402;否则,执行步骤408。
考虑到对对象采样点进行单点峰值抵消操作时,也会对对象采样点周围的其他采样点的幅值造成影响,因此本步骤中针对更新后的待处理信号。并且,这里的采样点是以其在待处理信号中所处的位置进行区分的,即使某一位置的采样点幅值发生变化,也只能够视为一个采样点,而不能够以两个采样点来看待。例如,某一采样点在步骤402至406中进行了单点峰值抵消处理,虽然其幅值发生了变化,但是由于该采样点在单点峰值抵消处理之前已经与当前比较门限进行了比较,因此该采样点就应该被视为是已被比较过的采样点。
在步骤408~409中,判断是否所有比较门限均已被使用,如果是,则结束本流程;否则,从未被使用的比较门限中选择一个比较门限作为当前比较门限,并返回执行步骤402。
以两个比较门限为例,在根据一个比较门限对所有的采样点执行比较处理后,还需要再从步骤402开始,根据第二个比较门限重新对所有的采样点进行比较,进而进行单点峰值抵消。需要说明的是,此时步骤402中的待处理信号也是前一次单点峰值抵消处理中更新后的待处理信号。
至此,结束本实施例中多载波信号处理流程。
从上述的描述可见,本实施例中对幅值高于比较门限的采样点进行单点峰值抵消处理,并在每次单点峰值抵消处理后均更新待处理信号,以便及时且全面地反映出峰值抵消给对象采样点和其他采样点带来的幅值影响,以便在对后续采样点进行处理时,很大程度上避免采样点幅值未被准确反映而导致的峰值再升,从而能够有效地提高PAPR抑制性能。
再者,本实施例在挑选下一待处理采样点时以更新后的待处理信号为准,这样若某一采样点的幅值原本高于当前比较门限,但是在其他采样点的单点峰值抵消处理后降低到小于或者等于当前比较门限,则无需在对该采样点进行单独峰值抵消,这样能够有效地降低多载波信号处理过程中的复杂度。
此外,本实施例还预先设置了多个比较门限。在具体执行过程中,逐一依照各个比较门限对待处理信号进行单点峰值抵消处理操作。若按照从高至低的趋势来选用比较门限,当然无需严格地按照降序顺序,此时幅值较大的采样点被先抵消掉。由于采样点的幅值越大,其峰值抵消处理对相邻采样点的影响也越大,越早对这样的采样点进行处理,则越能够尽快降低其他采样点所受到的影响,从而在后续过程中需要进行单点峰值抵消处理的采样点越少,因此能够进一步降低多载波信号处理过程中的复杂度。
以上为仅进行一次迭代的情况,本实施例中还可以采用多次迭代的方式来进一步提高PAPR抑制性能。此时,在步骤401之前增加将迭代次数初始化为0的操作。步骤402变为:将待处理信号中未在本次迭代中被比较的一个采样点作为对象采样点。步骤407变为:判断更新后的待处理信号中是否存在本次迭代中未与当前比较门限比较的采样点,如果是,则返回执行步骤402;否则,执行步骤408。步骤408和409变为:判断是否所有比较门限均已在本次迭代中被使用,如果是,则执行步骤410;否则,选择一个本次迭代中未被使用的比较门限作为当前比较门限,并返回执行步骤402。步骤410和步骤411是新增的两个步骤,具体来说,步骤410为:判断迭代次数是否达到预先设置的峰值抵消最大迭代次数,如果是,则结束本流程,否则,执行步骤411;步骤411为:将迭代次数加1,根据更新后的待处理信号的平均功率更新比较门限和目标门限,以便下一次迭代使用,然后从更新后的比较门限中选择当前比较门限,并返回执行步骤402。由于本实施例中的目标门限和比较门限均以预设的PAPR门限以及待处理信号的平均功率为依据进行设置的,而在一次迭代完成后,幅值较大的采样点被抵消掉,待处理信号的平均功率势必有所降低。因此,在PAPR门限不变的情况下,下一次迭代的目标门限和比较门限就可以根据本次迭代后待处理信号的平均功率进行更新。具体操作时,下次迭代的目标门限为本次迭代的目标门限乘以更新后待处理信号的平均功率,下次迭代的比较门限为本次迭代的目标门限乘以更新后待处理信号的平均功率。
在这种多次迭代的情况下,每次迭代均针对前次迭代更新后的待处理信号,分别针对各个比较门限对待处理信号进行单点峰值抵消处理,这样能够进一步降低待处理信号中超过比较门限的采样点的幅值,从而在前次迭代的基础上再次提高PAPR抑制性能;而且,本实施例中在每次迭代之后均利用最新待处理的平均功率对目标门限和比较门限进行更新,达到根据实际信号变化调整两种门限的目的,使得后续处理中这两种门限的数值更接近实际需要,从而提高PAPR抑制性能;再一方面,在多次迭代的情况下,由于在先的迭代中能够将大部分超过比较门限的采样点的幅值降低到允许范围之内,那么,从概率的角度看,在后迭代需要处理的采样点将越来越少,因此迭代次数的增加不会带来明显的复杂度升高。此外,无论是一次迭代还是多次迭代,除了最初在生成作为待处理信号的时域OFDM信号时需执行一次快速傅立叶逆变换(IFFT)操作之外,在进行单点峰值抵消时均无需进行任何额外的快速傅立叶变换(FFT)操作或者IFFT操作,因此在获得相同的PAPR抑制性能时,本实施例中的多载波信号处理方案比现有的迭代限幅滤波方式具有更低的复杂度。
以时域OFDM信号 为例,假设迭代次数为2,目标门限为1dB,对应的幅值为1.12,比较门限为{2.5dB,1dB},对应的幅值为{1.33,1.12},在降序处理比较门限的情况下,每次迭代均先处理第一比较门限1.33,再处理第二比较门限1.12;并且抵消函数为加矩形截断窗长为3的sinc函数,其函数的时域表达式为g(n)={1,0.64,0,0,0,0,0,0.64},过采样因子J为2;再者,在采用指针位置查询表来实现随机选取对象采样点的方式下,该指针位置查询表INDEX={3,1,7,5,0,6,2,4}。
按照本实施例的方法,第一次迭代的操作包括:
1、首先确定当前比较门限为第一比较门限1.33;
2、读取指针位置查询表,得到对象采样点为待处理信号中序号为3的采样点s(3)=-0.52。由于|s(3)|=0.52<1.33,则不对该采样点进行处理,转向下一采样点;
3、读取指针位置查询表,得到对象采样点为|s(1)|=1.34>1.33,那么按照步骤404中的公式3得到该采样点的限幅比率为 ,圆周移位后的抵消函数为g(n-1)=[0.64,1,0.64,0,0,0,0,0],则抵消信号为α1g(n-1)=[0.1408,0.22,0.1408,0,0,0,0,0]。在原始的待处理信号上叠加该抵消信号后,得到更新后的待处理信号为:s=[-1.1592,-1.12,-1.0292,-0.52,-0.53,0.41,-1.45,0.55]。此时可以看出,在经过一个采样点的单点幅抵消操作后,除对象采样点之外,其他采样点中也存在幅值降低的情况,例如序号为0的采样点和序号为2的采样点,这样后续处理中这些采样点的峰值抵消结果将会更加理想,从而能够有效地降低PAPR。
4、按照s(7)、s(5)、s(0)、s(6)、s(2)、s(4)的顺序分别与第一比较门限1.33进行比较,对大于比较门限的点进行单点峰值抵消并更新待处理信号,那么,最后得到的待处理信号为:s=[-1.1592,-1.12,-1.0292,-0.52,-0.53,0.6212,-1.12,0.7612]。完成第一比较门限的处理。
5、针对第二比较门限重复前述处理步骤,得到待处理信号为:s=[-1.12,-1.0949,-1.0292,-0.52,-0.53,0.6212,-1.12,0.7863],其平均功率为E{|s(n)|2}=0.79。
6、根据此时待处理信号的平均功率将目标门限更新为1.12×0.79=1.00,第一比较门限更新为:1.33×0.79=1.18,第二比较门限更新为:1.12×0.79=1.00。
此后,在第二次迭代中,在当前比较门限为第一比较门限时,由于不存在幅值大于第一比较门限的采样点,则不对任何一个采样点进行单点峰值抵消;在当前比较门限为第二比较门限时,仍然按照前述s(3)、s(1)、s(7)、s(5)、s(0)、s(6)、s(2)、s(4)的顺序将各采样点作为对象采样点,得到最终的待处理信号为[-1,-0.9621,-0.9685,-0.52,-0.53,0.6980,-1,0.9010]。此时待处理信号的峰值功率为1,平均功率为0.7142,PAPR为1.46dB。可见,与第一次迭代开始之前待处理信号的PAPR3.23dB相比,此时的PAPR已经有了大幅度的降低。
图5示出了本实施例中多载波信号处理装置的结构图。参见图5,该装置包括:峰值抵消模块和存储模块,其中的峰值抵消模块包括控制子模块和计算子模块。
具体来说,本实施例中的控制子模块用于读取存储模块所保存的预设的比较门限和目标门限,当确定存在未被使用的比较门限时,选择当前比较门限,并将目标门限和所选择的当前比较门限发送给计算子模块,通知计算子模块启动单点峰值抵消处理,并且接收来自于计算子模块的更新后的待处理信号。计算子模块按照图4中的步骤402至步骤407进行操作,即在控制子模块的通知下,根据目标门限,对待处理信号中幅值高于当前比较门限的采样点逐一进行单点峰值抵消和用单点幅抵消结果更新待处理信号,在完成所有采样点的处理后,将更新后的待处理信号发送给控制子模块。
进一步,当采用多次迭代的方式时,存储模块还用于保存预先设置的峰值抵消最大迭代次数。相应地,控制子模块还用于从存储模块中读取峰值抵消最大迭代次数,在迭代次数未达到该峰值抵消最大迭代次数时,将迭代次数加1,并利用来自于计算子模块的更新后的待处理信号的平均功率更新预设的比较门限和目标门限,将更新结果作为下次迭代的比较门限和目标门限,发送给存储模块,并继续从前述读取存储模块所保存的比较门限和目标门限的操作开始执行。
以上多载波信号处理装置的各个部分均可以位于待处理信号的发射端内部。
实施例2
前述实施例1的多载波信号处理过程主要在信号的发射端执行,这一过程中,出于降低PAPR的目的,对一些采样点进行了至少一次的单点峰值抵消处理,这些采样点的幅值变化就会在待处理信号中产生PAPR抑制噪声,有可能对接收端接收信号的误码率(BER)造成负面影响。为了降低PAPR抑制噪声带来的影响,从而降低接收端的BER,本实施例在实施例1的基础上,还在接收端增加限幅噪声重构(Clipping Noise Reconstruction,CNR)处理。具体来说,在待处理信号中不存在未被峰值抵消的采样点时,对传输至接收端的所述待处理信号的时域信号进行限幅噪声重构处理。这里不存在未被峰值抵消的采样点是指待处理信号中的所有采样点均经历了前述实施例1的过程。例如,针对各个比较门限,完成了对所有的采样点进行了处理;或者,在多次迭代的情况下,待处理信号中的所有采样点均已完成所有迭代中的处理。而需要说明的是,本发明中在采样点的幅值被利用实施例1中操作进行了降低时,认为该采样点本身被执行了单点幅值抵消;而采样点在各次迭代中均与比较门限比较过,则认为该采样点已经历峰值抵消过程。
图6示出了本实施例中CNR处理的方法流程图。参见图6,该方法包括:
在步骤601中,接收端接收到来自发送端的待处理信号的时域信号后,对接收到的时域信号进行FFT处理,经过解调和调制映射后再进行IFFT处理,得到时域内的估计信号。
由于发射端发出的信号为时域信号,而接收端通常是以频域信号为对象进行解调和调制映射,因此本步骤中首先将接收到的时域信号进行TTF处理,转换为频域信号。另一方面,由于接收端进行CNR处理又需要针对时域信号,因此在解调和调制映射后进行IFFT处理,得到时域信号,作为对发射端所发出信号的估计信号。
在步骤602中,按照发射端进行峰值抵消处理时所采用的参数对估计信号进行单点峰值抵消处理。
由于抵消噪声为根据单点峰值抵消处理前、后的信号获得,因此本步骤中需要对估计信号进行处理,对其中幅值较大的采样点进行抵消,得到抵消处理后的估计单点峰值抵消结果。
具体来说,本步骤中接收端对估计信号执行与发射端完全相同的操作,并且所采用的目标门限、比较门限、迭代次数、选择对象采样点的方式等参数也与发射端相同。例如,发射端按照前述步骤401至409的一次迭代方式对待处理信号进行单点峰值抵消,则这里接收端按照完全相同的目标门限和比较门限步骤对该估计信号执行一次迭代的单点峰值抵消;或者,当发射端采用多次迭代方式时,接收端也对估计信号进行多次迭代的单点峰值抵消。
在步骤603中,根据估计信号和该序列的单点峰值抵消处理结果,重构出峰值抵消噪声信号。
本步骤中,预先根据单点峰值抵消的初始目标门限,确定衰减因子,以反映发射端的峰值抵消处理对待处理信号影响的统计特性;而后,再利用单点峰值抵消结果减去衰减因子与估计信号的乘积,得到峰值抵消噪声信号。假设衰减因子为β,估计信号为
,单点峰值抵消处理结果为
,峰值抵消噪声信号为
,则
在步骤604中,从接收端接收到的时域信号中消除重构的峰值抵消噪声信号,得到噪声消除结果。
假设接收端在步骤601中接收到的时域信号为R,则噪声消除结果
为:
从上述的描述可见,本实施例中根据估计信号在单点峰值抵消前后的变化,确定出由峰值抵消处理而带来的噪声,再将该噪声从接收到的时域信号中消除,这样接收端接收到的信号中的噪声所占的比例有所下降,则接收端在对经过CNR的信号进行分析时,能够获得正确结果的概率大大提高,因此接收端接收信号的BER能够得到有效地改善。
当然,还可以按照预先设置的限幅噪声重构处理最大迭代次数重复执行上述的步骤601至604,从第二次迭代之后将噪声去除结果作为步骤601中的时域信号,以便更有效地消除峰值抵消噪声信号,从而进一步提高接收端接收信号的BER。具体来说,在执行过程中,在限幅噪声重构处理迭代次数未达到预先设置的限幅噪声重构处理最大迭代次数时,将噪声消除结果作为传输至接收端的待处理信号的时域信号,将限幅噪声重构处理迭代次数加1,并返回执行对传输至接收端的待处理信号的时域信号进行FFT处理的操作。
采用上述CNR处理方式时,本实施例中多载波信号处理装置只需在实施例1中图5的基础上增加位于接收端的CNR处理模块,用于对接收到的信号进行CNR处理。图7示出了本实施例多载波处理装置中CNR处理模块的结构示意图。参见图7,该CNR处理模块包括:信号准备子模块、峰值抵消子模块以及噪声去除子模块。其中的信号准备子模块用于对接收到的时域信号进行FFT、解调、调制以及IFFT处理,得到时域内的估计信号;峰值抵消子模块用于按照待处理信号在发射端进行峰值抵消处理时所采用的参数对该估计信号进行单点峰值抵消处理;噪声去除子模块用于根据估计信号和该序列的单点峰值抵消处理结果,重构出峰值抵消噪声信号,从接收端接收到的时域信号中消除重构的峰值抵消噪声信号,得到噪声消除结果。
若采用多次迭代的CNR方式,则所述噪声去除子模块进一步在限幅噪声重构处理迭代次数未达到预先设置的限幅噪声重构处理最大迭代次数时,将限幅噪声重构处理迭代次数加1,将噪声去除结果作为传输至接收端的待处理信号的时域信号发送给信号准备子模块,通知该信号准备子模块返回执行对传输至接收端的待处理信号的时域信号进行FFT处理的操作。
图8示出了现有峰值抵消与实施例1中方案的PAPR比较仿真图。这里的仿真参数为:子载波数目为1024、调制方式为16正交幅度调制(QAM)、过采样倍数J=4,并且仅执行一次迭代。参见图8,其横轴为以dB为单位的PAPR参考值PAPR0、纵轴为实际PAPA大于PAPR0的概率Pr;并且带有方框的线条代表目标门限为5dB时利用现有峰值抵消处理后得到的PAPR,带有圆圈的线条代表目标门限为5dB时利用实施例1方案处理后得到的PAPR,带有三角形的线条代表目标门限为3dB时利用现有峰值抵消处理后得到的PAPR,带有菱形的线条代表目标门限为3dB时利用实施例1方案处理后得到的PAPR。由图8可见,现有峰值抵消方法中,PAPR抑制性能随着目标门限的降低而恶化,在目标门限为3dB时,PAPR已得不到有效的抑制,即现有峰值抵消方法失效;而应用本发明实施例1中的方案时,即使在目标门限为3dB,其PAPR抑制性能依然比较理想。
图9示出了现有迭代限幅方法与实施例1中方案的PAPR比较仿真图。这里假设时域截断窗长为51,目标门限为5dB。在图9中,横轴为以dB为单位的PAPR参考值PAPR0、纵轴为实际PAPA大于PAPR0的概率Pr;并且带有菱形的线条代表迭代次数为1时现有迭代限幅滤波对应的PAPR,带有方框的线条代表迭代次数为5时现有迭代限幅滤波对应的PAPR,带有三角形的线条代表迭代次数为9时现有迭代限幅滤波对应的PAPR,带有圆圈的线条代表迭代次数为1时实施例1对应的PAPR,带有星形框的线条代表迭代次数为3时实施例1对应的PAPR。由图9可见,在PAPR概率为10-3处,实施例1中的方案获得的PAPR抑制性能明显优于现有迭代限幅滤波,并且实施例1中的方案在经过一次迭代时获得的PAPR抑制性能相当于现有迭代限幅滤波中经过五次迭代获得的PAPR抑制性能;而本实施例中的方案在经过三次迭代后则相当于现有迭代限幅滤波经过九迭代后获得的PAPR抑制性能。再有,经过计算可知,实施例1中的方案在经过三次迭代时对应的计算复杂度仅为现有九次迭代限幅滤波的10%左右。可见,迭代次数相同的情况下,本发明实施例1中的方案获得的PAPR抑制性能远远好于现有的迭代限幅滤波;而在PAPR抑制性能相近的情况下,本发明实施例1中的方案对应的复杂度远远低于现有的迭代限幅滤波。
图10为本发明实施例1和实施例2方案的BER性能比较仿真图,这里的时域截断窗长为51,目标门限为3dB,迭代次数为1。参见图10,其横轴为以dB为单位的信噪比(Eb/N0),纵轴为BER;并且,带有圆圈的线条代表实施例1中方案对应的BER,带有菱形的线条代表实施例2中方案对应的BER。由图10可知,当实施例2在接收端进行CNR处理后,对应的BER更低。
图11示出了应用实施例1中峰值抵消方法时不同比较门限数目下PAPR性能比较仿真图,横轴为以dB为单位的PAPR参考值PAPR0、纵轴为实际PAPA大于PAPR0的概率Pr;并且带有菱形的线条代表迭代次数为1、一个比较门限时实施例1对应的PAPR,带有圆圈的线条代表迭代次数为1、两个比较门限时实施例1对应的PAPR,带有方框的线条代表迭代次数为1、三个比较门限时实施例1对应的PAPR,带有星形框的线条代表迭代次数为1、四个比较门限时实施例1对应的PAPR。这里,当仅存在一个比较门限时,该比较门限值为7dB;当存在两个比较门限时,其取值为[7dB,3dB];当存在三个比较门限时,其取值为[7dB,5dB,3dB];当存在四个比较门限时,其取值为[7dB,5dB,4dB,3dB]。由图11中的仿真图可知,比较门限数目越多,PAPR抑制性能越好,但是,当比较门限数目达到某值,例如:3个比较门限,此后PAPR抑制性能的提高程度就会变缓。经计算可知,上述四种情况下的复杂度分别为:32000、36044、40140、44236次乘加运算。这样,在实际应用中,就可以综合考虑对PAPR抑制性能和复杂度的实际要求,来选择比较门限的数目。例如,当存在3个比较门限时,PAPR抑制性能比两个比较门限的情况有明显提高,而复杂度并未有明显的增加,因此可以将3个比较门限作为优选。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。