CN108426942B - 数字锁相解调的滤波器的实现方法及装置 - Google Patents

数字锁相解调的滤波器的实现方法及装置 Download PDF

Info

Publication number
CN108426942B
CN108426942B CN201810113634.6A CN201810113634A CN108426942B CN 108426942 B CN108426942 B CN 108426942B CN 201810113634 A CN201810113634 A CN 201810113634A CN 108426942 B CN108426942 B CN 108426942B
Authority
CN
China
Prior art keywords
channel
digital
signal
frequency
analog
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.)
Active
Application number
CN201810113634.6A
Other languages
English (en)
Other versions
CN108426942A (zh
Inventor
李勇
刘泽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Jiaotong University
Original Assignee
Beijing Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Jiaotong University filed Critical Beijing Jiaotong University
Priority to CN201810113634.6A priority Critical patent/CN108426942B/zh
Publication of CN108426942A publication Critical patent/CN108426942A/zh
Application granted granted Critical
Publication of CN108426942B publication Critical patent/CN108426942B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/72Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
    • G01N27/82Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B61RAILWAYS
    • B61KAUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
    • B61K9/00Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
    • B61K9/08Measuring installations for surveying permanent way
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03LAUTOMATIC CONTROL, STARTING, SYNCHRONISATION, OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
    • H03L7/00Automatic control of frequency or phase; Synchronisation
    • H03L7/06Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
    • H03L7/08Details of the phase-locked loop
    • H03L7/085Details of the phase-locked loop concerning mainly the frequency- or phase-detection arrangement including the filtering or amplification of its output signal
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0202Two or more dimensional filters; Filters for complex signals
    • H03H2017/0205Kalman filters

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mechanical Engineering (AREA)
  • Biochemistry (AREA)
  • Electrochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)

Abstract

本发明提供了一种数字锁相解调的滤波器的实现方法及装置。该方法包括:在数字锁相解调的滤波器的两个通道分别输入数字信号,在每个通道对每一时刻输入的数字信号分别与该通道的正弦参考信号或者余弦参考信号相乘,获得每个通道的频移结果,在每个通道通过卡尔曼滤波器对一个周期点数的频移结果进行累加求平均,将计算出的平均值作为该通道当前时刻的观测值;通过卡尔曼滤波器利用每个通道当前时刻的观测值对该通道的当前时刻的预测值进行校正更新,再利用每个通道更新后的预测值求出所述数字信号的幅值和相位。本发明提供的数字锁相解调的滤波器改进设计方法既保证了解调的精度高又保证了解调的快速性,且稳定性高。

Description

