CN108920419B - 基于单参数的分数阶Fourier逆变换分级扫描方法 - Google Patents
基于单参数的分数阶Fourier逆变换分级扫描方法 Download PDFInfo
- Publication number
- CN108920419B CN108920419B CN201810636384.4A CN201810636384A CN108920419B CN 108920419 B CN108920419 B CN 108920419B CN 201810636384 A CN201810636384 A CN 201810636384A CN 108920419 B CN108920419 B CN 108920419B
- Authority
- CN
- China
- Prior art keywords
- scanning
- processing
- order
- wfrft
- formula
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06E—OPTICAL COMPUTING DEVICES; COMPUTING DEVICES USING OTHER RADIATIONS WITH SIMILAR PROPERTIES
- G06E3/00—Devices not provided for in group G06E1/00, e.g. for processing analogue or hybrid data
- G06E3/001—Analogue devices in which mathematical operations are carried out with the aid of optical or electro-optical elements
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2626—Arrangements specific to the transmitter only
- H04L27/2627—Modulators
- H04L27/2634—Inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators in combination with other circuits for modulation
- H04L27/2636—Inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators in combination with other circuits for modulation with FFT or DFT modulators, e.g. standard single-carrier frequency-division multiple access [SC-FDMA] transmitter or DFT spread orthogonal frequency division multiplexing [DFT-SOFDM]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Discrete Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Nonlinear Science (AREA)
- Signal Processing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Algebra (AREA)
- Optics & Photonics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Complex Calculations (AREA)
Abstract
基于单参数的分数阶Fourier逆变换分级扫描方法,本发明适用于采用单参数四项加权分数傅里叶变换SP‑4‑WFRFT的通信系统。针对基于SP‑4‑WFRFT的数字通信系统特性,在变换阶数未知的条件下,为快速地恢复数据信号,数据信号可是调制信号、基带序列、扩频序列等,或估计变换阶数,在逆变换过程中,建立一种分级扫描方法。主要通过定义分级终值常量,利用分级级数变量、扫描变量对变换阶数进行分级,并在每一级进行SP‑4‑WFRFT逆变换,从而得到所需的待解调数据或基带数据。本发明的优点是:SP‑4‑WFRFT逆变换过程中,变换阶数的分级扫描具有优先级,且可依据用户需求灵活地设置。
Description
技术领域
本发明涉及信号处理、数字通信技术领域,尤其是涉及基于单参数的分数阶Fourier(傅里叶)逆变换分级扫描方法。
背景技术
在信号处理领域,随着DFT算法的提出,傅里叶变换占据了主导地位。然而,傅里叶变换也存在局限性,比如对非平稳信号的分析与处理,为了更好地描述信号特性,出现了分数阶傅里叶变换(Fractional Fourier Transform,FRFT),目前FRFT主要有经典类FRFT(Chirp-type Fractional Fourier Transform,CFRFT)和加权类FRFT(WeightedFractional Fourier Transform,WFRFT),而加权类FRFT又主要有单参数四项加权分数傅里叶变换(Single-Parameter 4-Weighted Fractional Fourier Transform,SP-4-WFRFT)和多参数分数傅里叶变换(Multiple-parameter Fractional Fourier Transform,MPFRFT)。分数阶Fourier变换的应用已从最初的光学领域扩展到其他领域,近年已在通信系统的信号检测、信道估计、同步算法等领域得到广泛应用,且其完全适用于现有的通信发射、接收系统,并不需要额外的装置和系统变化。
由于分数阶Fourier变换的变换阶数敏感度较高,而为了能够正确地恢复原始数据,对于逆变换的变换阶数要求较为苛刻,尤其是在通信系统应用中存在误差或影响条件时。为此,对于分数阶Fourier逆变换过程中,变换阶数的扫描方法研究也成为了研究的重点。
发明内容
本发明的目的是提供基于单参数的分数阶Fourier逆变换分级扫描方法,适用于采用单参数四项加权分数傅里叶变换(SP-4-WFRFT)的通信系统。本发明针对基于SP-4-WFRFT的数字通信系统特性,在变换阶数未知的条件下,为快速地恢复数据信号(数据信号可是调制信号、基带序列、扩频序列等)或估计变换阶数,在逆变换过程中,建立一种分级扫描方法。
采用的技术方案是:
基于单参数的分数阶Fourier逆变换分级扫描方法:
在接收端,设定信号经过前端处理后的基带信号为:
r(n)=Fα(s(n))+λ0(n) (1)。
其中,λ0(n)为高斯白噪声项,s(n)为基带信号,α为变换阶数,Fα(·)为α阶的SP-4-WFRFT变换,具体过程为:
Fα(s(n))=ω0(α)s(n)+ω1(α)S(n)+ω2(α)s(-n)+ω3(α)S(-n) (2)。
式(2)中,s(n)、S(n)、s(-n)和S(-n)这4种“态函数”是s(n)分别作0、1、2、3次傅里叶变换得到的结果。加权系数ωl(α)的定义为式(3)。
结合DFT的公式,式(2)可等效为:
由于傅里叶变换具有周期为4的特性,为此随着α的变化,加权系数ωl(α)也呈现周期为4的变化,Fα(·)也满足周期特性,且α的取值主周期为[04)。
考虑r(n)为SP-4-WFRFT变换后的信号,利用SP-4-WFRFT变换的旋转可加性,为恢复数据信号或估计变换阶数α需对接收信号进行SP-4-WFRFT逆变换,也等效为-β阶的SP-4-WFRFT变换:
式中,当β与α相等时,接收信号可得:
r”(n)=F0(S(n))+λ'0(n)=S(n)+λ'0(n) (6)。
由于SP-4-WFRFT变换的保范性,高斯白噪声经过变换前后的统计特性不变,故可认为λ0(n)与λ'0(n)有相同的影响,为此逆变换处理得到的信号r'(n)中除了噪声影响外,即可得到数据信号。
然而,在变换阶数α未知的条件下,欲满足式(6)成立,即:β=α的要求,则通过以[04)为主周期对β扫描的方式,达到对α扫描的目的。
设定逆变换的变换阶数与正变换的变换阶数间差值为Δα=α-β,结合变换阶数扫描间隔与误码率的关系,Δα越小则逆变换处理后信号的误码率越小,也即r'(n)与真实信号的误差越小。然而,Δα越小,在0~4区间需要扫描的次数就越多,例如:当Δα=0.01时,需要扫描的次数为400次。当Δα=0.0078125时,需要扫描的次数为512次。
考虑逆变换过程中,在变换阶数α未知的条件下,扫描次数较大、扫描速度较慢等问题,针对基于SP-4-WFRFT的数字通信系统特性,为快速地恢复数据信号(数据信号可是调制信号、基带序列、扩频序列等)或估计变换阶数α,建立基于单参数的分数阶傅里叶逆变换分级扫描方法,原理如下:
定义分级终值常量M,其是分级的最大限定,M是用户依据精度需求自定义的一个正整数,数值大于等于1。进而定义分级级数变量m,m∈[1,2ΛM],随着m从1到M递增。定义变量j为扫描变量。从第一级到第M级对变换阶数进行扫描,每一级的扫描次数Lm为:
从第一级到第M级对变换阶数进行扫描,每一级的扫描优先级为m,每一级的扫描敏感度为21-m,遍历扫描后总体的扫描敏感度为21-M,遍历扫描后总体的扫描次数为式(10)。
在信号处理领域,为了更好地描述信号特性,出现了分数阶傅里叶变换(Fractional Fourier Transform,FRFT),其应用已从最初的光学领域扩展到其他领域,近年已在通信系统的信号检测、信道估计、同步算法等领域得到广泛应用。由于分数阶傅里叶变换的变换阶数敏感度较高,而为了能够正确地恢复原始数据,对于逆变换的变换阶数要求较为苛刻,尤其是在通信系统应用中存在误差或影响条件时。为此,对于分数阶傅里叶逆变换过程中,变换阶数的扫描方法研究也成为了研究的重点。本发明适用于采用单参数四项加权分数傅里叶变换(SP-4-WFRFT)的通信系统。针对基于SP-4-WFRFT的数字通信系统特性,在变换阶数未知的条件下,为快速地恢复数据信号(数据信号可是调制信号、基带序列、扩频序列等)或估计变换阶数,在逆变换过程中,建立一种分级扫描方法。主要通过定义分级终值常量,利用分级级数变量、扫描变量对变换阶数进行分级,并在每一级进行SP-4-WFRFT逆变换,从而得到所需的待解调数据或基带数据。本发明的优点是:SP-4-WFRFT逆变换过程中,变换阶数的分级扫描具有优先级,且可依据用户需求灵活地设置。
其优点在于:
主要通过定义分级终值常量,利用分级级数变量、扫描变量对变换阶数进行分级,并在每一级进行SP-4-WFRFT逆变换,从而得到所需的待解调数据或基带数据。本发明的优点是:SP-4-WFRFT逆变换过程中,变换阶数的分级扫描具有优先级,且可依据用户需求灵活地设置。
附图说明
图1是本发明方法的处理流程图。
具体实施方式
基于单参数的分数阶Fourier逆变换分级扫描方法:
Step1:用户依据所需的扫描精度21-M,设置分级终值常量M。
Step2:输入前端处理后的基带信号:
r(n)=Fα(s(n))+λ0(n) (11)。
式中r(n)、R(n)、r(-n)和R(-n)这4种“态函数”是r(n)分别作0、1、2、3次傅里叶变换的结果。加权系数ωl(α)为:
Step9:进行第三级扫描处理参数计算,即m=3,优先级继续递减,利用式(7)计算得到扫描次数Lm=2m=8。并计算扫描变量j的遍历范围为j∈[1 Lm-1]=[1 4]。进一步利用式(9)的分级规则,从而得到:
Step10:进行第三级扫描处理,利用对前端处理后的基带信号r(n)进行Lm=8次SP-4-WFRFT逆变换处理,第j次处理如式(16)~(17),从而得到8组结果。其中“态函数”和加权系数规则与Step 4相同。
Step12:进行第m级扫描参数计算,随着m的增大优先级递减,利用式(7)计算得到扫描次数Lm=2m。并计算扫描变量j的遍历范围为j∈[1 Lm-1]。进一步利用式(9)的分级规则,从而得到:
Step13:进行第m级扫描处理,利用对前端处理后的基带信号r(n)进行2m次SP-4-WFRFT逆变换处理,第j次处理如式(19)~(20),从而得到2m组结果。其中“态函数”和加权系数规则与Step 4相同。
当扫描的或与α近似相等时,在此设与α近似相等,则SP-4-WFRFT逆变换处理可得式(21),又由于λ0(n)与λ'0(n)有相同的影响,为此SP-4-WFRFT逆变换处理得到的信号中除了噪声影响外,即可得到所需数据信号。
Step15:重复Step 12~Step 14,直至m=M后终止。
Claims (1)
1.基于单参数的分数阶Fourier逆变换分级扫描方法,其特征在于包括下列步骤:
定义分级终值常量M,其是分级的最大限定,M是用户依据精度需求自定义的一个正整数,数值大于等于1;进而定义分级级数变量m,m∈[1,2ΛM],随着m从1到M递增;定义变量j为扫描变量;从第一级到第M级对变换阶数进行扫描,每一级的扫描次数Lm为:
从第一级到第M级对变换阶数进行扫描,每一级的扫描优先级为m,每一级的扫描敏感度为21-m,遍历扫描后总体的扫描敏感度为21-M,遍历扫描后总体的扫描次数为式(10);
当Step1:用户依据所需的扫描精度21-M,设置分级终值常量M;
Step2:输入前端处理后的基带信号:
r(n)=Fα(s(n))+λ0(n) (11);
其中,λ0(n)为高斯白噪声项,s(n)为基带信号,α为变换阶数,Fα(·)为α阶的SP-4-WFRFT变换;
式中r(n)、R(n)、r(-n)和R(-n)这4种“态函数”是r(n)分别作0、1、2、3次傅里叶变换的结果;加权系数ωl(α)为:
Step9:进行第三级扫描处理参数计算,即m=3,优先级继续递减,利用式(7)计算得到扫描次数Lm=2m=8;并计算扫描变量j的遍历范围为j∈[1 Lm-1]=[14];进一步利用式(9)的分级规则,从而得到:
Step10:进行第三级扫描处理,利用对前端处理后的基带信号r(n)进行Lm=8次SP-4-WFRFT逆变换处理,第j次处理如式(16)~(17),从而得到8组结果;其中“态函数”和加权系数规则与Step 4相同;
Step12:进行第m级扫描参数计算,随着m的增大优先级递减,利用式(7)计算得到扫描次数Lm=2m;并计算扫描变量j的遍历范围为j∈[1 Lm-1];进一步利用式(9)的分级规则,从而得到:
Step13:进行第m级扫描处理,利用对前端处理后的基带信号r(n)进行2m次SP-4-WFRFT逆变换处理,第j次处理如式(19)~(20),从而得到2m组结果;其中“态函数”和加权系数规则与Step4相同;
当扫描的或与α近似相等时,在此设与α近似相等,则SP-4-WFRFT逆变换处理可得式(21),又由于λ0(n)与λ'0(n)有相同的影响,为此SP-4-WFRFT逆变换处理得到的信号中除了噪声影响外,即可得到所需数据信号;
Step15:重复Step12~Step14,直至m=M后终止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810636384.4A CN108920419B (zh) | 2018-06-20 | 2018-06-20 | 基于单参数的分数阶Fourier逆变换分级扫描方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810636384.4A CN108920419B (zh) | 2018-06-20 | 2018-06-20 | 基于单参数的分数阶Fourier逆变换分级扫描方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108920419A CN108920419A (zh) | 2018-11-30 |
CN108920419B true CN108920419B (zh) | 2022-03-01 |
Family
ID=64420128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810636384.4A Active CN108920419B (zh) | 2018-06-20 | 2018-06-20 | 基于单参数的分数阶Fourier逆变换分级扫描方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108920419B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109391292B (zh) * | 2018-12-20 | 2021-07-16 | 哈尔滨工业大学 | 加权分数傅里叶变换域双时隙分集与复用的协同传输方法 |
CN110138539B (zh) * | 2019-05-24 | 2022-06-07 | 沈阳理工大学 | 一种基于混沌映射参数池的wfrft安全通信方法 |
CN110933011B (zh) * | 2019-12-11 | 2022-06-10 | 沈阳理工大学 | 一种wfrft信号的高分辨盲扫描方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101771530A (zh) * | 2010-01-26 | 2010-07-07 | 哈尔滨工业大学 | 一种利用一般信号序列隐藏加权信号的隐蔽通信方法 |
CN102201914A (zh) * | 2011-07-14 | 2011-09-28 | 哈尔滨工业大学 | 一种基于多维单/多参数四项加权分数傅里叶变换的保密通信方法 |
CN103441822A (zh) * | 2013-09-16 | 2013-12-11 | 哈尔滨工业大学 | 基于加权分数傅立叶变换及能量扩展变换的信号调制方法 |
CN103516658A (zh) * | 2013-09-13 | 2014-01-15 | 天津理工大学 | 一种降低光正交频分复用系统中高峰均值功率比的方法 |
CN103746803A (zh) * | 2014-01-17 | 2014-04-23 | 哈尔滨工业大学 | 基于加权分数傅里叶变换的多序列联合扩频保密通信方法 |
CN105530217A (zh) * | 2016-01-26 | 2016-04-27 | 哈尔滨工业大学 | 基于加权分数傅里叶变换的gfdm系统的信号发射和接收方法 |
CN106102045A (zh) * | 2016-05-25 | 2016-11-09 | 赵京磊 | 一种采用分数阶傅里叶变换的保密通信方法 |
CN107707501A (zh) * | 2017-10-13 | 2018-02-16 | 哈尔滨工业大学 | 基于时间交织的多矢量wfrft混合载波并行传输方法 |
-
2018
- 2018-06-20 CN CN201810636384.4A patent/CN108920419B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101771530A (zh) * | 2010-01-26 | 2010-07-07 | 哈尔滨工业大学 | 一种利用一般信号序列隐藏加权信号的隐蔽通信方法 |
CN102201914A (zh) * | 2011-07-14 | 2011-09-28 | 哈尔滨工业大学 | 一种基于多维单/多参数四项加权分数傅里叶变换的保密通信方法 |
CN103516658A (zh) * | 2013-09-13 | 2014-01-15 | 天津理工大学 | 一种降低光正交频分复用系统中高峰均值功率比的方法 |
CN103441822A (zh) * | 2013-09-16 | 2013-12-11 | 哈尔滨工业大学 | 基于加权分数傅立叶变换及能量扩展变换的信号调制方法 |
CN103746803A (zh) * | 2014-01-17 | 2014-04-23 | 哈尔滨工业大学 | 基于加权分数傅里叶变换的多序列联合扩频保密通信方法 |
CN105530217A (zh) * | 2016-01-26 | 2016-04-27 | 哈尔滨工业大学 | 基于加权分数傅里叶变换的gfdm系统的信号发射和接收方法 |
CN106102045A (zh) * | 2016-05-25 | 2016-11-09 | 赵京磊 | 一种采用分数阶傅里叶变换的保密通信方法 |
CN107707501A (zh) * | 2017-10-13 | 2018-02-16 | 哈尔滨工业大学 | 基于时间交织的多矢量wfrft混合载波并行传输方法 |
Non-Patent Citations (3)
Title |
---|
The multiple-parameter fractional Fourier transform;J. Lang 等;《Science China Information Sciences》;20080831;第1-15页 * |
加权类分数傅立叶变换及其在通信系统中的应用;梅林;《博士学位论文》;20130320;第1-10页 * |
加权类分数傅里叶变换信号的抗截获性能分析;冯红;《硕士学位论文》;20130320;第1-12页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108920419A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108920419B (zh) | 基于单参数的分数阶Fourier逆变换分级扫描方法 | |
US7133808B2 (en) | System and method for analyzing data streams | |
Carrillo et al. | Iterative hard thresholding for compressed sensing with partially known support | |
Yang | Random sampling and reconstruction in multiply generated shift-invariant spaces | |
WO2017088587A1 (zh) | 一种数据处理方法及装置 | |
Artyushenko et al. | The demodulation signal under the influence of additive and multiplicative non-Gaussian noise | |
Tanajura et al. | A sequential data assimilation method based on the properties of a diffusion-type process | |
CN115118334B (zh) | 卫星通信帧头捕获方法、装置、通信设备和存储介质 | |
CN108846407B (zh) | 基于独立成分高序不确定脑网络的磁共振影像分类方法 | |
CN108196221B (zh) | 一种基于多基线干涉仪角度模糊区间的去野值方法 | |
CN113705335B (zh) | 一种低频稀疏线谱信号时频超分辨率提取方法 | |
Solomon et al. | Signal denoising using Ramanujan periodic transform | |
Hage et al. | Sparse signal recovery via kalman-filter-based ℓ1 minimization | |
CN113972951A (zh) | 一种低信噪比高动态信号载波多普勒频偏的梳状搜索方法 | |
Yang et al. | A greedy pursuit algorithm for arbitrary block sparse signal recovery | |
Carrillo et al. | Bayesian compressed sensing using generalized Cauchy priors | |
CN107506779B (zh) | 一种植物茎干含水量的估算方法及系统 | |
CN111462479A (zh) | 一种基于傅里叶-递归神经网络的交通流预测方法 | |
Artyushenko et al. | Analysis of Moment and Cumulant Description of non-Gaussian Random Processes, Signals and Noise | |
CN115015731B (zh) | 器件击穿电压确定方法、装置、芯片、电子设备及介质 | |
Yang et al. | Arbitrary block-sparse signal reconstruction based on incomplete single measurement vector | |
CN110826212B (zh) | 一种基于直方图筛选及最小二乘拟合的序列周期估计方法 | |
Obidin et al. | Wavelets and adaptive thresholding | |
APANTAKU et al. | A MODIFIED GENERALIZED CHAIN RATIO IN REGRESSION ESTIMATOR | |
Niedźwiecki et al. | Adaptive notch filters based on combined parametric and nonparametric approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |