你的浏览器版本过低,可能导致网站不能正常访问!
为了你能正常使用网站功能,请使用这些浏览器。

基于STM32利用TIMER测量信号周期及占空比的精度经验分享

[复制链接]
攻城狮Melo 发布时间:2024-2-29 14:30
有人使用STM32G4系列的通用型TIMER基于捕获功能对外来信号进行周期及占空比的测量。他用TIM3产生频率、占空比可调的PWM输出做为被测信号。TIM4用来进行频率测量,工作在复位从模式,被测信号接到其通道2的输入脚。然后经内部边沿检测和内部滤波电路后,兵分两路分别连接到IC1和IC2。显然IC2使用直接输入模式,IC1使用间接输入模式,如下图功能框图示意。
5 _3 k; a/ J3 q1 l5 @( v+ C
0 T& }: g. h' T- g6 [
微信图片_20240229143004.png
: R' r& t4 m& f5 |: G" g  b
3 d# H; q+ d. B6 l" ~9 _
但他发现被测信号频率较高时,误差就明显变大了。当然,他也知道,任何测量肯定是有极限的。他现在就是想知道,能否基于现有方案将测量极限拉高点。比方说,他现在测量20KHz信号时就明显误差过大,导致测算结果难以采用。是否可以将可靠的测量结果提升到25KHz或更高呢。: o. j, O( d# R

+ w) S# @! S% U6 }/ f4 R6 j  [
微信图片_20240229143000.png
: S% ]& z9 p$ p
" }' w- l/ c7 {9 \6 ]4 K, c
我们不妨一起看看这个问题。依然保持相同的测量方案,TIM4工作在复位从模式,TIM3输出的被测信号连接到TIM4的CH2,然后兵分两路连接都TIM4的IC1与IC2。* B% K& m' f2 J( S* ^

