程序中如何實(shí)現(xiàn)濾波、峰值檢測(cè)和心率計(jì)算

繼續(xù)講解心率算法！從大西北回到南京后，本人就一直在完善之前寫的程序，終于在十月初完成一版了。本文將給出算法程序的全部?jī)?nèi)容！當(dāng)然給出的只是浮點(diǎn)程序，如果要轉(zhuǎn)為單片機(jī)中的C程序，首先要進(jìn)行定點(diǎn)處理。先來看看程序中如何實(shí)現(xiàn)濾波、峰值檢測(cè)和心率計(jì)算的！只有先完成性能仿真，符合了要求，才能進(jìn)行下一階段的工作。接著上篇文章的內(nèi)容繼續(xù)講解，首先來看看算法程序！算法程序的主要功能是能夠從PPG數(shù)據(jù)中得到準(zhǔn)確的心率值！文章末尾給出數(shù)據(jù)的下載鏈接！

% PPG_25hz_dsp_6v2

% data：2021-10-8

%%% 程序說明

% 讀取csv文件！

% 濾波采用差分算法！

% 采用峰值比較算法3

clc;clear all;close all;

%---------------- program main --------------------------%

%%%%%%%%%%%%% 數(shù)據(jù)讀取

%%%%%%% 采集數(shù)據(jù)

%%% 文件名

PPG_filename = '20210814-2229.csv';

% PPG數(shù)據(jù) 文件已經(jīng)上傳到網(wǎng)盤！

%%% 讀入數(shù)據(jù)

PPG_data = readmatrix(PPG_filename);

%%% 心率數(shù)據(jù)

% 采樣頻率 25Hz

Fs = 25 ;

heart_rate_data = PPG_data(:,2);

% time 959second

%%%%%%%%%%%%%%%%%% 單片機(jī)實(shí)現(xiàn)部分

%%%%%%%%%%% 數(shù)據(jù)濾波處理

%%% 差分

% hr_data1 = diff(heart_rate_data);

%%% 五點(diǎn)差分濾波 二階導(dǎo)數(shù)

責(zé)任編輯：haq