一种频率分段的通道均衡方法
【技术领域】
[0001] 本发明属于阵列信号技术，涉及对一定带宽内通道不匹配的均衡方法。
【背景技术】
[0002] 在现代雷达和现代通信中，双通道或多通道阵列信号处理的应用越来越广泛。例 如在雷达信号处理中，传统的雷达旁瓣对消技术（SLC)，以及自适应波束形成技术（ADBF)， 就是利用两个或多个通道来对消来自与目标不同角度的干扰；空时自适应处理技术 (STAP)则可利用多个通道来对消空时耦合的各类杂波。然而通道之间频率特性的差异将造 成对消的剩余，从而造成信干噪比（SINR)或信杂噪比（SCNR)的下降。通道均衡能够校正 通道间频率特性的差异，改善通道对消性能，从而提高各种先进处理方法的SINR或SCNR。
[0003] 关于通道均衡，不少学者对之进行了研究，并发表了相关论文，如王峰等人的《基 于傅里叶变换的雷达通道均衡算法性能分析及改进》，李建平的《宽带数字阵列雷达通道均 衡算法研究》，冉涌的《宽带数字阵列雷达通道校正技术》，以及赵立艳的《数字阵列雷达中 通道均衡技术的研究》等。但以上论文提及的方法都是基于如何构建时域滤波器系数的方 法，这种时域处理方法需通过逻辑时序单元实现，很难叠加到其它处理过程（如脉压过程） 同时实现。另一方面，即使已知均衡滤波器的频率响应，直接在脉压处理时应用需要极大的 存储空间，限制了频域均衡的工程应用。
【发明内容】
[0004] 发明创造的目的
[0005] 本发明的目的是：提出一种频率分段的通道均衡方法，以满足多通道阵列信号处 理的需要。
[0006] 技术方案
[0007] 本发明的技术方案是：一种频率分段的通道均衡方法，基于已知的通道频域均衡 滤波器频率特性e = [e [1]，e [2]，…，e [P] ]τ，其中e [P]代表从-fs/2至fs/2内P个不同 频点的均衡系数，P = 1…P，fs为采样频率；已采集的参考通道与待均衡通道的基带复信号 xdP xA，xM= [xM[l], xM[2],…，xM[L]]%参考通道的复信号，xM[l]是参考通道第1个复采 样，1 = 1···ΜχΑ= [χΑ[1]，χΑ[2]，···，χΑα]]τ为待均衡通道的复信号，x M[l]是参考通道第 1个复采样，1 = 1···? ;信号所占带宽为BW，其中[· ]τ表示转置；其特征在于，通道均衡的 步骤如下：
[0008] 第1步、将频域均衡滤波器转换为N点滤波器：选取N为大于等于L的最小二的 幂，求频域均衡滤波器的N点频率特性<，
[0010] 其中FFT (·，Ν)表示N点快速傅里叶变换，IFFT (·，Ν)表示N点快速傅里叶逆变 换；
[0011] 第2步、求出分段均衡的各段滤波器系数：
[0012] 2. 1、根据需要将N点分为K段，每段M个频点，即N = MXK，将4变为MXK的矩 阵蜮
[0014] 其中reshape ( ·，Μ，Κ)表示将向量变为MXK的矩阵，并将4写成如下的块矩阵：
[0016] 其中丨为MX (Ffl)的块矩阵，代表Ε|的1~F1-I列，Β,°/χ巧-以）为 MX (Wl)的块矩阵，代表祀的匕~F2列，为MX (K-F2)的块矩阵，代表Eg的 K-F2+l~K列，匕分别表示K个频段中的F ^ F2个子频段，F广F 2个子频段覆盖信号带 宽BW ;
[0017] 2. 2、对进行处理，得到频率分段均衡系数eK:
[0018] 2. 2. 1、首先，对求平均，
[0020] 其中，mean( ·)为对列求平均，
[0021] 2. 2. 2、随后，将ef的1~F1-I以及F2+l~K个元素置为1，得到各段滤波器系数 为e K，写为下式：
[0023] 其中，为频率分段处理后需要存储的频域均衡系数。
[0024] 2. 3、根据分段均衡系数％求得实际均衡时所采用N点频域均衡系数e N:
[0026] 其中，[· ]τ为矩阵转置；
[0027] 第3步、用eN对待均衡通道的基带复信号χΑ进行均衡处理，得到均衡后的复信号 序列X?;
[0028] 3. 1、将待均衡通道的基带复信号xA转换到频域；
[0029] Xan= FFT (x a, N)-------------------------[7]
[0030] 3. 2、在频域用eN对待均衡通道的频域数据进行均衡处理并转换到时域；
[0031] XANE= IFFT(XAN*eN)-------------------------[8]
[0032] 3. 3、对N点复信号序列xANE取前1*L个数据，得到均衡后的复信号序列x AE
[0033] xAE= xANE(l :L)-------------------------[9] 〇
[0034] 发明创造的优点
[0035] 本发明能够减少通道均衡所需的存储空间，改善通道间对消性能，从而提升阵列 处理性能，满足现代通信系统以及相控阵雷达系统等对多通道阵列信号处理的需要。
【具体实施方式】
[0036] 下面对本发明做进一步详细说明。一种频率分段的通道均衡方法，基于已知的通 道频域均衡滤波器频率特性e = [e [1]，e [2]，…，e [P] ]τ，其中e [p]，p = 1···Ρ代表从_fs/2 至fs/2内P个不同频点的均衡系数，fs为采样频率；已采集的参考通道与待均衡通道的基 带复?目号Xm和X A，xM= [X M[l]，xm[2]，…，Xm[L] ]%参考通道的复?目号，X M[l]，1 = 1…L是 参考通道第1个复采样；&=[^[1]，&[2]，一，&1]^为待均衡通道的复信号，^[1]，1 = 1…L是参考通道第1个复采样；信号所占带宽BW，其中[· ]τ表示转置；其特征在于，通道 均衡的步骤如下：
[0037] 1、将频域均衡滤波器转换为N点滤波器：为便于L点时域信号χΜ[1]的频域滤波 处理，选取N为大于等于L的最小二的幂，求滤波器的N点频率特性eg，
[0039] 其中FFT (·，N)表示N点快速傅里叶变换，IFFT (·，N)表示快速傅里叶逆变换；
[0040] 2、求出分段均衡的各段滤波器系数，将实际使用的系数由N点减少为（F2-Ffl) 点，以减少后续处理的运算量：
[0041] 2. 1、根据需要将N点分为K段，每段M个频点，即N = MXK，将 < 变为MXK的矩 阵
[0043] 其中reshape(*，M，K)表示将向量变为MXK的矩阵，并将写成如下的块矩阵：
[0045] 其中誇】为MX (F「l)的块矩阵，代表Ef.的1~F1-I列，为 MX (F2-FM)的块矩阵，代表E总的FfF2列，咐-F2)为MX (K-F2)的块矩阵，代表Eg的 K-F2+l~K列，匕分别表示K个频段中的F ^ F2个子频段，F广F 2个子频段覆盖信号带 宽BW ;1~F1-I及K-F2+l~K个子频段为带外部分；
[0046] 2. 2、对Ef进行处理，得到频率分段均衡系数eK:
[0047] 2. 2. 1、首先，对的求平均，
[0049] 其中，mean( ·)为对列求平均，
[0050] 2. 2. 2、随后，将ef的1~F1-I以及F2