stm32l475+MAX30102心率血氧测量+发送数据到串口调试助手

STM32L475与MAX30102实现心率血氧监测

最新推荐文章于 2025-12-06 15:55:34 发布

原创最新推荐文章于 2025-12-06 15:55:34 发布 · 463 阅读

6 ·

CC 4.0 BY-SA版权

文章标签：

#stm32 #嵌入式硬件 #单片机

该文章已生成可运行项目，

一、MAX30102心率血氧传感器介绍
在健康监测领域，脉搏血氧仪和心率传感器是至关重要的部件。MAX30102作为一款集成度高、性能出色的传感器，能够精确测量人体的心率和血氧饱和度，广泛应用于可穿戴设备、医疗监护仪器等产品中。了解并掌握MAX30102的使用方法，对于开发具有创新性的健康监测设备具有重要意义。

二、MAX30102心率血氧传感器工作原理介绍

（一）光电容积脉搏波原理（PPG）

MAX30102利用光电容积脉搏波原理来检测脉搏信号。它内置了红外LED和红光LED，通过向人体组织发射光线，并检测反射或透射回来的光强度变化。当心脏跳动时，血管内的血液量会发生周期性变化，这导致对光线的吸收和散射程度也随之改变。传感器接收到的光强度变化信号经过处理后，可得到脉搏波信号。

（二）血氧饱和度计算原理

血氧饱和度（SpO2）是指血液中氧合血红蛋白（HbO2）占总血红蛋白（Hb）的百分比。MAX30102通过同时检测红外光和红光的吸收情况，利用特定的算法来计算血氧饱和度。由于HbO2和Hb对红外光和红光的吸收特性不同，根据朗伯 - 比尔定律，通过测量两种光的吸收比例，可推算出血氧饱和度的值。

三、stm32cubemx配置

四、代码

main.c

#include "stdint.h"
/* USER CODE END Includes */

/* Private typedef -----------------------------------------------------------*/
/* USER CODE BEGIN PTD */

/* USER CODE END PTD */

/* Private define ------------------------------------------------------------*/
/* USER CODE BEGIN PD */

/* USER CODE END PD */

/* Private macro -------------------------------------------------------------*/
/* USER CODE BEGIN PM */

/* USER CODE END PM */

/* Private variables ---------------------------------------------------------*/

/* USER CODE BEGIN PV */
extern int heart;
extern float SpO2;
/* USER CODE END PV */

/* Private function prototypes -----------------------------------------------*/
void SystemClock_Config(void);
/* USER CODE BEGIN PFP */

/* USER CODE END PFP */

/* Private user code ---------------------------------------------------------*/
/* USER CODE BEGIN 0 */
 int fputc(int ch,FILE *f)
 {
	 HAL_UART_Transmit(&huart1,(uint8_t *)&ch,1,1000);
	 return ch;
 }
/* USER CODE END 0 */

/**
  * @brief  The application entry point.
  * @retval int
  */
int main(void)
{

  /* USER CODE BEGIN 1 */

  /* USER CODE END 1 */

  /* MCU Configuration--------------------------------------------------------*/

  /* Reset of all peripherals, Initializes the Flash interface and the Systick. */
  HAL_Init(); 

  /* USER CODE BEGIN Init */

  /* USER CODE END Init */

  /* Configure the system clock */
  SystemClock_Config();

  /* USER CODE BEGIN SysInit */

  /* USER CODE END SysInit */

  /* Initialize all configured peripherals */
  MX_GPIO_Init();
  MX_USART1_UART_Init();
  MX_I2C2_Init();
  /* USER CODE BEGIN 2 */
  Max30102_reset();     //����max30102
	MAX30102_Config();    //����max30102
	printf("MAX30102初始化完成\n");
  /* USER CODE END 2 */

  /* Infinite loop */
  /* USER CODE BEGIN WHILE */
  while (1)
  {
		blood_Loop();
    /* USER CODE END WHILE */

    /* USER CODE BEGIN 3 */
  }
  /* USER CODE END 3 */
}

/**
  * @brief System Clock Configuration
  * @retval None
  */
