1. 项目概述：基于STM32的Proteus仿真测频计
在电子测量领域，频率测量是最基础也最关键的参数之一。传统硬件测频方案往往需要购置专用仪器或复杂的外围电路，而今天我要分享的这套基于STM32的纯仿真方案，仅用200行代码就实现了1%精度的多波形频率测量。这个项目的独特之处在于完全基于Proteus仿真环境开发，既避免了硬件焊接的麻烦，又能完整验证所有核心功能。
这个测频计支持测量方波、正弦波和锯齿波等多种波形，测量范围覆盖100Hz-100kHz。核心原理是利用STM32的定时器输入捕获功能，通过计算信号边沿时间差来推算频率。整个开发过程涉及定时器配置、中断处理、信号调理等多个关键技术点，接下来我将逐一拆解每个环节的实现细节。
2. 硬件仿真环境搭建
2.1 Proteus工程配置
在Proteus中新建工程时，关键是要选择正确的单片机型号。这里我们使用STM32F103C8，这是最经典的Cortex-M3内核单片机，具有丰富的外设资源。工程创建后需要做以下配置：
- 设置晶振频率为8MHz（后续通过代码倍频到72MHz）
- 添加虚拟示波器用于观察输入信号
- 放置信号发生器作为测试源
- 配置PA0引脚为TIM2_CH1输入通道
注意：Proteus默认的电源电压是5V，而STM32的IO口虽然标称3.3V，但实际可以耐受5V输入。如果使用真实硬件，建议添加电平转换电路保护MCU。
2.2 信号调理电路设计
不同波形需要不同的预处理方式：
- 方波：直接接入PA0引脚即可
- 正弦波：需要添加施密特触发器整形电路
- 锯齿波：建议并联104电容平滑波形
在Proteus中可以使用"Schmitt Trigger"元件实现波形整形，具体参数设置：
- 正向阈值电压：1.6V
- 负向阈值电压：1.4V
- 电源电压：5V
3. 软件实现原理
3.1 定时器配置详解
定时器是测频的核心，配置TIM2需要关注以下参数：
c复制TIM_TimeBaseInitTypeDef TIM_InitStructure;
TIM_InitStructure.TIM_Prescaler = 72-1; // 72MHz/72=1MHz
TIM_InitStructure.TIM_CounterMode = TIM_CounterMode_Up;
TIM_InitStructure.TIM_Period = 0xFFFF; // 16位最大值
TIM_InitStructure.TIM_ClockDivision = 0;
TIM_TimeBaseInit(TIM2, &TIM_InitStructure);
这段配置将定时器时钟降为1MHz（每个计数代表1μs），采用向上计数模式。Prescaler值设为71是因为STM32的预分频器实际值为设置值+1。
3.2 输入捕获配置
输入捕获的配置决定了如何检测信号边沿：
c复制TIM_ICInitTypeDef TIM_ICInitStructure;
TIM_ICInitStructure.TIM_Channel = TIM_Channel_1;
TIM_ICInitStructure.TIM_ICPolarity = TIM_ICPolarity_Rising; // 上升沿触发
TIM_ICInitStructure.TIM_ICSelection = TIM_ICSelection_DirectTI;
TIM_ICInitStructure.TIM_ICPrescaler = TIM_ICPSC_DIV1; // 不分频
TIM_ICInitStructure.TIM_ICFilter = 0x0; // 不滤波
TIM_ICInit(TIM2, &TIM_ICInitStructure);
这里选择直接输入模式（DirectTI），关闭了数字滤波器（Filter=0），因为仿真环境下信号比较干净。实际硬件应用中，建议设置Filter=0x4，可以滤除短于4个时钟周期的干扰。
4. 核心算法实现
4.1 中断服务函数
捕获中断是测量的核心，这里使用全局变量记录时间差：
c复制volatile uint32_t last_capture = 0;
volatile uint32_t period = 0;
void TIM2_IRQHandler(void) {
if(TIM_GetITStatus(TIM2, TIM_IT_CC1) == SET) {
uint16_t current_capture = TIM_GetCapture1(TIM2);
period = current_capture - last_capture;
last_capture = current_capture;
TIM_ClearITPendingBit(TIM2, TIM_IT_CC1);
}
}
volatile关键字确保编译器不会优化掉这些变量，这在嵌入式开发中尤为重要。中断服务函数要尽可能简短，避免影响下一次捕获。
4.2 频率计算与误差处理
基础频率计算公式很简单：
c复制float get_freq(void) {
if(period == 0) return 0;
return 1000000.0f / period; // 1MHz时钟换算
}
但在实际应用中需要考虑以下特殊情况：
- 信号频率超过500kHz时，周期测量模式误差增大
- 输入信号存在抖动或噪声
- 占空比不对称的波形
改进方案是采用多次测量取平均的方法：
c复制#define SAMPLE_COUNT 5
float get_accurate_freq(void) {
uint32_t samples[SAMPLE_COUNT];
for(int i=0; i<SAMPLE_COUNT; i++) {
while(period == 0); // 等待新数据
samples[i] = period;
period = 0;
}
// 去掉最大最小值后取平均
uint32_t sum = 0;
uint32_t min = 0xFFFFFFFF, max = 0;
for(int i=0; i<SAMPLE_COUNT; i++) {
if(samples[i] < min) min = samples[i];
if(samples[i] > max) max = samples[i];
sum += samples[i];
}
sum = sum - min - max;
return 1000000.0f / (sum / (SAMPLE_COUNT-2));
}
5. 精度优化技巧
5.1 时钟源选择
STM32的时钟源选择直接影响测量精度：
- 内部HSI时钟：精度±1%，不适合精密测量
- 外部HSE时钟：精度可达±0.005%，推荐使用
在代码中启用外部时钟并配置PLL：
c复制RCC_DeInit();
RCC_HSEConfig(RCC_HSE_ON);
while(RCC_GetFlagStatus(RCC_FLAG_HSERDY) == RESET);
RCC_PLLConfig(RCC_PLLSource_HSE_Div1, RCC_PLLMul_9); // 8MHz*9=72MHz
RCC_PLLCmd(ENABLE);
while(RCC_GetFlagStatus(RCC_FLAG_PLLRDY) == RESET);
RCC_SYSCLKConfig(RCC_SYSCLKSource_PLLCLK);
5.2 中断优先级管理
为确保定时器捕获的实时性，需要合理设置中断优先级：
c复制NVIC_InitTypeDef NVIC_InitStructure;
NVIC_InitStructure.NVIC_IRQChannel = TIM2_IRQn;
NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 0; // 最高抢占优先级
NVIC_InitStructure.NVIC_IRQChannelSubPriority = 0;
NVIC_InitStructure.NVIC_IRQChannelCmd = ENABLE;
NVIC_Init(&NVIC_InitStructure);
同时，在捕获中断中临时禁用其他高优先级中断：
c复制void TIM2_IRQHandler(void) {
__disable_irq(); // 禁用全局中断
if(TIM_GetITStatus(TIM2, TIM_IT_CC1) == SET) {
// 捕获处理代码
}
__enable_irq(); // 恢复全局中断
}
6. 实测结果分析
6.1 不同波形测试数据
使用Proteus内置信号发生器测试得到以下结果：
| 波形类型 | 设定频率 | 测量频率 | 相对误差 |
|---|---|---|---|
| 方波 | 1kHz | 1001Hz | +0.1% |
| 正弦波 | 10kHz | 9989Hz | -0.11% |
| 锯齿波 | 50Hz | 49.98Hz | -0.04% |
| 方波 | 100kHz | 100.3kHz | +0.3% |
6.2 误差来源分析
主要误差来源包括：
- 定时器时钟精度（外部晶振决定）
- 中断响应延迟（约12个时钟周期）
- 信号边沿抖动
- 软件计算时的浮点舍入
对于100kHz以上的信号，建议改用测周期数法：在固定时间窗口（如1ms）内统计脉冲数量，这种方法在高频时更准确。
7. 项目扩展方向
7.1 添加LCD实时显示
可以连接一个16x2字符LCD来实时显示频率值：
c复制// LCD初始化
void LCD_Init(void) {
// 初始化代码...
}
// 显示频率
void LCD_DisplayFreq(float freq) {
char buf[16];
if(freq >= 1000) {
sprintf(buf, "Freq:%7.3fkHz", freq/1000);
} else {
sprintf(buf, "Freq:%7.3fHz ", freq);
}
LCD_WriteString(buf);
}
7.2 自动量程切换
实现自动量程切换可以扩大测量范围：
c复制void AutoRangeAdjust(float freq) {
if(freq > 50000) {
// 切换到高频测量模式
TIM_ICInitStructure.TIM_ICPolarity = TIM_ICPolarity_BothEdges;
TIM_ICInit(TIM2, &TIM_ICInitStructure);
} else {
// 常规测量模式
TIM_ICInitStructure.TIM_ICPolarity = TIM_ICPolarity_Rising;
TIM_ICInit(TIM2, &TIM_ICInitStructure);
}
}
7.3 添加串口输出功能
通过串口将测量结果发送到电脑：
c复制void USART_Config(void) {
// USART初始化代码...
}
void SendFreqToPC(float freq) {
char msg[32];
sprintf(msg, "Frequency: %.2f Hz\r\n", freq);
USART_SendString(msg);
}
8. 常见问题与解决方案
8.1 测量结果不稳定
可能原因及解决方法：
- 信号幅度不足 → 确保信号峰峰值大于2V
- 未使用施密特触发器 → 添加整形电路
- 中断优先级设置不当 → 调整NVIC优先级
8.2 高频测量误差大
优化方案：
- 改用输入捕获的双边沿模式
- 缩短中断服务函数执行时间
- 使用DMA传输捕获数据
8.3 Proteus仿真异常
常见仿真问题：
- 信号源设置错误 → 检查幅值和偏置
- 单片机未加载hex文件 → 重新编译并加载
- 仿真速度太慢 → 调整仿真步长
9. 关键代码优化技巧
9.1 使用硬件除法器
STM32F103具有硬件除法器，可以加速频率计算：
c复制float get_freq_optimized(void) {
if(period == 0) return 0;
return 1000000.0f / (float)(period);
}
9.2 减少浮点运算
对于只需要整数结果的场合，可以改用定点数运算：
c复制uint32_t get_freq_int(void) {
if(period == 0) return 0;
return 1000000 / period;
}
9.3 循环展开优化
在多次采样求平均时，可以手动展开循环：
c复制float get_accurate_freq_opt(void) {
uint32_t s1, s2, s3, s4, s5;
while(period == 0); s1 = period; period = 0;
while(period == 0); s2 = period; period = 0;
while(period == 0); s3 = period; period = 0;
while(period == 0); s4 = period; period = 0;
while(period == 0); s5 = period; period = 0;
// 排序并去掉最大最小值
// ...
}
10. 项目总结与心得
通过这个项目，我深刻体会到STM32定时器功能的强大之处。仅仅利用基本的输入捕获功能，配合适当的软件算法，就能实现相当精度的频率测量。在开发过程中，有几点特别值得注意：
- 时钟配置是精度基础，务必使用外部晶振
- 中断服务函数要尽可能精简
- 多次测量取平均能显著提高稳定性
- 不同波形需要不同的预处理方式
这个仿真方案虽然省去了硬件调试的麻烦，但所有核心原理和代码都可以直接移植到真实硬件上。我建议有兴趣的开发者可以先用Proteus验证思路，然后再动手制作实体频率计，这样能大大降低开发风险。