CN110297199A - 一种基于全相位fft的铯光泵磁力仪频率测量方法及*** - Google Patents
一种基于全相位fft的铯光泵磁力仪频率测量方法及*** Download PDFInfo
- Publication number
- CN110297199A CN110297199A CN201910606439.1A CN201910606439A CN110297199A CN 110297199 A CN110297199 A CN 110297199A CN 201910606439 A CN201910606439 A CN 201910606439A CN 110297199 A CN110297199 A CN 110297199A
- Authority
- CN
- China
- Prior art keywords
- phase
- fft
- window
- real
- imaginary parts
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/175—Spectrum analysis; Fourier analysis by delay means, e.g. tapped delay lines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/26—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux using optical pumping
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Mathematical Physics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
一种基于全相位FFT的铯光泵磁力仪频率测量方法及***,将铯光泵磁力仪输出的磁共振信号经放大整形成方波信号,再将方波信号以及方波信号延迟一个固定的时间后分别通过ADC进行采样,得到离散的输入信号序列及输入信号延时序列;控制FPGA分别对两个序列使用汉宁窗进行全相位预处理;调用FFT IP核对进行了全相位预处理后的序列分别做FFT变换,输出数据的实部和虚部,在一个采样周期内选取出最大幅值处的实部和虚部;通过CORDIC算法根据最大幅值处的实部和虚部推算出对应的相位值,将求出的两个相位值传入单片机进行频率计算,从而得出频率测量结果。本发明可以有效地防止频谱泄漏,并可以极大地节省FPGA的资源,提高效率,且具有高精度和良好的稳定性。
Description
技术领域
本发明涉及铯光泵磁力仪领域,更具体地说,涉及一种基于全相位FFT的铯光泵磁力仪频率测量方法及***。
背景技术
铯光泵磁力仪高分辨率高采样率的优点使其能够应用于航空磁测等复杂环境下的磁场测量场合,它是以铯原子能级在外磁场中的塞曼***为理论基础,通过光抽运效应和光磁共振技术结合而制成的仪器。由于铯原子的磁共振频率值与被测外磁场值成正比,所以只要能够测出铯原子的共振信号频率即可获得外磁场的值。常见的频率测量方法有直接测频法、直接测周期法、等精度测频法、多周期同步法等,但其无法同时满足航空磁测所需的高分辨率、高精度且能连续测量的要求。
发明内容
本发明要解决的技术问题在于,针对现有技术常见的频率测量方法无法同时满足航空磁测所需的高分辨率、高精度且能连续测量的要求的技术缺陷,本发明提出了一种利用全相位FFT对铯光泵磁力仪的磁共振信号进行频率测量的数字方及***法。
根据本发明的其中一方面,本发明解决其技术问题所采用的技术方案是:构造一种基于全相位FFT的铯光泵磁力仪频率测量方法,包含如下步骤:
S1、将铯光泵磁力仪输出的磁共振信号经放大整形成方波信号,再将方波信号以及方波信号延迟一个固定的时间后分别通过ADC进行采样,得到离散的输入信号序列及输入信号延时序列,完成模数转换过程;
S2、控制FPGA分别对输入信号序列及延时序列分别使用汉宁窗进行全相位预处理;
S3、调用FFT IP核对进行了全相位预处理后的序列分别做FFT变换,输出数据的实部和虚部,然后根据输出的缩放因子对所述数据的实部和虚部进行符号位拓展和移位得到真实的实部和虚部值;
S4、由真实的实部和虚部推算原始数据的幅值,在一个采样周期内选取出最大幅值处的实部和虚部;
S5、通过CORDIC算法根据最大幅值处的实部和虚部推算出对应的相位值,将求出的两个相位值传入单片机进行频率计算,从而得出频率测量结果。
进一步地,在本发明的基于全相位FFT的铯光泵磁力仪频率测量方法中,步骤S2中,所述使用汉宁窗进行全相位预处理时,所采用的N点汉宁窗序列与归一化后的2N-1点汉宁窗由MATLAB产生并转化为mif文件存储在FPGA的单端口ROM模块中,供FPGA调用,N为FFT点数。
进一步地,在本发明的基于全相位FFT的铯光泵磁力仪频率测量方法中,所述使用汉宁窗进行全相位预处理具体包括如下步骤:
S21、对N点汉宁窗自身求卷积,得到2N-1点卷积窗;
S22、将得到的2N-1点卷积窗相加求和;
S23、将2N-1点卷积窗相加求和后进行归一化处理,得到2N-1点汉宁窗;
S24、将得到归一化后的2N-1点汉宁窗依次与输入序列各项相乘,从而进行汉宁加窗,加窗时,第1项乘积与第N+1项乘积相加,第2项乘积与第2N-1项乘积相加,…,第N-1项乘积与第2N-1项乘积相加。
进一步地,在本发明的基于全相位FFT的铯光泵磁力仪频率测量方法中,在步骤S3中,所述调用FFT IP核是指调用Quartus II软件自带的FFT IP核。
进一步地,在本发明的基于全相位FFT的铯光泵磁力仪频率测量方法中,步骤S5中,单片机进行频率计算具体是在接收到两个所述相位值后,依据下述公式计算出最终的频率值:
其中为两个所述相位值的差值,Δt为步骤S2中所述固定的时间。
根据本发明的另一方面,本发明为解决其技术问题,还提供了一种基于全相位FFT的铯光泵磁力仪频率测量***,包含如下模块:
序列获取模块,用于将铯光泵磁力仪输出的磁共振信号经放大整形成方波信号,再将方波信号以及方波信号延迟一个固定的时间后分别通过ADC进行采样,得到离散的输入信号序列及输入信号延时序列,完成模数转换过程;
预处理模块,用于控制FPGA分别对输入信号序列及延时序列分别使用汉宁窗进行全相位预处理;
FFT变换模块,用于调用FFT IP核对进行了全相位预处理后的序列分别做FFT变换,输出数据的实部和虚部,然后根据输出的缩放因子对所述数据的实部和虚部进行符号位拓展和移位得到真实的实部和虚部值;
极值确定模块,用于由真实的实部和虚部推算原始数据的幅值,在一个采样周期内选取出最大幅值处的实部和虚部;
结果计算模块,用于通过CORDIC算法根据最大幅值处的实部和虚部推算出对应的相位值,将求出的两个相位值传入单片机进行频率计算,从而得出频率测量结果。
进一步地,在本发明的基于全相位FFT的铯光泵磁力仪频率测量***中,预处理模块中,所述使用汉宁窗进行全相位预处理时,所采用的N点汉宁窗序列与归一化后的2N-1点汉宁窗由MATLAB产生并转化为mif文件存储在FPGA的单端口ROM模块中,供FPGA调用,N为FFT点数。
进一步地,在本发明的基于全相位FFT的铯光泵磁力仪频率测量***中,所述使用汉宁窗进行全相位预处理具体由下述单元实现:
自身求卷积单元,用于对N点汉宁窗自身求卷积,得到2N-1点卷积窗;
卷积窗求和单元,用于将得到的2N-1点卷积窗相加求和;
归一化处理单元,用于将2N-1点卷积窗相加求和后进行归一化处理,得到2N-1点汉宁窗;
加窗处理单元,用于将得到归一化后的2N-1点汉宁窗依次与输入序列各项相乘,从而进行汉宁加窗,加窗时,第1项乘积与第N+1项乘积相加,第2项乘积与第2N-1项乘积相加,…,第N-1项乘积与第2N-1项乘积相加。
进一步地,在本发明的基于全相位FFT的铯光泵磁力仪频率测量***中,在FFT变换模块中,所述调用FFT IP核是指调用Quartus II软件自带的FFT IP核。
进一步地,在本发明的基于全相位FFT的铯光泵磁力仪频率测量***中,结果计算模块中,单片机进行频率计算具体是在接收到两个所述相位值后,依据下述公式计算出最终的频率值:
其中为两个所述相位值的差值,Δt为步骤S2中所述固定的时间。
实施本发明的基于全相位FFT的铯光泵磁力仪频率测量方法及***,具有以下有益效果:
本发明使用了全相位方法对输入序列进行预处理,可以有效地防止频谱泄漏,其主谱线所对应地相位值即为初始相位;调用Quartus II软件自带的FFT IP核可以极大地节省FPGA的资源,提高效率;采用快速傅里叶变换的数字测频法,使其误差不受闸门时间、频率计数等影响,具有高精度和良好的稳定性。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明的基于全相位FFT的铯光泵磁力仪频率测量方法的总体架构图;
图2是本发明的基于全相位FFT的铯光泵磁力仪频率测量方法一实施例的流程图;
图3是汉宁加窗的流程图;
图4是调用FFT ip核后生成的bsf文件图;
图5是基于全相位FFT的铯光泵磁力仪频率测量***一实施例的原理图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
参考图1以及图2,一种基于全相位FFT的铯光泵磁力仪频率测量方法,包含如下步骤:
S1、将铯光泵磁力仪输出的磁共振信号经放大整形成方波信号,再将方波信号以及方波信号延迟一个固定的时间后分别通过ADC进行采样,得到离散的输入信号序列及输入信号延时序列,完成模数转换过程。
S2、控制FPGA分别对输入信号序列及延时序列分别使用汉宁窗进行全相位预处理。步骤S2中,所述使用汉宁窗进行全相位预处理时,所采用的N点汉宁窗序列与归一化后的2N-1点汉宁窗由MATLAB产生并转化为mif文件存储在FPGA的单端口ROM模块中,供FPGA调用,N为FFT点数。
参考图3,使用汉宁窗进行全相位预处理具体包括如下步骤:
S21、对N点汉宁窗自身求卷积,得到2N-1点卷积窗;
S22、将得到的2N-1点卷积窗相加求和;卷积窗实际上就是一组数据,相加就是把所有数据加起来得到一个总的和值,然后在后续步骤用每个数据除以和值就得到了归一化之后的数据,这些都是本领域常用技术手段,这里不具体展开说明。
S23、将2N-1点卷积窗相加求和后进行归一化处理,得到2N-1点汉宁窗;
S24、将得到归一化后的2N-1点汉宁窗依次与输入序列各项相乘,从而进行汉宁加窗,加窗时,第1项乘积与第N+1项乘积相加,第2项乘积与第2N-1项乘积相加,…,第N-1项乘积与第2N-1项乘积相加。
S3、调用FFT IP核对进行了全相位预处理后的序列分别做FFT变换,输出数据的实部和虚部,然后根据输出的缩放因子对所述数据的实部和虚部进行符号位拓展和移位得到真实的实部和虚部值。在本实施例中,调用FFT IP核是指调用Quartus II软件自带的FFTIP核。
参考图4,其为调用FFT ip核后生成的bsf文件图source_real为输出数据的实部,source_imag为输出数据的虚部,source_exp为数据的缩放因子,FFT ip核的输入输出数据和缩放因子都是有符号数,实际运算的结果是实部、虚部与缩放因子的结合,如果缩放因子为负,实部和虚部要左移相应的位数,如果为正则右移。例如:缩放因子宽度为6,其二进制值为101011,该值为有符号数,表示-21,那么实部和虚部要左移21位才是最终结果。
S4、由真实的实部和虚部推算原始数据的幅值,在一个采样周期内选取出最大幅值处的实部和虚部。S5、通过CORDIC算法根据最大幅值处的实部和虚部推算出对应的相位值,将求出的两个相位值传入单片机进行频率计算,从而得出频率测量结果。
单片机进行频率计算具体是在接收到两个所述相位值后,依据下述公式计算出最终的频率值:
其中为两个所述相位值的差值,Δt为步骤S2中所述固定的时间。
参考图5,在本实施例的基于全相位FFT的铯光泵磁力仪频率测量***中,其包含如下模块:序列获取模块51、预处理模块52、FFT变换模块53、极值确定模块54以及结果计算模块55。
序列获取模块51用于将铯光泵磁力仪输出的磁共振信号经放大整形成方波信号,再将方波信号以及方波信号延迟一个固定的时间后分别通过ADC进行采样,得到离散的输入信号序列及输入信号延时序列,完成模数转换过程。预处理模块52用于控制FPGA分别对输入信号序列及延时序列分别使用汉宁窗进行全相位预处理。FFT变换模块53用于调用FFTIP核对进行了全相位预处理后的序列分别做FFT变换,输出数据的实部和虚部。极值确定模块54用于由输出的实部和虚部推算原始数据的幅值,在一个采样周期内选取出最大幅值处的实部和虚部。结果计算模块55用于通过CORDIC算法根据最大幅值处的实部和虚部推算出对应的相位值,将求出的两个相位值传入单片机进行频率计算,从而得出频率测量结果。
预处理模块52中,所述使用汉宁窗进行全相位预处理时,所采用的N点汉宁窗序列与归一化后的2N-1点汉宁窗由MATLAB产生并转化为mif文件存储在FPGA的单端口ROM模块中,供FPGA调用,N为FFT点数。
所述使用汉宁窗进行全相位预处理具体由下述单元实现:
自身求卷积单元,用于对N点汉宁窗自身求卷积,得到2N-1点卷积窗;
卷积窗求和单元,用于将得到的2N-1点卷积窗相加求和;
归一化处理单元,用于将2N-1点卷积窗相加求和后进行归一化处理,得到2N-1点汉宁窗;
加窗处理单元,用于将得到归一化后的2N-1点汉宁窗依次与输入序列各项相乘,从而进行汉宁加窗,加窗时,第1项乘积与第N+1项乘积相加,第2项乘积与第2N-1项乘积相加,…,第N-1项乘积与第2N-1项乘积相加。
在FFT变换模块53中,所述调用FFT IP核是指调用Quartus II软件自带的FFT IP核。
结果计算模块55中,单片机进行频率计算具体是在接收到两个所述相位值后,依据下述公式计算出最终的频率值:
其中为两个所述相位值的差值,Δt为步骤S2中所述固定的时间。
与目前广泛使用的频率测量方法相比,本发明的优点是:
1、全相位FFT能够有效地抑制频谱泄漏,使其谱线可读性强,不需要任何校正措施即可得到输入序列地初始相位值;
2、采用数字测频法可以在保证连续测量的基础上减小误差,提高测量精度;
3、其结构简单,由于FPGA进行数据处理较困难,故发送至单片机进行最后的频率计算,可以有效地节约资源,提高速率。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。
Claims (10)
1.一种基于全相位FFT的铯光泵磁力仪频率测量方法,其特征在于,包含如下步骤:
S1、将铯光泵磁力仪输出的磁共振信号经放大整形成方波信号,再将方波信号以及方波信号延迟一个固定的时间后分别通过ADC进行采样,得到离散的输入信号序列及输入信号延时序列,完成模数转换过程;
S2、控制FPGA分别对输入信号序列及延时序列分别使用汉宁窗进行全相位预处理;
S3、调用FFTIP核对进行了全相位预处理后的序列分别做FFT变换,输出数据的实部和虚部,然后根据输出的缩放因子对所述数据的实部和虚部进行符号位拓展和移位得到真实的实部和虚部值;
S4、由真实的实部和虚部推算原始数据的幅值,在一个采样周期内选取出最大幅值处的实部和虚部;
S5、通过CORDIC算法根据最大幅值处的实部和虚部推算出对应的相位值,将求出的两个相位值传入单片机进行频率计算,从而得出频率测量结果。
2.根据权利要求1所述基于全相位FFT的铯光泵磁力仪频率测量方法,其特征在于,步骤S2中,所述使用汉宁窗进行全相位预处理时,所采用的N点汉宁窗序列与归一化后的2N-1点汉宁窗由MATLAB产生并转化为mif文件存储在FPGA的单端口ROM模块中,供FPGA调用,N为FFT点数。
3.根据权利要求2所述基于全相位FFT的铯光泵磁力仪频率测量方法,其特征在于,所述使用汉宁窗进行全相位预处理具体包括如下步骤:
S21、对N点汉宁窗自身求卷积,得到2N-1点卷积窗;
S22、将得到的2N-1点卷积窗相加求和;
S23、将2N-1点卷积窗相加求和后进行归一化处理,得到2N-1点汉宁窗;
S24、将得到归一化后的2N-1点汉宁窗依次与输入序列各项相乘,从而进行汉宁加窗,加窗时,第1项乘积与第N+1项乘积相加,第2项乘积与第2N-1项乘积相加,…,第N-1项乘积与第2N-1项乘积相加。
4.根据权利要求1所述基于全相位FFT的铯光泵磁力仪频率测量方法,其特征在于,在步骤S3中,所述调用FFT IP核是指调用Quartus II软件自带的FFTIP核。
5.根据权利要求1所述基于全相位FFT的铯光泵磁力仪频率测量方法,其特征在于,步骤S5中,单片机进行频率计算具体是在接收到两个所述相位值后,依据下述公式计算出最终的频率值:
其中为两个所述相位值的差值,Δt为步骤S2中所述固定的时间。
6.一种基于全相位FFT的铯光泵磁力仪频率测量***,其特征在于,包含如下模块:序列获取模块,用于将铯光泵磁力仪输出的磁共振信号经放大整形成方波信号,再将方波信号以及方波信号延迟一个固定的时间后分别通过ADC进行采样,得到离散的输入信号序列及输入信号延时序列,完成模数转换过程;
预处理模块,用于控制FPGA分别对输入信号序列及延时序列分别使用汉宁窗进行全相位预处理;
FFT变换模块,用于调用FFT IP核对进行了全相位预处理后的序列分别做FFT变换,输出数据的实部和虚部,然后根据输出的缩放因子对所述数据的实部和虚部进行符号位拓展和移位得到真实的实部和虚部值;
极值确定模块,用于由真实的实部和虚部推算原始数据的幅值,在一个采样周期内选取出最大幅值处的实部和虚部;
结果计算模块,用于通过CORDIC算法根据最大幅值处的实部和虚部推算出对应的相位值,将求出的两个相位值传入单片机进行频率计算,从而得出频率测量结果。
7.根据权利要求6所述基于全相位FFT的铯光泵磁力仪频率测量***,其特征在于,预处理模块中,所述使用汉宁窗进行全相位预处理时,所采用的N点汉宁窗序列与归一化后的2N-1点汉宁窗由MATLAB产生并转化为mif文件存储在FPGA的单端口ROM模块中,供FPGA调用,N为FFT点数。
8.根据权利要求7所述基于全相位FFT的铯光泵磁力仪频率测量***,其特征在于,所述使用汉宁窗进行全相位预处理具体由下述单元实现:
自身求卷积单元,用于对N点汉宁窗自身求卷积,得到2N-1点卷积窗;
卷积窗求和单元,用于将得到的2N-1点卷积窗相加求和;
归一化处理单元,用于将2N-1点卷积窗相加求和后进行归一化处理,得到2N-1点汉宁窗;
加窗处理单元,用于将得到归一化后的2N-1点汉宁窗依次与输入序列各项相乘,从而进行汉宁加窗,加窗时,第1项乘积与第N+1项乘积相加,第2项乘积与第2N-1项乘积相加,…,第N-1项乘积与第2N-1项乘积相加。
9.根据权利要求6所述基于全相位FFT的铯光泵磁力仪频率测量***,其特征在于,在FFT变换模块中,所述调用FFT IP核是指调用Quartus II软件自带的FFTIP核。
10.根据权利要求6所述基于全相位FFT的铯光泵磁力仪频率测量***,其特征在于,结果计算模块中,单片机进行频率计算具体是在接收到两个所述相位值后,依据下述公式计算出最终的频率值:
其中为两个所述相位值的差值,Δt为步骤S2中所述固定的时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910606439.1A CN110297199B (zh) | 2019-07-05 | 2019-07-05 | 一种基于全相位fft的铯光泵磁力仪频率测量方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910606439.1A CN110297199B (zh) | 2019-07-05 | 2019-07-05 | 一种基于全相位fft的铯光泵磁力仪频率测量方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110297199A true CN110297199A (zh) | 2019-10-01 |
CN110297199B CN110297199B (zh) | 2021-08-27 |
Family
ID=68030516
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910606439.1A Active CN110297199B (zh) | 2019-07-05 | 2019-07-05 | 一种基于全相位fft的铯光泵磁力仪频率测量方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110297199B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111597498A (zh) * | 2020-05-18 | 2020-08-28 | 电子科技大学 | 一种基于大点数fft电路的频谱获取方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6876186B1 (en) * | 2003-08-29 | 2005-04-05 | Xilinx, Inc. | Measurement of circuit delay |
CN104375111A (zh) * | 2014-11-16 | 2015-02-25 | 甘肃省机械科学研究院 | 对密集频谱进行快速高精度细化校正的方法 |
CN106226590A (zh) * | 2016-07-19 | 2016-12-14 | 国网河北省电力公司电力科学研究院 | 一种电力***同步相量测量方法 |
CN106814249A (zh) * | 2017-02-28 | 2017-06-09 | 中国科学院国家空间科学中心 | 一种用于光泵磁力仪频率计的频率测量方法 |
CN106885998A (zh) * | 2017-03-08 | 2017-06-23 | 中国地质大学(武汉) | 一种提高铯光泵磁共振信号测频精度和速度的方法和电路 |
CN206710575U (zh) * | 2017-03-08 | 2017-12-05 | 中国地质大学(武汉) | 一种提高铯光泵磁共振信号测频精度和速度的电路 |
-
2019
- 2019-07-05 CN CN201910606439.1A patent/CN110297199B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6876186B1 (en) * | 2003-08-29 | 2005-04-05 | Xilinx, Inc. | Measurement of circuit delay |
CN104375111A (zh) * | 2014-11-16 | 2015-02-25 | 甘肃省机械科学研究院 | 对密集频谱进行快速高精度细化校正的方法 |
CN106226590A (zh) * | 2016-07-19 | 2016-12-14 | 国网河北省电力公司电力科学研究院 | 一种电力***同步相量测量方法 |
CN106814249A (zh) * | 2017-02-28 | 2017-06-09 | 中国科学院国家空间科学中心 | 一种用于光泵磁力仪频率计的频率测量方法 |
CN106885998A (zh) * | 2017-03-08 | 2017-06-23 | 中国地质大学(武汉) | 一种提高铯光泵磁共振信号测频精度和速度的方法和电路 |
CN206710575U (zh) * | 2017-03-08 | 2017-12-05 | 中国地质大学(武汉) | 一种提高铯光泵磁共振信号测频精度和速度的电路 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111597498A (zh) * | 2020-05-18 | 2020-08-28 | 电子科技大学 | 一种基于大点数fft电路的频谱获取方法 |
CN111597498B (zh) * | 2020-05-18 | 2023-04-18 | 电子科技大学 | 一种基于大点数fft电路的频谱获取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110297199B (zh) | 2021-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
US3881097A (en) | Fully digital spectrum analyzer using time compression and discrete fourier transform techniques | |
WO2018188228A1 (zh) | 高精度频率测量***及方法 | |
CN103257271B (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN104007316B (zh) | 一种欠采样速率下的高精度频率测量方法及其测量仪 | |
CN103197141A (zh) | 一种测量电力***信号频率及谐波参数的方法 | |
CN103207319A (zh) | 数字化变电站电力信号非同步采样条件下的谐波测量方法 | |
CN109444515B (zh) | 一种基于sdft算法的无功、不平衡与谐波检测方法 | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
Sapozhnikov et al. | Advantages of using the probabilistic form of information representation in information-control systems | |
CN103983957B (zh) | 一种多普勒偏移测量方法及其装置 | |
CN102998523A (zh) | 一种用于电能计量的谐波功率计算方法 | |
CN105911341B (zh) | 一种谐波无功功率的测量方法 | |
CN108918965A (zh) | 多通道信号相位、幅度高精度测量方法 | |
CN106018960B (zh) | 一种基于压缩传感的同步相量测量方法 | |
CN106546949A (zh) | 一种基于频率预估计的双阵元正弦信号来波方向估计方法 | |
CN102680945A (zh) | 基于fpga的多普勒调频率估计方法 | |
CN110068729B (zh) | 一种信号相量计算方法 | |
CN110297199A (zh) | 一种基于全相位fft的铯光泵磁力仪频率测量方法及*** | |
CN106645780A (zh) | 一种基于dsp的转速检测方法及*** | |
CN104333348B (zh) | 一种高阶数字滤波***及方法 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN104266732B (zh) | 双台位流量计检定装置脉冲计数精度补偿方法 | |
CN105182072A (zh) | 一种基于虚拟滤波因子的谐波检测方法与*** | |
CN102520246B (zh) | 定频率相量提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |