dogukan barkan / Mbed 2 deprecated HSP_SpO2

SpO2 algorithm used in MAX30101 onboard MAX32620HSP. Algorithm is not perfect but it works

Dependencies:   mbed MAX14720 USBDevice

algorithm/algorithm.h@0:af4d13000e95, 2019-03-28 (annotated)

Committer:
douqan93
Date:
Thu Mar 28 08:07:54 2019 +0000
Revision:
0:af4d13000e95
This program for MAX32620HSP's onboard SpO2 sensor (MAX30101). Algorithm is not perfect but it works.

Who changed what in which revision?

UserRevisionLine numberNew contents of line
douqan93 0:af4d13000e95 1 /** \file algorithm.h ******************************************************
douqan93 0:af4d13000e95 2 *
douqan93 0:af4d13000e95 3 * Project: MAXREFDES117#
douqan93 0:af4d13000e95 4 * Filename: algorithm.c
douqan93 0:af4d13000e95 5 * Description: This module is the heart rate/SpO2 calculation algorithm header file
douqan93 0:af4d13000e95 6 *
douqan93 0:af4d13000e95 7 * Revision History:
douqan93 0:af4d13000e95 8 *\n 1-18-2016 Rev 01.00 SK Initial release.
douqan93 0:af4d13000e95 9 *\n
douqan93 0:af4d13000e95 10 *
douqan93 0:af4d13000e95 11 * --------------------------------------------------------------------
douqan93 0:af4d13000e95 12 *
douqan93 0:af4d13000e95 13 * This code follows the following naming conventions:
douqan93 0:af4d13000e95 14 *
douqan93 0:af4d13000e95 15 *\n char ch_pmod_value
douqan93 0:af4d13000e95 16 *\n char (array) s_pmod_s_string[16]
douqan93 0:af4d13000e95 17 *\n float f_pmod_value
douqan93 0:af4d13000e95 18 *\n int32_t n_pmod_value
douqan93 0:af4d13000e95 19 *\n int32_t (array) an_pmod_value[16]
douqan93 0:af4d13000e95 20 *\n int16_t w_pmod_value
douqan93 0:af4d13000e95 21 *\n int16_t (array) aw_pmod_value[16]
douqan93 0:af4d13000e95 22 *\n uint16_t uw_pmod_value
douqan93 0:af4d13000e95 23 *\n uint16_t (array) auw_pmod_value[16]
douqan93 0:af4d13000e95 24 *\n uint8_t uch_pmod_value
douqan93 0:af4d13000e95 25 *\n uint8_t (array) auch_pmod_buffer[16]
douqan93 0:af4d13000e95 26 *\n uint32_t un_pmod_value
douqan93 0:af4d13000e95 27 *\n int32_t * pn_pmod_value
douqan93 0:af4d13000e95 28 *
douqan93 0:af4d13000e95 29 * ------------------------------------------------------------------------- */
douqan93 0:af4d13000e95 30 /*******************************************************************************
douqan93 0:af4d13000e95 31 * Copyright (C) 2015 Maxim Integrated Products, Inc., All Rights Reserved.
douqan93 0:af4d13000e95 32 *
douqan93 0:af4d13000e95 33 * Permission is hereby granted, free of charge, to any person obtaining a
douqan93 0:af4d13000e95 34 * copy of this software and associated documentation files (the "Software"),
douqan93 0:af4d13000e95 35 * to deal in the Software without restriction, including without limitation
douqan93 0:af4d13000e95 36 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
douqan93 0:af4d13000e95 37 * and/or sell copies of the Software, and to permit persons to whom the
douqan93 0:af4d13000e95 38 * Software is furnished to do so, subject to the following conditions:
douqan93 0:af4d13000e95 39 *
douqan93 0:af4d13000e95 40 * The above copyright notice and this permission notice shall be included
douqan93 0:af4d13000e95 41 * in all copies or substantial portions of the Software.
douqan93 0:af4d13000e95 42 *
douqan93 0:af4d13000e95 43 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
douqan93 0:af4d13000e95 44 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
douqan93 0:af4d13000e95 45 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
douqan93 0:af4d13000e95 46 * IN NO EVENT SHALL MAXIM INTEGRATED BE LIABLE FOR ANY CLAIM, DAMAGES
douqan93 0:af4d13000e95 47 * OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
douqan93 0:af4d13000e95 48 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
douqan93 0:af4d13000e95 49 * OTHER DEALINGS IN THE SOFTWARE.
douqan93 0:af4d13000e95 50 *
douqan93 0:af4d13000e95 51 * Except as contained in this notice, the name of Maxim Integrated
douqan93 0:af4d13000e95 52 * Products, Inc. shall not be used except as stated in the Maxim Integrated
douqan93 0:af4d13000e95 53 * Products, Inc. Branding Policy.
douqan93 0:af4d13000e95 54 *
douqan93 0:af4d13000e95 55 * The mere transfer of this software does not imply any licenses
douqan93 0:af4d13000e95 56 * of trade secrets, proprietary technology, copyrights, patents,
douqan93 0:af4d13000e95 57 * trademarks, maskwork rights, or any other form of intellectual
douqan93 0:af4d13000e95 58 * property whatsoever. Maxim Integrated Products, Inc. retains all
douqan93 0:af4d13000e95 59 * ownership rights.
douqan93 0:af4d13000e95 60 *******************************************************************************
douqan93 0:af4d13000e95 61 */
douqan93 0:af4d13000e95 62 #ifndef ALGORITHM_H_
douqan93 0:af4d13000e95 63 #define ALGORITHM_H_
douqan93 0:af4d13000e95 64
douqan93 0:af4d13000e95 65 #include "mbed.h"
douqan93 0:af4d13000e95 66
douqan93 0:af4d13000e95 67 #define true 1
douqan93 0:af4d13000e95 68 #define false 0
douqan93 0:af4d13000e95 69 #define FS 100
douqan93 0:af4d13000e95 70 #define BUFFER_SIZE (FS* 5)
douqan93 0:af4d13000e95 71 #define HR_FIFO_SIZE 7
douqan93 0:af4d13000e95 72 #define MA4_SIZE 4 // DO NOT CHANGE
douqan93 0:af4d13000e95 73 #define HAMMING_SIZE 5// DO NOT CHANGE
douqan93 0:af4d13000e95 74 #define min(x,y) ((x) < (y) ? (x) : (y))
douqan93 0:af4d13000e95 75
douqan93 0:af4d13000e95 76 const uint16_t auw_hamm[31]={ 41, 276, 512, 276, 41 }; //Hamm= long16(512* hamming(5)');
douqan93 0:af4d13000e95 77 //SPO2table is computed as -45.060*ratioAverage* ratioAverage + 30.354 *ratioAverage + 94.845 ;
douqan93 0:af4d13000e95 78 const uint8_t uch_spo2_table[184]={ 95, 95, 95, 96, 96, 96, 97, 97, 97, 97, 97, 98, 98, 98, 98, 98, 99, 99, 99, 99,
douqan93 0:af4d13000e95 79 99, 99, 99, 99, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100,
douqan93 0:af4d13000e95 80 100, 100, 100, 100, 99, 99, 99, 99, 99, 99, 99, 99, 98, 98, 98, 98, 98, 98, 97, 97,
douqan93 0:af4d13000e95 81 97, 97, 96, 96, 96, 96, 95, 95, 95, 94, 94, 94, 93, 93, 93, 92, 92, 92, 91, 91,
douqan93 0:af4d13000e95 82 90, 90, 89, 89, 89, 88, 88, 87, 87, 86, 86, 85, 85, 84, 84, 83, 82, 82, 81, 81,
douqan93 0:af4d13000e95 83 80, 80, 79, 78, 78, 77, 76, 76, 75, 74, 74, 73, 72, 72, 71, 70, 69, 69, 68, 67,
douqan93 0:af4d13000e95 84 66, 66, 65, 64, 63, 62, 62, 61, 60, 59, 58, 57, 56, 56, 55, 54, 53, 52, 51, 50,
douqan93 0:af4d13000e95 85 49, 48, 47, 46, 45, 44, 43, 42, 41, 40, 39, 38, 37, 36, 35, 34, 33, 31, 30, 29,
douqan93 0:af4d13000e95 86 28, 27, 26, 25, 23, 22, 21, 20, 19, 17, 16, 15, 14, 12, 11, 10, 9, 7, 6, 5,
douqan93 0:af4d13000e95 87 3, 2, 1 } ;
douqan93 0:af4d13000e95 88 static int32_t an_dx[ BUFFER_SIZE-MA4_SIZE]; // delta
douqan93 0:af4d13000e95 89 static int32_t an_x[ BUFFER_SIZE]; //ir
douqan93 0:af4d13000e95 90 static int32_t an_y[ BUFFER_SIZE]; //red
douqan93 0:af4d13000e95 91
douqan93 0:af4d13000e95 92
douqan93 0:af4d13000e95 93 void maxim_heart_rate_and_oxygen_saturation(uint32_t *pun_ir_buffer , int32_t n_ir_buffer_length, uint32_t *pun_red_buffer , int32_t *pn_spo2, int8_t *pch_spo2_valid , int32_t *pn_heart_rate , int8_t *pch_hr_valid);
douqan93 0:af4d13000e95 94 void maxim_find_peaks( int32_t *pn_locs, int32_t *pn_npks, int32_t *pn_x, int32_t n_size, int32_t n_min_height, int32_t n_min_distance, int32_t n_max_num );
douqan93 0:af4d13000e95 95 void maxim_peaks_above_min_height( int32_t *pn_locs, int32_t *pn_npks, int32_t *pn_x, int32_t n_size, int32_t n_min_height );
douqan93 0:af4d13000e95 96 void maxim_remove_close_peaks( int32_t *pn_locs, int32_t *pn_npks, int32_t *pn_x, int32_t n_min_distance );
douqan93 0:af4d13000e95 97 void maxim_sort_ascend( int32_t *pn_x, int32_t n_size );
douqan93 0:af4d13000e95 98 void maxim_sort_indices_descend( int32_t *pn_x, int32_t *pn_indx, int32_t n_size);
douqan93 0:af4d13000e95 99
douqan93 0:af4d13000e95 100 #endif /* ALGORITHM_H_ */