数字锁相解调的滤波器的实现方法及装置
技术领域
本发明涉及数字锁相解调技术领域,特别涉及一种数字锁相解调的滤波器的实现方法及装置。
背景技术
为适应高速钢轨缺陷的在线检测,需提高解调的速度以快速解调出高频变化的感应电压。但是在实际检测中,检测信号往往非常微弱且受到模拟电路中的白噪声以及钢轨牵引回流的干扰。电磁检测缺陷的准确度受到强噪声的影响而下降,为此需要从噪声中提取有效信号以提高缺陷的识别能力。
数字锁相解调是一种对交变信号进行相敏检波的方式。数字锁相解调方法利用与被测信号有相同频率和相位关系的参考信号作为比较基准,只对被测信号本身以及与参考信号同频(或者倍频)、同相的噪声分量有响应,具有大幅度提高信噪比的优点。数字锁相解调方法的稳定性和精度很高,但在速度上,由于传统的低通滤波器的计算量大,解调速度慢。
现有技术中的数字锁相解调的低通滤波器设计主要分为两类:IIR滤波器和FIR滤波器。在对现有设计分析后,发明人发现现有的数字锁相解调的低通滤波器设计至少存在如下不足:虽然使用IIR滤波器解调的精度很高,但是收敛的速度很慢,以致影响最终解调的精度或速度;而使用FIR滤波器速度快,但是容易受到噪声和谐波的干扰,抗扰动能力差。
因此,目前迫切需要一种精度高、检测速度快、成本低的在线钢轨缺陷检测装置。
发明内容
本发明的实施例提供了一种数字锁相解调的滤波器的实现方法及装置,精度高、检测速度快、成本低的在线钢轨缺陷检测装置。
为了实现上述目的,本发明采取了如下技术方案。
根据本发明的一个方面,提供了一种数字锁相解调的滤波器的实现方法,包括:
在数字锁相解调的滤波器的两个通道分别输入数字信号,所述两个通道包括通道1和通道2;
在每个通道对每一时刻输入的数字信号分别与该通道的正弦参考信号或者余弦参考信号相乘,获得每个通道的频移结果,所述正弦参考信号与所述余弦参考信号正交;
在每个通道通过卡尔曼滤波器对一个周期点数的频移结果进行累加求平均,将计算出的平均值作为该通道当前时刻的观测值;
通过卡尔曼滤波器利用每个通道当前时刻的观测值对该通道的当前时刻的预测值进行校正更新,再利用每个通道更新后的预测值求出所述数字信号的幅值和相位。
进一步地,所述的在数字锁相解调的滤波器的两个通道分别输入数字信号之前,还包括:
根据输入的数字信号的激励频率和采样频率,计算出相互正交的正弦参考信号和余弦参考信号,所述通道1和所述通道2分别使用所述正弦参考信号和所述余弦参考信号中的一种信号,将N个相互正交的正弦参考信号和余弦参考信号预先存储在内存当中,所述N是以采样频率fs采样目标信号一个周期所需要的点数,所述参考信号与所述数字信号的激励信号的频率相同。
进一步地,所述的在数字锁相解调的滤波器的两个通道分别输入数字信号之前,还包括:
模数转换器以采样频率fs将接收到的当前时刻的模拟差分信号转换成当前时刻的数字信号S[k],当前时刻的数字信号S[k]表示为:
S[k]=A0sin[Ω0k+θ0]
A0为数字信号中频率与激励频率f0相同的幅值,θ0为数字信号的相位,Ω0为数字信号的角频率,Ω0=2π/N,其中N是以采样频率为fs对频率为f0的差分信号采集一个周期所需要的点数,N=fs/f0
进一步地,所述的在每个通道对每一时刻输入的数字信号分别与该通道的正弦参考信号或者余弦参考信号相乘,获得每个通道的频移结果,包括:
在通道1和通道2,分别对每一时刻输入的数字信号通过乘法器进行频移,当前时刻输入的数字信号的频移结果表示如下:
SPSD=S[k]·Sr[k]=(R[k] I[k])T
R[k]和I[k]分别是当前时刻通道1和通道2的频移结果。
进一步地,所述的在每个通道通过卡尔曼滤波器对一个周期点数的频移结果进行累加求平均,将计算出的平均值作为该通道当前时刻的观测值,包括:
在每个通道利用卡尔曼滤波器对前N个时刻的频移结果进行累加后求平均,将计算出的平均值作为该通道当前时刻的观测值Y[k],所述卡尔曼滤波器的观测方程表示如下:
Figure BDA0001570022850000031
式中,
SPSD[k]=(R[k] I[k])T
R[k]和I[k]分别是通道1和通道2第k时刻的频移结果,V[k]是当前时刻的观测噪声向量。
进一步地,所述的通过卡尔曼滤波器利用每个通道当前时刻的观测值对该通道的预测值进行校正更新,再利用每个通道更新后的预测值求出所述输入的数字信号的幅值和相位,包括:
在每个通道,利用当前时刻的观测值Y[k]来校正当前时刻的卡尔曼滤波器的预测值X[k|k-1],得到该通道当前时刻的最优估计值,通道1和通道2当前时刻的最优估计值分别为R′[k]和I′[k];
每个通道当前时刻的卡尔曼滤波器的预测值X[k|k-1]等于上一时刻该通道的最优估计值,卡尔曼滤波器的预测方程的表示如下:
X[k|k-1]=(R′[k-1] I′[k-1])T
R′[k-1]和I′[k-1]分别是上一时刻通道1和通道2的最优估计值;
所述输入的数字信号的幅值和相位的幅值A0和相位θ0的计算公式如下:
Figure BDA0001570022850000041
Figure BDA0001570022850000042
根据本发明的另一个方面,提供了一种基于所述的数字锁相解调的滤波器的实现方法的在线钢轨缺陷检测的装置,包括:依次连接的数模转换器、功率放大器、线圈传感器、差分放大器、模拟带通滤波器、模数转换器和单片机;
所述数模装换器,用于将数字信号转换成激励模拟信号,将所述激励模拟信号传输给功率放大器;
所述功率放大器,用于将所述激励模拟信号放大,将放大后的激励模拟信号传输给所述线圈传感器的激励线圈;
所述线圈传感器,用于包括中间线圈和前后线圈,所述中间线圈为激励线圈,所述前后线圈为检测线圈,所述前后线圈利用所述中间线圈激发的交变磁场感应检测钢轨缺陷,所述前后线圈产生的感应信号传输给差分放大器;
差分放大器,用于将所述两个检测线圈的感应信号进行差分放大,产生模拟的差分信号,将模拟的差分信号传输给模拟带通滤波器;
模拟带通滤波器,用于滤除所述模拟的差分信号的低频和高频干扰,将处理后的差分信号传输给所述模数转换器;
模数转换器,用于将接收到的差分信号转换成当前时刻的数字信号,将当前时刻的数字信号传输给单片机;
单片机,用于对接收到的数字信号进行基于卡尔曼滤波的数字锁相解调处理,得到数字信号的幅值和相位。
进一步地,所述的装置还包括:
WiFi模块,用于接入互联网络,与所述单片机电路连接,传输单片机解调得到的幅值和相位至远程服务器,以方便用户观看。
进一步地,所述的装置还包括:
显示模块,用于与所述单片机电路连接,显示单片机解调得到的幅值和相位。
由上述本发明的实施例提供的技术方案可以看出,本发明提供的数字锁相解调的滤波器改进设计方法既保证了解调的精度高又保证了解调的快速性,且稳定性高。本发明提供的在线钢轨缺陷检测装置,成本低,操作简单,测量精度高。
本发明附加的方面和优点将在下面的描述中部分给出,这些将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种基于数字锁相解调的滤波器的在线钢轨缺陷检测装置的结构示意图;
图2为本发明实施例提供的一种线圈传感器检测钢轨缺陷示意图;
图3为本发明实施例提供的一种基于卡尔曼滤波数字锁相解调的流程图。
具体实施方式
下面详细描述本发明的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本发明,而不能解释为对本发明的限制。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解,当我们称元件被“连接”或“耦接”到另一元件时,它可以直接连接或耦接到其他元件,或者也可以存在中间元件。此外,这里使用的“连接”或“耦接”可以包括无线连接或耦接。这里使用的措辞“和/或”包括一个或更多个相关联的列出项的任一单元和全部组合。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。
为便于对本发明实施例的理解,下面将结合附图以几个具体实施例为例做进一步的解释说明,且各个实施例并不构成对本发明实施例的限定。
本发明实施例提供的一种基于数字锁相解调的滤波器的在线钢轨缺陷检测的装置的结构图如图1所示,包括:依次连接的数模转换器、功率放大器、线圈传感器、低噪声差分放大器、模拟带通滤波器、模数转换器、单片机、显示模块以及WiFi模块。所述装置被安装在列车底部,并设置在钢轨踏面正上方一定距离处。
所述数模装换器用于将数字信号转换成激励模拟信号,将所述激励模拟信号传输给功率放大器;
所述功率放大器,用于将所述激励模拟信号放大,将放大后的激励模拟信号通入所述线圈传感器的激励线圈;
所述线圈传感器,用于包括中间线圈和前后线圈,所述中间线圈为激励线圈,所述前后线圈为检测线圈,所述前后线圈利用所述中间线圈激发的交变磁场感应检测钢轨缺陷,所述前后线圈产生的感应信号传输给差分放大器。图2为本发明实施例提供的一种线圈传感器检测钢轨缺陷示意图。其中间线圈为激励线圈,而前后线圈为检测线圈。
低噪声差分放大器,用于将所述两个检测线圈的感应信号进行差分放大,产生模拟的差分信号,将模拟的差分信号传输给模拟带通滤波器;
模拟带通滤波器,用于对差分信号宽泛地滤除低频和高频干扰,将滤除低频和高频干扰后的差分信号传输给模数转换器。
模数转换器,用于将接收到的差分信号转换成当前时刻的数字信号,将当前时刻的数字信号传输给单片机;
单片机,用于对数字信号进行基于卡尔曼滤波的数字锁相解调处理,得到数字信号的幅值和相位;
WiFi模块,用于接入互联网络,与所述单片机电路连接,并传输单片机解调得到的幅值和相位至远程服务器,以方便用户观看。
显示模块,用于与单片机电路连接,显示单片机解调得到的幅值和相位。
图3为本发明实施例提供的一种单片机中执行的基于卡尔曼滤波的数字锁相解调处理的流程图,具体处理过程包括:
根据输入信号的激励频率和采样频率,计算出通道1和通道2使用的相互正交的正弦参考信号和余弦参考信号,并将相互正交的正弦参考信号和余弦参考信号预先存储在内存当中。在实际应用中,可以通道1使用正弦参考信号,通道2使用与正弦参考信号正交的余弦参考信号;或者,通道1使用余弦参考信号,通道2使用正弦参考信号。
相互正交的正余弦参考信号是通过目标信号的同频采样计算得到的,预先存储N个正弦参考信号和余弦参考信号的值在内存,所述N是以采样间隔为1/fs采样目标信号一个周期所需要的点数,不同目标信号频率的N不同。所述参考信号与激励信号的频率相同。
正余弦参考信号的向量表示如下:
Sr[k]=(Srs[k] Src[k])T=(sin[Ω0k] cos[Ω0k])T
Srs[k]和Src[k]分别为当前时刻通道1和通道2的参考信号;
采样当前k个时刻输入的数字信号,如果当前k时刻不大于或等于采样点数N时继续进行信号采样,否则往下执行相敏检测。
模数转换器以采样频率fs将接收到的当前时刻的模拟差分信号转换成当前时刻的数字信号S[k],当前时刻的数字信号S[k]表示为:
S[k]=A0sin[Ω0k+θ0],
式(1)中,A0为数字信号中频率与激励频率f0相同的幅值,θ0为数字信号的相位,Ω0为数字信号的角频率,Ω0=2π/N,其中N是以采样频率为fs对频率为f0的差分信号采集一个周期所需要的点数,N=fs/f0
在每个通道(通道1或者通道2)对每一时刻输入的数字信号,分别与该通道的正弦参考信号或者余弦参考信号相乘,获得每个通道的频移结果。每一个通道对一个周期点数的频移结果进行累加求平均,将计算出的平均值作为该通道当前时刻的观测值。卡尔曼滤波器利用每个通道当前时刻的观测值对该通道的预测值进行校正更新,再利用每个通道更新后的预测值迭代地求出输入的数字信号的幅值和相位。
在通道1和通道2,分别对每一时刻输入的数字信号通过乘法器进行频移,当前时刻输入的数字信号的频移结果表示如下:
SPSD=S[k]·Sr[k]=(R[k] I[k])T
R[k]和I[k]分别是当前时刻通道1和通道2的频移结果;
在每个通道利用卡尔曼滤波器对前N个时刻的频移结果进行累加后求平均,将计算出的平均值作为该通道当前时刻的观测值Y[k],所述卡尔曼滤波器的观测方程表示如下:
Figure BDA0001570022850000091
式中,
SPSD[k]=(R[k] I[k])T
R[k]和I[k]分别是通道1和通道2第k时刻的频移结果,V[k]是当前时刻的观测噪声向量。
在每个通道,通过式(4)得到的当前时刻的观测值Y[k]来校正当前时刻的卡尔曼滤波器的预测值X[k|k-1],得到该通道当前时刻的最优估计值,通道1和通道2当前时刻的最优估计值分别为R′[k]和I′[k]。
每个通道当前时刻的卡尔曼滤波器的预测值X[k|k-1]等于上一时刻该通道的最优估计值,卡尔曼滤波器的预测方程的表示如下:
X[k|k-1]=(R′[k-1] I′[k-1])T
R′[k-1]和I′[k-1]分别是上一时刻双通道的最优估计值;
所述基于卡尔曼滤波数字锁相解调方法求得的所述输入的数字信号的幅值和相位,分别表示为:
Figure BDA0001570022850000101
卡尔曼滤波器的预测方程的误差方差表示如下:
P[k|k-1]=P[k-1]+qE,
式中,q是噪声的方差系数;
更新卡尔曼增益,更新方程表示如下:
K=P[k|k-1](P[k|k-1]+rE)-1
式中,r是测量误差的方差系数,其值满足rE=E(V[k]V T[k])T
当前k时刻最优估计更新方程表示如下:
X[k]=X[k|k-1]+K(Y[k]-X[k|k-1]),
P[k]=(E-K)P[k-1],
式中,X[k]=(R′[k]I′[k])T表示当前k时刻的最优估计,P[k]表示当前k时刻的最优方差估计;
所述噪声方差系数q和测量误差方差系数r是先验证值,需根据检测信号的噪声大小程度调整。
综上所述,本发明实施例提供的数字锁相解调的滤波器的实现方法,既保证信号解调的快速性,又提高了信号解调精度,且稳定性高,能够适应高速检测钢轨缺陷。
本发明实施例还公开了一种基于数字锁相解调的滤波器的用于在线钢轨缺陷快速检测的装置,所述装置通过无线WiFi链接服务器,其成本低、操作简单、测量精度高且能够装配在高速运行的列车上。
本领域普通技术人员可以理解:附图只是一个实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置或系统实施例而言,由于其基本相似于方法实施例,所以描述得比较简单,相关之处参见方法实施例的部分说明即可。以上所描述的装置及系统实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。

