繼續(xù)講解心率算法！從大西北回到南京后，本人就一直在完善之前寫的程序，終于在十月初完成一版了。本文將給出算法程序的全部內(nèi)容！當(dāng)然給出的只是浮點程序，如果要轉(zhuǎn)為單片機中的C程序，首先要進行定點處理。先來看看程序中如何實現(xià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

%%%%%%%%%%%%%%%%%% 單片機實現(xiàn)部分

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

%%% 差分

% hr_data1 = diff(heart_rate_data);

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

責(zé)任編輯：haq