void SystemClock_Config(void)
{
  RCC_OscInitTypeDef RCC_OscInitStruct = {0};
  RCC_ClkInitTypeDef RCC_ClkInitStruct = {0};

  /** Configure the main internal regulator output voltage
  */
  if (HAL_PWREx_ControlVoltageScaling(PWR_REGULATOR_VOLTAGE_SCALE1) != HAL_OK)
  {
    Error_Handler();
  }

  /** Initializes the RCC Oscillators according to the specified parameters
  * in the RCC_OscInitTypeDef structure.
  */
  RCC_OscInitStruct.OscillatorType = RCC_OSCILLATORTYPE_HSE;
  RCC_OscInitStruct.HSEState = RCC_HSE_ON;
  RCC_OscInitStruct.PLL.PLLState = RCC_PLL_ON;
  RCC_OscInitStruct.PLL.PLLSource = RCC_PLLSOURCE_HSE;
  RCC_OscInitStruct.PLL.PLLM = 1;
  RCC_OscInitStruct.PLL.PLLN = 18;
  RCC_OscInitStruct.PLL.PLLP = RCC_PLLP_DIV7;
  RCC_OscInitStruct.PLL.PLLQ = RCC_PLLQ_DIV2;
  RCC_OscInitStruct.PLL.PLLR = RCC_PLLR_DIV2;
  if (HAL_RCC_OscConfig(&RCC_OscInitStruct) != HAL_OK)
  {
    Error_Handler();
  }

  /** Initializes the CPU, AHB and APB buses clocks
  */
  RCC_ClkInitStruct.ClockType = RCC_CLOCKTYPE_HCLK|RCC_CLOCKTYPE_SYSCLK
                              |RCC_CLOCKTYPE_PCLK1|RCC_CLOCKTYPE_PCLK2;
  RCC_ClkInitStruct.SYSCLKSource = RCC_SYSCLKSOURCE_PLLCLK;
  RCC_ClkInitStruct.AHBCLKDivider = RCC_SYSCLK_DIV1;
  RCC_ClkInitStruct.APB1CLKDivider = RCC_HCLK_DIV4;
  RCC_ClkInitStruct.APB2CLKDivider = RCC_HCLK_DIV2;

  if (HAL_RCC_ClockConfig(&RCC_ClkInitStruct, FLASH_LATENCY_4) != HAL_OK)
  {
    Error_Handler();
  }
}

/* USER CODE BEGIN 4 */

/* USER CODE END 4 */

/**
  * @brief  This function is executed in case of error occurrence.
  * @retval None
  */
void Error_Handler(void)
{
  /* USER CODE BEGIN Error_Handler_Debug */
  /* User can add his own implementation to report the HAL error return state */
  __disable_irq();
  while (1)
  {
  }
  /* USER CODE END Error_Handler_Debug */
}

#ifdef  USE_FULL_ASSERT
/**
  * @brief  Reports the name of the source file and the source line number
  *         where the assert_param error has occurred.
  * @param  file: pointer to the source file name
  * @param  line: assert_param error line source number
  * @retval None
  */
void assert_failed(uint8_t *file, uint32_t line)
{
  /* USER CODE BEGIN 6 */
  /* User can add his own implementation to report the file name and line number,
     ex: printf("Wrong parameters value: file %s on line %d\r\n", file, line) */
  /* USER CODE END 6 */
}
#endif /* USE_FULL_ASSERT */

max30102.c

/**
 * ************************************************************************
 * 
 * @file MAX30102.c
 * @author zxr
 * @brief 
 * 
 * ************************************************************************
 * @copyright Copyright (c) 2024 zxr 
 * ************************************************************************
 */
#include "MAX30102.h"

uint16_t fifo_red;  //¶¨ÒåFIFOÖÐµÄºì¹âÊý¾Ý
uint16_t fifo_ir;   //¶¨ÒåFIFOÖÐµÄºìÍâ¹âÊý¾Ý

/**
 * ************************************************************************
 * @brief ÏòMAX30102¼Ä´æÆ÷Ð´ÈëÒ»¸öÖµ
 * 
 * @param[in] addr  ¼Ä´æÆ÷µØÖ·
 * @param[in] data  ´«ÊäÊý¾Ý
 * 
 * @return 
 * ************************************************************************
 */
uint8_t max30102_write_reg(uint8_t addr, uint8_t data)
{
  HAL_I2C_Mem_Write(&hi2c2, MAX30102_Device_address,	addr,	1,	&data,1,HAL_MAX_DELAY);
  return 1;
}


/**
 * ************************************************************************
 * @brief ¶ÁÈ¡MAX30102¼Ä´æÆ÷µÄÒ»¸öÖµ
 * 
 * @param[in] addr  ¼Ä´æÆ÷µØÖ·
 * 
 * @return 
 * ************************************************************************
 */
uint8_t max30102_read_reg(uint8_t addr )
{
  uint8_t data=0;
  HAL_I2C_Mem_Read(&hi2c2, MAX30102_Device_address, addr, 1, &data, 1, HAL_MAX_DELAY);
  return data;
}


/**
 * ************************************************************************
 * @brief MAX30102´«¸ÐÆ÷¸´Î»
 * 
 * 
 * @return 
 * ************************************************************************
 */
uint8_t Max30102_reset(void)
{
	if(max30102_write_reg(REG_MODE_CONFIG, 0x40))
        return 1;
    else
        return 0;    
}

/**
 * ************************************************************************
 * @brief MAX30102´«¸ÐÆ÷Ä£Ê½ÅäÖÃ
 * 
 * 
 * ************************************************************************
 */
