用C语言实战：最小公倍数在嵌入式编程和单片机开发中的一个具体应用案例

用C语言实战：最小公倍数在嵌入式编程中的精准调度应用

在嵌入式系统开发中，时间管理是核心挑战之一。想象一下这样的场景：你的STM32单片机需要同时处理两个周期性任务——一个传感器每3毫秒采集一次数据，另一个显示屏每5毫秒刷新一次内容。如何确保这两个任务能够和谐共处，不会互相干扰？这就是最小公倍数算法大显身手的地方。

不同于教科书上的纯数学示例，我们将深入探讨如何用C语言实现最小公倍数计算，并将其应用于真实的嵌入式调度问题。通过这个案例，你会看到算法不仅仅是抽象的概念，而是解决实际工程问题的利器。本文面向有一定C语言基础的嵌入式初学者，将从实际问题出发，逐步构建解决方案。

1. 理解嵌入式系统中的定时调度需求

在单片机编程中，定时器中断是实现周期性任务的常用手段。假设我们有以下两个任务：

• 任务A：环境传感器数据采集，每3ms执行一次
• 任务B：OLED显示屏刷新，每5ms执行一次

当这两个任务独立运行时，系统需要确定它们何时会同时触发，以便合理分配CPU资源，避免冲突。这正是最小公倍数能够解决的问题——3和5的最小公倍数是15，意味着两个任务每15ms会同时触发一次。

典型的问题场景包括：

• 避免中断服务程序(ISR)执行时间重叠
• 协调多个外设的访问时序
• 实现复合功能的同步点
• 优化CPU利用率

在STM32等ARM Cortex-M系列单片机中，系统滴答定时器(SysTick)通常作为时间基准，配合多个通用定时器(TIM)可以实现复杂的调度需求。理解这些定时器的协调机制，是嵌入式开发者的必备技能。

2. 最小公倍数算法的C语言实现

针对嵌入式环境的特殊要求，我们需要高效、可靠的最小公倍数算法实现。以下是三种经过优化的方法：

2.1 基于递增的暴力搜索法

/** * 暴力搜索法求最小公倍数 * @param a 第一个数 * @param b 第二个数 * @return 最小公倍数 */ uint32_t lcm_brute_force(uint32_t a, uint32_t b) { uint32_t max = (a > b) ? a : b; while(1) { if((max % a == 0) && (max % b == 0)) { return max; } max++; } }

这种方法直接明了，但效率较低，特别是在大数情况下。它的时间复杂度为O(n)，在资源受限的单片机上可能不是最佳选择。

2.2 利用最大公约数的数学方法

更高效的方式是利用数学关系：LCM(a,b) = (a×b)/GCD(a,b)。我们需要先实现欧几里得算法求最大公约数：

