CN104158777A - 联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法 - Google Patents
联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法 Download PDFInfo
- Publication number
- CN104158777A CN104158777A CN201410432148.2A CN201410432148A CN104158777A CN 104158777 A CN104158777 A CN 104158777A CN 201410432148 A CN201410432148 A CN 201410432148A CN 104158777 A CN104158777 A CN 104158777A
- Authority
- CN
- China
- Prior art keywords
- signal
- interference
- vector
- ofdm
- beam forming
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 77
- 230000001629 suppression Effects 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 51
- 238000005070 sampling Methods 0.000 claims abstract description 27
- 230000000295 complement effect Effects 0.000 claims abstract description 18
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 18
- 230000008030 elimination Effects 0.000 claims description 6
- 238000003379 elimination reaction Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 238000010606 normalization Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 abstract description 10
- 238000012545 processing Methods 0.000 abstract description 2
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 238000004891 communication Methods 0.000 description 11
- 238000001228 spectrum Methods 0.000 description 11
- 238000004088 simulation Methods 0.000 description 6
- 230000000116 mitigating effect Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 101000654471 Mus musculus NAD-dependent protein deacetylase sirtuin-1 Proteins 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000790 scattering method Methods 0.000 description 1
- 101150026756 sir1 gene Proteins 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Noise Elimination (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
一种联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法:将来自阵列天线的射频信号经过射频前端转换为模拟基带信号,现通过A/D转换为数字基带信号,将多个数字基带信号通过正交投影算法计算得到干扰信号正交补空间的投影矩阵,并作用于接收信号矢量消除测距仪脉冲信号的干扰,进一步通过循环波束成型方法获取阵列天线的波束成型权值矢量,用其对消除脉冲干扰后的信号矢量进行波束成型提取OFDM直射径信号;移除循环前缀后,通过FFT运算转换为频域信号,进一步通过4倍频域下采样,再通过信道估计信道特性,并通过均衡器进行信道均衡,再依次通过解调器、解交织器及译码器得到发送比特序列的估计值。本发明可有效克服测距仪脉冲及散射多径信号的干扰。
Description
技术领域
本发明涉及一种测距仪脉冲干扰抑制方法。特别是涉及一种联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法。
背景技术
为保障民航新一代空中交通管理系统安全、可靠、高效的运行,国际民航组织(ICAO)提出了两种地空数据链候选技术方案:L频段数字航空通信系统1(L-DACS1)与L频段数字航空通信系统2(L-DACS2),其中,L-DACS1系统采用多载波正交频分复用(OFDM)传输方案,L-DACS2系统采用单载波高斯最小移频键控(GMSK)传输方案。相对于L-DACS2系统,L-DACS1系统具有更高的频谱效率、更高的传输容量、更适合频率选择性衰落信道传输的优点,因此L-DACS1系统获得民用航空界的广泛关注。与此同时,为解决民用航空甚高频通信频率资源匮乏的问题,2007年世界无线电大会(WRC)批准未来民航地空数据链系统部署在航空无线电导航频段,且以内嵌方式工作于测距仪(DME)频道间(DME系统的频道间隔为1MHz,L-DACS系统部署在相邻DME频道的中央,占用带宽0.5MHz)。由于测距仪信号频谱与L-DACS1系统的频谱存在部分交叠,且DME发射机以大功率、突发脉冲方式工作,因此不可避免出现DME发射的脉冲信号干扰L-DACS1系统OFDM接收机的问题。
文献M Schnell,S Brandes,S Gligorevic,Interference Mitigation for Broadband L-DACS,2008 Digital Avionics Systems Conference(DASC 2008),pp.1–12,St.Paul,MN,USA,October2008和文献S Brandes,M Schnell,Mitigation of Pulsed Interference in OFDM Based Systems,2008 International OFDM Workshop(InOWo 2008),pp.193–197,Hamburg,Germany,August2008首先建立给出测距仪脉冲干扰的信号模型,并仿真研究脉冲熄灭法对OFDM接收机链路差错性能的影响,仿真研究表明脉冲熄灭法可一定程度消除测距仪脉冲干扰,但脉冲熄灭会导致OFDM接收机产生子载波间干扰(ICI);为应用脉冲熄灭法消除测距仪脉冲干扰,文献UEpple,M Schnell,German Aerospace Center(DLR),Overview of Interference Situation andMitigation Techniques for LDACS1,2011 IEEE/AIAA 30th Digital Avionics Systems Conference(DASC),pp:4C5-1-4C5-12,2011,提出两种测距仪脉冲干扰信号出现位置的检测方法,并仿真研究了所提出方法的检测性能;针对脉冲熄灭法产生OFDM接收机子载波间干扰问题,文献S Brandes,M Schnell,Compensation of the Impact of Interference Mitigation by PulseBlanking in OFDM Systems,2009 IEEE Global Telecommunications Conference(GLOBECOM2009),pp.1–6,Honolulu,USA,Nov.2009,提出了迭代ICI干扰重构并补偿的方法,并仿真研究ICI干扰补偿方法的性能,研究发现ICI干扰补偿方法的运算复杂度高;针对OFDM系统存在的非相关随机脉冲干扰,文献G Caire,T Y.Al-Naffouri,and A.Narayanan,Impulse noisecancellation in OFDM:an application of compressed sensing,IEEE Int.Symp.on Inform.Theory,Canada,Jul.2008,将压缩感知理论应用OFDM系统,提出基于凸优化的脉冲干扰消除方法;针对电力线通信系统存在的块稀疏非相关脉冲干扰,文献L Lampe.Bursty impulse noisedetection by compressed sensing,2011 IEEE International Symposium on Power LineCommunications and Its Applications(ISPLC),On page(s):29–34,给出块基于块稀疏压缩感知脉冲干扰消除方法。
OFDM接收机脉冲干扰消除的方法主要包括:脉冲熄灭方法、脉冲限幅方法、联合脉冲熄灭与ICI干扰补偿方法、压缩感知脉冲干扰重构与消除方法,以下分别叙述这些技术的缺陷与不足。
脉冲熄灭方法的基本思想:接收机根据接收信号的幅值,首先确定脉冲干扰信号出现的位置,然后,将存在脉冲干扰的信号样值设置为零,以消除脉冲干扰的影响。在实际应用中,该方法时存在以下两个方面的问题:脉冲干扰信号出现的位置不易确定、脉冲熄灭后导致ICI干扰。在确定脉冲干扰信号出现的位置时,通常采用门限比较的方法,当接收信号幅值超过设定的门限时,接收机认为这些信号样值包含脉冲干扰,由于OFDM信号自身峰均比(PAPR)较高,且经过多径信道传播,接收机通常不易确定脉冲熄灭门限,这将导致接收机比特差错性能的损失。另一个方面,接收机采用脉冲熄灭后,使得接收信号产生子载波间干扰(ICI),子载波间干扰最终也将导致接收机比特差错性能的恶化。
脉冲限幅方法的基本思想:接收机根据接收信号确定OFDM信号的峰值幅度,然后将接收信号中超过OFDM信号峰值幅度的样值均设置为峰值幅度值,以降低脉冲干扰信号的影响。脉冲限幅方法也存在两个方面的问题:脉冲限幅值不易确定、脉冲限幅后信号仍残留脉冲干扰。由于OFDM信号自身峰均比(PAPR)较高,且经过多径信道传播,接收机通常不易确定脉冲限幅值,脉冲限幅值设置不合理将导致接收机比特差错性能的损失;此外,脉冲限幅方法并不能完全消除脉冲干扰,限幅后的信号仍包含部分脉冲干扰,这些残留的脉冲干扰将导致接收机比特差错性能的恶化。
联合脉冲熄灭与ICI干扰补偿方法的基本思想:首先利用脉冲熄灭方法消除脉冲干扰,然后进一步重构子载波间干扰(ICI),最后对脉冲熄灭产生的ICI干扰进行补偿,从而降低脉冲熄灭法产生ICI干扰的影响,提高链路传输的可靠性。该方法存在两个方面的问题,第一脉冲熄灭门限不易确定,第二迭代ICI干扰重构及补偿的运算复杂度较高。
针对OFDM系统存在的稀疏、非相关脉冲干扰,还可使用压缩感知重构并消除脉冲干扰,压缩感知方法的适用范围:脉冲干扰信号具有稀疏特性,脉冲干扰信号为非相关随机干扰。而在L频段数字航空通信系统1(L-DACS1)中,测距仪产生脉冲干扰为强相关脉冲干扰,因此压缩感知方法难于直接应用于OFDM系统测距仪脉冲干扰抑制。
发明内容
本发明所要解决的技术问题是,提供一种可有效克服测距仪脉冲及散射多径信号的干扰,提高L频段数字航空通信系统的链路传输可靠性的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法。
本发明所采用的技术方案是:一种联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,其特征在于,包括:将来自阵列天线的射频信号经过射频前端转换为模拟基带信号,模拟基带信号通过A/D转换为数字基带信号,将多个接收通道输出的数字基带信号通过正交投影算法计算得到干扰信号正交补空间的投影矩阵,将投影矩阵作用于接收信号矢量消除测距仪脉冲信号的干扰,消除脉冲干扰后的信号矢量进一步通过循环波束成型方法获取阵列天线的波束成型权值矢量,利用得到权值矢量对消除脉冲干扰后的信号矢量进行波束成型提取OFDM直射径信号;波束成型器输出的OFDM直射径信号在移除循环前缀后,通过FFT运算转换为频域信号,并进一步通过4倍频域下采样,下采样器输出的信号通过信道估计估计信道特性,并通过均衡器进行信道均衡,均衡器输出信号依次通过解调器、解交织器及译码器得到发送比特序列的估计值。
在模拟基带信号转换为数字信号时,为避免测距仪脉冲信号采样产生频谱混叠干扰OFDM信号的接收,接收机使用了过四倍采样。
所述的正交投影算法计算得到干扰信号正交补空间的投影矩阵,包括:
假设接收信号使用均匀线阵,阵元通道数目为M,阵元间隔为λ/2,则阵列天线接收信号矢量x(n)表示为:
其中,n代表信号采样序号,s(n)代表OFDM信号,aθ代表OFDM信号对应的导向矢量,wk(n)代表第k个干扰信号,dk代表第k个干扰信号对应的导向矢量,K代表干扰源的总数,n(n)代表信道输入的复高斯白噪声矢量;
将(1)式进一步表示为:
x(n)=s(n)+w(n)+n(n) (2)
其中,s(n)=s(n)aθ,
阵列天线接收信号矢量x(n)的协方差矩阵定义为:
Rxx=E{x(n)·xH(n)} (3)
考虑到OFDM信号、测距仪脉冲信号、噪声信号统计独立,则(3)式化简为:
Rxx=Rs+Rw+Rn (4)
其中,Rs代表OFDM信号的协方差矩阵;Rw代表干扰信号的协方差矩阵;Rn代表噪声信号的协方差矩阵;
由于测距仪脉冲信号的强度远大于OFDM信号及噪声信号,因此(4)式进一步化简为:
Rxx=Rw+Rv (5)
其中,Rv=Rs+Rn;
对矩阵Rxx进行特征值分解得到:
其中,代表协方差矩阵Rxx的非零特征值,ei,i=1,2...M代表Rxx的非零特征值对应的特征向量,代表Rxx的非零特征值构成的对角矩阵,U代表Rxx的特征值对应特征列向量构成的酉矩阵,满足UUH=I,定义K个脉冲干扰信号对应的特征值为λi,其中i=1,2,..K,由于测距仪脉冲干扰信号强度远大于OFDM信号与噪声信号,因此λi的取值显著大于特征值所述的K个脉冲干扰信号对应的特征值λi称为主特征值,主特征值λi对应的特征向量张成干扰信号子空间,记为Sw=span{e1,e2,...,eK},其中,ei,i=1,2...K代表协方差矩阵Rxx的K个主特征值λi对应的特征向量;
得到干扰信号子空间后,计算得到干扰信号的正交补空间投影矩阵:
然后,将投影矩阵作用于接收信号矢量以消除测距仪脉冲信号的干扰。
所述的将投影矩阵作用于接收信号矢量以消除测距仪脉冲信号的干扰,包括:
将接收信号矢量x(n)向干扰信号的正交补空间投影得到:
考虑到最后(8)式化简为:
(9)式表明:通过将接收信号矢量向干扰信号正交补空间进行投影可消除信道输入的测距仪脉冲干扰。
所述的通过循环波束成型方法获取阵列天线的波束成型权值矢量,包括:
首先引入以下参量,Ns代表一个OFDM符号的采样点数,Ng代表OFDM符号循环前缀的采样点数,Nu代表除循环前缀外符号的采样点数,且Ns=Nu+Ng;
针对正交投影脉冲干扰消除后输出的信号矢量中,定义两个矢量z(j,i)与z(j,i+Nu):
z(j,i)=[z1[(j-1)·Ns+i],...,zm[(j-1)·Ns+i],...,zM[(j-1)·Ns+i]]T (10)
z(j,i+Nu)=[z1[(j-1)Ns+Nu+i],...,zm[(j-1)Ns+Nu+i],...,zM[(j-1)Ns+Nu+i]]T·ej2πα (11)
其中,zm[(j-1)·Ns+i]其代表M个阵元通道中的第m个通道的第j个OFDM符号的第i个采样样值,其中,i=1,...,Ng,j=1,...,N,m=1,...,M;z(j,i)代表所有M个阵元通道的第j个OFDM符号的第i个采样值构成的信号矢量;z(j,i+Nu)代表在修正接收信号多普勒频偏后所有M个通道的第j个OFDM符号的第i个采样值构成的信号矢量,其中,α代表接收信号的归一化多普勒频偏,α=Δf/FS;Δf代表接收信号多普勒频偏;FS代表OFDM系统子载波间隔;
假设接收机已经完全建立符号定时同步,则接收机循环波束成型的目标:寻找最优权矢量ω使期望信号ωHz(j,i)与参考信号cHz(j,i+Nu)之间的相关值最大化,所述的相关值最大化的数学描述为:
(12)式中的ω和c满足ωHω=cHc=1
其中,i=1,2,...,Ng,j=1,...,N,N代表一帧中OFDM符号总数,进一步引入定义则(12)式进一步表示为:
(13)式中的ω和c满足ωHω=cHc=1
其中,c代表辅助矢量;(13)式进一步表示为在||ω||=||c||=1约束下,ωHRccHRHω最大化的问题,对(13)式采用拉格朗日乘子法求微分并令等式为零可化简得到:
ωHRc=ξmax (14)
其中,ξmax代表R最大奇异值,当最优权矢量ω和辅助矢量c取值为R的最大奇异值ξmax对应的左右奇异值向量时,(13)式取最大值,最优权矢量ω是OFDM信号对应的导向矢量aθ的直接估计,即ω∝aθ。
所述的信道估计是采用频域LS估计和频域插值的方法获得。
本发明的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,可有效克服测距仪脉冲及散射多径信号的干扰,提高L频段数字航空通信系统的链路传输的可靠性。具有如下效果:
1、本方法不存在脉冲熄灭门限设置问题,不存在子载波间干扰(ICI)的问题,可直接利用接收信号矢量构造干扰信号正交补空间的投影矩阵进而作用于接收信号,从而有效抑制测距仪脉冲信号。
2、本方法不存在脉冲限幅门限设置问题,不存在残留脉冲干扰问题,具有更好的链路差错传输性能。
3、本方法不存在脉冲熄灭门限设置问题,也无需迭代重构ICI干扰,因此本发明的运算复杂度较低。
4、本方法不仅可消除随机脉冲干扰,同时巧妙利用OFDM循环前缀的对称特性有效提取出期望信号。本发明应用范围更广,且干扰抑制效果更好。
附图说明
图1是本发明联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法的构成框图;
图2a是OFDM发射信号功率谱图;
图2b是滤波后残留干扰信号功率谱图;
图2c是接收机接收信号功率谱图;
图2d是抑制脉冲干扰后的信号功率谱图;
图3是测距仪脉冲干扰抑制前后时域信号波形比较图(单个DME干扰,信干比=-10dB,噪声功率=0);
图4是循环波束成型的波束图(400次蒙特卡罗试验);
图5是系统差错性能曲线(QPSK调制)。
具体实施方式
下面结合实施例和附图对本发明的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法做出详细说明。
为解决L频段数字航空通信系统1(L-DACS1)以内嵌方式部署在L频段而产生的测距仪强脉冲干扰L频段数字航空通信系统OFDM接收机的问题。本发明提出了的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,该方法利用航空器在航路飞行阶段,航空信道主要以直射信道、散射方式传播,且直射径与散射径电波入射角度扩展较大的特性,提出了联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,接收机首先通过信号正交投影方法消除高强度测距仪脉冲干扰,然后利用OFDM信号循环前缀的对称特性,使用循环波束成型方法提取OFDM直射径信号。利用本发明可有效克服测距仪脉冲及散射多径信号的干扰,提高L频段数字航空通信系统的链路传输的可靠性。
如图1所示,本发明的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,包括:将来自阵列天线的射频信号经过射频前端转换为模拟基带信号,模拟基带信号通过A/D转换为数字基带信号,在模拟基带信号转换为数字信号时,为避免测距仪脉冲信号采样产生频谱混叠干扰OFDM(正交频分复)信号的接收,接收机使用了过四倍采样(L-DACS1接收机标准采样频率为625kHz,四倍过采样频率为2.5MHz)。将多个接收通道输出的数字基带信号通过正交投影算法计算得到干扰信号正交补空间的投影矩阵,将投影矩阵作用于接收信号矢量消除测距仪脉冲信号的干扰,消除脉冲干扰后的信号矢量进一步通过循环波束成型方法获取阵列天线的波束成型权值矢量,利用得到权值矢量对消除脉冲干扰后的信号矢量进行波束成型提取OFDM直射径信号;波束成型器输出的OFDM直射径信号在移除循环前缀后,通过FFT运算转换为频域信号,并进一步通过4倍频域下采样,下采样器输出的信号通过信道估计估计信道特性,并通过均衡器进行信道均衡,均衡器输出信号依次通过解调器、解交织器及译码器得到发送比特序列的估计值。所述的信道估计是采用频域LS估计和频域插值的方法获得。和常规OFDM接收方法相比,本专利提出的接收方案主要体现以下两个方面:信号正交投影脉冲干扰抑制及循环波束成型OFDM直射径信号提取。
所述的正交投影算法计算得到干扰信号正交补空间的投影矩阵,在信号正交投影器中,利用阵列天线输出信号矢量计算接收信号矢量的协方差矩阵,对接收信号矢量的协方差矩阵进行特征值分解后,确定大特征值的个数并确定大特征值对应的特征向量,随后利用大特征值对应的特征向量张成干扰信号子空间,由该干扰子空间构造出正交补空间的投影矩阵,将该矩阵作用于接收信号矢量即可消除测距仪脉冲干扰。具体包括:
假设接收信号使用均匀线阵,阵元通道数目为M,阵元间隔为λ/2,则阵列天线接收信号矢量x(n)表示为:
其中,n代表信号采样序号,s(n)代表OFDM信号,aθ代表OFDM信号对应的导向矢量,wk(n)代表第k个干扰信号,dk代表第k个干扰信号对应的导向矢量,K代表干扰源的总数,n(n)代表信道输入的复高斯白噪声矢量;
将(1)式进一步表示为:
x(n)=s(n)+w(n)+n(n) (2)
其中,s(n)=s(n)aθ,
阵列天线接收信号矢量x(n)的协方差矩阵定义为:
Rxx=E{x(n)·xH(n)} (3)
考虑到OFDM信号、测距仪脉冲信号、噪声信号统计独立,则(3)式化简为:
Rxx=Rs+Rw+Rn (4)
其中,Rs代表OFDM信号的协方差矩阵;Rw代表干扰信号的协方差矩阵;Rn代表噪声信号的协方差矩阵;
在本研究问题中,由于测距仪脉冲信号的强度远大于OFDM信号及噪声信号,因此(4)式进一步化简为:
Rxx=Rw+Rv (5)
其中,Rv=Rs+Rn;
对矩阵Rxx进行特征值分解得到:
其中,代表协方差矩阵Rxx的非零特征值,ei,i=1,2...M代表Rxx的非零特征值对应的特征向量,代表Rxx的非零特征值构成的对角矩阵,U代表Rxx的特征值对应特征列向量构成的酉矩阵,满足UUH=I,定义K个脉冲干扰信号对应的特征值为λi,其中i=1,2,..K,由于测距仪脉冲干扰信号强度远大于OFDM信号与噪声信号,因此λi的取值显著大于特征值所述的K个脉冲干扰信号对应的特征值λi称为主特征值,主特征值λi对应的特征向量张成干扰信号子空间,记为Sw=span{e1,e2,...,eK},其中,ei,i=1,2...K代表协方差矩阵Rxx的K个主特征值λi对应的特征向量;
得到干扰信号子空间后,计算得到干扰信号的正交补空间投影矩阵:
然后,将投影矩阵作用于接收信号矢量以消除测距仪脉冲信号的干扰。
所述的将投影矩阵作用于接收信号矢量以消除测距仪脉冲信号的干扰,包括:
将接收信号矢量x(n)向干扰信号的正交补空间投影得到:
考虑到最后(8)式化简为:
(9)式表明:通过将接收信号矢量向干扰信号正交补空间进行投影可消除信道输入的测距仪脉冲干扰。
所述的通过循环波束成型方法获取阵列天线的波束成型权值矢量,在循环波束成型器中,利用OFDM信号循环前缀的对称特性,以期望信号与参考信号相关值最大化为目标函数寻找最优权矢量。同时,以OFDM信号循环前缀信号矢量跟自身频移信号矢量构造循环相关矩阵。最优权矢量最终确定为循环相关矩阵的最大奇异值对应的左奇异值矢量,将该矢量作用于已消除测距仪脉冲干扰的信号矢量并进一步通过循环波束成型的方法提取出OFDM直射径信号。具体包括:
为方便叙述OFDM接收机循环波束成型方法的工作原理,首先引入以下参量,Ns代表一个OFDM符号的采样点数,Ng代表OFDM符号循环前缀的采样点数,Nu代表除循环前缀外符号的采样点数,且Ns=Nu+Ng;
针对正交投影脉冲干扰消除后输出的信号矢量中,定义两个矢量z(j,i)与z(j,i+Nu):
z(j,i)=[z1[(j-1)·Ns+i],...,zm[(j-1)·Ns+i],...,zM[(j-1)·Ns+i]]T (10)
z(j,i+Nu)=[z1[(j-1)Ns+Nu+i],...,zm[(j-1)Ns+Nu+i],...,zM[(j-1)Ns+Nu+i]]T·ej2πα (11)
其中,zm[(j-1)·Ns+i]其代表M个阵元通道中的第m个通道的第j个OFDM符号的第i个采样样值,其中,i=1,...,Ng,j=1,...,N,m=1,...,M;z(j,i)代表所有M个阵元通道的第j个OFDM符号的第i个采样值构成的信号矢量;z(j,i+Nu)代表在修正接收信号多普勒频偏后所有M个通道的第j个OFDM符号的第i个采样值构成的信号矢量,其中,α代表接收信号的归一化多普勒频偏,α=Δf/FS;Δf代表接收信号多普勒频偏;FS代表OFDM系统子载波间隔;
假设接收机已经完全建立符号定时同步,则接收机循环波束成型的目标:寻找最优权矢量ω使期望信号ωHz(j,i)与参考信号cHz(j,i+Nu)之间的相关值最大化,所述的相关值最大化的数学描述为:
(12)式中的ω和c满足ωHω=cHc=1
其中,i=1,2,...,Ng,j=1,...,N,N代表一帧中OFDM符号总数,进一步引入定义则(12)式进一步表示为:
(13)式中的ω和c满足ωHω=cHc=1
其中,c代表辅助矢量;(13)式进一步表示为在||ω||=||c||=1约束下,ωHRccHRHω最大化的问题,对(13)式采用拉格朗日乘子法求微分并令等式为零可化简得到:
ωHRc=ξmax (14)
其中,ξmax代表R最大奇异值,当最优权矢量ω和辅助矢量c取值为R的最大奇异值ξmax对应的左右奇异值向量时,(13)式取最大值,最优权矢量ω是OFDM信号对应的导向矢量aθ的直接估计,即ω∝aθ。
利用得到权值矢量对消除脉冲干扰后的信号矢量进行波束成型,即可提取OFDM直射径信号。
循环波束成型器输出的信号由循环前缀移除器移除循环前缀后,通过FFT变换器转换为频域信号,并进一步通过4倍频域下采样,下采样器输出的信号送入信道估计器,信道估计器通过频域LS估计和频域插值算法估计信道特性,进而由信道均衡器通过迫零均衡算法对信道均衡处理,均衡器输出的信号送入信号解调器解调输出,又通过解交织器完成解交织处理,最后由信道译码器恢复出原始发送信息比特的估计值。
下面结合附图说明本发明联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法的效果。
图2a、图2b、图2c、图2d显示给出了脉冲干扰抑制前后信号功率谱的比较曲线。每个图的横坐标代表频率(MHz),纵坐标代表信号功率(dBw)。其中,图2a显示OFDM发射信号的功率谱(发射信号功率为1),由图2a可观测到:OFDM信号主要频率分量位于-250KHz~+250KHz,OFDM通频带内信号功率谱取值为-30dBw。图2b显示测距仪脉冲信号经等效抗混叠滤波器后残留干扰信号的功率谱(测距仪信号载波偏置500KHz,信干比SIR=-10dB),由图2b可观测到:经过接收机抗混叠滤波器后残留干扰信号的主要频率成分位于+250KHz左右,且残留干扰信号仍具有较强的信号功率。图2c显示给出接收机接收信号的功率谱(第1接收通道,测距仪载波偏置+500KHz,SIR=-10dB,噪声功率为0),由图2c可观测到:在-250KHz~+150KHz频率范围内,OFDM信号功率谱取值为-30dBw,而在+150KHz~+250KHz频带范围内明显可观测到残留的脉冲干扰信号的频率分量,且干扰信号最高强度取值为-10dBw。图2d显示给出脉冲干扰抑制后信号的功率谱(第1输出通道),由图2d可观测到:在-250KHz~+250KHz的频率范围内,OFDM信号功率谱取值为-30dBw,而在150KHz~250KHz频率范围内观测不到测距仪脉冲干扰信号分量。图2a~图2d比较表明:正交投影方法可显著抑制测距仪脉冲干扰。
图3显示给出正交投影算法抑制测距仪干扰前后时域波形比较,图中横坐标代表时间,纵坐标表示信号幅值(实部),图3中标有“-.-”的曲线代表脉冲干扰抑前时域信号波形(第1输入通道),标有“—”的曲线代表脉冲干扰抑制后时域信号波形(第1输出通道)。曲线比较表明:利用阵列天线正交投影方法可显著抑制测距仪脉冲干扰。
图4显示给出了阵列天线循环波束成型算法的波束图,图4横坐标代表信号的到达角(度),纵坐标表示归一化的波束增益(dB)。仿真试验中,OFDM信号直射径来向20度,散射路径来向80度,信噪比为6dB,DME干扰信号来向60度,信干比为-20dB,DME载波偏置为500KHz。图4可观测到:(1)阵列天线的主波束指向OFDM直射路径方向;(2)阵列天线在测距仪脉冲信号来向方向60度形成较深零陷,信号衰减达到-45dB;(3)多次仿真显示阵列天线波束成型后主瓣突出,零陷方向稳定。
图5显示给出L-DACS1系统的比特差错性能曲线,图5中横坐标代表信噪比,纵坐标表示比特差错概率。图5中标有“▲”的曲线代表无DME干扰时系统的比特差错性能曲线;标有“■”曲线代表存在单个DME干扰时的比特差错性能曲线(DME载波偏置500KHz,信干比SIR=-3dB);标有“●”的曲线代表存在两个DME干扰时的比特差错性能曲线(DME载波偏置分别为500KHz与-500KHz,SIR1=-10dB,SIR2=-5dB)。曲线比较表明:(1)存在单个测距仪脉冲干扰时,本发明提出的方法可完成消除测距仪脉冲干扰,系统的比特差错性能与不存在测距仪脉冲干扰时系统的比特差错性能曲线基本相同;(2)当系统存在两个测距仪脉冲干扰时,所提出的方法仍可获得满意的比特差错性能。
Claims (6)
1.一种联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,其特征在于,包括:将来自阵列天线的射频信号经过射频前端转换为模拟基带信号,模拟基带信号通过A/D转换为数字基带信号,将多个接收通道输出的数字基带信号通过正交投影算法计算得到干扰信号正交补空间的投影矩阵,将投影矩阵作用于接收信号矢量消除测距仪脉冲信号的干扰,消除脉冲干扰后的信号矢量进一步通过循环波束成型方法获取阵列天线的波束成型权值矢量,利用得到权值矢量对消除脉冲干扰后的信号矢量进行波束成型提取OFDM直射径信号;波束成型器输出的OFDM直射径信号在移除循环前缀后,通过FFT运算转换为频域信号,并进一步通过4倍频域下采样,下采样器输出的信号通过信道估计估计信道特性,并通过均衡器进行信道均衡,均衡器输出信号依次通过解调器、解交织器及译码器得到发送比特序列的估计值。
2.根据权利要求1所述的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,其特征在于,在模拟基带信号转换为数字信号时,为避免测距仪脉冲信号采样产生频谱混叠干扰OFDM信号的接收,接收机使用了过四倍采样。
3.根据权利要求1所述的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,其特征在于,所述的正交投影算法计算得到干扰信号正交补空间的投影矩阵,包括:
假设接收信号使用均匀线阵,阵元通道数目为M,阵元间隔为λ/2,则阵列天线接收信号矢量x(n)表示为:
其中,n代表信号采样序号,s(n)代表OFDM信号,aθ代表OFDM信号对应的导向矢量,wk(n)代表第k个干扰信号,dk代表第k个干扰信号对应的导向矢量,K代表干扰源的总数,n(n)代表信道输入的复高斯白噪声矢量;
将(1)式进一步表示为:
x(n)=s(n)+w(n)+n(n) (2)
其中,s(n)=s(n)aθ,
阵列天线接收信号矢量x(n)的协方差矩阵定义为:
Rxx=E{x(n)·xH(n)} (3)
考虑到OFDM信号、测距仪脉冲信号、噪声信号统计独立,则(3)式化简为:
Rxx=Rs+Rw+Rn (4)
其中,Rs代表OFDM信号的协方差矩阵;Rw代表干扰信号的协方差矩阵;Rn代表噪声信号的协方差矩阵;
由于测距仪脉冲信号的强度远大于OFDM信号及噪声信号,因此(4)式进一步化简为:
Rxx=Rw+Rv (5)
其中,Rv=Rs+Rn;
对矩阵Rxx进行特征值分解得到:
其中,代表协方差矩阵Rxx的非零特征值,ei,i=1,2...M代表Rxx的非零特征值对应的特征向量,代表Rxx的非零特征值构成的对角矩阵,U代表Rxx的特征值对应特征列向量构成的酉矩阵,满足UUH=I,定义K个脉冲干扰信号对应的特征值为λi,其中i=1,2,..K,由于测距仪脉冲干扰信号强度远大于OFDM信号与噪声信号,因此λi的取值显著大于特征值所述的K个脉冲干扰信号对应的特征值λi称为主特征值,主特征值λi对应的特征向量张成干扰信号子空间,记为Sw=span{e1,e2,...,eK},其中,ei,i=1,2...K代表协方差矩阵Rxx的K个主特征值λi对应的特征向量;
得到干扰信号子空间后,计算得到干扰信号的正交补空间投影矩阵:
然后,将投影矩阵作用于接收信号矢量以消除测距仪脉冲信号的干扰。
4.根据权利要求1或3所述的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,其特征在于,所述的将投影矩阵作用于接收信号矢量以消除测距仪脉冲信号的干扰,包括:
将接收信号矢量x(n)向干扰信号的正交补空间投影得到:
考虑到最后(8)式化简为:
(9)式表明:通过将接收信号矢量向干扰信号正交补空间进行投影可消除信道输入的测距仪脉冲干扰。
5.根据权利要求1所述的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,其特征在于,所述的通过循环波束成型方法获取阵列天线的波束成型权值矢量,包括:
首先引入以下参量,Ns代表一个OFDM符号的采样点数,Ng代表OFDM符号循环前缀的采样点数,Nu代表除循环前缀外符号的采样点数,且Ns=Nu+Ng;
针对正交投影脉冲干扰消除后输出的信号矢量中,定义两个矢量z(j,i)与z(j,i+Nu):
z(j,i)=[z1[(j-1)·Ns+i],...,zm[(j-1)·Ns+i],...,zM[(j-1)·Ns+i]]T (10)
z(j,i+Nu)=[z1[(j-1)Ns+Nu+i],...,zm[(j-1)Ns+Nu+i],...,zM[(j-1)Ns+Nu+i]]T·ej2πα (11)
其中,zm[(j-1)·Ns+i]其代表M个阵元通道中的第m个通道的第j个OFDM符号的第i个采样样值,其中,i=1,...,Ng,j=1,...,N,m=1,...,M;z(j,i)代表所有M个阵元通道的第j个OFDM符号的第i个采样值构成的信号矢量;z(j,i+Nu)代表在修正接收信号多普勒频偏后所有M个通道的第j个OFDM符号的第i个采样值构成的信号矢量,其中,α代表接收信号的归一化多普勒频偏,α=Δf/FS;Δf代表接收信号多普勒频偏;FS代表OFDM系统子载波间隔;
假设接收机已经完全建立符号定时同步,则接收机循环波束成型的目标:寻找最优权矢量ω使期望信号ωHz(j,i)与参考信号cHz(j,i+Nu)之间的相关值最大化,所述的相关值最大化的数学描述为:
(12)式中的ω和c满足ωHω=cHc=1
其中,i=1,2,...,Ng,j=1,...,N,N代表一帧中OFDM符号总数,进一步引入定义则(12)式进一步表示为:
(13)式中的ω和c满足ωHω=cHc=1
其中,c代表辅助矢量;(13)式进一步表示为在||ω||=||c||=1约束下,ωHRccHRHω最大化的问题,对(13)式采用拉格朗日乘子法求微分并令等式为零可化简得到:
ωHRc=ξmax (14)
其中,ξmax代表R最大奇异值,当最优权矢量ω和辅助矢量c取值为R的最大奇异值ξmax对应的左右奇异值向量时,(13)式取最大值,最优权矢量ω是OFDM信号对应的导向矢量aθ的直接估计,即ω∝aθ。
6.根据权利要求1所述的联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法,其特征在于,所述的信道估计是采用频域LS估计和频域插值的方法获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410432148.2A CN104158777B (zh) | 2014-08-28 | 2014-08-28 | 联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410432148.2A CN104158777B (zh) | 2014-08-28 | 2014-08-28 | 联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104158777A true CN104158777A (zh) | 2014-11-19 |
CN104158777B CN104158777B (zh) | 2017-09-08 |
Family
ID=51884185
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410432148.2A Expired - Fee Related CN104158777B (zh) | 2014-08-28 | 2014-08-28 | 联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104158777B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104601514A (zh) * | 2015-01-27 | 2015-05-06 | 中国民航大学 | 联合子空间投影与clean的干扰抑制方法 |
CN105207966A (zh) * | 2015-08-10 | 2015-12-30 | 中国民航大学 | 一种基于空频编码的压缩感知脉冲干扰抑制系统 |
CN105227258A (zh) * | 2015-11-03 | 2016-01-06 | 重庆大学 | 基于高阶统计量的l-dacs1系统自适应干扰消除方法 |
CN105791198A (zh) * | 2015-01-09 | 2016-07-20 | 想象技术有限公司 | 脉冲噪声抑制 |
CN111082837A (zh) * | 2019-11-25 | 2020-04-28 | 宁波大学 | 基于递归压缩感知的电力线通信系统脉冲噪声抑制方法 |
RU2724923C2 (ru) * | 2018-08-10 | 2020-06-26 | Российская Федерация, от имени которой выступает Министерство промышленности и торговли Российской Федерации (Минпромторг) | Способ скрытного мониторинга радиомолчащих объектов |
CN116054854A (zh) * | 2023-01-13 | 2023-05-02 | 电子科技大学 | 一种低复杂度相控阵自干扰数字域抑制方法 |
CN116112323A (zh) * | 2021-11-10 | 2023-05-12 | 大唐移动通信设备有限公司 | 干扰抑制方法、装置、设备及存储介质 |
CN118400082A (zh) * | 2024-06-20 | 2024-07-26 | 上海金卓科技有限公司 | 一种无线信号同步方法、装置、设备及存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101807979A (zh) * | 2010-03-12 | 2010-08-18 | 北京航空航天大学 | 一种无线通信中基于单天线功率约束的多用户预编码方法 |
CN102694628A (zh) * | 2012-05-07 | 2012-09-26 | 西安电子科技大学 | 基于多用户mimo协同中继系统的干扰抑制方法 |
CN102830387A (zh) * | 2012-08-23 | 2012-12-19 | 北京理工大学 | 一种基于数据预处理的协方差矩阵正交化波束形成方法 |
TWI388138B (zh) * | 2009-06-30 | 2013-03-01 | Jeich Mar | 結合環形陣列天線與多載波直接展頻序列倍頻展頻技術於雙態雷達估測低訊雜比目標來向之接收裝置與估測方法 |
CN103634038A (zh) * | 2013-12-09 | 2014-03-12 | 南京邮电大学 | 基于多天线的联合doa估计和波束形成的多径信号接收方法 |
US8737498B2 (en) * | 2011-05-09 | 2014-05-27 | Electronics And Telecommunications Research Institute | Secondary transmitter included in cognitive radio communication system and communication method of the secondary transmitter |
CN103942449A (zh) * | 2014-05-05 | 2014-07-23 | 北京理工大学 | 一种基于信源个数估计的特征干扰相消波束形成方法 |
-
2014
- 2014-08-28 CN CN201410432148.2A patent/CN104158777B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI388138B (zh) * | 2009-06-30 | 2013-03-01 | Jeich Mar | 結合環形陣列天線與多載波直接展頻序列倍頻展頻技術於雙態雷達估測低訊雜比目標來向之接收裝置與估測方法 |
CN101807979A (zh) * | 2010-03-12 | 2010-08-18 | 北京航空航天大学 | 一种无线通信中基于单天线功率约束的多用户预编码方法 |
US8737498B2 (en) * | 2011-05-09 | 2014-05-27 | Electronics And Telecommunications Research Institute | Secondary transmitter included in cognitive radio communication system and communication method of the secondary transmitter |
CN102694628A (zh) * | 2012-05-07 | 2012-09-26 | 西安电子科技大学 | 基于多用户mimo协同中继系统的干扰抑制方法 |
CN102830387A (zh) * | 2012-08-23 | 2012-12-19 | 北京理工大学 | 一种基于数据预处理的协方差矩阵正交化波束形成方法 |
CN103634038A (zh) * | 2013-12-09 | 2014-03-12 | 南京邮电大学 | 基于多天线的联合doa估计和波束形成的多径信号接收方法 |
CN103942449A (zh) * | 2014-05-05 | 2014-07-23 | 北京理工大学 | 一种基于信源个数估计的特征干扰相消波束形成方法 |
Non-Patent Citations (2)
Title |
---|
CHAN KYU KIM,ET.AL: "《ADAPTIVE BEAMFORMING ALGORITHM FOR OFDM SYSTEMS WITH ANTENNA ARRAYS》", 《IEEE TRANSACTIONS ON CONSUMER ELECTRONICS》 * |
刘海涛等: "《联合正交变换与信号交织的测距仪脉冲干扰抑制方法》", 《航空学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105791198B (zh) * | 2015-01-09 | 2020-09-22 | 想象技术有限公司 | 脉冲噪声抑制方法和ofdm接收机 |
CN105791198A (zh) * | 2015-01-09 | 2016-07-20 | 想象技术有限公司 | 脉冲噪声抑制 |
CN104601514B (zh) * | 2015-01-27 | 2017-09-22 | 中国民航大学 | 联合子空间投影与clean的干扰抑制方法 |
CN104601514A (zh) * | 2015-01-27 | 2015-05-06 | 中国民航大学 | 联合子空间投影与clean的干扰抑制方法 |
CN105207966A (zh) * | 2015-08-10 | 2015-12-30 | 中国民航大学 | 一种基于空频编码的压缩感知脉冲干扰抑制系统 |
CN105227258A (zh) * | 2015-11-03 | 2016-01-06 | 重庆大学 | 基于高阶统计量的l-dacs1系统自适应干扰消除方法 |
CN105227258B (zh) * | 2015-11-03 | 2017-09-15 | 重庆大学 | 基于高阶统计量的l‑dacs1系统自适应干扰消除方法 |
RU2724923C2 (ru) * | 2018-08-10 | 2020-06-26 | Российская Федерация, от имени которой выступает Министерство промышленности и торговли Российской Федерации (Минпромторг) | Способ скрытного мониторинга радиомолчащих объектов |
CN111082837A (zh) * | 2019-11-25 | 2020-04-28 | 宁波大学 | 基于递归压缩感知的电力线通信系统脉冲噪声抑制方法 |
CN111082837B (zh) * | 2019-11-25 | 2021-04-23 | 宁波大学 | 基于递归压缩感知的电力线通信系统脉冲噪声抑制方法 |
CN116112323A (zh) * | 2021-11-10 | 2023-05-12 | 大唐移动通信设备有限公司 | 干扰抑制方法、装置、设备及存储介质 |
CN116112323B (zh) * | 2021-11-10 | 2024-06-07 | 大唐移动通信设备有限公司 | 干扰抑制方法、装置、设备及存储介质 |
CN116054854A (zh) * | 2023-01-13 | 2023-05-02 | 电子科技大学 | 一种低复杂度相控阵自干扰数字域抑制方法 |
CN118400082A (zh) * | 2024-06-20 | 2024-07-26 | 上海金卓科技有限公司 | 一种无线信号同步方法、装置、设备及存储介质 |
CN118400082B (zh) * | 2024-06-20 | 2024-09-13 | 上海金卓科技有限公司 | 一种无线信号同步方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN104158777B (zh) | 2017-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104158777A (zh) | 联合正交投影与循环波束成型的测距仪脉冲干扰抑制方法 | |
CN103117781B (zh) | 一种复杂电磁环境下的天线阵列校准方法及其装置 | |
CN104092634B (zh) | 联合doa估计与主波束成型的测距仪脉冲干扰抑制方法 | |
CN101807977B (zh) | 基于波形特征的空时盲自适应抗干扰方法 | |
US7961824B2 (en) | Receiver apparatus | |
WO2011080644A2 (en) | Radio channel analyzer to determine doppler shifts across multiple frequencies of a wideband signal | |
CN103414677B (zh) | 联合正交变换与信号交织的正交频分复用传输方法 | |
CN102291363A (zh) | 一种用于ofdm系统的信道估计及数据检测方法 | |
CN105227258B (zh) | 基于高阶统计量的l‑dacs1系统自适应干扰消除方法 | |
CN106789760A (zh) | 一种抗干扰的单载波信道估计与均衡方法与装置 | |
CN104539562A (zh) | 基于多输入多输出正交频分复用的宽带短波信道估计方法 | |
CN104601514B (zh) | 联合子空间投影与clean的干扰抑制方法 | |
CN103560991A (zh) | 适用于正交频分复用接收机抑制测距仪脉冲干扰的方法 | |
Kandaurov et al. | Signal-Code Constructions for Wideband Signals Based on M-Cpfsk and Non-Binary Error-Correcting Codes | |
Gajewski | Design of OFDM-based radio communication systems for coast-to-sea and coast-to-air propagation environments | |
Zieliński et al. | Wireless OTFS-based Integrated Sensing and Communication for Moving Vehicle Detection | |
CN107147597B (zh) | 毫米波通信系统中分数间隔频域均衡硬件实现方法 | |
Doosti-Aref et al. | Efficient cooperative multicarrier underwater acoustic communication over the Persian Gulf channel | |
Karpovich et al. | Field tests of a random-padded OTFSM waveform in a joint sensing and communication system | |
CN109274486A (zh) | 基于多径分离的频分双工系统中互易信道增益构建方法 | |
CN100493052C (zh) | 一种弱能量并行pn序列进行信道估计的方法 | |
CN113438047A (zh) | 一种基于高阶统计特性的窄带物联网干扰抑制方法 | |
Zhou et al. | Impact of DME Interference on LDACS Channel Estimation in En-route Phase | |
Liu | BPSK/BOC Modulation Signal System for GPS Satellite Navigation Signals | |
CN103227766A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170908 |
|
CF01 | Termination of patent right due to non-payment of annual fee |