void MAX30102_Config(void)
{
	max30102_write_reg(REG_INTR_ENABLE_1,0xc0); //INTR setting
	max30102_write_reg(REG_INTR_ENABLE_2,0x00);//
	max30102_write_reg(REG_FIFO_WR_PTR,0x00);//FIFO_WR_PTR[4:0]
	max30102_write_reg(REG_OVF_COUNTER,0x00);//OVF_COUNTER[4:0]
	max30102_write_reg(REG_FIFO_RD_PTR,0x00);//FIFO_RD_PTR[4:0]
	
	max30102_write_reg(REG_FIFO_CONFIG,0x0f);//sample avg = 1, fifo rollover=false, fifo almost full = 17
	max30102_write_reg(REG_MODE_CONFIG,0x03);//0x02 for Red only, 0x03 for SpO2 mode 0x07 multimode LED
	max30102_write_reg(REG_SPO2_CONFIG,0x27);	// SPO2_ADC range = 4096nA, SPO2 sample rate (50 Hz), LED pulseWidth (400uS)  
	max30102_write_reg(REG_LED1_PA,0x32);//Choose value for ~ 10mA for LED1
	max30102_write_reg(REG_LED2_PA,0x32);// Choose value for ~ 10mA for LED2
	max30102_write_reg(REG_PILOT_PA,0x7f);// Choose value for ~ 25mA for Pilot LED
}

/**
 * ************************************************************************
 * @brief ¶ÁÈ¡FIFO¼Ä´æÆ÷µÄÊý¾Ý
 * 
 * 
 * ************************************************************************
 */
void max30102_read_fifo(void)
{
  uint16_t un_temp;
  fifo_red=0;
  fifo_ir=0;
  uint8_t ach_i2c_data[6];
  
  //read and clear status register
  max30102_read_reg(REG_INTR_STATUS_1);
  max30102_read_reg(REG_INTR_STATUS_2);
  
  ach_i2c_data[0]=REG_FIFO_DATA;
	
	HAL_I2C_Mem_Read(&hi2c2,MAX30102_Device_address,REG_FIFO_DATA,1,ach_i2c_data,6,HAL_MAX_DELAY);
	
  un_temp=ach_i2c_data[0];
  un_temp<<=14;
  fifo_red+=un_temp;
  un_temp=ach_i2c_data[1];
  un_temp<<=6;
  fifo_red+=un_temp;
  un_temp=ach_i2c_data[2];
	un_temp>>=2;
  fifo_red+=un_temp;
  
  un_temp=ach_i2c_data[3];
  un_temp<<=14;
  fifo_ir+=un_temp;
  un_temp=ach_i2c_data[4];
  un_temp<<=6;
  fifo_ir+=un_temp;
  un_temp=ach_i2c_data[5];
	un_temp>>=2;
  fifo_ir+=un_temp;
	
	if(fifo_ir<=10000)
	{
		fifo_ir=0;
	}
	if(fifo_red<=10000)
	{
		fifo_red=0;
	}
}

max30102.h

/**
 * ************************************************************************
 * 
 * @file MAX30102.h
 * @author zxr
 * @brief 
 * 
 * ************************************************************************
 * @copyright Copyright (c) 2024 zxr 
 * ************************************************************************
 */
#ifndef _MAX30102_H
#define _MAX30102_H

#include "main.h"                  // Device header
#include "i2c.h"
#include "stdbool.h"

#define MAX30102_Device_address 			0xAE

//register addresses
#define REG_INTR_STATUS_1 	                0x00
#define REG_INTR_STATUS_2 	                0x01
#define REG_INTR_ENABLE_1 	                0x02
#define REG_INTR_ENABLE_2 	                0x03
#define REG_FIFO_WR_PTR 		            0x04
#define REG_OVF_COUNTER 		            0x05
#define REG_FIFO_RD_PTR 		            0x06
#define REG_FIFO_DATA 			            0x07
#define REG_FIFO_CONFIG 		            0x08
#define REG_MODE_CONFIG 		            0x09
#define REG_SPO2_CONFIG 		            0x0A
#define REG_LED1_PA 				        0x0C
#define REG_LED2_PA 				        0x0D
#define REG_PILOT_PA 				        0x10
#define REG_MULTI_LED_CTRL1                 0x11
#define REG_MULTI_LED_CTRL2                 0x12
#define REG_TEMP_INTR 			            0x1F
#define REG_TEMP_FRAC 			            0x20
#define REG_TEMP_CONFIG 		            0x21
#define REG_PROX_INT_THRESH                 0x30
#define REG_REV_ID 					        0xFE
#define REG_PART_ID 				        0xFF

#define SAMPLES_PER_SECOND 					100	//¼ì²âÆµÂÊ

uint8_t Max30102_reset(void);
void MAX30102_Config(void);
void max30102_read_fifo(void);

uint8_t max30102_write_reg(uint8_t addr, uint8_t data);
uint8_t max30102_read_reg(uint8_t addr );

#endif

algorithm.c

/**
 * ************************************************************************
 * 
 * @file algorithm.c
 * @author zxr
 * @brief 
 * 
 * ************************************************************************
 * @copyright Copyright (c) 2024 zxr 
 * ************************************************************************
 */