/** * 欧几里得算法求最大公约数 * @param a 第一个数 * @param b 第二个数 * @return 最大公约数 */ uint32_t gcd(uint32_t a, uint32_t b) { while(b != 0) { uint32_t temp = b; b = a % b; a = temp; } return a; } /** * 通过GCD计算LCM * @param a 第一个数 * @param b 第二个数 * @return 最小公倍数 */ uint32_t lcm_via_gcd(uint32_t a, uint32_t b) { return (a / gcd(a, b)) * b; // 先除后乘避免溢出 }

这种方法的时间复杂度主要取决于GCD计算，通常是O(log(min(a,b)))，效率显著高于暴力搜索。

2.3 优化后的倍数检查法

结合前两种方法的优点，我们可以实现一个折中方案：

/** * 优化的倍数检查法 * @param a 第一个数 * @param b 第二个数 * @return 最小公倍数 */ uint32_t lcm_optimized(uint32_t a, uint32_t b) { uint32_t multiple = a; while(multiple % b != 0) { multiple += a; // 每次增加a的值，而不是1 } return multiple; }

这种方法比纯暴力搜索更快，但不如GCD方法高效。在实际应用中，应根据具体场景选择最合适的实现。

性能对比表：

3. 在STM32中实现多定时器协调

现在我们将最小公倍数算法应用到实际的STM32定时器调度中。假设使用TIM2和TIM3两个定时器，分别配置为3ms和5ms触发中断。

3.1 硬件定时器配置

首先配置定时器的基础参数：

// TIM2配置为3ms中断 void TIM2_Config(void) { RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM2, ENABLE); TIM_TimeBaseInitTypeDef TIM_TimeBaseStructure; TIM_TimeBaseStructure.TIM_Period = 2999; // 自动重装载值 TIM_TimeBaseStructure.TIM_Prescaler = 71; // 预分频值 TIM_TimeBaseStructure.TIM_ClockDivision = 0; TIM_TimeBaseStructure.TIM_CounterMode = TIM_CounterMode_Up; TIM_TimeBaseInit(TIM2, &TIM_TimeBaseStructure); TIM_ITConfig(TIM2, TIM_IT_Update, ENABLE); TIM_Cmd(TIM2, ENABLE); } // TIM3配置为5ms中断 void TIM3_Config(void) { RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM3, ENABLE); TIM_TimeBaseInitTypeDef TIM_TimeBaseStructure; TIM_TimeBaseStructure.TIM_Period = 4999; // 自动重装载值 TIM_TimeBaseStructure.TIM_Prescaler = 71; // 相同的预分频值 TIM_TimeBaseStructure.TIM_ClockDivision = 0; TIM_TimeBaseStructure.TIM_CounterMode = TIM_CounterMode_Up; TIM_TimeBaseInit(TIM3, &TIM_TimeBaseStructure); TIM_ITConfig(TIM3, TIM_IT_Update, ENABLE); TIM_Cmd(TIM3, ENABLE); }

注意：预分频值71是基于72MHz系统时钟计算的，使计数器每1μs递增一次。

3.2 中断服务程序中的调度逻辑

在中断服务程序中，我们需要实现基于最小公倍数的协同调度：

volatile uint32_t timer2_count = 0; volatile uint32_t timer3_count = 0; volatile uint32_t sync_point = 0; void TIM2_IRQHandler(void) { if(TIM_GetITStatus(TIM2, TIM_IT_Update) != RESET) { TIM_ClearITPendingBit(TIM2, TIM_IT_Update); timer2_count++; // 执行3ms任务 Sensor_Read(); // 检查同步点 if((timer2_count % 5 == 0) && (timer3_count % 3 == 0)) { sync_point = 1; } } } void TIM3_IRQHandler(void) { if(TIM_GetITStatus(TIM3, TIM_IT_Update) != RESET) { TIM_ClearITPendingBit(TIM3, TIM_IT_Update); timer3_count++; // 执行5ms任务 Display_Refresh(); // 检查同步点 if((timer2_count % 5 == 0) && (timer3_count % 3 == 0)) { sync_point = 1; } } }

3.3 主程序中的协同处理

在主循环中，我们可以利用同步点执行需要协调的操作：

int main(void) { SystemInit(); TIM2_Config(); TIM3_Config(); NVIC_Config(); // 配置中断优先级 // 计算理论同步周期 uint32_t lcm_period = lcm_via_gcd(3, 5); // 15ms while(1) { if(sync_point) { sync_point = 0; // 执行需要同步的操作 Process_Synchronized_Tasks(); // 可以在此添加调试信息 Debug_Print("Sync at %lu ms", lcm_period * (timer2_count / 5)); } // 其他后台任务 Background_Task(); } }

4. 进阶优化与错误处理

在实际工程应用中，我们需要考虑更多细节来确保系统的稳定性和可靠性。

4.1 定时器漂移补偿

即使使用精密的晶体振荡器，长时间运行后定时器也可能出现微小漂移。我们可以实现一个简单的补偿机制：

#define MAX_DRIFT_NS 1000 // 允许的最大漂移(纳秒) void Adjust_Timer_Drift(uint32_t expected, uint32_t actual) { int32_t drift = (int32_t)(actual - expected); if(abs(drift) > MAX_DRIFT_NS) { // 调整定时器预分频值或自动重装载值 uint32_t new_period = TIM_GetAutoreload(TIM2) + (drift / 1000); TIM_SetAutoreload(TIM2, new_period); } }

4.2 中断优先级管理

正确处理中断优先级对系统稳定性至关重要：

void NVIC_Config(void) { NVIC_InitTypeDef NVIC_InitStructure; // TIM2中断配置 NVIC_InitStructure.NVIC_IRQChannel = TIM2_IRQn; NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 1; NVIC_InitStructure.NVIC_IRQChannelSubPriority = 1; NVIC_InitStructure.NVIC_IRQChannelCmd = ENABLE; NVIC_Init(&NVIC_InitStructure); // TIM3中断配置 NVIC_InitStructure.NVIC_IRQChannel = TIM3_IRQn; NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 2; NVIC_InitStructure.NVIC_IRQChannelSubPriority = 1; NVIC_Init(&NVIC_InitStructure); }

4.3 资源冲突预防

当多个任务需要访问共享资源时，最小公倍数同步点可以作为天然的互斥点：

void Process_Synchronized_Tasks(void) { // 在同步点安全地访问共享资源 static uint32_t shared_data; // 使用临界区保护 __disable_irq(); // 读写共享数据 shared_data = Sensor_Get_Value(); Display_Set_Data(shared_data); __enable_irq(); // 执行其他需要同步的操作 Log_System_Status(); }

4.4 动态周期调整

某些应用可能需要动态调整任务周期，这时我们需要重新计算最小公倍数：

void Dynamic_Period_Adjust(uint32_t new_period_A, uint32_t new_period_B) { // 计算新的最小公倍数 uint32_t new_lcm = lcm_via_gcd(new_period_A, new_period_B); // 调整定时器配置 TIM_SetAutoreload(TIM2, (new_period_A * 1000) - 1); TIM_SetAutoreload(TIM3, (new_period_B * 1000) - 1); // 重置计数器 timer2_count = 0; timer3_count = 0; sync_point = 0; }