离散余弦变换（Discrete Cosine Transform, DCT），信号去噪

介绍

离散余弦变换（Discrete Cosine Transform, DCT）是一种常用的信号处理工具，特别是在数据压缩、图像处理和模式识别等领域中。DCT的基本思想是将信号从空间域或时间域转换到频率域，以揭示信号中的频率成分。与离散傅里叶变换（DFT）相比，DCT具有一些独特的性质，使其在特定应用场景中更加有效。

DCT的背景与动机

DCT的开发最初是为了克服傅里叶变换在处理边界条件和处理实值信号时的一些局限性。傅里叶变换和DFT将信号表示为一系列复数指数函数的叠加，其中包括正弦和余弦成分。这种表示虽然功能强大，但也引入了复数运算和相位信息，这在某些应用中并不必要。

DCT将信号表示为一组余弦函数的叠加，所有这些余弦函数都具有不同的频率，但没有复数部分，因此DCT只处理实值数据。对于许多信号，特别是图像和音频信号，DCT能够集中能量在低频部分，从而在压缩和特征提取过程中表现出色。

DCT的定义

离散余弦变换有多种形式，最常用的为DCT-II。假设我们有一个长度为$N$的序列$x[n]$，其中$n=0,1,2,...,N-1$。DCT-II 的定义如下：
$$X[k]=\sum _{n=0}^{N-1}x[n]\cdot cos[\frac{\pi }{N}(n+\frac{1}{2})k]$$
其中$k=0,1,2,...,N-1$,$X[k]$是信号在频率 𝑘处的DCT系数。
DCT的逆变换（IDCT）用于将频域信号转换回时间或空间域，其定义为：
$$x[n]=\sum _{n=0}^{N-1}x[n]\cdot cos[\frac{\pi }{N}(n+\frac{1}{2})k]$$
DCT的其他形式包括DCT-I、DCT-III和DCT-IV，它们在边界条件和具体应用上有所不同，但DCT-II是最常用的版本，特别是在图像压缩中。

DCT的性质

DCT具有一些重要的性质，使其在信号处理中的应用非常广泛：

能量压缩（Energy Compaction）：
DCT能够将信号的大部分能量集中在较少的低频系数上，这对于信号压缩特别有利。在许多应用中，DCT系数的前几项就能很好地近似原始信号，而高频系数则可以被忽略或进一步压缩。

实值变换（Real-Valued Transform）：
DCT只处理实数信号，没有复数部分，这简化了计算，并使得DCT在处理实值数据（如图像和音频）时非常有效。

对称性：
DCT的对称性使其在处理边界条件时表现得更好，特别是在处理周期信号或以某种方式对称的信号时。

可逆性：
DCT是可逆的，意味着我们可以通过逆DCT（IDCT）从频域信号恢复原始信号，前提是没有对DCT系数进行过度的削减。

DCT的应用

DCT在各种应用中发挥着关键作用，以下是一些主要的应用领域：

图像压缩：
在JPEG图像压缩标准中，DCT是核心技术。JPEG压缩将图像分成8x8的块，对每个块应用DCT，然后量化和压缩DCT系数，去除人眼不易察觉的高频成分，从而达到压缩的目的。

音频压缩：
在音频压缩（如MP3、AAC）中，DCT用于将时间域信号转换到频域，以便更有效地压缩。音频信号通常具有较强的低频成分，DCT能够很好地将这些能量集中在少数几个系数中。

特征提取：
在模式识别、语音识别等领域，DCT用于提取信号的频率特征。例如，在语音识别中，DCT是计算MFCC（Mel Frequency Cepstral Coefficients，梅尔频率倒谱系数）的关键步骤。

数据压缩与去噪：
DCT也用于数据压缩和去噪。通过只保留DCT系数中的显著成分，可以去除噪声或压缩数据，同时尽量保留原始信号的重要特征。

DCT与其他变换的比较

DCT与DFT（离散傅里叶变换）和DWT（离散小波变换）是信号处理中常用的三种变换，各有优缺点：

DCT vs. DFT：
DCT与DFT相比，DCT只涉及实数计算，没有相位信息，因此计算上更简单且能量集中效果更好。但DFT更适合处理周期性信号和需要复杂频率分析的场景。

DCT vs. DWT：
DCT是线性变换，处理时间复杂度低；而DWT具有多分辨率分析能力，适合处理具有突变或非平稳特性的信号。DWT在图像处理中的应用也非常广泛，但DCT在压缩效率和计算复杂度上具有优势。