#include "algorithm.h"
#include "stm32l4xx.h"                  // Device header



#define XPI         (3.1415926535897932384626433832795)				//¶¨ÒåÔ²ÖÜÂÊÖµ£¬±£Áô31Î»
#define XENTRY      (100)
#define XINCL       (XPI/2/XENTRY)									//ÓÃÓÚÕýÏÒº¯ÊýµÄ¾«¶È¿ØÖÆ
#define PI 			3.1415926535897932384626433832795028841971    	//¶¨ÒåÔ²ÖÜÂÊÖµ£¬±£Áô40Î»


/**
 * ************************************************************************
 * @brief ¾²Ì¬ÕýÏÒÖµ¶ÔÓ¦±í
 * 
 * 	Ê¹ÓÃÔ¤ÏÈ¼ÆËãºÃµÄÊýÖµ¿ÉÒÔ½ÚÊ¡¼ÆËãÊ±¼ä
 * ************************************************************************
 */
static const double XSinTbl[] = 
{
	0.00000000000000000  , 0.015707317311820675 , 
	0.031410759078128292 , 0.047106450709642665 , 
	0.062790519529313374 , 0.078459095727844944 , 
	0.094108313318514325 , 0.10973431109104528  , 
	0.12533323356430426  , 0.14090123193758267  ,
	0.15643446504023087  , 0.17192910027940955  , 
	0.18738131458572463  , 0.20278729535651249  , 
	0.21814324139654256  , 0.23344536385590542  , 
	0.24868988716485479  , 0.26387304996537292  , 
	0.27899110603922928  , 0.29404032523230400  ,
	0.30901699437494740  , 0.32391741819814940  , 
	0.33873792024529142  , 0.35347484377925714  , 
	0.36812455268467797  , 0.38268343236508978  , 
	0.39714789063478062  , 0.41151435860510882  , 
	0.42577929156507272  , 0.43993916985591514  ,
	0.45399049973954680  , 0.46792981426057340  , 
	0.48175367410171532  , 0.49545866843240760  , 
	0.50904141575037132  , 0.52249856471594880  , 
	0.53582679497899666  , 0.54902281799813180  , 
	0.56208337785213058  , 0.57500525204327857  ,
	0.58778525229247314  , 0.60042022532588402  , 
	0.61290705365297649  , 0.62524265633570519  , 
	0.63742398974868975  , 0.64944804833018377  , 
	0.66131186532365183  , 0.67301251350977331  , 
	0.68454710592868873  , 0.69591279659231442  ,
	0.70710678118654757  , 0.71812629776318881  , 
	0.72896862742141155  , 0.73963109497860968  , 
	0.75011106963045959  , 0.76040596560003104  , 
	0.77051324277578925  , 0.78043040733832969  , 
	0.79015501237569041  , 0.79968465848709058  ,
	0.80901699437494745  , 0.81814971742502351  , 
	0.82708057427456183  , 0.83580736136827027  , 
	0.84432792550201508  , 0.85264016435409218  , 
	0.86074202700394364  , 0.86863151443819120  , 
	0.87630668004386369  , 0.88376563008869347  ,
	0.89100652418836779  , 0.89802757576061565  , 
	0.90482705246601958  , 0.91140327663544529  , 
	0.91775462568398114  , 0.92387953251128674  , 
	0.92977648588825146  , 0.93544403082986738  , 
	0.94088076895422557  , 0.94608535882754530  ,
	0.95105651629515353  , 0.95579301479833012  , 
	0.96029368567694307  , 0.96455741845779808  , 
	0.96858316112863108  , 0.97236992039767667  , 
	0.97591676193874743  , 0.97922281062176575  , 
	0.98228725072868872  , 0.98510932615477398  ,
	0.98768834059513777  , 0.99002365771655754  , 
	0.99211470131447788  , 0.99396095545517971  , 
	0.99556196460308000  , 0.99691733373312796  , 
	0.99802672842827156  , 0.99888987496197001  , 
	0.99950656036573160  , 0.99987663248166059  ,
	1.00000000000000000  
};


/**
 * ************************************************************************
 * @brief ÏòÏÂÈ¡Õûº¯Êý
 * 
 * @param[in] x  ÐèÒªÈ¡ÕûµÄ¸¡µãÊý²ÎÊý
 * 
 * @return 
 * ************************************************************************
 */
double my_floor(double x)
{
	double y=x;
    if( (*( ( (int *) &y)+1) & 0x80000000)  != 0) //»òÕßif(x<0)
        return (float)((int)x)-1;
    else
        return (float)((int)x);
}


/**
 * ************************************************************************
 * @brief È¡Óàº¯Êý
 * 
 * @param[in] x  ²ÎÊý1
 * @param[in] y  ²ÎÊý2
 * @note	±ÜÃâÁË¶Ô¸¡µãÊýµÄÖ±½Ó³ý·¨ÔËËã£¬´Ó¶øÌá¸ßÁËÐ§ÂÊ
 * @return 
 * ************************************************************************
 */
