1.本发明属于谐波测量技术领域，具体涉及一种同步跟踪的谐波测量装置及方法。
背景技术：2.在电力系统的测量中，对交流电压(或电流)的谐波测量是一个十分重要的测量项目，尤其是针对目前正热的新能源领域汽车电子、光伏等行业，谐波测量逐步成为国标要求中的产品出厂检测项，生产测量中对模块化扩展与测量精度提出了更高的要求。
3.目前，应用最多的谐波测量方法有两种，一种为加窗插值法，一种为硬件锁相环倍频采样法。
4.加窗插值法的优点是：对硬件要求比较低，可在低采样率时通过插值保证足够计算点数，有效防止谐波测量时点数不够带来的栅栏效应而引起的频率泄露误差。加窗插值法的缺点是：虽然可以通过hanning窗等滤波算法修正误差，但由于是模拟信号插值，插值相位终究有误差，测量精度会受到影响，此种方式测量精度普遍不高且特殊频率下会出现误差偏大的风险。
5.硬件锁相法是通过硬件锁相电路对基波进行倍频，其优点是：准确度较高，随时跟踪基波频率进行倍频采样，响应快。硬件锁相法的缺点是：受硬件制约，对锁相环倍频的稳定性与准确度要求较高，若想与加窗插值法达到相同基波频率的谐波测量能力，电路复杂且成本高，且每次扩展均需要扩展硬件，若进行多路不同步谐波测量，则需要有多路独立锁相环电路，硬件冗余严重。其次，由于谐波测量往往伴随其它数据测量，硬件锁相环倍频直接产生采样频率，控制采样速度，导致采样率低于正常采样，影响其它数值精度，包括有效值rms、谐波分量hrms等。
6.所以，如何获得更精准的谐波测量，更标准的模块和更便捷的扩展，是本发明要解决的问题。
技术实现要素：7.为了解决现有技术中存在的技术问题，本发明提供一种同步跟踪的谐波测量装置及方法。本发明所采用的技术方案如下：
8.一种同步跟踪的谐波测量装置，包括：主控单元、同步采样跟踪单元和模拟信号采样单元，所述的主控单元采用mcu，所述的采样跟踪单元采用fpga，所述的模拟信号采样单元为若干个测量卡；主控单元中设置有第一双口ram模块和fft计算模块，同步采样跟踪单元中设置有第二双口ram模块、时钟模块a、时钟模块b、数据转存模块、频率跟踪模块、adc计算模块、频率计算模块和adc控制器，频率计算模块和adc控制器分别与测量卡进行数据交互；频率计算模块从测量卡读取信号过零点信息，计算当前频率；adc控制器接收跟踪处理后的时钟信号，并转化为cnv控制adc转换速度，通过spi总线与adc进行数据传输；频率跟踪模块实时同步跟踪频率计算模块所得频率，结合当前谐波倍频数换算并控制时钟模块a和时钟模块b的时钟输出；时钟模块a和时钟模块b根据输入的cnv调节参量k对时钟进行任意
倍数分频，得到目标频率时钟，输入adc控制器控制adc转换频率；adc控制器向adc计算模块传输数据，adc计算模块除处理原始波点外同时计算均方根值和峰值，adc计算模块将数据发送至数据转存模块，第二双口ram模块通过fsmc数据总线与mcu建立连接进行通信，通过寻址进行ram区数据读取，进行adc采样数据传输；第一双口ram模块通过第二双口ram模块读取数据转存模块中的数据，fft计算模块对读取的adc采样数据进行傅里叶变换并进行谐波数据处理。
9.一种同步跟踪的谐波测量方法，应用前述的一种同步跟踪的谐波测量装置，包括以下步骤：
10.①
捕获模拟信号采样单元中过零点信息，采样同步跟踪单元通过频率计算模块计算得到被采集信号的频率；
11.②
在采样同步跟踪单元中计算目标采样频率fn，fn＝f
×
ct，f为被采集信号频率，ct为单周期进行采样的次数；
12.③
在采样同步跟踪单元中计算fpga最大采样频率f
nmax
，f
nmax
＝fn×
p
max
，p
max
为不高于采样能力的最大整数倍；
13.④
计算cnv调节参量k值，使用k值控制时钟，k值根据被测频率及目标频率进行计算，k＝f
nmax
÷
f0，f0等于时钟模块的最大支持频率除以计数器的满值；
14.⑤
根据k值切换adc转换速度，接到k值的时钟模块a分频，生成目标cnv并输出给adc控制器，adc控制器根据该信号改变采样cnv，采样过程中，频率跟踪模块再次重新计算k值，并下发给时钟模块b，当本周期达到采样点数时，切换adc控制器的时钟模块b，时钟模块a和时钟模块b交错使用，在采样频率发生变化时，adc快速切换无相位延迟；
15.⑥
根据采样频率f
nmax
控制cnv信号进行adc转换，得到整周期数组s(i)，sf(i)＝s(p
max
×
i)，分别取sf(0)＝s(0)，sf(1)＝s(3)，sf(2)＝s(6)，sf(i)＝s(p
max
×
i)，
……
，i指第n个数，得到新数组sf(ct)，得到原始采样点，满足fft运算条件；
16.⑦
将sf(ct)送入主控单元中，运用fft模块进行fft快速傅里叶变换，得到谐波数据。
17.本发明的有益效果：
18.1、精度高、易计算：保证数据基于基波频率整数倍采样，可使用fft快速傅里叶变换进行计算，方便快捷，准确度高，保证了测量精度。
19.2、低成本、易扩展：多路谐波测量不需额外增加硬件，增加fpga软核节省成本与空间，解决现有产品难以升级的问题。
20.3、模块化程度高，系统架构清晰，可扩展性强，保持高采样率，测量精度不会受锁相环带来的降频影响。
附图说明
21.图1是本发明实施例的原理拓扑图；
22.图2是本发明实施例的测量步骤流程图。
具体实施方式
23.下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施
例是本发明的一部分实施例，而不是全部的实施例。
24.如图1所示，是本发明实施例的原理拓扑图，图1中的箭头表示数据传输。本发明使用模拟信号同步跟踪采样技术，灵活控制采样频率，达到保证精度的同时实现快速谐波测量的目的，系统总共分为三部分：
25.模块一：主控单元mcu。主控单元mcu中设置有第一双口ram模块(m1)和fft计算模块(m2)，第一双口ram模块(m1)和fft计算模块(m2)均为软件功能模块，第一双口ram模块(m1)通过fsmc总线与模块二中的第二双口ram模块(f1)进行数据交互。
26.fft计算模块(m2)主要用来进行数据处理，将m1采集到的adc采样数据进行傅里叶变换并进行谐波数据处理。
27.这里使用arm官方库函数进行计算：
28.(1)傅里叶变换使用函数arm_cfft_f32(&arm_cfft_f32_len1024,fft_inputbuf,0,1)来对输入数据进行傅里叶变换。它的输入参数含义如下：
29.arm_cfft_f32_len1024为傅里叶变换结构体，1024是要运算的点数。在进行其他点数的计算时，比如说32个点的fft，就可以用arm_cfft_sr_f32_len32。fft_inputbuf为傅里叶变换需要处理的数据的首地址，第三个参数0是正变换，1是反变换，第四个参数默认为1。
30.经过傅里叶变换后的结果仍然为复数，虚部和实部的比可以计算出频率点的相位，这个在这里不进行考虑，直接对复数取模。
31.(2)arm_cfft_radix4_f32进行傅里叶变换，得到变换后的数据实部与虚部。
32.(3)arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,fft_length)用于计算模值，执行取模操作，后可求出幅值、相位等参数。
33.其中，ft_inputbuf源数据，形式为复数，fft_outputbuf取完模后的数据，形式为实数，fft_length是取模的点数。
34.内部计算公式为：
[0035][0036]
模块二：同步采样跟踪单元(fpga)。以fpga作为硬件支撑，在其中设置软件功能模块，包括：第二双口ram模块(f1)、时钟模块a[100m](f2)、时钟模块b[100m](f3)、数据转存模块(f4)、频率跟踪模块(f5)、adc计算模块(f6)、频率计算模块(f7)和adc控制器(f8)。
[0037]
频率计算模块(f7)和adc控制器(f8)分别与若干个测量卡进行数据交互，频率计算模块(f7)从测量卡读取信号过零点信息，计算当前频率；adc控制器(f8)接收跟踪处理后的时钟信号，并转化为cnv(convst adc开始转换信号)控制adc采样速度，通过spi总线与adc进行数据传输。
[0038]
频率跟踪模块(f5)：实时同步跟踪频率计算模块(f7)所得频率，结合当前谐波倍频数换算并控制时钟模块a和b(f2、f3)的时钟输出。
[0039]
第二双口ram模块(f1)通过fsmc数据总线与mcu建立连接进行通信，通过寻址进行ram区数据读取，这里主要进行adc采样数据传输。
[0040]
时钟模块a[100m](f2)和时钟模块b[100m](f3)根据输入的cnv调节参量k对时钟进行任意倍数分频，得到目标频率时钟，输入adc控制器(f8)从而控制adc采样频率。f2、f3两组时钟模块交替使用，当一组正在使用时，另一组进行配置、初始化新的时钟，当频率更
改时快速切换并输出控制信号，避免了由于时钟初始化引入的延时，此处100m为最大支持100mhz频率的时钟。
[0041]
adc控制器(f8)向adc计算模块(f6)传输数据，adc计算模块(f6)除处理原始波点外还需计算均方根值和峰值等相关数据。adc计算模块(f6)将数据发送至数据转存模块(f4)，数据转存模块(f4)是ram缓存区，转存adc计算模块(f6)数据，等待mcu读取。
[0042]
模块三：模拟信号采样单元，包括n个测量卡。
[0043]
测量卡是现有的产品，测量卡主要构成有采样电阻、运放与adc电路，主要功能是对接入的信号进行采样，读回的数据共有两部分，一部分通过fpga中的adc控制器(f8)spi总线将数据读出，读回的为固定个数采样点的ad字。另一部分由频率计算模块(f7)读回，读回的是过零点信号，通过过零点的周期数与过零点计时进行频率值计算。
[0044]
由于使用fft(快速傅里叶变换，fast fourier transform)进行谐波计算，需要保证采样个数为谐波周期的整数倍，而要精准跟踪被测频率，并及时控制模拟信号采样部分adc的转换，需要本发明中各个模块紧密配合。
[0045]
如图2所示，是本发明实施例的测量步骤流程图。一种同步跟踪的谐波测量方法，应用前述的一种同步跟踪的谐波测量装置，包括以下步骤：
[0046]
①
捕获模拟信号采样单元中过零点信息，采样同步跟踪单元(f5)通过频率计算模块(f7)计算得到被采集信号的频率。
[0047]
②
在采样同步跟踪单元(f5)中计算目标采样频率fn。
[0048]
例如：被采集信号频率f为50hz的信号，计算100次谐波，单周期进行ct次采样，为保证谐波数据准确度这里ct选择1024次，所以目标采样频率应为：
[0049]fn
＝f
×
ct＝50
×
1024＝51.2khz。
[0050]
③
在采样同步跟踪单元(f5)中计算fpga最大采样频率。
[0051]
例：若模拟信号采样单元最大采样能力为200khz，若直接用51.2khz采样会损失谐波数据精度，本发明可将采样频率增加至不高于采样能力的fn最大整数倍频率，故当fn低于200khz时，fn最大整数倍频率f
nmax
为：
[0052]fnmax
＝fn×
p
max
(最大整数倍)＝51.2*3＝153.6khz。
[0053]
优于单纯锁相所提供的频率。
[0054]
④
计算cnv调节参量k值，这里使用k值来控制时钟，k值则是根据被测频率及目标频率进行计算。
[0055]
通过控制时钟模块(f2)和(f3)给adc控制器(f8)发送时钟信号，保证adc转换快速进行，很好的利用了fpga的实时性与并发性，保证信号采样同步。fpga中的寄存器和片上内存(bram)是属于各自的控制逻辑的，无需不必要的仲裁和缓存且fpga是单时钟周期信息处理，所有操作基本瞬时完成，结合以上两点保证了fpga的实时性与同步性，所以通过fpga控制可保障的采样与数据处理同步性。
[0056]
本发明使用100m时钟模块，使用32位计数器进行计时，满值为2
32
＝4294967296，即计数频率为：
[0057]
f0＝100000000(100mhz)
÷232
＝0.02328306hz，将计算好的f
nmax
带入可得到k值，系数k＝f
nmax
÷
f0＝153600
÷
0.02328306＝6597071。
[0058]
⑤
根据k值切换adc转换速度。
[0059]
首先，采样过程中，计算k值k＝f
nmax
÷
f0，其中的f0是稳定不变的，f
nmax
是随信号变化的，和被测信号f有关，切换时钟发生器为时钟模块a(f2)，并下发k值至时钟模块a(f2)，接到k值的时钟模块a分频，生成目标cnv并输出给adc控制器(f8)，adc控制器(f8)此时根据该信号改变采样cnv，采样过程中，频率跟踪模块(f5)再次重新计算k值，并下发给时钟模块b(f3)，当本周期达到采样点数时，切换adc控制器(f8)的时钟模块b(f3)。采样cnv是指convst adc开始转换信号，是一个脉冲信号，连接adc的convst引脚，控制adc转换速度，在使用时钟模块b(f3)采样过程中，频率跟踪模块(f5)再次重新计算k值，等待时钟切换，这样ab时钟交错使用方式，保证在采样频率发生变化时，adc可快速切换，无相位延迟。
[0060]
⑥
得到原始采样点：
[0061]
根据采样频率f
nmax
控制cnv信号进行adc转换，可得到整周期数组s(i)，此时，单周期得到的点数为p
max
×
ct＝3
×
1024＝3072，分别取sf(0)＝s(0)，sf(1)＝s(3)，sf(2)＝s(6)，sf(i)＝s(p
max
×
i)，
……
(i指第n个数)，得到新数组sf(ct)(ct＝1024)，该数组为整周期50hz 1024次采样，满足fft运算条件。本发明实施例中，采样信号进行fft变换只需要使用51.2khz采样率单周期进行1024次采样即可，但是，为了保证均方根等值计算精度，本发明加大了采样频率，所以使用153.6k进行采样即3倍的采样频率，最后得到3072个采样点，所以，在进行fft计算前，需要从3072个采样点中均匀筛选出1024个采样点进行fft运算。
[0062]
⑦
谐波数据计算：
[0063]
最后将sf(ct)送入mcu中运用fft模块(m2)进行fft快速傅里叶变换，可得到谐波数据。
技术特征：1.一种同步跟踪的谐波测量装置，包括：主控单元、同步采样跟踪单元和模拟信号采样单元，所述的主控单元采用mcu，所述的采样跟踪单元采用fpga，所述的模拟信号采样单元为若干个测量卡；其特征在于，主控单元中设置有第一双口ram模块和fft计算模块，同步采样跟踪单元中设置有第二双口ram模块、时钟模块a、时钟模块b、数据转存模块、频率跟踪模块、adc计算模块、频率计算模块和adc控制器，频率计算模块和adc控制器分别与测量卡进行数据交互；频率计算模块从测量卡读取信号过零点信息，计算当前频率；adc控制器接收跟踪处理后的时钟信号，并转化为cnv控制adc转换速度，通过spi总线与adc进行数据传输；频率跟踪模块实时同步跟踪频率计算模块所得频率，结合当前谐波倍频数换算并控制时钟模块a和时钟模块b的时钟输出；时钟模块a和时钟模块b根据输入的cnv调节参量k对时钟进行任意倍数分频，得到目标频率时钟，输入adc控制器控制adc转换频率；adc控制器向adc计算模块传输数据，adc计算模块除处理原始波点外同时计算均方根值和峰值，adc计算模块将数据发送至数据转存模块，第二双口ram模块通过fsmc数据总线与mcu建立连接进行通信，通过寻址进行ram区数据读取，进行adc采样数据传输；第一双口ram模块通过第二双口ram模块读取数据转存模块中的数据，fft计算模块对读取的adc采样数据进行傅里叶变换并进行谐波数据处理。2.根据权利要求1所述的一种同步跟踪的谐波测量装置，其特征在于，时钟模块a和时钟模块b两组时钟模块交替使用，当一组正在使用时，另一组进行配置、初始化新的时钟，当频率更改时快速切换并输出控制信号。3.一种同步跟踪的谐波测量方法，其特征在于，应用如权利要求2所述的一种同步跟踪的谐波测量装置，包括以下步骤：
①
捕获模拟信号采样单元中过零点信息，采样同步跟踪单元通过频率计算模块计算得到被采集信号的频率；
②
在采样同步跟踪单元中计算目标采样频率f
n
，f
n
＝f
×
ct，f为被采集信号频率，ct为单周期进行采样的次数；
③
在采样同步跟踪单元中计算fpga最大采样频率f
nmax
，f
nmax
＝f
n
×
p
max
，p
max
为不高于采样能力的最大整数倍；
④
计算cnv调节参量k值，使用k值控制时钟，k值根据被测频率及目标频率进行计算，k＝f
nmax
÷
f0，f0等于时钟模块的最大支持频率除以计数器的满值；
⑤
根据k值切换adc转换速度，接到k值的时钟模块a分频，生成目标cnv并输出给adc控制器，adc控制器根据该信号改变采样cnv，采样过程中，频率跟踪模块再次重新计算k值，并下发给时钟模块b，当本周期达到采样点数时，切换adc控制器的时钟模块b，时钟模块a和时钟模块b交错使用，在采样频率发生变化时，adc快速切换无相位延迟；
⑥
根据采样频率f
max
控制cnv信号进行adc转换，得到整周期数组s(i)，sf(i)＝s(p
max
×
i)，分别取sf(0)＝s(0)，sf(1)＝s(3)，sf(2)＝s(6)，sf(i)＝s(p
max
×
i)，
……
，i指第n个数，得到新数组sf(ct)，得到原始采样点，满足fft运算条件；
⑦
将sf(ct)送入主控单元中，运用fft模块进行fft快速傅里叶变换，得到谐波数据。
技术总结本发明属于谐波测量技术领域，具体涉及一种同步跟踪的谐波测量装置及方法，谐波测量装置包括：主控单元、同步采样跟踪单元和模拟信号采样单元，主控单元中设置有第一双口RAM模块和FFT计算模块，同步采样跟踪单元中设置有第二双口RAM模块、时钟模块A、时钟模块B、数据转存模块、频率跟踪模块、ADC计算模块、频率计算模块和ADC控制器，时钟模块A和时钟模块B根据输入的CNV调节参量K对时钟进行任意倍数分频，得到目标频率时钟，输入ADC控制器控制ADC采样频率；FFT计算模块对读取的ADC采样数据进行傅里叶变换并进行谐波数据处理。本发明精度高、易计算，低成本、易扩展，模块化程度高。模块化程度高。模块化程度高。
技术研发人员：许晓东 白洪超
受保护的技术使用者：青岛艾诺智能仪器有限公司
技术研发日：2022.07.25
技术公布日：2022/11/1
