Junichi Katsu
心拍・酸素飽和度モニタモジュール MAXREFDES117から取得した心拍の値でmicro:bit上のLEDに表示したハートを鼓動させるプログラムです。
Fork of
microbit-Heartbeat-serial
by 
Junichi Katsu
algorithm/algorithm.h@1:83ace7df2c63, 2018-04-16 (annotated)
- Committer:
- JKsoft_main
- Date:
- Mon Apr 16 16:52:36 2018 +0000
- Revision:
- 1:83ace7df2c63
??
Who changed what in which revision?
| User | Revision | Line number | New contents of line |
|---|---|---|---|
| JKsoft_main | 1:83ace7df2c63 | 1 | /** \file algorithm.h ****************************************************** |
| JKsoft_main | 1:83ace7df2c63 | 2 | * |
| JKsoft_main | 1:83ace7df2c63 | 3 | * Project: MAXREFDES117# |
| JKsoft_main | 1:83ace7df2c63 | 4 | * Filename: algorithm.c |
| JKsoft_main | 1:83ace7df2c63 | 5 | * Description: This module is the heart rate/SpO2 calculation algorithm header file |
| JKsoft_main | 1:83ace7df2c63 | 6 | * |
| JKsoft_main | 1:83ace7df2c63 | 7 | * Revision History: |
| JKsoft_main | 1:83ace7df2c63 | 8 | *\n 1-18-2016 Rev 01.00 SK Initial release. |
| JKsoft_main | 1:83ace7df2c63 | 9 | *\n |
| JKsoft_main | 1:83ace7df2c63 | 10 | * |
| JKsoft_main | 1:83ace7df2c63 | 11 | * -------------------------------------------------------------------- |
| JKsoft_main | 1:83ace7df2c63 | 12 | * |
| JKsoft_main | 1:83ace7df2c63 | 13 | * This code follows the following naming conventions: |
| JKsoft_main | 1:83ace7df2c63 | 14 | * |
| JKsoft_main | 1:83ace7df2c63 | 15 | *\n char ch_pmod_value |
| JKsoft_main | 1:83ace7df2c63 | 16 | *\n char (array) s_pmod_s_string[16] |
| JKsoft_main | 1:83ace7df2c63 | 17 | *\n float f_pmod_value |
| JKsoft_main | 1:83ace7df2c63 | 18 | *\n int32_t n_pmod_value |
| JKsoft_main | 1:83ace7df2c63 | 19 | *\n int32_t (array) an_pmod_value[16] |
| JKsoft_main | 1:83ace7df2c63 | 20 | *\n int16_t w_pmod_value |
| JKsoft_main | 1:83ace7df2c63 | 21 | *\n int16_t (array) aw_pmod_value[16] |
| JKsoft_main | 1:83ace7df2c63 | 22 | *\n uint16_t uw_pmod_value |
| JKsoft_main | 1:83ace7df2c63 | 23 | *\n uint16_t (array) auw_pmod_value[16] |
| JKsoft_main | 1:83ace7df2c63 | 24 | *\n uint8_t uch_pmod_value |
| JKsoft_main | 1:83ace7df2c63 | 25 | *\n uint8_t (array) auch_pmod_buffer[16] |
| JKsoft_main | 1:83ace7df2c63 | 26 | *\n uint32_t un_pmod_value |
| JKsoft_main | 1:83ace7df2c63 | 27 | *\n int32_t * pn_pmod_value |
| JKsoft_main | 1:83ace7df2c63 | 28 | * |
| JKsoft_main | 1:83ace7df2c63 | 29 | * ------------------------------------------------------------------------- */ |
| JKsoft_main | 1:83ace7df2c63 | 30 | /******************************************************************************* |
| JKsoft_main | 1:83ace7df2c63 | 31 | * Copyright (C) 2015 Maxim Integrated Products, Inc., All Rights Reserved. |
| JKsoft_main | 1:83ace7df2c63 | 32 | * |
| JKsoft_main | 1:83ace7df2c63 | 33 | * Permission is hereby granted, free of charge, to any person obtaining a |
| JKsoft_main | 1:83ace7df2c63 | 34 | * copy of this software and associated documentation files (the "Software"), |
| JKsoft_main | 1:83ace7df2c63 | 35 | * to deal in the Software without restriction, including without limitation |
| JKsoft_main | 1:83ace7df2c63 | 36 | * the rights to use, copy, modify, merge, publish, distribute, sublicense, |
| JKsoft_main | 1:83ace7df2c63 | 37 | * and/or sell copies of the Software, and to permit persons to whom the |
| JKsoft_main | 1:83ace7df2c63 | 38 | * Software is furnished to do so, subject to the following conditions: |
| JKsoft_main | 1:83ace7df2c63 | 39 | * |
| JKsoft_main | 1:83ace7df2c63 | 40 | * The above copyright notice and this permission notice shall be included |
| JKsoft_main | 1:83ace7df2c63 | 41 | * in all copies or substantial portions of the Software. |
| JKsoft_main | 1:83ace7df2c63 | 42 | * |
| JKsoft_main | 1:83ace7df2c63 | 43 | * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS |
| JKsoft_main | 1:83ace7df2c63 | 44 | * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF |
| JKsoft_main | 1:83ace7df2c63 | 45 | * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. |
| JKsoft_main | 1:83ace7df2c63 | 46 | * IN NO EVENT SHALL MAXIM INTEGRATED BE LIABLE FOR ANY CLAIM, DAMAGES |
| JKsoft_main | 1:83ace7df2c63 | 47 | * OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, |
| JKsoft_main | 1:83ace7df2c63 | 48 | * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR |
| JKsoft_main | 1:83ace7df2c63 | 49 | * OTHER DEALINGS IN THE SOFTWARE. |
| JKsoft_main | 1:83ace7df2c63 | 50 | * |
| JKsoft_main | 1:83ace7df2c63 | 51 | * Except as contained in this notice, the name of Maxim Integrated |
| JKsoft_main | 1:83ace7df2c63 | 52 | * Products, Inc. shall not be used except as stated in the Maxim Integrated |
| JKsoft_main | 1:83ace7df2c63 | 53 | * Products, Inc. Branding Policy. |
| JKsoft_main | 1:83ace7df2c63 | 54 | * |
| JKsoft_main | 1:83ace7df2c63 | 55 | * The mere transfer of this software does not imply any licenses |
| JKsoft_main | 1:83ace7df2c63 | 56 | * of trade secrets, proprietary technology, copyrights, patents, |
| JKsoft_main | 1:83ace7df2c63 | 57 | * trademarks, maskwork rights, or any other form of intellectual |
| JKsoft_main | 1:83ace7df2c63 | 58 | * property whatsoever. Maxim Integrated Products, Inc. retains all |
| JKsoft_main | 1:83ace7df2c63 | 59 | * ownership rights. |
| JKsoft_main | 1:83ace7df2c63 | 60 | ******************************************************************************* |
| JKsoft_main | 1:83ace7df2c63 | 61 | */ |
| JKsoft_main | 1:83ace7df2c63 | 62 | #ifndef ALGORITHM_H_ |
| JKsoft_main | 1:83ace7df2c63 | 63 | #define ALGORITHM_H_ |
| JKsoft_main | 1:83ace7df2c63 | 64 | |
| JKsoft_main | 1:83ace7df2c63 | 65 | #include "mbed.h" |
| JKsoft_main | 1:83ace7df2c63 | 66 | |
| JKsoft_main | 1:83ace7df2c63 | 67 | #define true 1 |
| JKsoft_main | 1:83ace7df2c63 | 68 | #define false 0 |
| JKsoft_main | 1:83ace7df2c63 | 69 | #define FS 2 |
| JKsoft_main | 1:83ace7df2c63 | 70 | #define BUFFER_SIZE (FS* 5) |
| JKsoft_main | 1:83ace7df2c63 | 71 | #define HR_FIFO_SIZE 7 |
| JKsoft_main | 1:83ace7df2c63 | 72 | #define MA4_SIZE 4 // DO NOT CHANGE |
| JKsoft_main | 1:83ace7df2c63 | 73 | #define HAMMING_SIZE 5// DO NOT CHANGE |
| JKsoft_main | 1:83ace7df2c63 | 74 | #define min(x,y) ((x) < (y) ? (x) : (y)) |
| JKsoft_main | 1:83ace7df2c63 | 75 | |
| JKsoft_main | 1:83ace7df2c63 | 76 | const uint16_t auw_hamm[31]={ 41, 276, 512, 276, 41 }; //Hamm= long16(512* hamming(5)'); |
| JKsoft_main | 1:83ace7df2c63 | 77 | //SPO2table is computed as -45.060*ratioAverage* ratioAverage + 30.354 *ratioAverage + 94.845 ; |
| JKsoft_main | 1:83ace7df2c63 | 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, |
| JKsoft_main | 1:83ace7df2c63 | 79 | 99, 99, 99, 99, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100, |
| JKsoft_main | 1:83ace7df2c63 | 80 | 100, 100, 100, 100, 99, 99, 99, 99, 99, 99, 99, 99, 98, 98, 98, 98, 98, 98, 97, 97, |
| JKsoft_main | 1:83ace7df2c63 | 81 | 97, 97, 96, 96, 96, 96, 95, 95, 95, 94, 94, 94, 93, 93, 93, 92, 92, 92, 91, 91, |
| JKsoft_main | 1:83ace7df2c63 | 82 | 90, 90, 89, 89, 89, 88, 88, 87, 87, 86, 86, 85, 85, 84, 84, 83, 82, 82, 81, 81, |
| JKsoft_main | 1:83ace7df2c63 | 83 | 80, 80, 79, 78, 78, 77, 76, 76, 75, 74, 74, 73, 72, 72, 71, 70, 69, 69, 68, 67, |
| JKsoft_main | 1:83ace7df2c63 | 84 | 66, 66, 65, 64, 63, 62, 62, 61, 60, 59, 58, 57, 56, 56, 55, 54, 53, 52, 51, 50, |
| JKsoft_main | 1:83ace7df2c63 | 85 | 49, 48, 47, 46, 45, 44, 43, 42, 41, 40, 39, 38, 37, 36, 35, 34, 33, 31, 30, 29, |
| JKsoft_main | 1:83ace7df2c63 | 86 | 28, 27, 26, 25, 23, 22, 21, 20, 19, 17, 16, 15, 14, 12, 11, 10, 9, 7, 6, 5, |
| JKsoft_main | 1:83ace7df2c63 | 87 | 3, 2, 1 } ; |
| JKsoft_main | 1:83ace7df2c63 | 88 | static int32_t an_dx[ BUFFER_SIZE-MA4_SIZE]; // delta |
| JKsoft_main | 1:83ace7df2c63 | 89 | static int32_t an_x[ BUFFER_SIZE]; //ir |
| JKsoft_main | 1:83ace7df2c63 | 90 | static int32_t an_y[ BUFFER_SIZE]; //red |
| JKsoft_main | 1:83ace7df2c63 | 91 | |
| JKsoft_main | 1:83ace7df2c63 | 92 | |
| JKsoft_main | 1:83ace7df2c63 | 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); |
| JKsoft_main | 1:83ace7df2c63 | 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 ); |
| JKsoft_main | 1:83ace7df2c63 | 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 ); |
| JKsoft_main | 1:83ace7df2c63 | 96 | void maxim_remove_close_peaks( int32_t *pn_locs, int32_t *pn_npks, int32_t *pn_x, int32_t n_min_distance ); |
| JKsoft_main | 1:83ace7df2c63 | 97 | void maxim_sort_ascend( int32_t *pn_x, int32_t n_size ); |
| JKsoft_main | 1:83ace7df2c63 | 98 | void maxim_sort_indices_descend( int32_t *pn_x, int32_t *pn_indx, int32_t n_size); |
| JKsoft_main | 1:83ace7df2c63 | 99 | |
| JKsoft_main | 1:83ace7df2c63 | 100 | #endif /* ALGORITHM_H_ */ |