单频振幅抑噪测量方法与流程

1.本发明涉及电子通信测控技术领域，具体涉及一种单频振幅抑噪测量方法。


背景技术：

2.因为在电子通信测控等应用中大量运用信号振幅调制，例如调幅通信应用就是通过调制某频率信号的振幅大小来传输通信数据，所以测量电子信号振幅相关应用非常广泛；因为在测量过程会混入各种噪声，所以实际测量数据必然存在误差；目前主要采用滤波、平均、以及加权等方式降低噪声的干扰；常见滤波方案需要事先知道信号所在的频带，然后设计滤波器的通频带，用抑制带外噪声的方法降低噪声干扰；还有傅立叶分解分析频谱的方法，就是把信号频带内所有频率分量全都测量出来，之后再寻找出其中最大的那个信号振幅，计算量巨大。


技术实现要素：

3.本发明为解决现有测量电子信号振幅的方法需要将信号频带内所有频率分量全部测量，然后寻找最大信号的振幅，导致存在计算量大的问题，提供一种单频振幅抑噪测量方法。
4.单频振幅抑噪测量方法，该方法通过单频振幅抑噪测量系统实现，所述单频振幅抑噪测量系统包括数据输入端，存贮器，处理器和输出端；
5.数据输入端输入待测信号保存在存贮器中；
6.所述存贮器保存从数据输入端获得的待测信号；在每一次测量过程中在逻辑上连续存贮n个数据，形成一个一维数组，用s[i]表示；其中s为数组名称，i为下标索引序号；下标索引序号i的最小值为1，最大值为n；
[0007]
所述处理器分析保存在存贮器中的数据，计算所述待测信号的振幅a；具体过程为：
[0008]
步骤一、测量待测信号区间；
[0009]
步骤a1、定义变量i，初值为1；定义变量j，初值为零；定义有n个元素的一维数组q；
[0010]
步骤a2、在数组s中取得以i，i+1，i+2，以及i+3为索引序号的四个元素s[i]，s[i+1]，s[i+2]，以及s[i+3]，计算这四个元素的和，保存在数组q中以i为下标索引序号的元素中；
[0011]
步骤a3、变量i＝i+1，如果变量i大于n，则执行频骤a4，否则，执行步骤a2；
[0012]
步骤a4、将变量i的值设置为1，然后将变量u以及变量d都设置成数组q中第1个元素；
[0013]
步骤a5、比较变量u与数组q中以i为索引序号的元索q[i]，如果变量u小于q[i]，则将变量u的值设置成为q[i]的值；
[0014]
步骤a6、比较变量d与数组q中以i为索引序号的元索q[i]，如果变量d大于q[i]，则将变量d的值设置成为q[i]；
[0015]
步骤a7、变量i＝i+1，如果变量i大于n，则执行步骤a8，否则，执行步骤a5；
[0016]
步骤a8、计算变量u和变量d的平均值c；
[0017]
步骤a9、定义变量w，如果数组q中第一个元素q[1]大于变量c，则将w的值设为1，否则将w的值设为0；
[0018]
步骤a10、将变量i的值设置为1；
[0019]
步骤a11、比较数组q中以i为索引序号的元素q[i]与变量c的大小，再比较数组q中以i+1为索引序号的元素q[i+1]与变量c的大小；如果q[i]大于变量c，并且q[i+1]小于变量c，并且w的值为1，则将变量d的值改为变量i之后执行步骤a13；否则，如果q[i]小于变量c，并且q[i+1]大于变量c，并且w的值为0，则将变量d的值改为变量i之后执行步骤a13；
[0020]
步骤a12、将变量i的值加1之后，如果变量i等于n，则执行步骤a13，否则，执行步骤11；
[0021]
步骤a13、比较数组q中以i为索引序号的元素q[i]与变量c的大小，再比较数组q中以i+1为索引序号的元素q[i+1]与变量c的大小；如果q[i]大于变量c，并且q[i+1]小于变量c，并且w的值为1，则将变量u的值改为变量i；否则，如果q[i]小于变量c，并且q[i+1]大于变量c，并且w的值为0，将变量u的值改为变量i；
[0022]
步骤a14、将变量i的值加1之后，如果变量i等于n，则执行步骤a15，否则，执行步骤a13；
[0023]
步骤a15、变量d保存数组s中测量区间的起点索引号，变量u保存终点索引号；定义变量g作为数据数量，设置g的值等于u-d+1；执行步骤二；
[0024]
步骤二、去除直流成分；具体步骤为：
[0025]
步骤b1、设置变量i的初值为d；定义变量l，初值为零；
[0026]
步骤b2、将变量l加上数组s以i为索引序号的元素s[i]；
[0027]
步骤b3、将变量i＝i+1，如果变量i大于g，则执行步骤b4，否则，执行步骤b2；
[0028]
步骤b4、将变量l除以g以后的数值，重新保存到变量l中；再将变量i的值设置为d；
[0029]
步骤b5、取数组s以i为索引序号的元素s[i]，将s[i]的值减去变量l之后的数值，重新保存到s[i]中；
[0030]
步骤b6、将变量i＝i+1，如果变量i大于g，则执行步骤b7，否则，执行步骤b5；
[0031]
步骤b7、数组s中保存的数据则为去除直流成分的信号数据；
[0032]
步骤三、对步骤二去除直流成分的信号数据计算平方均值m与绝对平均值e；
[0033]
步骤c1、将变量i的值设为d；重新定义变量m，初值也为零；定义变量e，初值为零；
[0034]
步骤c2、取得数组s以i为索引序号的元素s[i]，将变量m加上s[i]的平方；将变量e加上s[i]的绝对值；
[0035]
步骤c3、i＝i+1，如果变量i大于g，则执行步骤c4，否则，执行步骤c2；
[0036]
步骤c4、将变量m除以g以后的数值，重新保存到变量m中；将变量e除以n以后的数值，重新保存到变量e中；
[0037]
步骤四、根据步骤三获得的计算原始信号的振幅a；
[0038]
步骤d1、定义变量x，数值等于变量e乘以圆周率π；
[0039]
步骤d2、定义变量y，数值等于x的平方减去6倍的变量m；
[0040]
步骤d3、定义变量z，数值等于变量y开平方；
[0041]
步骤d4、变量x加上变量z之后，再除以3等于抑制噪声之后原始信号的振幅a，将所述振幅a从输出端输出。
[0042]
本发明的有益效果：本发明所述的测量方法可以在混有随机噪声的信号中测量出最强那个信号振幅的方法；具体过程是通过分析经过数字量化后的信号数据，先去除直流成分，再计算数据的方均值和绝对平均值，最后即可计算获得混合信号中强度最大那个频率分量的振幅。
[0043]
本发明与滤波器方案相比，因为测量过程与信号频率无关，所以不需要设计通频带；还可以抑制传统滤波方案无法解决的带内噪声；本方法计算量小，与傅立叶方案相比，本发明因无须计算出每种频率分量，所以可在快速测量场合应用；数字振幅调制方式大量应用在现代各种通信网络、4g、5g、wifi、蓝牙、以及数字广播电视等通信系统中，在发送端振幅调制的载波信号，被接收端通过测量信号的瞬时振幅而还原通信数据；调幅广播的节目声音信号经过振幅调制载波，无线传输到收音机之后被检波器测量振幅从而还原节目内容；运用本发明所述的方法可以为这些电子系统提高抗干扰能力。
附图说明
[0044]
图1为本发明所述的单频振幅抑噪测量方法原理框图。
具体实施方式
[0045]
结合图1说明本实施方式，单频振幅抑噪测量方法，该方法通过单频振幅抑噪测量系统实现，所述系统包括数据输入端，存贮器，处理器以及输出端；
[0046]
数据输入端输入待测信号，若待测信号为模拟信号，则需要用模数转换器转换之后再送给数据输入端；要求对输入待测信号以固定频率f采样得到数字量数据；必须保证采样频率f大于上限数值h的8倍以上；要求待测信号在单次测量的期间，振幅保持不变；采样时长的最小值为待测信号的一个周期。
[0047]
所述待测信号为原始信号混入随机噪声信号；原始信号为单一频率的余弦信号，具体频率用f表示，f的具体数值未知，但是知道f必定低于一个最高频率上限h；另外还已知原始信号的振幅大于噪声信号振幅；在单次测量的过程中，原始信号的振幅为a保持不变；本系统的目标是在待测信号中抑制噪声信号的干扰测量出原始信号的振幅a。
[0048]
在测量精度方面，采样时长接近待测信号整数倍周期时精度较高；因此，如果事先不知道信号的大约频率，可通过分析过零点，极大值，以及极小值等特征点，粗测得到信号大约频率，然后设计最合适的采样时长就能获得更高的测量精度；因为加大采样数量有利于提高测量精度，所以当采样数量较大时，即使采样时长并没有接近待测信号的整数倍也能获得较高的测量精度。
[0049]
所述存贮器2保存从数据输入端获得的待测信号数据；在每一次测量过程中需要一批数据有n个数值，在逻辑上连续存贮，形成一个一维数组的数据结构，用s[i]表示；其中s为数组名称，i为下标索引序号；下标索引序号i的最小值为1，最大值为n；
[0050]
所述处理器3分析保存在存贮器2中的数据，分析计算出信号振幅的具体方式分成测量信号区间，去除直流成分,计算方均m与绝对平均值e，以及计算原始信号的振幅a几个步骤；对于没有直流分量的纯交流信号，则可省去测量信号区间，以及去除直流成分这两个
步骤；
[0051]
一、测量信号区间；
[0052]
a1、定义变量i，初值为1；定义变量j，初值也为零；定义变量中，初值也为零；定义有n个元素的一维数组q；
[0053]
a2、在数组s中取得以i，i+1，i+2，以及i+3为索引序号的四个元素s[i]，s[i+1]，s[i+2]，以及s[i+3]，计算这四个元素的和，保存在数组q中以i为下标索引序号的元素中；
[0054]
a3、把变量i加1之后，如果变量i大于n，则执行频骤a4，否则执行频骤a2；
[0055]
a4、把变量i设为1，然后把变量u以及变量d都设置成数组q中第1个元素；即：u＝d＝q[0]；
[0056]
a5、比较变量u与数组q中以i为索引序号的元索q[i]，如果变量u小于q[i]，那么就把变量u的值设置成为q[i]；
[0057]
a6、比较变量d与数组q中以i为索引序号的元索q[i]，如果变量d大于q[i]，那么就把变量d的值设置成为q[i]；
[0058]
a7、把变量i加1之后，如果变量i大于n，则执行步骤a8，否则执行步骤a5；
[0059]
a8、计算变量u和变量d的平均值c；
[0060]
a9、定义一个变量w，如果数组q中第一个元素q[1]大于变量c，则把w的值设为1，否则把w的值设为0；
[0061]
a10、把变量i设为1；
[0062]
a11、比较数组q中以i为索引序号的元素q[i]与变量c的大小，再比较比较数组q中以i+1为索引序号的元素q[i+1]与变量c的大小；如果q[i]大于变量c，并且q[i+1]小于变量c，并且w的值为1，那么把把变量d的值改为变量i之后执行步骤a13；否则如果q[i]小于变量c，并且q[i+1]大于变量c，并且w的值为0，那么把把变量d的值改为变量i之后执行步骤a13；
[0063]
a12、把变量i的值加1之后，如果变量i等于n，则执行下一步，否则执行步骤a11；
[0064]
a13、比较数组q中以i为索引序号的元素q[i]与变量c的大小，再比较比较数组q中以i+1为索引序号的元素q[i+1]与变量c的大小；如果q[i]大于变量c，并且q[i+1]小于变量c，并且w的值为1，那么把把变量u的值改为变量i；否则如果q[i]小于变量c，并且q[i+1]大于变量c，并且w的值为0，那么把把变量u的值改为变量i；
[0065]
a14、把变量i的值加1之后，如果变量i等于n，则执行下一步，否则执行步骤a13；
[0066]
a15、此时，变量d保存数组s中测量区间的起点索引号，变量u保存终点索引号；定义变量g作为数据数量，设置g的值等于u-d+1；
[0067]
二、去除直流成分；
[0068]
传统计算信号中直流成分的方案是计算所有数据，未能考虑信号数据的周期区间问题；因为只有在信号完整周期区间的计算才准确，所以传统方式的精度较低；本发明在从第一大步骤“测量信号区间”获得了信号数据的周期区间，所以计算直流成分精度更高；具体去除直流成分的步骤为；
[0069]
b1、定义变量i，初值为d；定义变量l，初值为零；
[0070]
b2、把变量l加上数组s以i为索引序号的元素s[i]；
[0071]
b3、把变量i加1之后，如果变量i大于g，则执行步骤b4，否则执行步骤b2；
[0072]
b4、把变量l除以g以后的数值，重新保存到变量l中；再把变量i的值设为d；
[0073]
b5、取得数组s以i为索引序号的元素s[i]，把s[i]减去变量l之后的数值，重新保存到s[i]中；
[0074]
b6、把变量i加1之后，如果变量i大于g，则执行步骤b7，否则执行步骤b5；
[0075]
b7、此时数组s中保存的数据就是已被去除直流成分的信号数据；
[0076]
三、计算方均m与绝对平均值e；
[0077]
c1、把变量i的值设为d；定义变量m，初值也为零；定义变量e，初值也为零；
[0078]
c2、取得数组s以i为索引序号的元素s[i]，把变量m加上s[i]的平方；把变量e加上s[i]的绝对值；
[0079]
c3、把变量i加1之后，如果变量i大于g，则执行步骤c4，否则执行步骤c2；
[0080]
c4、把变量m除以g以后的数值，重新保存到变量m中；把变量e除以n以后的数值，重新保存到变量e中；
[0081]
四、计算原始信号的振幅a；
[0082]
d1、定义变量x，数值等于变量e乘以圆周率π；
[0083]
d2、定义变量y，数值等于x的平方减去6倍的变量m；
[0084]
d3、定义变量z，数值等于变量y开平方；
[0085]
d4、变量x加上变量z之后，再除以3就等于抑制噪声之后原始信号的振幅a，得到原始信号的振幅a从输出端输出。
[0086]
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0087]
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。