double my_fmod(double x, double y)
{
	double temp, ret;

	if (y == 0.0)
		return 0.0;
	temp = my_floor(x/y);
   ret = x - temp * y;
	if ((x < 0.0) != (y < 0.0))
		ret = ret - y;
	return ret;
}


/**
 * ************************************************************************
 * @brief ÕýÏÒº¯Êý
 * 
 * @param[in] x  ½Ç¶ÈÖµ
 * 
 * @return 
 * 
 * @note Í¨¹ý²é±íºÍÌ©ÀÕÕ¹¿ªÊ½¼ÆËãÕýÏÒÖµ£¬Ïà¶ÔÓÚÖ±½Óµ÷ÓÃ±ê×¼¿âº¯Êý£¬
 * 		 ¿ÉÄÜ»áÎþÉüÒ»Ð©¾«¶È£¬µ«ÔÚÄ³Ð©Ç¶ÈëÊ½ÏµÍ³ÖÐ¿ÉÄÜ¸ü¼Ó¸ßÐ§
 * ************************************************************************
 */
double XSin( double x )
{
    int s = 0 , n;
    double dx , sx , cx;
    if( x < 0 )
        s = 1 , x = -x;
    x = my_fmod( x , 2 * XPI );
    if( x > XPI )
        s = !s , x -= XPI;
    if( x > XPI / 2 )
        x = XPI - x;
    n = (int)( x / XINCL );
    dx = x - n * XINCL;
    if( dx > XINCL / 2 )
        ++n , dx -= XINCL;
    sx = XSinTbl[n];
    cx = XSinTbl[XENTRY-n];
    x = sx + dx*cx - (dx*dx)*sx/2
        - (dx*dx*dx)*cx/6 
        + (dx*dx*dx*dx)*sx/24;

    return s ? -x : x;
}


/**
 * ************************************************************************
 * @brief ÓàÏÒº¯Êý
 * 
 * @param[in] x  ½Ç¶ÈÖµ
 * 
 * @return 
 * ************************************************************************
 */
double XCos( double x )
{
    return XSin( x + XPI/2 );
}


/**
 * ************************************************************************
 * @brief ¿ªÆ½·½
 * 
 * @param[in] a  ²ÎÊý
 * 
 * @return 
 * ************************************************************************
 */
int qsqrt(int a)
{
	uint32_t rem = 0, root = 0, divisor = 0;
	uint16_t i;
	for(i=0; i<16; i++)
	{
		root <<= 1;
		rem = ((rem << 2) + (a>>30));
		a <<= 2;
		divisor = (root << 1) + 1;
		if(divisor <= rem)
		{
			rem -= divisor;
			root++;
		}
	}
	return root;
}


/**
 * ************************************************************************
 * @brief Á½¸ö¸´ÊýÏà³Ë
 * 
 * @param[in] a  ¸´Êý1
 * @param[in] b  ¸´Êý2
 * @note ³Ë»ýµÄÊµ²¿ÎªÁ½¸ö¸´ÊýÊµ²¿µÄ³Ë»ý¼õÈ¥Ðé²¿µÄ³Ë»ý
 * 		 ³Ë»ýµÄÐé²¿ÎªÁ½¸ö¸´ÊýÊµ²¿µÄ³Ë»ý¼ÓÉÏÐé²¿µÄ³Ë»ý
 * @return 
 * ************************************************************************
 */
struct compx EE(struct compx a,struct compx b)
{
	struct compx c;
	c.real=a.real*b.real-a.imag*b.imag;
	c.imag=a.real*b.imag+a.imag*b.real;
	return(c);
}


/**
 * ************************************************************************
 * @brief ¶ÔÊäÈëµÄ¸´Êý×é½øÐÐ¿ìËÙ¸µÀïÒ¶±ä»»£¨FFT£©
 * 
 * @param[in] xin  ¸´Êý×é
 * 
 * ************************************************************************
 */
