CN1163038C - A phase stagger method to reduce the peak-to-average ratio of multi-carrier signals - Google Patents
A phase stagger method to reduce the peak-to-average ratio of multi-carrier signals Download PDFInfo
- Publication number
- CN1163038C CN1163038C CNB011054433A CN01105443A CN1163038C CN 1163038 C CN1163038 C CN 1163038C CN B011054433 A CNB011054433 A CN B011054433A CN 01105443 A CN01105443 A CN 01105443A CN 1163038 C CN1163038 C CN 1163038C
- Authority
- CN
- China
- Prior art keywords
- carrier
- phase
- initial phase
- envelope
- discrete
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
Description
技术领域technical field
本发明涉及一种降低多载波信号峰平比的相位错开方法。The invention relates to a phase stagger method for reducing the peak-to-average ratio of multi-carrier signals.
背景技术Background technique
目前,降低多载波合成信号的峰平比(峰值功率对平均功率之比)的主要方法有概率削波、信号补偿和相位错开法。其中,概率削波实现比较简单,但是由于削波带来的信号畸变往往会给系统引人杂散。信号补偿虽然不会引起信号的畸变,但是补偿信号的引人必然加大信号的发射功率。而相位错开避免了上述的缺陷,而被广泛采用。At present, the main methods to reduce the peak-to-average ratio (the ratio of peak power to average power) of multi-carrier composite signals include probability clipping, signal compensation and phase staggering. Among them, the implementation of probability clipping is relatively simple, but the signal distortion caused by clipping often introduces spurs to the system. Although signal compensation will not cause signal distortion, the introduction of the compensation signal will inevitably increase the transmission power of the signal. The phase stagger avoids the above-mentioned defects and is widely used.
相位错开方法大致可以分为两类,一类是适用于等幅等频率间隔的多载波信号合成,另一类是适用于幅度不等频率间隔不均匀的信号合成。Phase staggering methods can be roughly divided into two categories, one is suitable for multi-carrier signal synthesis with equal amplitude and frequency interval, and the other is suitable for signal synthesis with different amplitude and frequency interval.
一般来说,要采用相位错开的方法来压缩信号包络,最容易想到的是用穷尽查找的方法找到一组最优相位组。但是,穷尽查找法由于运算量太大而只适合于载波数比较少,并且实时性要求不高的应用中。为解决载波数较多时的相位错开问题,从大量的实验得出了经验公式(1)、经验公式(2)和半经验公式(3)。Generally speaking, to use the method of phase staggering to compress the signal envelope, the easiest thing to think of is to use the exhaustive search method to find a group of optimal phase groups. However, the exhaustive search method is only suitable for applications with a relatively small number of carriers and low real-time requirements due to the large amount of calculation. In order to solve the phase stagger problem when the number of carriers is large, empirical formula (1), empirical formula (2) and semi-empirical formula (3) are obtained from a large number of experiments.
θk=π(k-1)(k-1)/N (1)θ k = π(k-1)(k-1)/N (1)
θk=πk·k/N (2)θ k = πk·k/N (2)
θk=π(k-1)(k-2)/N (3)θ k = π(k-1)(k-2)/N (3)
以上各式中,θk为用于错开的各个载波的初始相位;N为载波数;k为各个载波的谐波数,对于频率分别为1、2、3、4...15其对应各个载波的k值分别为1、2、3、4...15。In the above formulas, θ k is the initial phase of each carrier used for staggering; N is the number of carriers; The k values of the carriers are 1, 2, 3, 4...15 respectively.
上述各种相位错开算法除了等幅等间隔的应用条件限制,要获得好的错开效果其载波数要求较大(通常要求在十几个载波以上)。In addition to the limitation of application conditions of equal amplitude and equal interval, the various phase stagger algorithms mentioned above require a large number of carriers (usually more than a dozen carriers) to obtain a good stagger effect.
相对于前者,变振幅、不等间隔多载波信号合成的相位错开则要复杂的多。由于振幅的变化和频率间隔的不定性,使得很难求得一个初始相位的解析解,因此通常采用的是迭代算法。目前比较有代表性的算法有:夏洛克(Schroeder)法,帕雷克(Patrick)法和时间频率域交换算法。这三种算法中,时间频率域交换算法由于效果明显,实现简单并且其涉及的主要运算都有快速算法而被广泛应用。Compared with the former, the phase staggering of multi-carrier signal synthesis with variable amplitude and unequal intervals is much more complicated. Due to the variation of the amplitude and the uncertainty of the frequency interval, it is difficult to obtain an analytical solution of the initial phase, so an iterative algorithm is usually used. At present, the representative algorithms are: Schroeder method, Patrick method and time-frequency domain exchange algorithm. Among the three algorithms, the time-frequency domain exchange algorithm is widely used because of its obvious effect, simple implementation and fast algorithms for the main operations involved.
如图1所示,为时间频率域交换算法的原理框图。该算法采用了倒推的思路,它从压缩包络峰值这一目的出发,来求解各个载波的初始相位值。首先设置各载波的初始相位,然后由初始相位和幅度构建傅立叶频谱,逆离散傅立叶变换(IDFT),得到时域波形,对时域波形采样,并对采样后的时域波形信号削波,并再进行离散傅立叶变换(DFT),计算初始相位。如此叠代多次,获得要求的各载波初始相位。在这里,削波的范围可以选择为原合成信号包络最大值的75%~95%。削得太少时不容易收敛,削得太多则所得到的相位错开效果不佳。该算法中,离散傅立叶变换用来求解各个载波的初始相位。As shown in Figure 1, it is a functional block diagram of the time-frequency domain exchange algorithm. The algorithm adopts the idea of backward deduction, and it starts from the purpose of compressing the peak value of the envelope to solve the initial phase value of each carrier. First, the initial phase of each carrier is set, and then the Fourier spectrum is constructed from the initial phase and amplitude, and the inverse discrete Fourier transform (IDFT) is used to obtain the time-domain waveform, the time-domain waveform is sampled, and the sampled time-domain waveform signal is clipped, and Discrete Fourier transform (DFT) is then performed to calculate the initial phase. This is repeated for many times to obtain the required initial phase of each carrier. Here, the range of clipping can be selected as 75%~95% of the maximum value of the envelope of the original synthesized signal. It is not easy to converge when you shave too little, and the phasing effect obtained by shaving too much is not good. In this algorithm, discrete Fourier transform is used to solve the initial phase of each carrier.
从以上的描述可以看出,时间频率域交换算法的主要运算量由三部分组成:离散傅立叶变换DFT、逆离散傅立叶变换IDFT和由初始相位和幅度构建傅立叶频谱。虽然DFT、IDFT都有相应的快速算法,但是对于采样点数较多时,运算量还是比较大。另外,以上算法中,由于各载波的初始相位是随机设置的,这样要通过增加叠代算法的次数来达到要求的各载波初始相位,因此增加了运算量,同时增加了系统的成本。From the above description, it can be seen that the main computational load of the time-frequency domain exchange algorithm consists of three parts: discrete Fourier transform DFT, inverse discrete Fourier transform IDFT and the construction of Fourier spectrum from the initial phase and amplitude. Although DFT and IDFT have corresponding fast algorithms, when the number of sampling points is large, the amount of calculation is still relatively large. In addition, in the above algorithm, since the initial phase of each carrier is randomly set, it is necessary to increase the number of iterative algorithms to achieve the required initial phase of each carrier, thus increasing the amount of calculation and increasing the cost of the system.
发明内容Contents of the invention
本发明的目的在于提供一种降低多载波信号峰平比的相位错开方法,该方法通过降低总的运算量使得系统成本降低。The purpose of the present invention is to provide a phase staggering method for reducing the peak-to-average ratio of multi-carrier signals, which reduces the system cost by reducing the total calculation amount.
本发明的目的是这样实现的,一种降低多载波信号峰平比的相位错开方法,包括下列步骤:设置各载波的初始相位,并设定迭代次数为N;利用上述的初始相位构建多载波的合成包络;对多载波合成包络抽样,得到离散包络;对离散包络的时域信号削波,获得削波后的样点;将该削波后的样点作为离散傅立叶变换的输入,进行离散傅立叶变换,得到傅立叶频谱;由傅立叶频谱计算,获得迭代次数N等于1的多载波的初始相位;重复上述过程,直至迭代次数为N时,所获得的初始相位即为所求的初始相位。The purpose of the present invention is achieved in this way, a phase stagger method for reducing the peak-to-average ratio of multi-carrier signals, comprising the following steps: setting the initial phase of each carrier, and setting the number of iterations as N; using the above-mentioned initial phase to construct multi-carrier Synthetic envelope of the multi-carrier composite envelope; the discrete envelope is obtained by sampling the multi-carrier synthetic envelope; the clipping of the time domain signal of the discrete envelope is obtained after clipping; the sample point after the clipping is used as the discrete Fourier transform Input, perform discrete Fourier transform to obtain the Fourier spectrum; calculate from the Fourier spectrum, obtain the initial phase of the multi-carrier with the number of iterations N equal to 1; repeat the above process until the number of iterations is N, the obtained initial phase is the desired initial phase.
由于采用了上述的方法,将已有的相位错开方法中所采用的傅立叶频谱构建和IDFT两个模块用对合成包络直接抽样来代替,这一改进可以减少一半的运算量,而不改变算法的效果,并采用粗细搜索方法来获得一个比较合适的初始相位,加快叠代运算的收敛速度,从而减少整个过程中的总的运算量,使系统实现的成本降低。Due to the adoption of the above method, the two modules of Fourier spectrum construction and IDFT used in the existing phase staggering method are replaced by direct sampling of the synthetic envelope. This improvement can reduce the amount of calculation by half without changing the algorithm. effect, and use the coarse and fine search method to obtain a more suitable initial phase, speed up the convergence speed of the iterative operation, thereby reducing the total amount of calculation in the whole process, and reducing the cost of system implementation.
下面结合实施例及附图对本发明作进一步的说明。The present invention will be further described below in conjunction with the embodiments and accompanying drawings.
附图说明Description of drawings
图1是现有的时间频率域交换算法原理框图;Fig. 1 is the functional block diagram of existing time-frequency domain exchange algorithm;
图2是本发明的相位错开方法的过程框图。Fig. 2 is a process block diagram of the phase shift method of the present invention.
具体实施方式Detailed ways
如图2所示,本发明,即:一种降低多载波信号峰平比的相位错开方法,包括下列步骤:As shown in Figure 2, the present invention, that is: a phase stagger method that reduces the peak-to-average ratio of multi-carrier signals, comprises the following steps:
第一步:对于待合并的任意k个载波:Step 1: For any k carriers to be combined:
xi=Ai(coswit+θi) (4)x i =A i (cosw i t+θ i ) (4)
其中,i=0,1,...,k-1。首先设置各载波的初始相位θi,并设定迭代次数为N。其中,各载波的初始相位的设置过程是:Wherein, i=0, 1, . . . , k-1. Firstly, the initial phase θ i of each carrier is set, and the number of iterations is set to N. Among them, the setting process of the initial phase of each carrier is:
a、使用随机数发生器产生M组随机的相位组:[θ0 j,θ1 j,...θk-1 j],j=0,1,...,M-1。然后,把这M个相位组分别代入公式(4)中,求出合成包络峰值最小的一组相位作为初搜索的结果。a. Use a random number generator to generate M groups of random phase groups: [θ 0 j , θ 1 j , . . . θ k-1 j ], j=0, 1, . . . , M-1. Then, these M phase groups are respectively substituted into the formula (4), and a group of phases with the smallest synthetic envelope peak value is obtained as the result of the initial search.
b、找出以上k个载波中,能量最大的一个(为叙述方便这里设为第m个载波)。并把以上初搜索所得的一组相位值,作为初始相位代入(4)式中。b. Find out the one with the largest energy among the above k carriers (here set as the m-th carrier for the convenience of description). And put a group of phase values obtained from the above initial search as the initial phase into (4) formula.
c、调整载波m的初始相位θm(调整幅度为2π/n,调整次数为n),并把调整后的相位θm l(l=0,1,n-1)分别代入(4)式中。c. Adjust the initial phase θ m of the carrier m (the adjustment range is 2π/n, and the number of adjustments is n), and substitute the adjusted phase θ m l (l=0, 1, n-1) into formula (4) middle.
d、将调整初始相位后的最大能量载波m与其它载波合成新的包络,并求出对应合成包络峰值最低的一组相位组。d. Synthesize the maximum energy carrier m after the initial phase adjustment and other carriers into a new envelope, and find a group of phase groups corresponding to the lowest peak value of the synthesized envelope.
第二步:利用上述细选出来的一组初始相位构建多载波的合成包络:y=∑yi。Step 2: Construct a composite envelope of multi-carriers using the group of initial phases finely selected above: y=∑yi.
第三步:对多载波合成包络抽样,得到离散包络:z=[z0,z1,…,zi]。Step 3: Sampling the multi-carrier synthesis envelope to obtain a discrete envelope: z=[z 0 , z 1 , . . . , z i ].
第四步:对离散包络的时域信号削波,获得削波后的样点:z′=[z0′,z1′,…,zi′]。Step 4: clipping the time-domain signal of the discrete envelope, and obtaining sample points after clipping: z′=[z 0 ′, z 1 ′, . . . , z i ′].
第五步:将削波后的样点作为离散傅立叶变换的输入,进行离散傅立叶变换,得到傅立叶频谱:y=[Y0,Y1,…,Yi]。Step 5: Use the sample points after clipping as the input of discrete Fourier transform, and perform discrete Fourier transform to obtain Fourier spectrum: y=[Y 0 , Y 1 , . . . , Y i ].
第六步:由傅立叶频谱,通过计算可获得迭代次数等于1的多载波的初始相位:[ψ0,ψ1,…,ψk-1]。Step 6: From the Fourier spectrum, the initial phase of the multi-carrier with the number of iterations equal to 1 can be obtained by calculation: [ψ 0 , ψ 1 , . . . , ψ k-1 ].
重复上述第二步到第六步,直至迭代次数为为N时,所获得的初始相位即为所求的初始相位。Repeat the above steps from the second to the sixth until the number of iterations is N, and the obtained initial phase is the desired initial phase.
Claims (1)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB011054433A CN1163038C (en) | 2001-02-27 | 2001-02-27 | A phase stagger method to reduce the peak-to-average ratio of multi-carrier signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB011054433A CN1163038C (en) | 2001-02-27 | 2001-02-27 | A phase stagger method to reduce the peak-to-average ratio of multi-carrier signals |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1371203A CN1371203A (en) | 2002-09-25 |
| CN1163038C true CN1163038C (en) | 2004-08-18 |
Family
ID=4654522
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB011054433A Expired - Fee Related CN1163038C (en) | 2001-02-27 | 2001-02-27 | A phase stagger method to reduce the peak-to-average ratio of multi-carrier signals |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1163038C (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8422434B2 (en) * | 2003-02-18 | 2013-04-16 | Qualcomm Incorporated | Peak-to-average power ratio management for multi-carrier modulation in wireless communication systems |
| CN101471722B (en) * | 2007-12-24 | 2013-01-16 | 中兴通讯股份有限公司 | Clipping method and apparatus for TD-SCDMA system containing multiple frequency points |
| CN101540043B (en) * | 2009-04-30 | 2012-05-23 | 南京理工大学 | Analytical Iterative Fast Spectrum Extrapolation Method for Single Image Restoration |
| CN102299882B (en) * | 2010-06-22 | 2014-04-02 | 华为技术有限公司 | Method, device and base station for determining initial phase |
-
2001
- 2001-02-27 CN CNB011054433A patent/CN1163038C/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN1371203A (en) | 2002-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1197312C (en) | Optimized synchronous prestructure | |
| KR101107448B1 (en) | Digital signal processing device and digital signal processing method | |
| CN1248439C (en) | Method for realizing frequency synchronization in orthogonal FDM system | |
| BR112012009375B1 (en) | SYSTEM CONFIGURED TO GENERATE A HIGH FREQUENCY COMPONENT FROM AN AUDIO SIGNAL, METHOD TO GENERATE A HIGH FREQUENCY COMPONENT FROM AN AUDIO SIGNAL AND METHOD TO DESIGN A HARMONIC TRANSPOSITOR | |
| EP4664454A2 (en) | Signal processing method and data processing method and apparatus | |
| CN106789828B (en) | A method of being fed back based on peak time tracking reduces FBMC-OQAM system peak-to-average power ratio | |
| CN1373930A (en) | Digital filter design method and apparatus for noise suppression by spectral substraction | |
| CN1163038C (en) | A phase stagger method to reduce the peak-to-average ratio of multi-carrier signals | |
| CN114881072B (en) | Fourier decomposition signal noise reduction method based on peak envelope spectrum | |
| CN103033653B (en) | There is the generation of the band limited noise of adjustable crest factor | |
| CN109525215A (en) | It is a kind of to compose the experience small wave converting method for determining sub-band boundary using kurtosis | |
| CN1677868A (en) | Dither system for a quantizing device | |
| Sarkar et al. | The interlaced chirp Z transform | |
| De Sena et al. | A fast Mellin transform with applications in DAFx | |
| Lu | Application of DFT filter bank to power frequency harmonic measurement | |
| CN1787413A (en) | Transmission tech. scheme for low peak equal ratio orthogonal frequency division multiplex | |
| CN117978607A (en) | A method for reducing PAPR using superimposed pilot OTFS signal using companding technology | |
| Lee et al. | Wireless communication system based on OFDM system using wavelet transform | |
| Fu et al. | Non‐orthogonal frequency division multiplexing based on sparse representation | |
| US11881829B2 (en) | Method for audio peak reduction using all-pass filter | |
| CN1822581A (en) | Interpolation module, interposer and method thereof | |
| CN100352176C (en) | Narrow-compressing synthesizing method for super wideband symbol pulse based on orthogonal wavelet form | |
| CN1127218C (en) | Long-and short-block switching method based on aware entropy in wavelet audio encode | |
| CN1795649A (en) | Reduced complexity time-frequency trained equalizer for discrete multi-tone based DSL systems | |
| Matsubara et al. | Improvement of wavelet-synchrosqueezing transform with time shifted angular frequency |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20040818 Termination date: 20170227 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |