Python科学计算——检包络与去包络

Python025

Python科学计算——检包络与去包络,第1张

检波 (detection):广义的检波通常称为 解调 ,是调制的逆过程,即从已调波提取调制信号的过程。狭义的检波是指从调幅波的包络提取调制信号的过程,这种检波方法也被称为 包络检波 希尔伯特变换 可以用作包络检波。

Hilbert 变换能在振幅保持不变的情况下将输入信号的相角偏移 90 度,简单地说就是能够将正弦波形转换为余弦波形:

相角偏移90度相当于复数平面上的点与虚数单位 1j 相乘,因此 Hilbert 变换的频率响应可以用如下公式表示:

Hilbert 变换可以用作 包络检波 。具体算法如下:

用频率扫描波可以测量滤波器的频率响应,也可以用它检测 Hilbert 变换用于包络检波的性能:

从上图可以看出,在高频和低频处包络计算出现较大的误差,而中频部分能很好地计算出包络的形状。

在 Hilbert 变换检测出包络的基础上,可以利用简单的去包络算法将包络从原始信号中去除而恢复载波信息,去包络算法用公式表示如下:

从上面的图可以看出,包络已经被很好的去掉,但是从时域图像,并不能完全确定包络被很好的去掉,我们需要从去包络前后信号的频率响应一探究竟:

从去包络前后信号的频率响应图可以看出,包络已经被很好的抑制,只剩下了单频载波信号。

从上述去包络前信号的频率响应图可以看出,其上下边带的幅度是相等的,当上下边带信号幅度不对等时,去包络算法效果会如何呢?

当边带信号幅度不对等时,包络的抑制效果就会变差,而且还会引入新的频率成份,这将会在一定程度上恶化信号。

FFT (Fast Fourier Transform, 快速傅里叶变换) 是离散傅里叶变换的快速算法,也是数字信号处理技术中经常会提到的一个概念。用快速傅里叶变换能将时域的数字信号转换为频域信号,转换为频域信号后我们可以很方便地分析出信号的频率成分。

当我们把双频信号FFT示例中的 fft_size 的值改为 2**12 时,这时,基频为 16Hz,不能被 1kHz整除,所以 1kHz 处发生了频谱泄露,而它能被 4kHz 整除,所以 4kHz 可以很好地被采样。

由于波形的前后不是连续的,出现波形跳变,而跳变处有着非常广泛的频谱,因此FFT的结果中出现了频谱泄漏。

为了减小FFT所截取的数据段前后的跳变,可以对数据先乘以一个窗函数,使得其前后数据能平滑过渡。常用的hanning窗函数的定义如下:

50Hz 正弦波与hann窗函数乘积之后的重复波形如下:

我们对频谱泄漏示例中的1kHz 和 4kHz 信号进行了 hann 窗函数处理,可以看出能量更加集中在 1kHz 和 4kHz,在一定程度上抑制了频谱泄漏。

以 1kHz 三角波为例,我们知道三角波信号中含有丰富的频率信息,它的傅里叶级数展开为:

当数字信号的频率随时间变化时,我们称之为扫频信号。以频率随时间线性变化的扫频信号为例,其数学形式如下:

其频率随时间线性变化,当我们在 [0,1] 的时间窗口对其进行采样时,其频率范围为 0~5kHz。当时间是连续时,扫频信号的频率也是连续的。但是在实际的处理中,是离散的点采样,因此时间是不连续的,这就使扫频信号的快速傅里叶变换问题退化为多点频信号快速傅里叶变换问题。其快速傅里叶变换得到的频谱图如下所示:

以 50Hz 正弦信号相位调制到 1kHz 的信号为例,其信号形式如下:

它的时域波形,频率响应和相位响应如下图所示:

以扫频信号为例,当我们要探究FFT中的能量守恒时,我们要回归到信号最初的形式: