通信信号变采样率频域实现方法.docVIP

通信信号变采样率频域实现方法 摘 要：提出一种通信信号变采样率的频域实现方法。首先，介绍了解析信号的概念及由实信号转换为解析信号的方法；然后，以2 倍升采样为例，通过理论分析证明了变采样率频域实现的可行性；接下来，针对整数倍升采样、整数倍降采样和分数倍变采样三种情况，分别介绍了相应的频域处理方法，并进行了仿真验证；最后，简单分析了算法的优点及适用领域。 关键词：变采样率；解析信号；频域 0 引言 在通信信号的处理过程中，经常要根据实际需要对采样率进行调整，即要求信号分析系统能够工作在多采样率状态。一般来说，采样率的调整都是在时域进行的，计算量比较大。虽然可以采用多相滤波器组来实现，以降低运算量提高运算速度，但对于大整数倍和分数倍变采样率，多相滤波器组的结构非常复杂，而且其结构完全是通过数学推导得到的，不能直观地反应变采样率的物理过程。与时域处理相比，通信信号的频域处理方法，如上下变频、频域滤波等，通常能够简化处理过程，降低算法的复杂度，并且直接与物理过程相对应，易于理解。 1 解析信号 对于一个通信信号可以用实函数形式或复函数形式表示，在进行信号处理之前往往要将实信号复化，变为它的解析形式。最简单的方法是直接用实信号s(t)作复信号z(t)的实部，并添加一个信号x(t) 作其虚部： z(t) = s(t) + jx(t) 。 （1） 通常的做法是保留其正频部分（并将幅度加倍，以使原信号的总能量保持不变），并剔除掉负频部分。由于实信号的频谱为共轭对称，剔除负频部分不会造成任何信息损失，也不会带来虚假信息。这样，其频谱就不存在共轭对称性，所对应的时域信号应为复信号。因此，通信中接收机接收的一般是正交的两路信号，即信号的实部和虚部。如果只接收了一路信号，那就要对它进行复化，这样才能得到通信信号的幅度和相位信息。剔除实信号的负频部分，在频域相当于把实函数的频谱通过一个阶跃式传输函数，其过程在时域可以表示为： （2） 式中，H[s(t)]称为s(t)的Hilbert变换，即： z(t) = s(t) + jH[s(t)]。 （3） 对于采样后的离散化数据，式（3）转换为： z(n) = s(n) + jH[s(n)]。 （4） 由于解析信号非常便于处理，因此在后面的讨论中，都首先将通信信号复化为解析信号然后再进行分析。 2 整数倍升采样 整数倍升采样的一般原理如图1 所示。其中L 为升采样的倍数，h(n) 为低通滤波器，其过程是首先在原序列相邻数据之间插入L-1 个0，然后在通过低通滤波滤除多余的频谱分量，就得到L 倍升采样后的序列。下面以2 倍内插为例讨论如何在频域实现这一过程。 图1 L 倍升采样原理 设原序列为{s(n)}，n=0,1,…,N-1，其复数表示为{z(n)}，即s(n)=Re{z(n)}，插入0 后的序列为{y(k)}，k=0,1,…,2N-1，不失一般性： （5） 设Z(m)为z(n)的DFT 变换： ，m=0,1,…,N-1， （6） Y(l)为y(k)的DFT 变换，即： ,l==0,1,…,2N-1， （7） 考察Y(l)： , (8) 当l≤N-1 时： Y(l) = Z(l)， （9） 当lN-1 时，令l=N+v，则： 。（10） 一般在处理之前都要去掉通信信号中的均值，因此=0，所以： ， （11） 由此可见，插入0 值后的频谱是原频谱延拓一个周期，而插零后的低通滤波，其作用就是滤除延拓的频谱。这一过程可以推广到L 倍升采样。L 倍升采样中，插入L-1 个0 后信号的频谱是原频谱延拓L-1 个周期，而滤波的过程就是将这L-1 个延拓频谱滤除。这一过程等价于在原有频谱Z(m)补(L-1)N 个0。综上所述，可以得到L 倍升采样的频域处理过程如下， ① 将输入实信号s(n)进行复化，得到解析信号z(n)； ② 求z(n)的频谱Z(m)； ③ 在Z(m)进行扩展，补(L-1)N 个0 得到频谱Y(l)； ④ 对Y(l)做LN 点的逆DFT 运算并取实部，即得到了L倍升采样后的时域信号y(k)。 图2是部分仿真结果。输入信号是频率为1kHz 的余弦信号，采样率为20 ks/s。为了便于比较，以样点数作为时域波形的横轴。 （a）原信号波形 (b) 2 倍升采样后信号波形 （c）3 倍升采样后信号波形 （d）4 倍升采样后信号波形 图2 整数倍升采样仿真结果 3 整数倍降采样 整数倍降采样的一般原理如图3 所示[3]。首先通过低通滤波器h(n)滤除多余频谱分量，避免降采样后发生混叠，然后进行