🏴☠️STM32Cubemx ADC+TIM+DMA欠采样采集高频信号
本文主要讲解ADC借助欠采样采集高频信号，比如使用100k左右的采样率去采集1M的信号。
所需工具：
- 开发板:STM32F103RCT6
- STM32CubeMX
- IDE: Keil-MDK
相关文章：
- STM32HAL ADC+TIM+DMA采集交流信号
- STM32cubeHAL ADC+TIM+DMA （二）
- stm32cubemx ADC+TIM+DMA超频采样_
- STM32ADC同步采样
😄原理简介
看过本文最一开始的“相关文章”中提到的文章后，会对信号采集有一定的概念。文章中使用的都是过采样，意思是采样率为待测信号频率的两倍及以上。这么说有些空，举例来说。如果想采集一个1K的正弦信号，一个周期采集4个点，那么采样率为4K。过程如下图：
题外话，没学过信号与系统的小朋友，可能会疑惑，这边采集4个点有啥用？根本看不出来是正弦，还原出来和三角波一样嘛！这就涉及到奈奎斯特采样定理的原理了，要好好学习《信号与系统》和《数字信号处理》，就能明白为什么这个正弦信号，一个周期采集2个点以上就可以复原出来。怎么复原呢？需要用到增采样的知识，我写过一篇matlab仿真增采样：Matlab增采样仿真
上图展示的采样率为4k，每隔0.25ms采集一个点。
增采样对采样率要求较高，如果要采集的信号是1M，一个周期4个点的话，那么我的采样率是不是就要4M了？有没有办法用低采样率，采集高频信号呢？
一个周期四个点嘛，我能否每个周期只采集一个点，然后拼在一起呢？每隔1.25ms采集一个点，这样采集的四个点拼接在一起，不也是一样的效果。
如图中所示，每个周期采集一个点，把所有点拼接在一起，就实现了每个周期采集4个点的效果。
采样率
=
1
1.25
∗
1
0
−
3
=
800
h
z
采样率=\frac{1}{1.25*10^{-3}}=800hz
采样率=1.25∗10−31=800hz
我们使用800hz的采样率去采集1K信号，达到了4K采样率的效果！这便是欠采样，用时间换采样率。
⚽例程1
工程建立
建立过程与STM32HAL ADC+TIM+DMA采集交流信号 基于cubemx_stm32 adc dma 采集交流_四臂西瓜的博客-CSDN博客一样，唯一不同之处是，采样率为800hz。文末会放置本文的完整工程。
采样率
=
时钟：
72
∗
1
0
6
分频：
1000
∗
90
=
800
h
z
采样率=\frac{时钟：72*10^{6}}{分频：1000*90}=800hz
采样率=分频：1000∗90时钟：72∗106=800hz
补充一个知识：我们平时说的预分频PSC是指上图中的Prescaler，此处为1000-1。ARR自动重装载值，是指图中的Counter Period（AutoReload Register），此处为90-1。
运行结果
信号发生器输出：1kHz，1V-2V幅值范围的正弦信号。
VOFA上观察到如下结果：
可以看到一个周期确实是4个点，此处形状和原理分析中、预想的三角形波形不一样，是因为相位没有从零度开始采集，如下图：
采集的点，连成线，就如图中的蓝色线一样，形成“锯齿”状。
这里就可以通过理想插值（增采样），来实现复原，如何增采样呢？可以关注我的“TCQ的电赛小站”，有时间会更新基于CMSIS-DSP的理想插值实现。
🏓例程2
工程建立
现在我们来思考一个问题，我们现在每个周期只采集4个点，如果想每个周期采集10个点，那么采样时间间隔就要从1.25ms变成1.1ms
换算成采样率为：
采样率
=
1
1.1
∗
1
0
−
3
=
909.0909...
h
z
采样率=\frac{1}{1.1*10^{-3}}=909.0909...hz
采样率=1.1∗10−31=909.0909...hz
会发现采样率并不是一个整数，这就麻烦了。如果我们现在想用定时器来出发ADC进行采样，那么分频如何设置？
分频
=
72
∗
1
0
6
909.0909
=
79200.00000792
分频=\frac{72*10^{6}}{909.0909}=79200.00000792
分频=909.090972∗106=79200.00000792
我们现在设置：
P
S
C
=
100
A
R
R
=
792
PSC=100\\ ARR=792
PSC=100ARR=792
运行结果
与例程1一样，信号发生器输出：1kHz，1V-2V幅值范围的正弦信号。
如果读者愿意，可以数一数，会发现一个周期正是10个点！
我们用
909.09
h
z
的采样率，去采集了一个
1
k
的信号，等效采样率
10
k
！！
我们用909.09hz的采样率，去采集了一个1k的信号，等效采样率10k！！
我们用909.09hz的采样率，去采集了一个1k的信号，等效采样率10k！！
🥊难点
是的，欠采样也不是万能的，家家有本难念的经，采样家族里，欠采样也有它的烦恼>_<😕。
在例程2中，我们试图每1.1ms采集一个点，换算成频率后，发现是不是一个整数，非常幸运的是，转换成分频后，可以找到合适的PSC和ARR使得采样率和预期值基本一致。
这个计算过程可以通过公式化简一下：
采样率
=
1
秒
采样时间间隔
采样率=\frac{1秒}{采样时间间隔}
采样率=采样时间间隔1秒
定时器分频 = P S C ∗ A R R = 定时器时钟 采样率 = 定时器时钟 1 秒 采样时间间隔 = 定时器时钟 ∗ 采样时间间隔 \begin{align} 定时器分频&=PSC*ARR=\frac{定时器时钟}{采样率}\\ &=\frac{定时器时钟}{\frac{1秒}{采样时间间隔}}\\ &=定时器时钟*采样时间间隔 \end{align} 定时器分频=PSC∗ARR=采样率定时器时钟=采样时间间隔1秒定时器时钟=定时器时钟∗采样时间间隔
现在有了这个公式，我们来重新计算下例程1、2的分频：
例程
1
：
P
S
C
∗
A
R
R
=
(
72
∗
1
0
6
)
∗
(
1.25
∗
1
0
−
3
)
=
90000
例程1：PSC*ARR=(72*10^{6})*(1.25*10^{-3})=90000
例程1：PSC∗ARR=(72∗106)∗(1.25∗10−3)=90000
例程 2 ： P S C ∗ A R R = ( 72 ∗ 1 0 6 ) ∗ ( 1.1 ∗ 1 0 − 3 ) = 79200 例程2：PSC*ARR=(72*10^{6})*(1.1*10^{-3})=79200 例程2：PSC∗ARR=(72∗106)∗(1.1∗10−3)=79200
挺好的呀，PSC*ARR都是比较好处理的整数。可是当频率高起来后，就不是特别好凑出来了，甚至于根本没法凑出来。比如以1.1us的时间间隔去采集1M信号，等效采样率为10M。
P
S
C
∗
A
R
R
=
(
72
∗
1
0
6
)
∗
(
1.1
∗
1
0
−
6
)
=
7.92
PSC*ARR=(72*10^{6})*(1.1*10^{-6})=7.92
PSC∗ARR=(72∗106)∗(1.1∗10−6)=7.92
看吧，PSC * ARR是小数，但是PSC和ARR必须是整数，没法凑出来🤔。
🥑练习
-
下载例程并且上电调试，观察结果。
-
尝试测量更高频率的信号，比如测量几兆的信号。
-
假如，PSC*ARR可以设置成小数，欠采样可以采集任意频率的信号吗？
提示，学习下ADC的“采样保持”。
🍉后记
本文章收录于：
唐承乾的电赛小站
本文为系列文章中的冰山一角，欢迎进入小站查看。
配套程序：文章来源：https://www.toymoban.com/news/detail-620165.html
配套工程 - Gitee.com文章来源地址https://www.toymoban.com/news/detail-620165.html
到了这里，关于☠️stm32Cubemx欠采样（等效采样）原理讲解与实现 采集高频信号的文章就介绍完了。如果您还想了解更多内容，请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章，希望大家以后多多支持TOY模板网！
