% 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:求绝对值
xlabel('频率(Hz)');
ylabel('幅值');
title('原始信号频域波形');
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:求绝对值
xlabel('频率(Hz)');
ylabel('幅值');
title('原始信号频域波形');