Claims (7)

1.一种数字锁相解调的滤波器的实现方法,其特征在于,包括:
在数字锁相解调的滤波器的两个通道分别输入数字信号,所述两个通道包括通道1和通道2;
在每个通道对每一时刻输入的数字信号分别与该通道的正弦参考信号或者余弦参考信号相乘,获得每个通道的频移结果,所述正弦参考信号与所述余弦参考信号正交;
在每个通道通过卡尔曼滤波器对一个周期点数的频移结果进行累加求平均,将计算出的平均值作为该通道当前时刻的观测值;
通过卡尔曼滤波器利用每个通道当前时刻的观测值对该通道的当前时刻的预测值进行校正更新,再利用每个通道更新后的预测值求出所述数字信号的幅值和相位;
所述的在每个通道通过卡尔曼滤波器对一个周期点数的频移结果进行累加求平均,将计算出的平均值作为该通道当前时刻的观测值,包括:
在每个通道利用卡尔曼滤波器对前N个时刻的频移结果进行累加后求平均,将计算出的平均值作为该通道当前时刻的观测值Y[k],所述卡尔曼滤波器的观测方程表示如下:
Figure FDA0002603631380000011
式中,
SPSD[k]=(R[k] I[k])T
R[k]和I[k]分别是通道1和通道2第k时刻的频移结果,V[k]是当前时刻的观测噪声向量,SPSD[k]是第k个时刻频移结果R[k]和I[k]组成的向量,SPSD[k-m]是第k个时刻往前数N个时刻的向量中的第m个时刻的向量;
所述的通过卡尔曼滤波器利用每个通道当前时刻的观测值对该通道的预测值进行校正更新,再利用每个通道更新后的预测值求出所述输入的数字信号的幅值和相位,包括:
在每个通道,利用当前时刻的观测值Y[k]来校正当前时刻的卡尔曼滤波器的预测值X[k|k-1],得到该通道当前时刻的最优估计值,通道1和通道2当前时刻的最优估计值分别为R′[k]和I′[k];
每个通道当前时刻的卡尔曼滤波器的预测值X[k|k-1]等于上一时刻该通道的最优估计值,卡尔曼滤波器的预测方程的表示如下:
X[k|k-1]=(R′[k-1] I′[k-1])T
R′[k-1]和I′[k-1]分别是上一时刻通道1和通道2的最优估计值;
所述输入的数字信号的幅值A0和相位θ0的计算公式如下:
Figure FDA0002603631380000021
Figure FDA0002603631380000022
2.根据权利要求1所述的方法,其特征在于,所述的在数字锁相解调的滤波器的两个通道分别输入数字信号之前,还包括:
根据输入的数字信号的激励频率和采样频率,计算出相互正交的正弦参考信号和余弦参考信号,所述通道1和所述通道2分别使用所述正弦参考信号和所述余弦参考信号中的一种信号,将N个相互正交的正弦参考信号和余弦参考信号预先存储在内存当中,所述N是以采样频率fs采样目标信号一个周期所需要的点数,所述正弦参考信号、余弦参考信号与所述数字信号的激励信号的频率相同。
3.根据权利要求2所述的方法,其特征在于,所述的在数字锁相解调的滤波器的两个通道分别输入数字信号之前,还包括:
模数转换器以采样频率fs将接收到的当前时刻的模拟差分信号转换成当前时刻的数字信号S[k],当前时刻的数字信号S[k]表示为:
S[k]=A0sin[Ω0k+θ0]
A0为数字信号中频率与激励频率f0相同的幅值,θ0为数字信号的相位,Ω0为数字信号的角频率,Ω0=2π/N,其中N是以采样频率为fs对频率为f0的差分信号采集一个周期所需要的点数,N=fs/f0
4.根据权利要求3所述的方法,其特征在于,所述的在每个通道对每一时刻输入的数字信号分别与该通道的正弦参考信号或者余弦参考信号相乘,获得每个通道的频移结果,包括:
在通道1和通道2,分别对每一时刻输入的数字信号通过乘法器进行频移,当前时刻输入的数字信号的频移结果表示如下:
SPSD[k]=S[k]·Sr[k]=(R[k] I[k])T
R[k]和I[k]分别是当前时刻通道1和通道2的频移结果;
Sr[k]=(Ssr[k] Scr[k])T
Ssr[k]Scr[k]分别是当前时刻的正弦参考信号和余弦参考信号,Sr[k]是当前时刻正弦参考信号和余弦参考信号组成的向量。
5.根据权利要求1至4任一项所述方法,其特征在于,设置在线钢轨缺陷检测的装置,包括:依次连接的数模转换器、功率放大器、线圈传感器、差分放大器、模拟带通滤波器、模数转换器和单片机;
所述数模转换器,用于将数字信号转换成激励模拟信号,将所述激励模拟信号传输给功率放大器;
所述功率放大器,用于将所述激励模拟信号放大,将放大后的激励模拟信号传输给所述线圈传感器的激励线圈;
所述线圈传感器,用于包括中间线圈和前后线圈,所述中间线圈为激励线圈,所述前后线圈为检测线圈,所述前后线圈利用所述中间线圈激发的交变磁场感应检测钢轨缺陷,所述前后线圈产生的感应信号传输给差分放大器;
差分放大器,用于将两个所述检测线圈的感应信号进行差分放大,产生模拟的差分信号,将模拟的差分信号传输给模拟带通滤波器;
模拟带通滤波器,用于滤除所述模拟的差分信号的低频和高频干扰,将处理后的差分信号传输给所述模数转换器;
模数转换器,用于将接收到的差分信号转换成当前时刻的数字信号,将当前时刻的数字信号传输给单片机;
单片机,用于通过权利要求1至4任一项所述的方法对接收到的数字信号进行基于卡尔曼滤波的数字锁相解调处理,得到数字信号的幅值和相位。
6.根据权利要求5所述的方法,其特征在于,所述的在线钢轨缺陷检测的装置还包括:
WiFi模块,用于接入互联网络,与所述单片机电路连接,传输单片机解调得到的幅值和相位至远程服务器,以方便用户观看。
7.根据权利要求6所述的方法,其特征在于,所述的在线钢轨缺陷检测的装置还包括:
显示模块,用于与所述单片机电路连接,显示单片机解调得到的幅值和相位。
CN201810113634.6A 2018-02-05 2018-02-05 数字锁相解调的滤波器的实现方法及装置 Active CN108426942B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810113634.6A CN108426942B (zh) 2018-02-05 2018-02-05 数字锁相解调的滤波器的实现方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810113634.6A CN108426942B (zh) 2018-02-05 2018-02-05 数字锁相解调的滤波器的实现方法及装置