void FFT(struct compx *xin)
{
	int f,m,nv2,nm1,i,k,l,j=0;
	struct compx u,w,t;

	nv2=FFT_N/2;                  //±äÖ·ÔËËã£¬¼´°Ñ×ÔÈ»Ë³Ðò±ä³Éµ¹Î»Ðò£¬²ÉÓÃÀ×µÂËã·¨
	nm1=FFT_N-1;  
	for(i=0;i<nm1;i++)        
	{
		if(i<j)                    //Èç¹ûi<j,¼´½øÐÐ±äÖ·
		{
			t=xin[j];           
			xin[j]=xin[i];
			xin[i]=t;
		}
		k=nv2;                    //ÇójµÄÏÂÒ»¸öµ¹Î»Ðò
		
		while(k<=j)               //Èç¹ûk<=j,±íÊ¾jµÄ×î¸ßÎ»Îª1   
		{           
			j=j-k;                 //°Ñ×î¸ßÎ»±ä³É0
			k=k/2;                 //k/2£¬±È½Ï´Î¸ßÎ»£¬ÒÀ´ÎÀàÍÆ£¬Öð¸ö±È½Ï£¬Ö±µ½Ä³¸öÎ»Îª0
		}
		
		j=j+k;                   //°Ñ0¸ÄÎª1
	}

	{  //FFTÔËËãºË£¬Ê¹ÓÃµûÐÎÔËËãÍê³ÉFFTÔËËã
		int le,lei,ip;                           
		f=FFT_N;
		for(l=1;(f=f/2)!=1;l++)                  //¼ÆËãlµÄÖµ£¬¼´¼ÆËãµûÐÎ¼¶Êý
			;
		for(m=1;m<=l;m++)                           // ¿ØÖÆµûÐÎ½á¼¶Êý
		{                                           //m±íÊ¾µÚm¼¶µûÐÎ£¬lÎªµûÐÎ¼¶×ÜÊýl=log£¨2£©N
			le=2<<(m-1);                            //leµûÐÎ½á¾àÀë£¬¼´µÚm¼¶µûÐÎµÄµûÐÎ½áÏà¾àleµã
			lei=le/2;                               //Í¬Ò»µûÐÎ½áÖÐ²Î¼ÓÔËËãµÄÁ½µãµÄ¾àÀë
			u.real=1.0;                             //uÎªµûÐÎ½áÔËËãÏµÊý£¬³õÊ¼ÖµÎª1
			u.imag=0.0;
			w.real=XCos(PI/lei);                     //wÎªÏµÊýÉÌ£¬¼´µ±Ç°ÏµÊýÓëÇ°Ò»¸öÏµÊýµÄÉÌ
			w.imag=-XSin(PI/lei);
			for(j=0;j<=lei-1;j++)                   //¿ØÖÆ¼ÆËã²»Í¬ÖÖµûÐÎ½á£¬¼´¼ÆËãÏµÊý²»Í¬µÄµûÐÎ½á
			{
				for(i=j;i<=FFT_N-1;i=i+le)            //¿ØÖÆÍ¬Ò»µûÐÎ½áÔËËã£¬¼´¼ÆËãÏµÊýÏàÍ¬µûÐÎ½á
				{
					ip=i+lei;                           //i£¬ip·Ö±ð±íÊ¾²Î¼ÓµûÐÎÔËËãµÄÁ½¸ö½Úµã
					t=EE(xin[ip],u);                    //µûÐÎÔËËã£¬Ïê¼û¹«Ê½
					xin[ip].real=xin[i].real-t.real;
					xin[ip].imag=xin[i].imag-t.imag;
					xin[i].real=xin[i].real+t.real;
					xin[i].imag=xin[i].imag+t.imag;
				}
				u=EE(u,w);                           //¸Ä±äÏµÊý£¬½øÐÐÏÂÒ»¸öµûÐÎÔËËã
			}
		}
	}
}


/**
 * ************************************************************************
 * @brief ÕÒµ½¾ßÓÐ×î´óÊµ²¿µÄÔªËØµÄË÷Òý
 * 
 * @param[in] data  Comment
 * @param[in] count  Comment
 * 
 * @return 
 * ************************************************************************
 */
int find_max_num_index(struct compx *data,int count)
{
	int i=START_INDEX;
	int max_num_index = i;
	float temp = data[i].real;
	for(i=START_INDEX;i<count;i++)
	{
		if(temp < data[i].real)
		{
			temp = data[i].real;
			max_num_index = i;
		}
	}
	return max_num_index; 
}


/**
 * ************************************************************************
 * @brief Ö±Á÷ÂË²¨Æ÷£¬È¥³ýÐÅºÅÖÐµÄÖ±Á÷³É·Ö£¬²¢Êä³öÂË²¨ºóµÄ½á¹û
 * 
 * @param[in] input  ÊäÈëÐÅºÅ
 * @param[in] df  Comment
 * 
 * @return 
 * ************************************************************************
 */
int dc_filter(int input,DC_FilterData * df) 
{

	float new_w  = input + df->w * df->a;
	int16_t result = 5*(new_w - df->w);
	df->w = new_w;
	
	return result;
}


/**
 * ************************************************************************
 * @brief ¶ÔÊäÈëÐÅºÅ½øÐÐ´ø¿íÏÞÖÆÂË²¨£¬ÏÞÖÆÊäÈëÐÅºÅµÄÆµÂÊ·¶Î§£¬²¢Êä³öÂË²¨ºóµÄ½á¹û
 * 
 * @param[in] input  ÊäÈëÐÅºÅ
 * @param[in] bw  Comment
 * 
 * @return 
 * ************************************************************************
 */
int bw_filter(int input,BW_FilterData * bw) 
{
    bw->v0 = bw->v1;
    
    bw->v1 = (1.241106190967544882e-2*input)+(0.97517787618064910582 * bw->v0);
    return bw->v0 + bw->v1;
}

