22年电赛F题信号调制度测量思路和常见问题分享，超详细。基于stm32f103

咱们的基本大体思路就是:通过对输入的调制信号通过单片机进行带通抽样采样后得到的频谱图，分析后由图像特性得到调制类型和调制度以及调制频率，再通过DDS发出对应的调制信号即可。

 

我们再细说:

1，采样频率的选择

由于我使用的stm32f103最大主频是72M,在不修改主频的情况下，满足ADC最快工作频率(14MHz)，6分频后为12MHz，考虑ADC最快采样周期12.5+1.5=14，最大采样频率只能达到0.8M左右，即使开双ADC也就区区1.6M，想要相对准确的的还原题目要求的30M，至少需要4*30=120MHz，看来是要上FPGA了，可是我们没有...

好在题目中有提到带通抽样定律，那什么是带通抽样，什么又是混频呢？我给挂个链接: 😁😁[带通抽样定律及其公式推导]

带通抽样理论告诉我们，如果我们以采样频率 ( f_s ) 对信号进行采样，采样后的信号频谱会在频域内以 ( f_s ) 为周期重复。为了防止在重复过程中发生频谱重叠，常见选择是 ( f_s ) 应该至少是信号带宽的两倍。为了确保这个条件的应用，使载波频率 ( f_c ) 落入合适的带通区域内，我们通常选择 ( f_s ) 使其是 ( f_c ) 的一个约数或倍数，这样可以避免频谱在变换过程中重叠。

观察题目中的频率步进频率，综合分析下来我们选择使用100kHz进行采样。

 

 

2 ，对频谱进行分析

2-1 频谱图像变化

 

频谱波峰总是出现在频谱的两边，是由实数信号的傅里叶变换共轭对称性、采样后的频域重复、以及对常见幅值谱和能量谱的关心所共同决定的。这些基本特性确保了信号的频谱在频率轴的正负部分对称出现，形成频谱两边看到的波峰。于是我们把图像进行移动位置，把他放在图像的中间，且要保证波峰之间的索引的差值不变，我们只需要对fftout数值的前前512个数据分别交换位置就可以了。(我的adc一组是采集了1024个数据）

我来画个图

 

前

 

后

 

这样是不是就好分析啦

 

 

2-2 设备没有频偏

假设你现在输入的是10MHz的AM调制信号，可能会得到一条直冲云霄的一条竖线，这时因为为了保证ADC采样的电压为正值，我们往往需要加上一个偏置电压，比如我是1v峰峰值的正弦信号加了1v的偏置，我们需要先把直流分量去掉，我这里进行的是均值滤波，将fftin数组的值都减去他们的均值即可。

如果你的设备足够好的话，打出来的频谱将得到中间高两边低3条竖线。

2-2-1 调制信号频率

调制信号的频率 = (主峰-次峰）的索引*(采样频率/1024) 。

2-2-2 调制类型判断

可以通过波峰数量判断调制类型，AM调制的波峰数量不会随着调制度变化，如果是FM调制会有大于3条的频谱线(调制度越大频谱尖峰数量越多），如果是CM(未调制）就只有一条线了

2-2-3再说调制度  

对于AM调制，由通信原理书上的公式可以改成，调制度 ma =k*(次峰幅值/主峰幅值），k值根据自己的情况设定，我的是1.5左右。

特别需要注意的是，如果发现你的峰值变化比较大难以测准调制度时，需要把fftin和fftout等数组的类型全部转化为浮点型，舍弃小数对幅值影响比较大。

 

对于FM调制，调制度 mf = k*(U/f）

U为解调后的波形幅度，f为解调后的波形频率

 

 

2-3 设备存在频偏

(比如我使用的设备，每10M信号大约会产生500k左右的频偏)

如果你和我一样，用的信号源过于古老，滤掉直流分量后原来的一条线都变成了2个峰，注意不是两条线，而是2个峰，这就导致可能会出现主峰边上的值也会大于2个次峰的情况(如上图我画出来的），所以寻找峰值的方法也得改，我这里加了个判断，只有第二个峰的与上个峰的索引差值大于5才认定为一个峰。

为了方便分析，我只分析了左半部分的频谱，也就是说只要按上述方法找到左半部分的3个峰，第二个和第三个峰的中间就是次峰的位置，次峰的幅值就取找到的第二大和第三大的均值。

主峰也是一样的方法，于是我定义了一个函数，可以找到给定区间内最大的2个值的索引，且两个值的差值必须大于5。

由于系统的变化不是线性的，在判读出给定信号后，可以根据范围自行设定值。

其余的思想和没有频谱的是一样的，只是数据处理的不同

 

2-4 就是调制类型的判断，在波峰的基础上再加上最大幅值的判断，我使用的方法比较简单粗暴，就是分析频谱发现我的FM主峰结果会一直大于AM，当然这得通过自己看频谱，或者把数据打到串口中观察。