Publications (2)

Publication Number Publication Date
CN108426942A CN108426942A (zh) 2018-08-21
CN108426942B true CN108426942B (zh) 2020-09-15

Family

ID=63156578

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810113634.6A Active CN108426942B (zh) 2018-02-05 2018-02-05 数字锁相解调的滤波器的实现方法及装置

Country Status (1)

Country Link
CN (1) CN108426942B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108507695A (zh) * 2017-02-27 2018-09-07 北京自动化控制设备研究所 一种基于锁相放大的原子气室温度检测方法
MX2018005076A (es) * 2018-04-25 2019-10-28 Inst Tecnologico Estudios Superiores Monterrey Sistema, método y aparato para la evaluación y seguimiento del rendimiento muscular con retroalimentación autoajustable.
CN109067695B (zh) * 2018-10-12 2020-11-17 苏州大学 基于级联卡尔曼滤波器减少噪声干扰的方法及系统
CN113395228B (zh) * 2021-06-10 2022-08-23 深圳清华大学研究院 一种针对2fsk信号的解调电路

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4791573A (en) * 1986-09-02 1988-12-13 Sanders Associates, Inc. State-deviation-estimation circuit employing a phase-locked-loop phase reference
WO1997003510A1 (en) * 1995-07-07 1997-01-30 Northern Telecom Limited Carrier recovery for digitally phase modulated signals, using a known sequence
CN102087362B (zh) * 2009-12-08 2014-04-16 北京邮电大学 信号跟踪的方法和装置
CN102403969A (zh) * 2011-11-03 2012-04-04 北京交通大学 一种锁相放大器及用于该锁相放大器的信号解调方法
CN204255900U (zh) * 2014-11-16 2015-04-08 吉林大学 一种钢件缺陷的电磁无损检测装置
CN105866239A (zh) * 2016-04-25 2016-08-17 江南大学 基于铁磁性试件的u型脉冲融合型传感器检测方法