algorithm.h

/**
 * ************************************************************************
 * 
 * @file algorithm.h
 * @author zxr
 * @brief 
 * 
 * ************************************************************************
 * @copyright Copyright (c) 2024 zxr 
 * ************************************************************************
 */
#ifndef __ALGORITHM_H
#define __ALGORITHM_H

#define FFT_N 			1024    //¶¨Òå¸µÀïÒ¶±ä»»µÄµãÊý
#define START_INDEX 	8  		//µÍÆµ¹ýÂËãÐÖµ

struct compx     	//¶¨ÒåÒ»¸ö¸´Êý½á¹¹
{
	float real;
	float imag;
};  

typedef struct		//¶¨ÒåÒ»¸öÖ±Á÷ÂË²¨Æ÷½á¹¹Ìå
{
	float w;
	int init;
	float a;
}DC_FilterData;		//ÓÃÓÚ´æ´¢Ö±Á÷ÂË²¨Æ÷µÄ²ÎÊý


typedef struct		//¶¨ÒåÒ»¸ö´ø¿íÂË²¨Æ÷½á¹¹Ìå
{
	double v0;
	double v1;
}BW_FilterData;		//ÓÃÓÚ´æ´¢´ø¿íÂË²¨Æ÷µÄ²ÎÊý



double my_floor(double x);

double my_fmod(double x, double y);

double XSin( double x );

double XCos( double x );

int qsqrt(int a);


struct compx EE(struct compx a,struct compx b);

void FFT(struct compx *xin);

int find_max_num_index(struct compx *data,int count);
int dc_filter(int input,DC_FilterData * df);
int bw_filter(int input,BW_FilterData * bw);


#endif

blood.c

/**
 * ************************************************************************
 * 
 * @file blood.c
 * @author zxr
 * @brief 
 * 
 * ************************************************************************
 * @copyright Copyright (c) 2024 zxr 
 * ************************************************************************
 */
#include "blood.h"
#include "usart.h"
#include "string.h"
#include "stdint.h"
#include "stdio.h"

int heart;		//¶¨ÒåÐÄÂÊ
double SpO2;		//¶¨ÒåÑªÑõ±¥ºÍ¶È

//µ÷ÓÃÍâ²¿±äÁ¿
extern uint16_t fifo_red;		//¶¨ÒåFIFOÖÐµÄºì¹âÊý¾Ý
extern uint16_t fifo_ir;		//¶¨ÒåFIFOÖÐµÄºìÍâ¹âÊý¾Ý

uint16_t g_fft_index = 0;         	 	//fftÊäÈëÊä³öÏÂ±ê
struct compx s1[FFT_N+16];           	//FFTÊäÈëºÍÊä³ö£º´ÓS[1]¿ªÊ¼´æ·Å£¬¸ù¾Ý´óÐ¡×Ô¼º¶¨Òå
struct compx s2[FFT_N+16];           	//FFTÊäÈëºÍÊä³ö£º´ÓS[1]¿ªÊ¼´æ·Å£¬¸ù¾Ý´óÐ¡×Ô¼º¶¨Òå



#define CORRECTED_VALUE			47   			//±ê¶¨ÑªÒºÑõÆøº¬Á¿

/**
 * ************************************************************************
 * @brief ¸üÐÂÑªÑõÊý¾Ý
 * @note ´Ó MAX30102 µÄ FIFO ÖÐ¶ÁÈ¡ºì¹âºÍºìÍâÊý¾Ý£¬²¢½«ËüÃÇ´æ´¢µ½Á½¸ö¸´ÊýÊý×és1ºÍs2ÖÐ,
 * 		 ÕâÐ©Êý¾ÝËæºó¿ÉÒÔÓÃÓÚ½øÐÐ¸µÀïÒ¶±ä»»µÈºóÐø´¦Àí
 * 
 * ************************************************************************
 */
void blood_data_update(void)  //ÊÖÖ¸·Åµ½ÉÏÃæ
{
	//±êÖ¾Î»±»Ê¹ÄÜÊ± ¶ÁÈ¡FIFO
	g_fft_index=0;
	while(g_fft_index < FFT_N)
	{
		while(max30102_read_reg(REG_INTR_STATUS_1)&0x40 )
		{
			//¶ÁÈ¡FIFO
			max30102_read_fifo();  //read from MAX30102 FIFO2
			//½«Êý¾ÝÐ´ÈëfftÊäÈë²¢Çå³ýÊä³ö
			if(g_fft_index < FFT_N)
			{
				//½«Êý¾ÝÐ´ÈëfftÊäÈë²¢Çå³ýÊä³ö
				s1[g_fft_index].real = fifo_red;
				s1[g_fft_index].imag= 0;
				s2[g_fft_index].real = fifo_ir;
				s2[g_fft_index].imag= 0;
				g_fft_index++;
//				g_fft_index+=100;
			}
		}
	}
}


