CN112346381A - 一种基于fpga的数字相敏检波方法及系统 - Google Patents
一种基于fpga的数字相敏检波方法及系统 Download PDFInfo
- Publication number
- CN112346381A CN112346381A CN202011181903.6A CN202011181903A CN112346381A CN 112346381 A CN112346381 A CN 112346381A CN 202011181903 A CN202011181903 A CN 202011181903A CN 112346381 A CN112346381 A CN 112346381A
- Authority
- CN
- China
- Prior art keywords
- signal
- square wave
- detected
- frequency
- wave signal
- 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
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
- G05B19/042—Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
- G05B19/0423—Input/output
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/21—Pc I-O input output
- G05B2219/21126—Signal processing, filter input
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Measuring Phase Differences (AREA)
Abstract
本发明涉及一种基于FPGA的数字相敏检波方法及系统,属于信号处理技术领域。本发明通过FPGA生成与待检测信号同频的正弦方波信号和余弦方波信号,将正弦方波信号和余弦方波信号作为参考信号,与待检测信号分别做相关运算,将两路相关运算的结果进行求模处理,并将求模结果乘以一个系数,以此作为待检测信号的幅值,从而实现从待检测信号提取出有用信号的幅值。本发明将方波作为检波时的参考信号,与正弦波作为参考信号相比,不仅简化了控制,还提高了计算效率,且操作简单、容易实现。
Description
技术领域
本发明涉及一种基于FPGA的数字相敏检波方法及系统,属于信号处理技术领域。
背景技术
在工程测量中,被测信号往往会被有源器件产生的高斯噪声、谐波噪声以及外界环境的随机噪声等各类噪声干扰或淹没,数字相敏检波算法能够很好地从强噪声中提取出有用信号的幅值、相位信息。因此,该方法被广泛应用于现有的信号检测系统。数字相敏检波主要采用高速ADC采集模拟正弦信号,与参考的数字正弦信号做相关,以数字方式实现相敏检波。数字相敏检波算法对谐波信号与有源器件引起的随机噪声具有很强的抑制能力,更高的精度,和更大的灵活性等优点,更适合微弱信号的检测。
在现有的数字相敏检波中往往利用与被测信号同频的正弦信号作为参考,通过对参考信号与被测信号进行信号处理,获得被测信号的幅度与相位信息。在FPGA实现中,正弦信号需要利用DA模块产生,操作复杂且不易控制,使得检波效率低。
发明内容
本发明的目的是提供一种基于FPGA的数字相敏检波方法及系统,以解决目前基于FPGA数字相敏检波存在的控制复杂、效率低的问题。
本发明为解决上述技术问题而提供一种基于FPGA的数字相敏检波方法,该方法包括以下步骤:
1)利用FPGA产生一个用于待检测信号同频的正弦方波信号和余弦方波信号;
2)将生成的正弦方波信号和余弦方波信号分别与待检测信号进行相关运算;
3)将相关运算结果进行求模处理,并将求模处理结果乘以频率相关系数,该结果为待检测信号的幅度。
本发明通过FPGA生成与待检测信号同频的正弦方波信号和余弦方波信号,将正弦方波信号和余弦方波信号作为参考信号,与待检测信号分别做相关运算,将两路相关运算的结果进行求模处理,并将求模结果乘以一个系数,以此作为待检测信号的幅值,从而实现从待检测信号提取出有用信号的幅值。本发明将方波作为检波时的参考信号,与正弦波作为参考信号相比,不仅简化了控制,还提高了计算效率,且操作简单、容易实现。
进一步地,所述的频率相关系数与待检测信号的频率相关,一个确定的频率对应一个确定的频率相关系数。
进一步地,为了准确得到频率相关系数,所述的频率相关系数通过实验仿真得到,仿真过程如下:
A.生成一个与待测信号频率相同的基准方波信号,其中待测信号为不同频率的已知幅值的标准信号;
B.将待测信号与基准方波信号作互相关计算结果,根据互相关结果确定方波信号提取的幅度,将方波信号提取的幅度与待测信号的幅值进行比较得到频率相关系数。
进一步地,所述的正弦方波信号和余弦方波信号分别为:square(2πfn/N+π/2)和square2πfn/N,
其中f为频率,N为一个周期的采样点数。
本发明还提供了一种基于FPGA的数字相敏检波系统,该检波系统采用FPGA,包括方波生成模块、相关运算模块和求模处理模块;
所述方波生成模块用于产生一个用于待检测信号同频的正弦方波信号和余弦方波信号;
所述相关运算模块用于将方波生成模块生成的正弦方波信号和余弦方波信号分别与待检测信号进行相关运算;
所述求模处理模块用于将相关运算模块的运算结果进行求模处理,并将求模处理结果乘以频率相关系数,该结果为待检测信号的幅度。
本发明通过FPGA生成与待检测信号同频的正弦方波信号和余弦方波信号,将正弦方波信号和余弦方波信号作为参考信号,与待检测信号分别做相关运算,将两路相关运算的结果进行求模处理,并将求模结果乘以一个系数,以此作为待检测信号的幅值,从而实现从待检测信号提取出有用信号的幅值。本发明将方波作为检波时的参考信号,与正弦波作为参考信号相比,不仅简化了控制,还提高了计算效率,且操作简单、容易实现。
进一步地,所述的频率相关系数与待检测信号的频率相关,一个确定的频率对应一个确定的频率相关系数。
进一步地,为了准确得到频率相关系数,所述的频率相关系数通过试实验仿真得到,仿真过程如下:
A.生成一个与待测信号频率相同的基准方波信号,其中待测信号为不同频率的已知幅值的标准信号;
B.将待测信号与基准方波信号作互相关计算结果,根据互相关结果确定方波信号提取的幅度,将方波信号提取的幅度与待测信号的幅值进行比较得到频率相关系数。
所述的正弦方波信号和余弦方波信号分别为:square(2πfn/N+π/2)和square2πfn/N,
其中f为频率,N为一个周期的采样点数。
附图说明
图1是本发明基于FPGA的数字相敏检波方法的原理示意图;
图2是本发明基于FPGA的数字相敏检波方法的流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步地说明。
方法实施例
本发明的相敏检波方法首先通过FPGA生成与待检测信号同频的正弦方波信号和余弦方波信号,然后将正弦方波信号和余弦方波信号作为参考信号,与待检测信号分别做相关运算,最后将两路相关运算的结果进行求模处理,并将求模结果乘以一个系数,以此作为待检测信号的幅值,从而实现从待检测信号提取出有用信号的幅值。该方法的实现原理如图1所示,流程如图2所示,具体实现过程如下。
1.利用FPGA生成方波信号。
目前的基于FPGA的相敏检波一般都是生成正弦信号作为参考信号,而利用FPGA生成正弦信号往往需要利用DA模块产生,DA模块操作复杂,不易控制。为此,本发明采用方波信号作为参考信号。与正弦波相比,FPGA在生成方波信号时只需要通过高电平控制即可,无需DA模块,大大简化了电路。本发明需要的方波信号与待检测信号的频率一致,由于方波信号可以看成是对正弦信或余弦信号整流得到。为此,本发明生成两路方波信号,一路是正弦方波信号,一路是余弦方波信号。
若待测信号为x(n)=A sin(2πn/N+ψ)+v(n),其中,A为待测信号幅度,N=f/fs(f为信号频率,fs为采样频率),ψ为相位,v(n)表示待测信号中携带的全部噪声。则生成的两路方波信号,信号分别为square(2πfn/N+π/2)和square2πfn/N。
2.进行相关处理。
经典的数字相敏检波原理推导过程:
将待测信号x(n)与正弦参考信号s(n)=cos(wn)+j sin(wn)进行互相关运算,则运算结果为:
方波与正弦波做相关处理在数学上不易计算,可以把方波进行傅里叶级数展开:方波可以写成多次谐波和的形式,而噪声与方波不相关,为此,本发明将生成的两个方波信号分别与待检测信号做相关运算,则可以消除待测信号中的噪声影响,得到只与有用信号有关的参数I和Q。
其中,RxΥ(0)对相关结果产生一定影响,但是由于高次谐波与基波做相关的复杂性,难以从理论上得出确定的影响关系,因此本发明需要结合matlab仿真,通过仿真结果确定一个比较合适的校正系数,从而方便实验应用。
3.求模处理。
由于高次谐波(例如三次谐波)与基波并非完全不相关的,即sin(wn)与sin((2i+1)wn)(i≥1)做相关不为零,本发明需要通过仿真计算一个系数k进行结果矫正。因此,本发明将得到的有关参数进行求模处理时,需将求模处理结果乘以一个频率相关系数,所得结果才为待检测信号的幅值。对本实施例而言,若相关运算得到的相关参数分别为I和Q,则求模结果为乘以一个频率相关系数k,最终得到的待检测信号的幅值A为其中频率相关系数k有待检测信号x(n)的频率有关,某一个确定频率对应着一个确定系数k。在实际操作过程中,可以先进行仿真,对仿真结果进行数值分析,以获得实验频率所对应的系数k。仿真采用MATLAB实现,过程如下:
选择一个标准信号y0=2sin(40πn+π/4)+v(n)作为待测信号,其中噪声信号v(n)为高斯白噪声。
生成一个基准方波信号,为了方便和现有技术比较,再生成一个基准正弦信号,其中基准方波信号和基准正弦信号分别为:
x1=square(40πn+π/2,50)+j*square(40πn,50)
x2=cos(40πn+π/2)+j*sin(40πn+π/2)
当信噪比为50db时,其仿真结果如表1所示表1
当信噪比从20db到80db变化时,其仿真结果如表2所示
表2
从以上两个表格可以看出,虽然方波提取出的信号幅度误差较大,但和正弦波一样,方波所得到的误差值比较稳定,由此可以推测,方波得到的幅度与正弦波的幅度只是相差一个系数s,通过s值可以校正方波所得误差。
根据实际操作中通常情况下产生的信噪比大小选择相应的仿真结果来确定k值,比如选取信噪比为50db时的数值,方波提取幅度结果大多数为2.75,标准信号幅度为2,则s=2.75/2=1.375。
在不好判断信噪比时,可以用取平均的方法确定k值,即对信噪比为30db,40db,50db,60db时方波提取幅度结果求平均(也可以采用多种信噪比条件下的结果,这里仅为举例,没有限定条件),再与标准信号幅度求商:s=[(2.7516+2.7468+2.7503+2.7501)/4]/2≈1.375。
当信号频率f从2Hz到100Hz变化,采样频率fs取10*f,采样点数保持200不变,信号幅度为1。其仿真结果如表3所示
表3
由表3可知,当信号频率变化时,方波提取到的幅度有较大的变化,而正弦波所得结果基本不变,因此做出判断,校正系数k与频率有关。
系统实施例
本发明基于FPGA的数字相敏检波系统采用FPGA实现,包括方波生成模块、相关运算模块和求模处理模块;方波生成模块用于产生一个用于待检测信号同频的正弦方波信号和余弦方波信号;关运算模块用于将方波生成模块生成的正弦方波信号和余弦方波信号分别与待检测信号进行相关运算;求模处理模块用于将相关运算模块的运算结果进行求模处理,并将求模处理结果乘以频率相关系数,该结果为待检测信号的幅度。上述各模块的具体实现手段已在方法的实施例中进行了详细说明,这里不再赘述。
因此,本发明与现有采用正弦波作为参考信号的相敏检波方法相比,能够大大简化实现电路,提高检波效率,且便于控制调整,实现简单。
Claims (8)
1.一种基于FPGA的数字相敏检波方法,其特征在于,该方法包括以下步骤:
1)利用FPGA产生一个用于待检测信号同频的正弦方波信号和余弦方波信号;
2)将生成的正弦方波信号和余弦方波信号分别与待检测信号进行相关运算;
3)将相关运算结果进行求模处理,并将求模处理结果乘以频率相关系数,该结果为待检测信号的幅度。
2.根据权利要求1所述的基于FPGA的数字相敏检波方法,其特征在于,所述的频率相关系数与待检测信号的频率相关,一个确定的频率对应一个确定的频率相关系数。
3.根据权利要求1或2所述的基于FPGA的数字相敏检波方法,其特征在于,所述的频率相关系数通过实验仿真得到,仿真过程如下:
A.生成一个与待测信号频率相同的基准方波信号,其中待测信号为不同频率的已知幅值的标准信号;
B.将待测信号与基准方波信号作互相关计算结果,根据互相关结果确定方波信号提取的幅度,将方波信号提取的幅度与待测信号的幅值进行比较得到频率相关系数。
4.根据权利要求1所述的基于PGA的数字相敏检波方法,其特征在于,所述的正弦方波信号和余弦方波信号分别为:square(2πfn/N+π/2)和square2πfn/N,其中f为频率,N为一个周期的采样点数。
5.一种基于FPGA的数字相敏检波系统,其特征在于,该检波系统采用FPGA,包括方波生成模块、相关运算模块和求模处理模块;
所述方波生成模块用于产生一个用于待检测信号同频的正弦方波信号和余弦方波信号;
所述相关运算模块用于将方波生成模块生成的正弦方波信号和余弦方波信号分别与待检测信号进行相关运算;
所述求模处理模块用于将相关运算模块的运算结果进行求模处理,并将求模处理结果乘以频率相关系数,该结果为待检测信号的幅度。
6.根据权利要求5所述的基于FPGA的数字相敏检波系统,其特征在于,所述的频率相关系数与待检测信号的频率相关,一个确定的频率对应一个确定的频率相关系数。
7.根据权利要求5或6所述的基于FPGA的数字相敏检波系统,其特征在于,所述的频率相关系数通过试实验仿真得到,仿真过程如下:
A.生成一个与待测信号频率相同的基准方波信号,其中待测信号为不同频率的已知幅值的标准信号;
B.将待测信号与基准方波信号作互相关计算结果,根据互相关结果确定方波信号提取的幅度,将方波信号提取的幅度与待测信号的幅值进行比较得到频率相关系数。
8.根据权利要求5所述的基于FPGA的数字相敏检波系统,其特征在于,所述的正弦方波信号和余弦方波信号分别为:square(2πfn/N+π/2)和square2πfn/N,其中f为频率,N为一个周期的采样点数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011181903.6A CN112346381B (zh) | 2020-10-29 | 2020-10-29 | 一种基于fpga的数字相敏检波方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011181903.6A CN112346381B (zh) | 2020-10-29 | 2020-10-29 | 一种基于fpga的数字相敏检波方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112346381A true CN112346381A (zh) | 2021-02-09 |
CN112346381B CN112346381B (zh) | 2022-07-08 |
Family
ID=74356529
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011181903.6A Active CN112346381B (zh) | 2020-10-29 | 2020-10-29 | 一种基于fpga的数字相敏检波方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112346381B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114422039A (zh) * | 2022-01-21 | 2022-04-29 | 中车大连电力牵引研发中心有限公司 | 一种能够去除信号中噪声的方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3944940A (en) * | 1974-09-06 | 1976-03-16 | Pertec Corporation | Versatile phase-locked loop for read data recovery |
CN102403969A (zh) * | 2011-11-03 | 2012-04-04 | 北京交通大学 | 一种锁相放大器及用于该锁相放大器的信号解调方法 |
CN102723921A (zh) * | 2012-04-28 | 2012-10-10 | 深圳职业技术学院 | 基于现场可编程门阵列的数字锁相放大实现方法及系统 |
KR101203041B1 (ko) * | 2011-10-24 | 2012-11-21 | 경희대학교 산학협력단 | 위상차를 이용한 정현파 진폭 측정 방법 및 그 장치 |
CN103840795A (zh) * | 2014-02-26 | 2014-06-04 | 中国科学院自动化研究所 | 一种基于dds芯片移相的正交检波电路 |
CN204788657U (zh) * | 2015-07-08 | 2015-11-18 | 北京光电技术研究所 | 基于红外辐射测量的激光功率计 |
CN105157830A (zh) * | 2015-07-08 | 2015-12-16 | 北京光电技术研究所 | 基于红外辐射测量的激光功率计 |
CN105572606A (zh) * | 2016-01-26 | 2016-05-11 | 上海交通大学 | 基于fpga的磁通门微小信号检测系统及方法 |
CN106374918A (zh) * | 2016-08-29 | 2017-02-01 | 北京航空航天大学 | 一种基于dsp和fpga的高速数字锁相放大器 |
CN107479089A (zh) * | 2017-09-08 | 2017-12-15 | 中国舰船研究设计中心 | 一种新型功率检波器 |
CN110971231A (zh) * | 2019-12-05 | 2020-04-07 | 上海交通大学 | 结合比率运算的数字锁相放大器 |
-
2020
- 2020-10-29 CN CN202011181903.6A patent/CN112346381B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3944940A (en) * | 1974-09-06 | 1976-03-16 | Pertec Corporation | Versatile phase-locked loop for read data recovery |
KR101203041B1 (ko) * | 2011-10-24 | 2012-11-21 | 경희대학교 산학협력단 | 위상차를 이용한 정현파 진폭 측정 방법 및 그 장치 |
CN102403969A (zh) * | 2011-11-03 | 2012-04-04 | 北京交通大学 | 一种锁相放大器及用于该锁相放大器的信号解调方法 |
CN102723921A (zh) * | 2012-04-28 | 2012-10-10 | 深圳职业技术学院 | 基于现场可编程门阵列的数字锁相放大实现方法及系统 |
CN103840795A (zh) * | 2014-02-26 | 2014-06-04 | 中国科学院自动化研究所 | 一种基于dds芯片移相的正交检波电路 |
CN204788657U (zh) * | 2015-07-08 | 2015-11-18 | 北京光电技术研究所 | 基于红外辐射测量的激光功率计 |
CN105157830A (zh) * | 2015-07-08 | 2015-12-16 | 北京光电技术研究所 | 基于红外辐射测量的激光功率计 |
CN105572606A (zh) * | 2016-01-26 | 2016-05-11 | 上海交通大学 | 基于fpga的磁通门微小信号检测系统及方法 |
CN106374918A (zh) * | 2016-08-29 | 2017-02-01 | 北京航空航天大学 | 一种基于dsp和fpga的高速数字锁相放大器 |
CN107479089A (zh) * | 2017-09-08 | 2017-12-15 | 中国舰船研究设计中心 | 一种新型功率检波器 |
CN110971231A (zh) * | 2019-12-05 | 2020-04-07 | 上海交通大学 | 结合比率运算的数字锁相放大器 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114422039A (zh) * | 2022-01-21 | 2022-04-29 | 中车大连电力牵引研发中心有限公司 | 一种能够去除信号中噪声的方法 |
CN114422039B (zh) * | 2022-01-21 | 2024-03-19 | 中车大连电力牵引研发中心有限公司 | 一种能够去除信号中噪声的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112346381B (zh) | 2022-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Belega et al. | Frequency estimation by two-or three-point interpolated Fourier algorithms based on cosine windows | |
CN102045036B (zh) | 一种数字锁相放大器 | |
CN112346381B (zh) | 一种基于fpga的数字相敏检波方法及系统 | |
Belega et al. | A high-performance procedure for effective number of bits estimation in analog-to-digital converters | |
CN107796977B (zh) | 一种三相电网电压参数检测方法及装置 | |
CN104215833B (zh) | 电力系统频率测量方法及装置 | |
CN101718816B (zh) | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN109541304A (zh) | 基于六项最小旁瓣窗插值的电网高次弱幅值谐波检测方法 | |
CN104965123A (zh) | 一种基于混沌振子的电力系统谐波、间谐波检测与估计新方法 | |
Wu et al. | Frequency estimation algorithm for ranging of millimeter wave LFMCW radar | |
Shao et al. | A frequency correcting method combining bilateral correction with weighted average for vortex flow sensor signal | |
Qin et al. | The order tracking of rotating machinery based on instantaneous frequency estimation | |
CN106768260A (zh) | 可抑制直流扰动的振动信号最大功率频率分量实时检测算法 | |
CN104407212B (zh) | 一种实时高精度相参脉冲信号测频方法 | |
CN107576842B (zh) | 一种宽频带同步采样方法 | |
Wen et al. | Comparative study of influence of noise on power frequency estimation of sine wave using interpolation FFT | |
Zhang et al. | Research on the signal separation method based on multi-sensor cross-correlation fusion algorithm | |
Lv et al. | Adaptive algorithm based on FFT for frequency estimation | |
Asnin et al. | DSP methods for dynamic estimation of frequency and magnitude parameters in power system transients | |
Belega et al. | Efficiency of the three-point interpolated DFT method on the normalized frequency estimation of a sine-wave | |
Shao et al. | Comparative of Influence of Noise on Frequency Estimation Provided by Three Points Interpolated Discrete Fourier Transform | |
CN112782451B (zh) | 一种基于时域的相位分析方法、装置及系统 | |
Lu et al. | FPGA Based Implementation of All-phase FFT Phase Difference Frequency Measurement | |
Belega et al. | Choice of the acquisition parameters for frequency estimation of a sine wave by interpolated DFT method | |
Belega et al. | Accurate sine-wave frequency estimation by means of an interpolated DTFT algorithm |
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 |