DCT的局限性

尽管DCT在许多方面表现出色，但它也有一些局限性：

边缘效应：
在处理边界数据时，DCT可能引入伪影，特别是在图像处理时，块之间的边界可能产生可见的块效应（Block Artifacts）。

全局性质：
DCT是一种全局变换，无法像小波变换那样提供多分辨率分析。因此，它对局部信号特征的捕捉能力有限。

应用范围：
DCT主要适用于实值信号和对称性较强的信号，对于复数信号或不对称信号，其应用效果可能不如其他变换。

本文代码

我们将展示了如何使用DCT去除信号中的噪声。代码包括以下几个步骤：生成含噪信号、应用DCT、使用阈值方法去除噪声、然后应用逆DCT恢复信号，并进行可视化对比

核心代码

% MATLAB Code for Denoising Using Discrete Cosine Transform (DCT)

% Step 1: Generate a Noisy Signal
Fs = 1000;                 % Sampling frequency
t = 0:1/Fs:1-1/Fs;         % Time vector (1 second duration)
freq1 = 50;                % Frequency of the first sinusoid
freq2 = 120;               % Frequency of the second sinusoid
signal = sin(2*pi*freq1*t) + sin(2*pi*freq2*t); % Clean signal

% Add Gaussian noise
SNR = 10;                  % Signal-to-noise ratio (in dB)
noisySignal = awgn(signal, SNR, 'measured');  % Noisy signal

% Step 2: Apply DCT to the Noisy Signal
dctCoeffs = dct(noisySignal);

% Step 3: Apply Thresholding to DCT Coefficients
% Thresholding is used to remove the noise components in the DCT domain
threshold = 0.2 * max(abs(dctCoeffs));  % Define a threshold based on max coefficient
dctCoeffs(abs(dctCoeffs) < threshold) = 0;  % Zero out coefficients below the threshold

% Step 4: Apply Inverse DCT to Recover the Denoised Signal

% Step 5: Visualize the Results
figure;

% Plot the original clean signal
subplot(4,1,1);
plot(t, signal);
title('Original Clean Signal');
xlabel('Time (s)');
ylabel('Amplitude');

% Plot the noisy signal
subplot(4,1,2);
plot(t, noisySignal);
title('Noisy Signal');
xlabel('Time (s)');
ylabel('Amplitude');

% Plot the DCT coefficients before and after thresholding
subplot(4,1,3);
plot(dctCoeffs);
title('DCT Coefficients after Thresholding');
xlabel('Coefficient Index');
ylabel('Magnitude');

% Plot the denoised signal
subplot(4,1,4);
plot(t, denoisedSignal);
title('Denoised Signal Using DCT');
xlabel('Time (s)');
ylabel('Amplitude');

% Step 6: Evaluate the Denoising Performance
% Compute the Mean Squared Error (MSE) between the original and denoised signal

% Compute the Signal-to-Noise Ratio (SNR) improvement
snrImprovement = 20*log10(norm(signal)/norm(signal-denoisedSignal));

% Display MSE and SNR improvement
fprintf('Mean Squared Error (MSE) between original and denoised signal: %.4f\n', mse);
fprintf('SNR Improvement: %.2f dB\n', snrImprovement);

代码说明

生成含噪信号：
生成了一个包含两个不同频率（50Hz 和 120Hz）正弦波的纯净信号。然后，通过添加高斯白噪声（AWGN）来创建一个信噪比为10dB的含噪信号。
应用DCT：
对含噪信号进行DCT变换，将其转换到频域。DCT将信号的主要能量集中在低频部分，因此，噪声通常表现为较小的高频系数。
阈值处理：
在DCT系数上应用阈值处理，以去除噪声。低于设定阈值的系数被置零，假设这些系数主要代表噪声分量。阈值设置为DCT系数中最大值的20%，这一比例可以根据具体应用场景调整。
逆DCT恢复信号：
应用逆DCT（IDCT）将处理后的DCT系数转换回时间域，恢复去噪后的信号。由于只保留了主要的信号成分，去除了噪声，因此恢复的信号质量应显著提高。
结果可视化：
代码通过多个子图展示了原始信号、含噪信号、阈值处理后的DCT系数以及去噪后的信号。这样可以直观地看到去噪过程的效果。
去噪效果评估：
使用均方误差（MSE）和信噪比（SNR）改进量来评估去噪效果。MSE越小，SNR提升越大，表明去噪效果越好。

效果

完整代码获取

关注下方卡片微信公众号，回复"DCT"获取完整代码