CN102624660B - 基于四项加权分数傅里叶变换的窄带干扰抑制的方法 - Google Patents
基于四项加权分数傅里叶变换的窄带干扰抑制的方法 Download PDFInfo
- Publication number
- CN102624660B CN102624660B CN201210120739.7A CN201210120739A CN102624660B CN 102624660 B CN102624660 B CN 102624660B CN 201210120739 A CN201210120739 A CN 201210120739A CN 102624660 B CN102624660 B CN 102624660B
- Authority
- CN
- China
- Prior art keywords
- signal
- point
- conversion
- fourier transform
- wfrft
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000001629 suppression Effects 0.000 title claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 10
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 238000006243 chemical reaction Methods 0.000 claims description 30
- 239000011159 matrix material Substances 0.000 claims description 21
- 230000011664 signaling Effects 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000004891 communication Methods 0.000 abstract description 4
- 238000005070 sampling Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
Landscapes
- Noise Elimination (AREA)
Abstract
基于四项加权分数傅里叶变换的窄带干扰抑制的方法,它涉及一种通信系统的干扰抑制方法,本发明为了解决现有的窄带抑制技术无法实现在分数域上直接进行窄带干扰抑制的问题。本发明的主要步骤为:α阶4-WFRFT变换;确定强干扰点的位置;系数替换:将第N-i+1点的变换前的信号点乘上相应的加权系数得到si,将已有的加权系数A=[A0,A1,A2,A3]经过一定的变换可得到相应的加权系数,利用置换矩阵H和矩阵A的乘积得到相应的加权系数A′=(HAH)H=AH,将加权系数A′和第i个点的被加权函数的乘积,作为新的分数域的点去代替受强干扰的点,变换系数替换的整个过程由公式si=AH[AH(AAH)-1sN-i+1]来表达;-α阶4-WFRFT变换。本方法用于在分数域上进行窄带干扰抑制。
Description
技术领域
本发明涉及一种通信系统的干扰抑制方法,具体的说是一种基于四项加权分数傅立叶变换的窄带抑制方法。
背景技术
窄带干扰是实际通信系统经常遇到的问题,有效地窄带干扰抑制方法能够改善宽带通信系统的性能。已有的窄带抑制技术分为时域抑制,变换域抑制和码辅助抑制。时域NBI抑制技术主要是利用窄带信号强相关性和宽带信号非强相关性进行抑制,其抗干扰性能能够达到最佳而被人们所重视,并得到广泛的研究。时域抑制技术实现简单,通过自适应算法产生最优权值,可有效抑制平稳窄带干扰,但多数方法需要长时间的迭代才能达到稳定,无法跟踪上快变的干扰,抗干扰性能仍有待改进。变换域进行干扰抑制的方法由于简单有效,并具有很多优点,被认为是极具潜力的抗干扰策略,由于在时域复杂的卷积(滤波)过程可以通过频域简单的乘积完成,而且频率还能够实现时域无法实现的理想滤波器,如矩形滤波器。总之,目前已有的窄带抑制技术无法实现在分数域上直接进行窄带干扰抑制的问题。目前还没有发现在四项加权分数傅立叶变换域抑制窄带干扰的技术出现。
发明内容
本发明的目的是提供一种基于四项加权分数傅里叶变换的窄带干扰抑制的方法,以解决现有的窄带抑制技术无法实现在分数域上直接进行窄带干扰抑制的问题。
本发明为解决上述技术问题采取的技术方案是:
本发明所述的基于四项加权分数傅里叶变换的窄带干扰抑制的方法的具体过程为:
步骤一、α阶4-WFRFT变换:对长度为N的数字信号X0={x1,x2,...,xN}进行α阶4-WFRFT变换,变换到分数域信号S={s1,s2,...sN},其中si表示第i个分数域采样点信号;
步骤二、确定强干扰点的位置:
设定分数域受到强干扰的点的位置为第i处,设i处信号为si+n,其中n是强干扰信号,
A0~A3是加权系数;
步骤三、系数替换:
步骤三(一)、利用矩阵广义逆分离出分数域信号点sN-i+1的变换前信号点 表示矩阵的转置,分离过程如下:
记A=[A0,A1,A2,A3],其中A0~A3是加权系数,由于A是行满秩矩阵,根据矩阵的广义逆可得,A-1=A+=AH(AAH)-1,其中AH表示A的共轭转置,
于是
步骤三(二)、根据sN-i+1和si的表达式,将第N-i+1点的变换前的信号点乘上相应的加权系数得到si,将已有的加权系数A=[A0,A1,A2,A3]经过一定的变换可得到相应的加权系数,定义如下的置换矩阵:
利用置换矩阵H和矩阵A的乘积得到相应的加权系数
步骤三(三)、将步骤三(二)得到的加权系数A′和第i个点的被加权函数的乘积,作为新的分数域的点去代替受强干扰的点;
上述变换系数替换的整个过程由下述公式表达:
si=AH[AH(AAH)-1sN-i+1];
步骤四、-α阶4-WFRFT变换:将消除干扰后的变换域信号,经过4-WFRFT逆变换,得到处理后的时域输出信号。
本发明的有益效果是:
本发明所述方法实现在分数域上直接进行窄带干扰抑制的问题,与传统的干扰置零法相比,本方法处理后的信号的最小均方误差明显减小。
本发明的方法采用以下仿真进行验证:
仿真条件:利用三角波信号,采样范围-5到5,采样周期0.1,三角波形如图3-a。经过α=0.6的4-WFRFT变换得到变换域的波形图,如图3-b;在采样点45的位置上有一个强干扰,如图3-c,干扰幅度是10。利用干扰置零方法的到的波形图是图3-d,这时的均方误差(MSE)是0.0018,而利用本发明的系数替换方法得到的波形是图3-e,这时的均方误差(MSE)是1.7713×10-4。由此可见,本发明方法均方误差相比利用干扰置零方法的均方误差大大降低。从波形匹配上分析,本发明的方法也更具有优势。
附图说明
图1是本发明所述的系数替换法原理示意图;
图2是本发明方法(系数替换法)的实现步骤框图;
图3-a至图3-e是受到窄带干扰后三角波在分数域(α=0.6)的系数替换和干扰置零的对比说明图:更具体的,图3-a是三角波时域信号波形图;图3-b是三角波在分数域α=0.6的波形图,采样点是总数101;图3-c是受窄带干扰的分数域(α=0.6)的波形图,在采样点45的位置上受到强干扰,干扰幅度是10;图3-d是利用干扰置零法得到的波形图,即采用干扰置零后的分数域(α=0.6)波形图;图3-e是经过系数替换法得到的波形图,即采用本发明所述系数替换后的分数域(α=0.6)波形图,从图形上看出受干扰点采用系数替换法替换,处理的效果更好。
具体实施方式
具体实施方式一:如图1~2所示,本实施方式所述的基于四项加权分数傅里叶变换的窄带干扰抑制的方法的具体过程为:
步骤一、α阶4-WFRFT变换:对长度为N的数字信号(X0={x1,x2,...,xN}),进行α阶4-WFRFT变换,变换到分数域信号S={s1,s2,...sN},其中si表示第i个分数域采样点信号;为了叙述方便,这里我们假设X1=(y1,y2,...,yN)为X0经过一次傅里叶变换后的信号,由傅里叶变换的性质,X1经过傅里叶变换后的信号X2=(xN,xN-1,...,x1),同样由傅里叶变换的性质,X2经过傅里叶变换后的信号X3=(yN,yN-1,...,y1);
步骤二、确定强干扰点的位置:
设定分数域受到强干扰的点的位置为第i处,设i处信号为(si+n),其中n是强干扰信号,
A0~A3是加权系数;
确定强干扰点的位置这个步骤基于现有技术即可实现的。
步骤三、系数替换:
步骤三(一)、利用矩阵广义逆分离出分数域信号点sN-i+1的变换前信号点 表示矩阵的转置,分离过程如下:
记A=[A0,A1,A2,A3],其中A0~A3是加权系数,由于A是行满秩矩阵,根据矩阵的广义逆可得,A-1=A+=AH(AAH)-1,其中AH表示A的共轭转置,
于是
步骤三(二)、通过sN-i+1和si的表达式可知,只要将第N-i+1点的变换前的信号点乘上相应的加权系数就可以得到si,这里已有的加权系数是A=[A0,A1,A2,A3],要经过一定的变换,得到相应的加权系数,这里我们定义如下的置换矩阵:
利用置换矩阵H和矩阵A的乘积得到相应的加权系数
步骤三(三)、将步骤三(二)得到的加权系数A′和第i个点的被加权函数的乘积,作为新的分数域的点去代替受强干扰的点;
上述变换系数替换的整个过程由下述公式表达:
si=AH[AH(AAH)-1sN-i+1];
至此,完成了变换系数替换的全过程;
步骤四、-α阶4-WFRFT变换:将消除干扰后的变换域信号,经过4-WFRFT逆变换,得到处理后的时域输出信号。
本发明方法先找到被干扰的信号准确的位置,例如第i个点受到了强干扰,这时可以利用四项加权分数傅里叶变换的被加权函数和加权系数的性质,第i个位置和第N-i+1个位置的被加权函数相同,但是与加权系数乘积的顺序不同,如图1所示。而利用第N-i+1位置的数据,经过变换系数,去替换受强干扰的第i个点。
本发明用到的数学工具是四项加权分数傅立叶变换,它是不同于经典傅里叶变换的一种新的数学变换方法,目前在四项加权分数傅里叶变换的文献中,未发现将其应用到抑制窄带干扰上的研究。
四项加权分数傅里叶变换(4-WFRFT)的定义式是:
Fs α[f](t)=A0(α)f(t)+A1(α)F(t)+A2(α)f(-t)+A3(α)F(-t)
其中Fs α[f](t)表示函数或者信号f(t)的四项加权分数傅里叶变换,A0~A3是加权系数。
具体可表示为:
其中,α的取值范围是[0,1],当α取0时,四项加权分数傅里叶变换的结果是原始信号f(t),是恒等变换,而当α取1时,四项加权分数傅里叶变换的结果是F(t),退化f(t)的傅里叶变换。其中f(t)和F(t)是互为傅里叶变换对,而且f(-t)是f(t)的关于原点对称的函数,即f(-t)是f(t)以原点为中心的反转,同样道理,F(-t)是F(t)以原点为中心的反转。傅里叶变换采用下列定义式:
四项加权分数傅里叶变换的定义式是连续的时间和频率混合表达,其变换结果通常为时频均不带限信号,不能直接在现有的电子系统实现。目前,已经有文献直接利用四项加权分数傅里叶变换处理离散信号。
Fs α[X0](n)=A0(α)X0(n)+A1(α)X1(n)+A2(α)X2(n)+A3(α)X3(n)
其中X0(n)~X3(n)是离散信号X0(n)的0-3阶的DFT,
加权分数傅里叶逆变换:
X0(n)=Fs -α[X0](n)=A0(-α)X0(n)+A1(-α)X1(n)+A2(-α)X2(n)+A3(-α)X3(n)。
实施例:
图1是本发明的原理示意图,信号的长度是N,采样点i是受强干扰的点,而本发明用采样点N-i+1的系数来替换采样点i,图2是本发明的原理实现步骤,基于图2举例如下:
设待检测信号为三角波信号X,采样周期0.1。经过变换域的信号设为S,受干扰的采样点位置是45,干扰幅度是10,变换阶数α=0.6;
一,对数字信号X,信号长度N=101,进行4-WFRFT变换,变换到分数域信号S;
二,估计出强干扰的点,根据已知条件,干扰的采样点是i=45;
三,利用系数替换法,si=AH[AH(AAH)-1sN-i+1],消除受干扰点的窄带干扰;
四,将消除干扰后的变换域信号,经过4-WFRFT逆变换,得到处理后的输出信号。
Claims (1)
1.一种基于四项加权分数傅里叶变换的窄带干扰抑制的方法,其特征在于:所述方法的具体过程为:
步骤一、α阶4-WFRFT变换:对长度为N的数字信号X0={x1,x2,...,xN}进行α阶4-WFRFT变换,变换到分数域信号S={s1,s2,...sN},其中si表示第i个分数域采样点信号;所述α阶4-WFRFT变换是指数字信号X0={x1,x2,...,xN}经过一次傅里叶变换后为信号X1=(y1,y2,...,yN),X1经过傅里叶变换后为信号X2=(xN,xN-1,...,x1),X2经过傅里叶变换后为信号X3=(yN,yN-1,...,y1);
步骤二、确定强干扰点的位置:
设定分数域受到强干扰的点的位置为第i处,设i处信号为si+n,其中n是强干扰信号,
A0~A3是加权系数;yi表示信号X1的第i个采样点信号,yN-i+1表示与yi对应的信号X3的第i个采样点信号;
步骤三、系数替换:
步骤三(一)、利用矩阵广义逆分离出分数域信号点sN-i+1的变换前信号点 表示矩阵的转置,分离过程如下:
记A=[A0,A1,A2,A3],其中A0~A3是加权系数,由于A是行满秩矩阵,根据矩阵的广义逆可得,A-1=A+=AH(AAH)-1,其中AH表示A的共轭转置,
于是
步骤三(二)、根据sN-i+1和si的表达式,将第N-i+1点的变换前的信号点乘上相应的加权系数得到si,将已有的加权系数A=[A0,A1,A2,A3]经过一定的变换可得到相应的加权系数,定义如下的置换矩阵:
利用置换矩阵H和矩阵A的乘积得到相应的加权系数
步骤三(三)、将步骤三(二)得到的加权系数A′和第i个点的被加权函数的乘积,作为新的分数域的点去代替受强干扰的点;
上述变换系数替换的整个过程由下述公式表达:
si=AH[AH(AAH)-1sN-i+1];
步骤四、-α阶4-WFRFT变换:将消除干扰后的变换域信号,经过4-WFRFT逆变换,得到处理后的时域输出信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210120739.7A CN102624660B (zh) | 2012-04-23 | 2012-04-23 | 基于四项加权分数傅里叶变换的窄带干扰抑制的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210120739.7A CN102624660B (zh) | 2012-04-23 | 2012-04-23 | 基于四项加权分数傅里叶变换的窄带干扰抑制的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102624660A CN102624660A (zh) | 2012-08-01 |
CN102624660B true CN102624660B (zh) | 2014-08-27 |
Family
ID=46564348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210120739.7A Expired - Fee Related CN102624660B (zh) | 2012-04-23 | 2012-04-23 | 基于四项加权分数傅里叶变换的窄带干扰抑制的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102624660B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102932311B (zh) * | 2012-11-26 | 2014-12-03 | 哈尔滨工业大学 | 基于四项加权分数傅里叶变换的带宽压缩方法 |
CN103067318A (zh) * | 2012-12-25 | 2013-04-24 | 哈尔滨工业大学 | 基于四项加权分数傅里叶变换加权矩阵的窄带干扰抑制方法 |
CN106559126A (zh) * | 2015-09-29 | 2017-04-05 | 华为技术有限公司 | 信号发送方法、信号接收方法、发射端及接收端 |
CN113630151B (zh) * | 2021-09-18 | 2022-07-19 | 哈尔滨工业大学 | 一种时频联合扩展传输方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101222326A (zh) * | 2008-01-18 | 2008-07-16 | 哈尔滨工业大学 | 一种基于四项加权分数傅立叶变换的通信保密和解密方法 |
CN101232488A (zh) * | 2008-02-25 | 2008-07-30 | 哈尔滨工业大学 | 能降低正交频分复用系统峰均功率比的通信方法 |
CN101237253A (zh) * | 2008-02-29 | 2008-08-06 | 哈尔滨工业大学 | 基于分数阶傅立叶变换的自适应干扰分离的信号收发装置 |
-
2012
- 2012-04-23 CN CN201210120739.7A patent/CN102624660B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101222326A (zh) * | 2008-01-18 | 2008-07-16 | 哈尔滨工业大学 | 一种基于四项加权分数傅立叶变换的通信保密和解密方法 |
CN101232488A (zh) * | 2008-02-25 | 2008-07-30 | 哈尔滨工业大学 | 能降低正交频分复用系统峰均功率比的通信方法 |
CN101237253A (zh) * | 2008-02-29 | 2008-08-06 | 哈尔滨工业大学 | 基于分数阶傅立叶变换的自适应干扰分离的信号收发装置 |
Non-Patent Citations (6)
Title |
---|
General Multifractional Fourier Transform Method Based on the Generalized Permutation Matrix Group;Qiwen Ran等;《IEEE TRANSACTIONS ON SIGNAL PROCESSING》;20050131;第53卷(第1期);第83~98页 * |
Qiwen Ran等.General Multifractional Fourier Transform Method Based on the Generalized Permutation Matrix Group.《IEEE TRANSACTIONS ON SIGNAL PROCESSING》.2005,第53卷(第1期),第83~98页. |
加权类分数傅立叶变换及其在通信系统中的应用;梅林;《中国博士学位论文全文数据库 信息科技辑》;中国学术期刊(光盘版)电子杂志社;20110815(第8期);第24~70页 * |
四项加权分数Fourier变换在通信系统中的应用研究;梅林等;《中国科学:信息科学》;《中国科学》杂志社;20100531;第40卷(第5期);第732~741页 * |
梅林.加权类分数傅立叶变换及其在通信系统中的应用.《中国博士学位论文全文数据库 信息科技辑》.中国学术期刊(光盘版)电子杂志社,2011,(第8期),第24~70页. |
梅林等.四项加权分数Fourier变换在通信系统中的应用研究.《中国科学:信息科学》.《中国科学》杂志社,2010,第40卷(第5期),第732~741页. |
Also Published As
Publication number | Publication date |
---|---|
CN102624660A (zh) | 2012-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103308804B (zh) | 基于快速k-s变换电能质量扰动信号时频参数提取方法 | |
CN102141403B (zh) | 基于小波阈值、中值滤波和均值滤波的实时混合去噪方法 | |
CN102915742B (zh) | 基于低秩与稀疏矩阵分解的单通道无监督语噪分离方法 | |
Zhang et al. | Reduced interference time-frequency representations and sparse reconstruction of undersampled data | |
CN104038181B (zh) | 一种基于nlms算法的自适应滤波器的构建方法 | |
CN102624660B (zh) | 基于四项加权分数傅里叶变换的窄带干扰抑制的方法 | |
CN106446868A (zh) | 一种基于emd与奇异值差分谱的侧信道信号特征提取方法 | |
CN103176946A (zh) | 一种面向块稀疏信号的稀疏分解去噪方法 | |
CN109031415B (zh) | 一种基于深度卷积神经网络的可控震源数据振铃压制方法 | |
CN113821978B (zh) | 基于改进步长lms自适应算法的行波检测方法和系统 | |
CN107729845A (zh) | 一种基于子空间特征值分解的实测频响函数降噪方法 | |
CN109813230A (zh) | 一种基于vmd的光栅莫尔条纹去噪方法 | |
CN102571034B (zh) | 基于随机循环矩阵的模拟压缩感知采样方法及系统 | |
CN104614767A (zh) | 基于分段延拓的时变地震子波相位校正方法 | |
CN102359815A (zh) | 一种基于小波分形组合的爆破振动信号特征提取方法 | |
CN104123462A (zh) | 实多项式求根实现均匀线阵的谱music方法 | |
CN108090270B (zh) | 一种基于形态学滤波和盲源分离的暂态振荡参数识别方法 | |
CN105740762A (zh) | 一种信噪分离的优化方法 | |
CN103078661A (zh) | 一种基于迭代门限的扩频系统干扰抑制方法 | |
CN103915102A (zh) | 一种lfm水声多途信号的噪声抑制方法 | |
CN103714825A (zh) | 基于听觉感知模型的多通道语音增强方法 | |
CN103176947A (zh) | 一种基于信号相关性的多通道信号去噪方法 | |
CN106501819A (zh) | 基于fpga的卫星接收机抗窄带干扰实现方法 | |
CN103577877B (zh) | 一种基于时频分析和bp神经网络的船舶运动预报方法 | |
CN101997788B (zh) | 一种信号恢复的优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140827 |