FFT能量归一化

简介: FFT能量归一化

前言

时域信号转换为频域信号时,涉及到幅度和能量的变化,目前大部分开源库在正变换和反变换时会忽略常数,因此当我们想将频域和时域信号归一化到统一尺度时(方便设置阈值),需要做归一化操作。


一、能量(功率)归一化有什么用?

答:添加功率归一化因子,目的在于使得不同调制方式(或者说对于所有映射方式)都能够取得相同的平均功率。

实际上,归一化是为了方便系统性能的比较,所以就要分清比较的模块是什么。比如,信道编码的增益问题,无论有无信道编码,比特能量是一样的,所以比较要以 Eb/No 为基准,而不是以进入信道前的符号能量 Es/No 为基准。再比如,在比较空时码系统和单天线系统中,还是以进入时空码编码前信号能量为基准,那么发送时的总能量一致,即时空码系统中各天线发射功率总和应和单天线系统发射功率相同。一般而言,归一化都在发射端处理。

二、归一化方式

1、单个频点幅度

X k = ∑ n = 0 N − 1 x n e − i 2 π k n / N X_k=\sum_{n=0}^{N-1}x_ne^{-i2\pi kn/N}Xk=n=0N1xnei2πkn/N

x n = 1 N ∑ k = 0 N − 1 X k e i 2 π n / N x_n=\frac{1}{N}\sum_{k=0}^{N-1}X_ke^{i2\pi n/N}xn=N1k=0N1Xkei2πn/N

根据傅里叶变换公式,使用 1 / N 1/N1/N 进行归一化,对于实数来说,有效频点为 N / 2 N/2N/2,因此使用 2 / N 2/N2/N

2、整帧能量

∑ n = 0 N − 1 ∣ x [ n ] ∣ 2 = 1 N ∑ k = 0 N − 1 ∣ X [ k ] ∣ 2 \sum_{n=0}^{N-1}|x[n]|^2=\frac{1}{N}\sum_{k=0}^{N-1}|X[k]|^2n=0N1x[n]2=N1k=0N1X[k]2

根据 DFT 变换的帕萨瓦尔定理,使用 1 / N 1/\sqrt N1/N

三、仿真测试

在进行 FFT、IFFT 时分别需除 sqrt(N),乘 sqrt(N),这样做的目的是使能量归一化,即使得时域和频域数据的能量一样。

MATLAB 测试:

Data_inf = randi([0,1],100,1)+1i*randi([0,1],100,1);
ifftData = ifft(Data_inf)*sqrt(100);
fftData = fft(Data_inf)./sqrt(100);
powerData = sum(abs(Data_inf).^2);
powerifft = sum(abs(ifftData).^2);
powerfft = sum(abs(fftData).^2);
fprintf('powerData = %d\n', powerData);
fprintf('powerifft = %d\n', powerifft);
fprintf('powerfft = %d\n', powerfft);

仿真结果:

powerData = 93
powerifft = 9.300000e+01
powerfft = 9.300000e+01

可知,数据能量都是 93,验证分析正确。

目录
相关文章
|
6月前
QAM 归一化因子
QAM 归一化因子
133 0
|
6月前
|
计算机视觉
傅里叶
傅里叶 “【5月更文挑战第23天】”
80 1
|
6月前
|
计算机视觉 Python
傅里叶滤波
傅里叶滤波
27 3
|
6月前
|
数据安全/隐私保护
matlab 曲线光滑,去毛刺,去离群值,数据滤波,高通滤波,低通滤波,带通滤波,带阻滤波
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
|
机器学习/深度学习
傅立叶变换之(一)——欧拉公式
傅立叶变换之(一)——欧拉公式
|
机器学习/深度学习 传感器 算法
【图像重建】基于离散余弦变换DCT图像重建含MSE、PSNR、压缩比附matlab代码
【图像重建】基于离散余弦变换DCT图像重建含MSE、PSNR、压缩比附matlab代码
|
机器学习/深度学习 传感器 算法
【滤波】基于平方根无迹卡尔曼滤波SR-UKF实现信号去噪附matlab代码
【滤波】基于平方根无迹卡尔曼滤波SR-UKF实现信号去噪附matlab代码
|
机器学习/深度学习 算法 计算机视觉
【MATLAB】离散余弦变换滤波算法(DCT)
之前介绍的所有滤波算法都是空间域滤波算法(即2D滤波算法)。离散余弦变换滤波算法(DCT)属于频率域滤波算法(即3D滤波算法)。
360 0
【MATLAB】离散余弦变换滤波算法(DCT)
|
计算机视觉
数字图像处理——图像变换(二维离散傅里叶正逆变换fft2/fft2、离散余弦正逆变换dct2/idct2、频谱正逆平移fftshift/ifftshift、幅度谱与相位谱)| 例题与分析
数字图像处理——图像变换(二维离散傅里叶正逆变换fft2/fft2、离散余弦正逆变换dct2/idct2、频谱正逆平移fftshift/ifftshift、幅度谱与相位谱)| 例题与分析
324 0
数字图像处理——图像变换(二维离散傅里叶正逆变换fft2/fft2、离散余弦正逆变换dct2/idct2、频谱正逆平移fftshift/ifftshift、幅度谱与相位谱)| 例题与分析
|
计算机视觉
数字图像处理实验(二)|图像变换{离散傅里叶变换fft2,离散余弦变换dct2、频谱平移fftshift}(附实验代码和截图)
数字图像处理实验(二)|图像变换{离散傅里叶变换fft2,离散余弦变换dct2、频谱平移fftshift}(附实验代码和截图)
410 0
数字图像处理实验(二)|图像变换{离散傅里叶变换fft2,离散余弦变换dct2、频谱平移fftshift}(附实验代码和截图)