! S5 v% T  s' M  C% R6 H: s. x7 [5 r显然,我们先要确定测量的计时起点。可以采用IC2的上沿捕获事件作为测量起点,即进入测试状态。也可以基于上沿触发信号产生定时器复位导致的更新事件加触发事件作为测量起点,我在下面就是使用后者来进行测量并组织相应代码。
' l/ o( n$ ?# e1 y4 j; c
/ r$ F7 I$ g2 |2 @
IC1的下降沿触发捕获,捕获到的计数器值存放于变量Value_1stCap【结合上图来看】。
* o3 g/ G, I+ R1 z% E; @' P# k3 ]2 d
IC2的上升沿作为TIM4的触发复位信号,且基于上升沿事件进行捕获, 捕获值存于Value_2ndCap。
3 W% \) [0 C" K# o
; f0 E: T$ V/ G  j% n0 u从发生复位事件到发生第2次捕获期间,对TIM4的更新事件次数进行统计,总的更新事件次数计为Total_Num_OvEvent。( t, d6 P$ j# p$ B

6 p% U# r7 v) I2 W* a% v从发生复位事件到发生第1次捕获期间,对TIM4的更新事件次数进行统计后,计为Front_Num_OvEvent。4 {- ]! V+ w7 B1 F

* \5 u1 U; i, |5 {" V现在使用STM32CubeMx进行配置。重点关注TIM4的配置。TIM4的时基和捕获配置参数如下:
- o0 x- Z/ o7 N, x
& g/ f& \) n1 ~6 [3 _ 微信图片_20240229142957.png : d( Y- m: q9 Z: o

( N; F/ U# w' W) \! X
结合上面配置我们不难看出,TIM4的溢出周期为20ms,TI2FP2作为TIM4的复位触发信号,上沿触发,同时IC2针对输入信号的上沿进行捕获。IC1针对输入信号的下沿进行捕获。- Q4 j& d5 v6 M# q' K
0 t% I) F3 d5 L0 f9 e7 O
至于TIM3的配置没啥特别的,就是产生PWM输出信号,在代码里将调整其PSC分频系数和ARR参数以及CCR参数,以改变其频率和占空比。选择其通道1做PWM输出。下面测试中,TIM3输出的PWM波形的占空比固定为40%。# l4 p; G" q* W$ O

8 |2 r+ W8 J. V! i% ]7 E 微信图片_20240229142953.png
' W! ~8 R0 V3 k  a" }+ u' H' c  D' \6 U7 y" O7 h
另外,STM32G4的系统主频配置为170MHz,使用HSE 时钟源。开启TIM4的中断响应使能。完成配置后创建工程。添加必要的用户初始及启动代码。
  P! U; i  K% y6 W% u. x- L- D
+ f% a  R! J3 e7 Q. q
  1. __IO uint32_t            Vaule_2ndCap = 0u;
    4 q# T  v! _3 K( D7 @) k# j
  2. __IO uint32_t            Vaule_1stCap = 0u;& U% o  f2 J, \1 C
  3. % x# b: b9 s, o# q! f- T7 R
  4. 1 @1 t6 F$ |- D9 f( I) P# X9 D/ A
  5. __IO    float          Signal_Cycle= 0.0f;
    * H$ X0 s5 x6 V3 K9 G: b# P6 _& Y
  6. __IO    float          Signal_Duty = 0.0f;% X" b& u$ R, e
  7. __IO    float          Signal_Freq = 0.0f;
    1 R/ I9 k8 [" P4 Q. D2 Y# ?

  8. 4 V# j7 p# S: T7 N2 b5 Y# ?. J
  9. 4 u/ b. G2 s4 [& }/ d
  10. __IO uint32_t Total_Num_OvEvent=0u;5 S6 R; I4 `, s9 x  U
  11. __IO uint32_t Front_Num_OvEvent=0u;
    ' R1 M0 L* L- l) W- q/ U

  12. % ?& k7 ?0 M! J6 L% ]: \* J
  13. , X- B/ C, {4 R  d- G3 W
  14. __IO uint32_t  Num_OvEvent =0u ;
    7 A/ o- v4 ?6 z* F" U
  15. __IO uint32_t  Measure_State = 0u;7 [; ]- R! Y0 z' Z/ y$ q: v
  16. 5 k" z: ?2 i! `" o: T" ^
  17. & U/ x7 u' S. {5 |# W( C: ^
  18. __IO uint32_t  Clk_Internal;1 y% N; ~3 G1 K) z) g9 T

  19. ! ^0 |) W) G) x( b+ w9 }2 s
  20. $ J9 c7 }3 B# j. t

  21. 4 p0 x+ t6 z& L% U# M

  22. 0 h3 g6 [5 H, a" ]8 [9 S4 Q; D
  23.     __HAL_TIM_CLEAR_IT(&htim4, TIM_IT_UPDATE);5 }  V$ S4 M1 R" k  T, C
  24.     __HAL_TIM_ENABLE_IT(&htim4, TIM_IT_UPDATE);  # _" P0 A. Y( d2 C# x

  25. . d( ]8 V! m9 p- Z

  26. 6 V0 \! h: ?8 Z/ f& e; D
  27.     HAL_TIM_IC_Start_IT(&htim4, TIM_CHANNEL_2); //enable IC1 interrupt of TIM4
    : j7 D0 V/ H8 @0 @) M/ n( K
  28.     HAL_TIM_IC_Start_IT(&htim4, TIM_CHANNEL_1); //enable IC2 interrupt of TIM46 X- g$ x; ~9 {4 T8 q
  29. % ?. b( L' Q* S( E

  30. " O& X6 w6 I( B! k& o
  31.   __HAL_TIM_CLEAR_IT(&htim4, TIM_IT_UPDATE);6 K6 ~- v4 Q, }9 M6 ~: g  Z
  32.   __HAL_TIM_CLEAR_IT(&htim4, TIM_IT_TRIGGER);//enable Update interrupt of TIM4
    # d. z- j) a1 [4 x  q2 s- |
  33.    
    - ^# m. w* u4 Z6 h6 z' f3 E3 T
  34.   Clk_Internal = HAL_RCC_GetHCLKFreq();  //170MHz for G4 Series9 Y0 ^! D. ]2 ?, S  A$ t
  35. ; E1 L, ]/ D; }; R
  36.   b. Y: ~. `- _  w2 J
  37.   Measure_State = 0x00;  //initial state of Measuring
    " r8 e4 U+ n& B5 ^) [! V# S1 Z  i; k
  38.   Num_OvEvent = 0x00;    //initial value of update count of TIM4
    - F. V6 C! r+ g2 [
  39. 9 J3 N( B8 _; k* f+ {
  40. 9 P' R) l7 F( Z+ I  j/ m

  41. 9 Y5 k3 u$ F$ n9 k; g

  42. 7 y* K8 ]% @3 T7 r# h" \; d/ a
  43.   HAL_TIM_PWM_Start(&htim3, TIM_CHANNEL_1);
复制代码

" H- k( r1 Z9 h  j$ ^, m5 Q6 A6 RTIM4事件的所有中断共用同一个中断矢量入口。为了尽可能地测得较高的信号频率,中断处理代码尽量要简洁些,避免过多、过深的函数调用。6 v" P+ }7 j$ X6 |
0 ]0 u8 {- g+ n' C  a# v4 {, R4 J
微信图片_20240229142949.png 0 i) G) x! ^: b+ \

; C) B& ^/ `% _9 ?7 I, |* L有关TIM4的捕获中断及更新中断的处理代码如下,包括IC1、IC2的捕获中断处理和更新中断处理代码,以及相关计算处理代码。
' i& }" G: o+ b7 ~' ^3 E
; Q' i  z! U% o2 }
  1. - m% ~4 a1 t5 p1 b/ e
  2. void HAL_TIM4_PWM_Measure_IRQ(TIM_HandleTypeDef *htim)
    ! F- G( u' K$ Q# [9 M) p1 ?
  3. {
    3 s% H$ v6 u! B/ y9 a3 |6 ~
  4. ' }8 O( v) U7 D. H
  5. if(__HAL_TIM_GET_FLAG(htim, TIM_FLAG_UPDATE) != RESET)
    0 ?& N0 G0 E% d4 f, r$ I7 e
  6.   {) [, F+ z7 U3 k$ i" T/ z5 V* A
  7.       __HAL_TIM_CLEAR_IT(htim, TIM_IT_UPDATE);8 d8 N0 @+ V( T6 p3 U& U4 Z

  8. : m# [. i: c, u
  9. if(__HAL_TIM_GET_FLAG(htim, TIM_FLAG_TRIGGER) != RESET), F' T: }* P( O1 P7 G0 U
  10.       {
    ) J7 P2 l8 e5 J4 ]" C# I
  11.           __HAL_TIM_CLEAR_IT(htim, TIM_IT_TRIGGER);1 G+ P" S7 h3 q

  12. % n! l4 q3 V% d' I
  13. if (Measure_State == 0). F( s* l7 {$ B% P

  14. $ q1 b9 F+ S! k# d! R6 _
  15.           {
    9 w6 C% U% H: J, v7 ^8 }
  16.             Num_OvEvent = 0;      // prepare count overflow events
    3 e3 M. _% ]! E1 E. w" H. g6 h
  17.             Measure_State = 0x01;   //start measurement
    9 F) q' w4 D! a: d  m7 l
  18.           }  E9 X: w' v# L/ u7 n
  19.          }
    0 l8 V" ]& z* v. U9 h
  20. $ k' h8 R  q% S1 P/ g2 g
  21. else# w2 S5 \( t& A  N' }
  22.     {
    , H7 Q+ K$ K( M( e; u' Q
  23. if(Measure_State != 0)  ' ?4 o( Y; v1 t
  24.       {
    - {: q/ [8 K; ~/ Z
  25.        Num_OvEvent++;/ m5 w3 @- c5 D% T
  26.       }  
    % U! w# m; o1 x% i! s1 U
  27.     }
    " K& o! _/ a3 b7 G: \3 U) O
  28.   }! w2 `  Z0 o0 R4 K; W1 F$ z

  29. $ L# d* G. Q' f6 I% Y0 ~
  30. /* Capture compare 1 event */- T' x/ T+ @7 s( W' \8 {2 y0 t: e; ^
  31. if(__HAL_TIM_GET_FLAG(htim, TIM_FLAG_CC1) != RESET)' P' \/ y" K) G4 D
  32.   {
    * ~2 D8 b: c" h# Q

  33. $ P5 X3 M9 W+ X
  34.         __HAL_TIM_CLEAR_IT(htim, TIM_IT_CC1); 1 v' r* q& G2 b2 ~( D8 i& o4 h
  35. / X. A  I9 G# c9 A. D

  36. 9 j1 C1 @/ I' H6 L; b- c2 I7 I) l" P. b
  37. if (Measure_State == 0x01)& a* H1 N: K$ d: p# i) |
  38.           {    . r; \0 u. L' Y
  39.             Front_Num_OvEvent = Num_OvEvent;
    . c) T! w2 d0 z% h: w) j
  40.             Measure_State =0x02;
    " j0 A3 t, T' {0 \8 p6 R5 f. v! |. m
  41.            }  / b" q3 y" _& o- @/ ]* c8 U; A

  42. 4 J+ C& M# ~5 [; [; [- W$ `
  43.   }
    2 K! e- j& E3 J1 b, O

  44. % T7 u- {. }3 [* R9 i: Z
  45. /* Capture compare 2 event */
    ! O+ P$ a9 E5 c$ v3 |$ l' g5 {
  46. if(__HAL_TIM_GET_FLAG(htim, TIM_FLAG_CC2) != RESET)
    + g! O) |' A# ]' w0 R+ N/ b0 p
  47. {8 F1 d1 F* K  o. X" l
  48.       __HAL_TIM_CLEAR_IT(htim, TIM_IT_CC2);6 q$ g8 y1 y  v
  49. 8 g$ _& A! g3 O
  50.   G) J6 C& ]) k1 _' [
  51. if(Measure_State == 0x02)
    : m8 a9 K1 P- X- a' A; Z! C; _; b
  52.         {6 B/ K& M5 b& r% j, z  G4 k- _
  53.           Measure_State =0x03;
    " m2 F8 d( N6 l- H
  54.           Total_Num_OvEvent = Num_OvEvent;3 p5 d0 }. W% G$ x2 j7 ^
  55.            HAL_TIM4_IC_CaptureCallback(htim);      //go to calculate pulse width and duty/ g, e4 {2 I9 w$ h+ F( Q2 g

  56. * K  z: ~, S- o/ g
  57.         }      
    " }) m) R" o8 l% |# S0 K
  58. 2 l. m* h+ A  w4 k+ W/ m6 {
  59. }
    7 h# v  d: j  M- @) b5 n5 b

  60. & G* v7 w$ ]; j* X( e5 k9 d: Q
  61. 8 T  t- l5 f) t- a/ K; z
  62. }
    0 A. H/ \8 t0 T1 M. y# h3 n+ ^) b
  63. # A/ r; f; }$ y5 V  U/ n8 p" V
  64. void HAL_TIM4_IC_CaptureCallback(TIM_HandleTypeDef *htim)9 F: s/ O8 l5 T0 _0 Y5 z
  65. {
    1 g. \) d' v' w* G% M* C8 ?  U

  66. : Y: F* p' V( z0 Q. i
  67. /* Get the Input Capture value */
    8 j6 q4 d, K! _2 F+ w) c
  68.     Vaule_1stCap = TIM4->CCR1;//HAL_TIM_ReadCapturedValue(htim, TIM_CHANNEL_1);
    ( B& [  J) ~, g, Y0 e! ?$ K
  69.     Vaule_2ndCap = TIM4->CCR2;//HAL_TIM_ReadCapturedValue(htim, TIM_CHANNEL_2);* Y$ ~; H3 t7 C, u3 Q
  70. 7 k6 l- p  ^& F" H) N
  71. - z0 `0 V7 T. f) F7 p4 L
  72.     Signal_Duty =(float)((Vaule_1stCap+(Front_Num_OvEvent*(TIM4_PERIOD+1)))*100u)/ \
    ) W8 D# N6 R0 e' z
  73.       (float)(Vaule_2ndCap+(Total_Num_OvEvent*(TIM4_PERIOD+1)));   
    * g  H# X1 d1 B; t! U) F
  74. & A( l( l6 k: u! o! I" C
  75. //  Clk_Internal = HAL_RCC_GetHCLKFreq();
    & \, ?# U/ f& w
  76. 0 I2 t8 P& Q: O( [9 z3 F  ?
  77.     Signal_Freq =(float)(Clk_Internal/(TIM4_PSC+1))/(float)((Vaule_2ndCap+(Total_Num_OvEvent*(TIM4_PERIOD+1))));
    : b" m  x8 ^$ ~$ A$ Z: u
  78.   Signal_Cycle = (float) (1.000f/Signal_Freq) ;
    ( @" D; f4 Q/ s0 Y- @2 C* c

  79. ! S9 O5 g' u$ H6 ^' d- q

  80. 5 i% d: ~* e) {* j0 [- m8 X" X
  81.       Measure_State = 0x00;
    : g( y4 X; z  C3 D
  82.       Num_OvEvent = 0x00;9 z3 |2 z# ~$ R$ H# a
  83.     __HAL_TIM_CLEAR_IT(htim, TIM_IT_CC1); 8 n% E& u1 M' [
  84.     __HAL_TIM_CLEAR_IT(htim, TIM_IT_CC2);
    0 T' @. g' C5 N3 }* m9 z  q
  85.     __HAL_TIM_CLEAR_IT(htim, TIM_IT_UPDATE);% C& ^  F4 g3 i9 t$ R
  86.     __HAL_TIM_CLEAR_IT(htim, TIM_IT_TRIGGER);
    8 t" I" T/ _2 @3 @) A4 u" H
  87.    
    5 b* z4 v- z# W7 q
  88. }
复制代码
, x( B, I, i" F, @5 _
现在基于上面的配置及用户代码进行测量验证。我们发现当TIM3输出的待测信号频率达到25KHz时,测得的频率及占空比就发生了明显的偏差。见下图标注问号的第4栏信息。
5 X( Q* O) Y+ a: X. m( a9 {- w. C 微信图片_20240229142947.png
4 u) l+ l7 Z. t" b* i) |9 m! Q* p* u# \
即当被测频率为25KHz,占空比为40%时,测得结果是下面的样子,明显偏差过大。# _9 N' `: o; d% X, T- f4 H
7 p% n" v; r, D0 {8 E; a2 F
微信图片_20240229142943.png " d/ }3 s/ e  O9 N7 w
如果说被测信号频率进一步提升的话,误差会变得更大。那么,这种情形是否有改善机会呢?即在当前的测试方案下,可准确测量的被测信号频率是否可以提高。. ~3 v* X1 _% _3 H* U- G
6 U5 y3 Z4 ^: n4 Q4 R
目前的中断处理代码应该说比较精简,没有什么可以优化的余地了。
) X# R8 M( x" k
5 V% h. v' ]" Z* }8 T聊到这里,有人可能发现了,我前面配置TIM4时,它是用来完成测量任务的,它的分频系数PSC为169。即TIM4的计数器计数时钟为1MHz,其计数分辨率为1us。不难理解,这个PSC系数应该会直接影响TIM4的计数精度,按理会影响到测量结果,尤其被测信号频率较高时。
3 n! R( T, P- k* Z. B4 j* ~; x% F
, }0 E; Z; Y$ n7 r$ j' q0 {
既然这样,我们将TIM4的PSC系数改为0,并适当调整其ARR值再来实施测量,看看结果会怎么样。见下图,TIM3输出的信号频率仍然是占空比为40% 、频率为25KHz的PWM信号。测量结果显著地明显改善,应该说此时结果是可以接受的,毕竟浮点运算也会带来些偏差。
8 L) W$ `& V  E' n: W$ Q
% \: ^; i+ P2 M; y* y+ ?1 \
微信图片_20240229142940.png
) m7 M7 _) p9 C) d* S4 i3 {/ @
/ {$ S! ?3 t. |3 J" z, o我们不妨在保持TIM3的PSC为0的条件下,将输出频率提升到50KHz、100KHz。继续看看测量结果,见下图:
3 q3 d7 d! A, @7 u8 z5 _
& _# K3 N' d  ?
微信图片_20240229142937.png
: j7 r: q' A( I* S/ I
; Q+ H: {) W5 T( P' k
从测试结果来看,当被测信号频率提升到50KHz,测量结果仍然很好,完全可以采用。即使当被测信号频率提升到100KHz时,测量结果虽发生了一些偏差,但此时的偏差相比PSC=169、被测信号为25KHz时的测试结果还要好得多。我把二者单列出来一起比较,见下图:3 z8 Z! ~7 s9 A! o
微信图片_20240229142932.png + i! A$ P& X8 K/ |" ]# H% }7 M

/ x8 Y. L. W4 Z7 ^8 o) u+ V
4 G0 Z6 |* ?3 k5 \0 o/ s6 }4 @
经过上面的讨论和验证,我们知道,在使用TIMER做信号的频率及占空比的测量时,当确定好测试方案后,为了尽可能地提升可以准确测试信号的频率,一方面代码要尽可能精简、优化,另一方面,因测试TIMER的分频系数会影响测试结果的精度及准确性,此时测试TIMER的分频系数要尽可能设置小、或不做分频,这点结合具体应用场景来定。$ j3 m& G8 W, W  q2 k9 \! `4 k
. s, N6 J' q) j; n3 `
关于上面的的话题,其实还可以有很多继续延伸的空间,因时间和篇幅问题,就先聊到这里。有兴趣的话,可以基于上面配置和分享的参考源码做些快速测试验证,也欢迎进一步地讨论交流。
: [; `+ ~8 M( K% `- q& {1 q$ c) j3 c/ c3 I8 P

8 f+ k5 G2 G9 ]9 U转载自; 茶话MCU
" h( v6 ?0 ^7 T" o5 U如有侵权请联系删除8 H0 R! d5 H7 B7 a+ ~# ]
# Y' y" c3 P# B4 Q- m5 |$ v
* E  E& y/ r& Y$ B( p
收藏 评论0 发布时间:2024-2-29 14:30

举报

0个回答
关于意法半导体
我们是谁
投资者关系
意法半导体可持续发展举措
创新和工艺
招聘信息
联系我们
联系ST分支机构
寻找销售人员和分销渠道
社区
媒体中心
活动与培训
隐私策略
隐私策略
Cookies管理
行使您的权利
关注我们
st-img 微信公众号
st-img 手机版