/**
 * ************************************************************************
 * @brief ÑªÒºÐÅÏ¢×ª»»
 * 
 * 
 * ************************************************************************
 */
void blood_data_translate(void)
{	
	double n_denom;
	uint16_t i;

	//Ö±Á÷ÂË²¨
	float dc_red =0; 
	float dc_ir =0;
	float ac_red =0; 
	float ac_ir =0;
	
	for (i=0 ; i<FFT_N ; i++ ) 
	{
		dc_red += s1[i].real ;
		dc_ir +=  s2[i].real ;
	}
		dc_red =dc_red/FFT_N ;
		dc_ir =dc_ir/FFT_N ;
	for (i=0 ; i<FFT_N ; i++ )  
	{
		s1[i].real =  s1[i].real - dc_red ; 
		s2[i].real =  s2[i].real - dc_ir ; 
	}

	//ÒÆ¶¯Æ½¾ùÂË²¨
	for(i = 1;i < FFT_N-1;i++) 
	{
		n_denom= ( s1[i-1].real + 2*s1[i].real + s1[i+1].real);
		s1[i].real=  n_denom/4.00; 
		
		n_denom= ( s2[i-1].real + 2*s2[i].real + s2[i+1].real);
		s2[i].real=  n_denom/4.00; 			
	}

	//°ËµãÆ½¾ùÂË²¨
	for(i = 0;i < FFT_N-8;i++) 
	{
		n_denom= ( s1[i].real+s1[i+1].real+ s1[i+2].real+ s1[i+3].real+ s1[i+4].real+ s1[i+5].real+ s1[i+6].real+ s1[i+7].real);
		s1[i].real=  n_denom/8.00; 
		
		n_denom= ( s2[i].real+s2[i+1].real+ s2[i+2].real+ s2[i+3].real+ s2[i+4].real+ s2[i+5].real+ s2[i+6].real+ s2[i+7].real);
		s2[i].real=  n_denom/8.00; 
		
	}

	//¿ªÊ¼±ä»»ÏÔÊ¾	
	g_fft_index = 0;	
	//¿ìËÙ¸µÀïÒ¶±ä»»
	FFT(s1);
	FFT(s2);
	
	for(i = 0;i < FFT_N;i++) 
	{
		s1[i].real=sqrtf(s1[i].real*s1[i].real+s1[i].imag*s1[i].imag);
		s1[i].real=sqrtf(s2[i].real*s2[i].real+s2[i].imag*s2[i].imag);
	}
	//¼ÆËã½»Á÷·ÖÁ¿
	for (i=1 ; i<FFT_N ; i++ ) 
	{
		ac_red += s1[i].real ;
		ac_ir +=  s2[i].real ;
	}
	
	for(i = 0;i < 50;i++) 
	{
		if(s1[i].real<=10)
			break;
	}
	
	//¶ÁÈ¡·åÖµµãµÄºá×ø±ê  ½á¹ûµÄÎïÀíÒâÒåÎª 
	int s1_max_index = find_max_num_index(s1, 60);
	int s2_max_index = find_max_num_index(s2, 60);

	//¼ì²éHbO2ºÍHbµÄ±ä»¯ÆµÂÊÊÇ·ñÒ»ÖÂ
	if(i>=45)
	{
		//ÐÄÂÊ¼ÆËã
		uint16_t Heart_Rate = 60.00 * SAMPLES_PER_SECOND * s1_max_index / FFT_N;
		heart = Heart_Rate;
		
		//ÑªÑõº¬Á¿¼ÆËã
		double R = (ac_ir*dc_red)/(ac_red*dc_ir);
		double sp02_num = -45.060*R*R+ 30.354 *R + 94.845;
		SpO2 = sp02_num;
		
		//×´Ì¬Õý³£
	}
	else //Êý¾Ý·¢ÉúÒì³£
	{
		heart = 0;
		SpO2 = 0;
	}
	//½áÊø±ä»»ÏÔÊ¾
}

/**
 * ************************************************************************
 * @brief ÐÄÂÊÑªÑõÑ»·º¯Êý
 * 
 * 
 * ************************************************************************
 */
void blood_Loop(void)
{
	//ÑªÒºÐÅÏ¢»ñÈ¡
	blood_data_update();
	//ÑªÒºÐÅÏ¢×ª»»
	blood_data_translate();

	SpO2 = (SpO2 > 99.99) ? 99.99:SpO2;  



  printf("ÐÄÂÊ%3d/min; ÑªÑõ%2d%%\n", heart, (int)SpO2);
}

blood.h

/**
 * ************************************************************************
 * 
 * @file blood.h
 * @author zxr
 * @brief 
 * 
 * ************************************************************************
 * @copyright Copyright (c) 2024 zxr 
 * ************************************************************************
 */
#ifndef _BLOOD_H
#define _BLOOD_H



#include "MAX30102.h"
#include "algorithm.h"
#include "math.h"


void blood_data_translate(void);
void blood_data_update(void);
void blood_Loop(void);

#endif

本文章已经生成可运行项目