PPG信号IIR低通滤波matlab实现

% 1. 导入数据
A=xlsread('baselinedata.xlsx'); %基线漂移明显的PPG信号

Fs=100;                         %采样频率，单位为Hz


% 2.IIR5Hz低通滤波，去高频噪声
A1=A(101:2148);
LA1=length(A1);                 %数据长度
tA1=(1:LA1)/Fs;                 %产生时间序列，A1的横坐标


% 2.1 对输入信号进行频率变换
A2=fft(A1,LA1);                 %计算傅立叶变换
ftA2=(2:LA1)*Fs/LA1;            %产生频率序列，A2的横坐标，以Hz为单位
A2=A2(2:length(A2));            %去掉fft变换后异常的第一个数据


figure('Name','原始信号时频图');  
subplot(2,1,1);
plot(tA1,A1); 
xlabel('时间(s)'); 
ylabel('幅值');
title('原始信号时域波形');


subplot(2,1,2);
plot(ftA2,abs(A2));              %求得傅里叶变换后的A2振幅，abs：求绝对值