CN105610748A - 一种频率分段的通道均衡方法 - Google Patents
一种频率分段的通道均衡方法 Download PDFInfo
- Publication number
- CN105610748A CN105610748A CN201410670589.6A CN201410670589A CN105610748A CN 105610748 A CN105610748 A CN 105610748A CN 201410670589 A CN201410670589 A CN 201410670589A CN 105610748 A CN105610748 A CN 105610748A
- Authority
- CN
- China
- Prior art keywords
- equalization
- frequency
- channel
- frequency domain
- complex signal
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明属于阵列信号技术,涉及对一定带宽内通道不匹配的均衡方法。通道均衡的步骤如下:第1步、将频域均衡滤波器转换为N点滤波器:选取N为大于等于L的最小二的幂,求频域均衡滤波器的N点频率特性第2步、求出分段均衡的各段滤波器系数;第3步、用eN对待均衡通道的基带复信号xA进行均衡处理,得到均衡后的复信号序列xAE。本发明能够减少通道均衡所需的存储空间,改善通道间对消性能,从而提升阵列处理性能,满足现代通信系统以及相控阵雷达系统等对多通道阵列信号处理的需要。
Description
技术领域
本发明属于阵列信号技术,涉及对一定带宽内通道不匹配的均衡方法。
背景技术
在现代雷达和现代通信中,双通道或多通道阵列信号处理的应用越来越广泛。例如在雷达信号处理中,传统的雷达旁瓣对消技术(SLC),以及自适应波束形成技术(ADBF),就是利用两个或多个通道来对消来自与目标不同角度的干扰;空时自适应处理技术(STAP)则可利用多个通道来对消空时耦合的各类杂波。然而通道之间频率特性的差异将造成对消的剩余,从而造成信干噪比(SINR)或信杂噪比(SCNR)的下降。通道均衡能够校正通道间频率特性的差异,改善通道对消性能,从而提高各种先进处理方法的SINR或SCNR。
关于通道均衡,不少学者对之进行了研究,并发表了相关论文,如王峰等人的《基于傅里叶变换的雷达通道均衡算法性能分析及改进》,李建平的《宽带数字阵列雷达通道均衡算法研究》,冉涌的《宽带数字阵列雷达通道校正技术》,以及赵立艳的《数字阵列雷达中通道均衡技术的研究》等。但以上论文提及的方法都是基于如何构建时域滤波器系数的方法,这种时域处理方法需通过逻辑时序单元实现,很难叠加到其它处理过程(如脉压过程)同时实现。另一方面,即使已知均衡滤波器的频率响应,直接在脉压处理时应用需要极大的存储空间,限制了频域均衡的工程应用。
发明内容
发明创造的目的
本发明的目的是:提出一种频率分段的通道均衡方法,以满足多通道阵列信号处理的需要。
技术方案
本发明的技术方案是:一种频率分段的通道均衡方法,基于已知的通道频域均衡滤波器频率特性e=[e[1],e[2],…,e[P]]T,其中e[p]代表从-fs/2至fs/2内P个不同频点的均衡系数,p=1…P,fs为采样频率;已采集的参考通道与待均衡通道的基带复信号xM和xA,xM=[xM[1],xM[2],…,xM[L]]T为参考通道的复信号,xM[l]是参考通道第l个复采样,l=1…L;xA=[xA[1],xA[2],…,xA[L]]T为待均衡通道的复信号,xM[l]是参考通道第l个复采样,l=1…L;信号所占带宽为BW,其中[·]T表示转置;其特征在于,通道均衡的步骤如下:
第1步、将频域均衡滤波器转换为N点滤波器:选取N为大于等于L的最小二的幂,求频域均衡滤波器的N点频率特性
其中FFT(·,N)表示N点快速傅里叶变换,IFFT(·,N)表示N点快速傅里叶逆变换;
第2步、求出分段均衡的各段滤波器系数:
2.1、根据需要将N点分为K段,每段M个频点,即N=M×K,将变为M×K的矩阵
其中reshape(·,M,K)表示将向量变为M×K的矩阵,并将写成如下的块矩阵:
其中为M×(F1-1)的块矩阵,代表的1~F1-1列,为M×(F2-F1+1)的块矩阵,代表的F1~F2列,为M×(K-F2)的块矩阵,代表的K-F2+1~K列,F1、F2分别表示K个频段中的F1、F2个子频段,F1~F2个子频段覆盖信号带宽BW;
2.2、对进行处理,得到频率分段均衡系数eK:
2.2.1、首先,对求平均,
其中,mean(·)为对列求平均,
2.2.2、随后,将的1~F1-1以及F2+1~K个元素置为1,得到各段滤波器系数为eK,写为下式:
其中,为频率分段处理后需要存储的频域均衡系数。
2.3、根据分段均衡系数eK求得实际均衡时所采用N点频域均衡系数eN:
其中,[·]T为矩阵转置;
第3步、用eN对待均衡通道的基带复信号xA进行均衡处理,得到均衡后的复信号序列xAE;
3.1、将待均衡通道的基带复信号xA转换到频域;
XAN=FFT(xA,N)┄┄┄┄┄┄┄┄┄┄┄┄[7]
3.2、在频域用eN对待均衡通道的频域数据进行均衡处理并转换到时域;
XANE=IFFT(XAN*eN)┄┄┄┄┄┄┄┄┄┄┄┄[8]
3.3、对N点复信号序列xANE取前1*L个数据,得到均衡后的复信号序列xAE
xAE=xANE(1:L)┄┄┄┄┄┄┄┄┄┄┄┄[9]。
发明创造的优点
本发明能够减少通道均衡所需的存储空间,改善通道间对消性能,从而提升阵列处理性能,满足现代通信系统以及相控阵雷达系统等对多通道阵列信号处理的需要。
具体实施方式
下面对本发明做进一步详细说明。一种频率分段的通道均衡方法,基于已知的通道频域均衡滤波器频率特性e=[e[1],e[2],…,e[P]]T,其中e[p],p=1…P代表从-fs/2至fs/2内P个不同频点的均衡系数,fs为采样频率;已采集的参考通道与待均衡通道的基带复信号xM和xA,xM=[xM[1],xM[2],…,xM[L]]T为参考通道的复信号,xM[l],l=1…L是参考通道第l个复采样;xA=[xA[1],xA[2],…,xA[L]]T为待均衡通道的复信号,xM[l],l=1…L是参考通道第l个复采样;信号所占带宽BW,其中[·]T表示转置;其特征在于,通道均衡的步骤如下:
1、将频域均衡滤波器转换为N点滤波器:为便于L点时域信号xM[l]的频域滤波处理,选取N为大于等于L的最小二的幂,求滤波器的N点频率特性
其中FFT(·,N)表示N点快速傅里叶变换,IFFT(·,N)表示快速傅里叶逆变换;
2、求出分段均衡的各段滤波器系数,将实际使用的系数由N点减少为(F2-F1+1)点,以减少后续处理的运算量:
2.1、根据需要将N点分为K段,每段M个频点,即N=M×K,将变为M×K的矩阵
其中reshape(·,M,K)表示将向量变为M×K的矩阵,并将写成如下的块矩阵:
其中为M×(F1-1)的块矩阵,代表的1~F1-1列,为M×(F2-F1+1)的块矩阵,代表的F1~F2列,为M×(K-F2)的块矩阵,代表的K-F2+1~K列,F1、F2分别表示K个频段中的F1、F2个子频段,F1~F2个子频段覆盖信号带宽BW;1~F1-1及K-F2+1~K个子频段为带外部分;
2.2、对进行处理,得到频率分段均衡系数eK:
2.2.1、首先,对求平均,
其中,mean(·)为对列求平均,
2.2.2、随后,将的1~F1-1以及F2+1~K个元素置为1,得到各段滤波器系数为eK,写为下式:
其中,为频率分段处理后需要存储的频域均衡系数,该系数为(F2-F1+1)点。
2.3、根据分段均衡系数eK求得实际均衡时所采用N点频域均衡系数eN:
其中,[·]T为矩阵转置;
此时,虽然eN是N点的滤波器系数,实际上需要存储的滤波器系数只有(F2-F1+1)点,减少了处理的存储量;
3、用eN对待均衡通道的基带复信号xA进行均衡处理,得到均衡后的复信号序列xAE;
3.1、将待均衡通道的基带复信号xA转换到频域,以便进行频域均衡处理;
XAN=FFT(xA,N)┄┄┄┄┄┄┄┄┄┄┄┄[7]
3.2、在频域用eN对待均衡通道的频域数据进行均衡处理并转换到时域;
XANE=IFFT(XAN*eN)┄┄┄┄┄┄┄┄┄┄┄┄[8]
3.3、对N点复信号序列xANE取前1□L个数据,得到均衡后的复信号序列xAE
xAE=xANE(1:L)┄┄┄┄┄┄┄┄┄┄┄┄[9]
实施例
本发明的一个实施例中,通道频域均衡滤波器频率特性e=[e[1],e[2],…,e[P]]T,P=256;fs=4MHz;参考通道与待均衡通道的基带复信号xM=[xM[1],xM[2],…,xM[L]]T,xA=[xA[1],xA[2],…,xA[L]]TL=1800;信号所占带宽BW=3MHz。
均衡过程如下:
1、取N=2048,求通道频域均衡滤波器的N点频率特性
2、求出分段均衡的各段滤波器系数,将实际使用的系数由2048点减少为24点,以减少后续处理的运算量:
2.1、将2048点分为32段,每段64个频点,将变为64×32的矩阵
并将写成如下的块矩阵,其中24个子频段(5~28)覆盖信号带宽范围3MHz:
2.2、对进行处理,得到eK=[11×4b1×2411×4]1×32,其中b1×24即为系统所需存储的频域均衡系数,该系数为仅为24点;
2.3、求得实际均衡时所采用N点频域均衡系数此时,虽然eN是64点的滤波器系数,实际上需要存储的滤波器系数只有24点,减少了处理的存储量;
3、用eN对待均衡通道的基带复信号xA进行均衡处理,得到均衡后的复信号序列xAE;
3.1、将待均衡通道的基带复信号xA转换到频域,以便进行频域均衡处理;
XAN=FFT(xA,2048)
3.2、在频域用eN对待均衡通道的频域数据进行均衡处理并转换到时域;
XANE=IFFT(XAN*eN)
3.3、对N点复信号序列xANE取前1*L个数据,得到均衡后的复信号序列xAE=xANE(1:1800)
均衡结果分析:利用本方法,用均衡后的复信号序列与参考通道的复信号序列求对消比CR,比均衡前改善了6.7dB;所需存储的数组为b1×24,只需存储24个复数即可,避免了存储2048个复数,大大减小了存储空间,提高了资源利用效率,同时使得只具备小存储空间的信号处理系统也能够实现良好的均衡效果。
Claims (1)
1.一种频率分段的通道均衡方法,基于已知的通道频域均衡滤波器频率特性e=[e[1],e[2],…,e[P]]T,其中e[p]代表从-fs/2至fs/2内P个不同频点的均衡系数,p=1…P,fs为采样频率;已采集的参考通道与待均衡通道的基带复信号xM和xA,xM=[xM[1],xM[2],…,xM[L]]T为参考通道的复信号,xM[l]是参考通道第l个复采样,l=1…L;xA=[xA[1],xA[2],…,xA[L]]T为待均衡通道的复信号,xM[l]是参考通道第l个复采样,l=1…L;信号所占带宽为BW,其中[·]T表示转置;其特征在于,通道均衡的步骤如下:
第1步、将频域均衡滤波器转换为N点滤波器:选取N为大于等于L的最小二的幂,求频域均衡滤波器的N点频率特性
其中FFT(·,N)表示N点快速傅里叶变换,IFFT(·,N)表示N点快速傅里叶逆变换;
第2步、求出分段均衡的各段滤波器系数:
2.1、根据需要将N点分为K段,每段M个频点,即N=M×K,将变为M×K的矩阵
其中reshape(·,M,K)表示将向量变为M×K的矩阵,并将写成如下的块矩阵:
其中为M×(F1-1)的块矩阵,代表的1~F1-1列,为M×(F2-F1+1)的块矩阵,代表的F1~F2列,为M×(K-F2)的块矩阵,代表的K-F2+1~K列,F1、F2分别表示K个频段中的F1、F2个子频段,F1~F2个子频段覆盖信号带宽BW;
2.2、对进行处理,得到频率分段均衡系数eK:
2.2.1、首先,对求平均,
其中,mean(·)为对列求平均,
2.2.2、随后,将的1~F1-1以及F2+1~K个元素置为1,得到各段滤波器系数为eK,写为下式:
其中,为频率分段处理后需要存储的频域均衡系数;
2.3、根据分段均衡系数eK求得实际均衡时所采用N点频域均衡系数eN:
其中,[·]T为矩阵转置;
第3步、用eN对待均衡通道的基带复信号xA进行均衡处理,得到均衡后的复信号序列xAE;
3.1、将待均衡通道的基带复信号xA转换到频域;
XAN=FFT(xA,N)┄┄┄┄┄┄┄┄┄┄┄┄[7]
3.2、在频域用eN对待均衡通道的频域数据进行均衡处理并转换到时域;
XANE=IFFT(XAN*eN)┄┄┄┄┄┄┄┄┄┄┄┄[8]
3.3、对N点复信号序列xANE取前1*L个数据,得到均衡后的复信号序列xAE
xAE=xANE(1:L)┄┄┄┄┄┄┄┄┄┄┄┄[9]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410670589.6A CN105610748B (zh) | 2014-11-20 | 2014-11-20 | 一种频率分段的通道均衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410670589.6A CN105610748B (zh) | 2014-11-20 | 2014-11-20 | 一种频率分段的通道均衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105610748A true CN105610748A (zh) | 2016-05-25 |
CN105610748B CN105610748B (zh) | 2018-11-16 |
Family
ID=55990290
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410670589.6A Active CN105610748B (zh) | 2014-11-20 | 2014-11-20 | 一种频率分段的通道均衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105610748B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107666453A (zh) * | 2016-07-28 | 2018-02-06 | 上海诺基亚贝尔股份有限公司 | 发射器和相应的方法 |
CN108132460A (zh) * | 2017-12-07 | 2018-06-08 | 西安电子工程研究所 | 基于频域通道均衡的脉冲压缩补偿算法 |
CN109387818A (zh) * | 2018-10-10 | 2019-02-26 | 电子科技大学 | 一种基于回波数据的双基地sar通道均衡方法 |
CN112214895A (zh) * | 2020-10-14 | 2021-01-12 | 中国船舶重工集团公司第七二四研究所 | 一种变换点数优化设计的傅里叶多波束形成方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050078772A1 (en) * | 2003-10-10 | 2005-04-14 | Hitachi Global Technologies Netherlands B.V. | Apparatus using a lengthened equalization target filter with a matched filter metric in a viterbi detector |
CN102404672A (zh) * | 2011-10-27 | 2012-04-04 | 苏州上声电子有限公司 | 数字化扬声器阵列系统的通道均衡与波束控制方法和装置 |
CN102724151A (zh) * | 2012-06-15 | 2012-10-10 | 武汉邮电科学研究院 | 通信系统、通信方法以及多通道自适应均衡器和ofdm信号解复用的方法 |
CN102883243A (zh) * | 2012-10-15 | 2013-01-16 | 苏州上声电子有限公司 | 在线迭代的声重放系统频响均衡方法和装置 |
CN103401825A (zh) * | 2013-08-14 | 2013-11-20 | 西安电子科技大学 | 基于块状导频的低复杂度单载波频域均衡方法 |
-
2014
- 2014-11-20 CN CN201410670589.6A patent/CN105610748B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050078772A1 (en) * | 2003-10-10 | 2005-04-14 | Hitachi Global Technologies Netherlands B.V. | Apparatus using a lengthened equalization target filter with a matched filter metric in a viterbi detector |
CN102404672A (zh) * | 2011-10-27 | 2012-04-04 | 苏州上声电子有限公司 | 数字化扬声器阵列系统的通道均衡与波束控制方法和装置 |
CN102724151A (zh) * | 2012-06-15 | 2012-10-10 | 武汉邮电科学研究院 | 通信系统、通信方法以及多通道自适应均衡器和ofdm信号解复用的方法 |
CN102883243A (zh) * | 2012-10-15 | 2013-01-16 | 苏州上声电子有限公司 | 在线迭代的声重放系统频响均衡方法和装置 |
CN103401825A (zh) * | 2013-08-14 | 2013-11-20 | 西安电子科技大学 | 基于块状导频的低复杂度单载波频域均衡方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107666453A (zh) * | 2016-07-28 | 2018-02-06 | 上海诺基亚贝尔股份有限公司 | 发射器和相应的方法 |
US20190173713A1 (en) * | 2016-07-28 | 2019-06-06 | Alcatel Lucent | Transmitter and corresponding method |
US10911287B2 (en) | 2016-07-28 | 2021-02-02 | Alcatel Lucent | Transmitter and corresponding method |
CN107666453B (zh) * | 2016-07-28 | 2021-05-18 | 上海诺基亚贝尔股份有限公司 | 发射器和相应的方法 |
CN108132460A (zh) * | 2017-12-07 | 2018-06-08 | 西安电子工程研究所 | 基于频域通道均衡的脉冲压缩补偿算法 |
CN109387818A (zh) * | 2018-10-10 | 2019-02-26 | 电子科技大学 | 一种基于回波数据的双基地sar通道均衡方法 |
CN112214895A (zh) * | 2020-10-14 | 2021-01-12 | 中国船舶重工集团公司第七二四研究所 | 一种变换点数优化设计的傅里叶多波束形成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105610748B (zh) | 2018-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108964731B (zh) | 基于快速卷积的无循环前缀滤波混合载波连续流传输方法 | |
CN105610748A (zh) | 一种频率分段的通道均衡方法 | |
CN100553249C (zh) | 用于wcdma直放站系统的在频域实现的干扰抵消方法 | |
CN110515038B (zh) | 一种基于无人机-阵列的自适应无源定位装置及实现方法 | |
CN106301498A (zh) | 子带处理方法及频空级联的宽带自适应波束获取方法 | |
CN101547173A (zh) | 基于平衡正交多小波变换的盲均衡方法 | |
CN101833949A (zh) | 一种用于消减噪声的有源噪声控制方法 | |
CN101685638A (zh) | 一种语音信号增强方法及装置 | |
CN104811867B (zh) | 基于阵列虚拟扩展的麦克风阵列空域滤波方法 | |
CN107037409B (zh) | 基于压缩感知的mimo雷达波形分离方法 | |
CN108365874B (zh) | 基于fdd大规模mimo贝叶斯压缩感知信道估计方法 | |
CN110320535A (zh) | 一种基于小波包变换与空域处理的卫星导航接收机多干扰抑制方法 | |
CN105679304A (zh) | 用于宽带主动噪声控制系统的可变带宽无延迟子带算法 | |
CN109188366B (zh) | 基于子带最大信噪比准则的宽带发射自适应波束形成方法 | |
CN104869086A (zh) | 基于二维压缩感知的mimo-ofdm通信系统下行信道估计方法、装置 | |
CN108845294B (zh) | 基于子带线性多约束最小方差准则的宽带发射自适应波束形成方法 | |
WO2020135875A1 (zh) | 一种基于维纳自适应的信道估计的方法及系统 | |
Lobov et al. | Investigation of the properties of an filterbank multicarrier modem in a broadband HF channel | |
CN105306117B (zh) | 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 | |
CN110166098B (zh) | 一种宽带唯相位发射自适应波束形成方法 | |
Abd El-Galil et al. | Efficient underwater acoustic communication with peak-to-average power ratio reduction and channel equalization | |
CN106646531B (zh) | 一种多星约束稳健空频抗干扰处理方法和装置 | |
CN108828536B (zh) | 基于二阶锥规划的宽带发射数字波束形成干扰设计方法 | |
CN112420068B (zh) | 一种基于Mel频率尺度分频的快速自适应波束形成方法 | |
CN107167809A (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 |