MATLAB已知PPG+ACC的叠加信号和ACC信号，如何分离出原始PPG信号？

Question

在网络上有查到用fft处理后逆ifft还原不同频率的信号，由于ppg和acc频率有重叠，但现在ACC信号是已知的，怎样还原出原始PPG信号？
MATLAB:

上图是MATLAB spectrogram plot出叠加频谱和ACC频谱，如何从叠加频谱中分离出ACC得到原始PPG频谱？

Answer 1

以下回答参考 皆我百晓生、券券喵儿 等免费微信小程序相关内容作答，并由本人整理回复。

要从MATLAB的spectrogram图中提取出原始PPG信号，可以使用以下步骤：

1. 首先，计算原始PPG信号的FFT（快速傅里叶变换）。
2. 然后，将FFT结果与ACC信号进行点对点比较。如果两个信号在某个特定频率处具有相似的幅度，则可以认为它们之间存在一个对应关系。

下面是实现这些步骤的MATLAB代码示例：

% 读取spectrogram数据
spectrogram_data = load('spectrogram.mat');
spectrogram = spectrogram_data.spectrogram;

% 计算原始PPG信号的FFT
original_ppg_fft = fft(spectrogram);

% 提取ACC信号的频带范围
acc_frequency = [700, 700];
acc_range = [0, 2]; % ACC信号的频率范围

% 检测ACC信号是否包含在原始PPG FFT中
for i = 1:length(acc_frequency)
    frequency_idx = find(abs(original_ppg_fft) > acc_range(i));
    
    if ~isempty(frequency_idx)
        % 如果找到了匹配的频率，说明找到了原始PPG信号
        original_ppg_signal = original_ppg_fft(frequency_idx);
        break;
    end
end

% 输出原始PPG信号
disp(['Original PPG signal: ', num2str(original_ppg_signal)]);

这段代码首先加载spectrogram数据，并将其转换为FFT格式。然后，它计算原始PPG信号的FFT，找到其中包含ACC信号的频率。最后，它输出这些频率对应的原始PPG信号值。

请注意，这个方法可能不适用于所有情况，因为ACC信号的频率范围可能覆盖了原始PPG信号的部分或全部频率范围。此外，这种方法依赖于FFT的结果与ACC信号的直接匹配，这在某些情况下可能会导致错误的识别。因此，在实际应用中，你可能需要根据实际情况调整参数以获得更好的效果。