Also Published As

Publication number Publication date
CN108426942A (zh) 2018-08-21

Similar Documents

Publication Publication Date Title
CN108426942B (zh) 数字锁相解调的滤波器的实现方法及装置
CN109495410B (zh) 高动态pcm/fm信号载波频率精确估计方法
CN106855628B (zh) 一种高动态卫星导航信号的快速捕获和跟踪系统和方法
JP2017526512A5 (zh)
CN107342713A (zh) 一种基于改进滑模观测器的永磁同步电机转子位置检测装置及方法
CN105676008B (zh) 一种数字式电场传感器
CN111884984A (zh) 快变载波多普勒频移捕获系统
CN111865865A (zh) 适用于高灵敏度星载ads-b接收机的频偏及相偏估计方法
CN104502701A (zh) 基于相位调制检测电力信号频率的方法和系统
CN106574848A (zh) 带有可变的载波频率和戈泽尔滤波的传感器装置
CN104502702A (zh) 检测电力信号的频率的方法和系统
CN115265691B (zh) 一种科氏流量计振动频率跟踪方法及系统
CN108344501A (zh) 一种信号相关性应用中共振辨识和消除方法及装置
CN106027454A (zh) 基于频率补偿的qam信号的频偏估计方法
US9461662B2 (en) System for analyzing the frequency of a signal, a method thereof and a system for measuring the relative phase between two input signals
WO2016062033A1 (zh) 频率补偿处理方法及装置
CN114301552A (zh) 一种数字调制信号测试方法及系统
US9094266B2 (en) Efficient de-spreading
CN115615467B (zh) 霍尔信号解调方法和装置
CN104618033A (zh) 一种多层自适应形态滤波重力信号噪声抑制方法
CN105158817A (zh) 深空探测器多普勒频率被动式测量方法
CN111580188B (zh) 一种大地电磁时间域标定方法及系统
CN116817739B (zh) 旋转变压器的误差补偿控制方法及设备
CN113721271B (zh) 一种北斗信号接收机接收信号快速并行捕获方法和系统
US11614344B1 (en) Device and method for determining an angular position of inductive position sensor

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