CN101806834B - 基于卡尔曼滤波器的信号实时时频谱仪 - Google Patents
基于卡尔曼滤波器的信号实时时频谱仪 Download PDFInfo
- Publication number
- CN101806834B CN101806834B CN2010101354382A CN201010135438A CN101806834B CN 101806834 B CN101806834 B CN 101806834B CN 2010101354382 A CN2010101354382 A CN 2010101354382A CN 201010135438 A CN201010135438 A CN 201010135438A CN 101806834 B CN101806834 B CN 101806834B
- Authority
- CN
- China
- Prior art keywords
- time
- module
- signal
- frequency
- parameter
- 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.)
- Expired - Fee Related
Links
- 238000001228 spectrum Methods 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 18
- 238000006243 chemical reaction Methods 0.000 claims abstract description 15
- 238000009499 grossing Methods 0.000 claims abstract description 10
- 238000012937 correction Methods 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims abstract description 6
- 238000007781 pre-processing Methods 0.000 claims description 17
- 238000009415 formwork Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 9
- 230000003278 mimic effect Effects 0.000 claims description 8
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 abstract description 3
- 238000011282 treatment Methods 0.000 abstract 2
- 238000004458 analytical method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000004304 visual acuity Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000009418 renovation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于卡尔曼滤波器的信号实时时频谱仪,包括传感器、A/D数据转换卡、处理器、控制面板、存储器和显示器。所述处理器包括依次连接的预处理模块、预测模块、修正模块和平滑模块。通过控制面板输入指令,将所述传感器获取的非平稳时变信号经过A/D数据转换卡进行采样和转换后,依次经过预处理模块、预测模块、修正模块和平滑模块对信号采样点处的参数进行预测、修正和平滑处理,将处理后的参数带入到时频谱函数后便可得到高分辨率的时频谱,最终在显示器上显示该时频谱。利用本发明频谱仪能够对信号进行实时、动态、连续的处理,可以得到具有高分辨率且又无交叉项干扰的时频谱,从而为后续的信号处理带来极大的便利。
Description
技术领域
本发明涉及一种智能动态测控系统,尤其涉及一种基于卡尔曼滤波器的信号实时时频谱生成系统。
背景技术
时频分析是通过设计时间和频率的联合函数来描述信号在不同时间和频率的能量密度或强度。时间和频率的这种联合函数简称为时频分布。利用时频分布来分析信号,能给出各个时刻的瞬时频率及其幅值,并且能够进行时频滤波和时变信号研究。对于时频分析作为分析时变非平稳信号的有力工具,成为现代信号处理研究的一个热点,时频分析作为一种新兴的信号处理方法,近年来受到越来越多的重视。在现代通讯、航空航天、装备制造、医学等领域有广阔的应用前景。
信号时频处理方法主要有短时傅里叶变换、Wigner-Ville分布、小波尺度分析等。其中,短时傅里叶变换实质是分段FFT,其前提是信号在每个分段内平稳或基本平稳。为了得到较高的频率分辨率,需要选择较大的窗长度,但是这却会影响到时域的分辨率,因此,短时傅里叶变换不能同时满足时域和频域分辨率。Wigner-Ville分布在当信号只包含单一的时频成分时,在时域和频域都能得到很高的分辨率,但在分析多频率成分信号时,由于该信号为二次型变换,不可避免地会出现交叉项干扰。小波尺度分析虽然能够在一定程度上改善上述干扰情况,但是小波尺度分析方法的时频分辨率会受到不确定性原理的限制,因此利用小波尺度分析方法无法同时在时域和频域都得到很好的分辨率。
因此,如何能够将含有多成分的信号明显的区分出来,并且在时域和频域都能得到较高的分辨率是目前时频分析领域中急需攻破的技术难题,解决该技术难题可以对后续的信号处理带来了很大的便利,具有重大的实用价值和指导意义。
发明内容
针对上述现有技术,本发明运用卡尔曼滤波原理,提供一种基于卡尔曼滤波器的信号实时时频谱仪,利用本发明频谱仪可以得到具有高分辨率且又无交叉项干扰的时频谱,从而为后续的信号处理带来极大的便利。
为了解决上述技术问题,本发明基于卡尔曼滤波器的信号实时时频谱仪予以实现的技术方案是:该时频谱仪包括传感器、A/D数据转换卡、处理器、控制面板、存储器和显示器。所述处理器包括依次连接的预处理模块、预测模块、修正模块和平滑模块。通过控制面板输入指令,将所述传感器获取的非平稳时变信号经过A/D数据转换卡进行采样和转换后,依次经过预处理模块、预测模块、修正模块和平滑模块对信号采样点处的参数进行预测、修正、平滑处理,通过运算后即可得到高分辨率的时频谱,最终在显示器上显示该时频谱图。
本发明基于卡尔曼滤波器的信号实时时频谱仪,其中,所述预处理模块、预测模块、修正模块和平滑模块分别是由运算放大器构成的模拟电路,所述运算放大器包括加法器、减法器、乘法器、矩阵转置单元和矩阵求逆单元。其中,所述预处理模块运用最大似然法对信号的白噪声方差σ2和参数的高斯随机噪声方差Qr进行估计,并根据具体信号来确定最佳阶数,在保证运算速度的前提下使拟合残差最小;所述预测模块用于根据上一采样点处得到的参数来预测下一个采样点处的参数;所述修正模块用于根据卡尔曼增益对上一环节得到的预测参数进行修正以提高参数的准确度;所述平滑模块用于对参数进行固定区间平滑。
与现有技术相比,本发明的有益效果是:
由于本发明是基于卡尔曼滤波技术获得时频谱,非常适合非平稳时变信号的处理,得到的时频谱分辨率高,并且没有交叉项干扰,能够将信号中不同成分清晰分辨出来。另外,由于卡尔曼滤波中方差的不断更新,因此所需要的存储空间更小,处理速度也大为提高,可以实现时频谱的实时生成。因而,本发明基于卡尔曼滤波器的信号实时时频谱仪能够对信号进行实时、动态、连续的处理。
附图说明
图1是本发明基于卡尔曼滤波器的信号实时时频谱仪的构成框图;
图2是传感器采集到的原始信号图;
图3是图2所示信号经本发明时频谱仪处理后得到的时频谱图;
图4是利用本发明时频谱仪生成高分辨率时频谱图的流程图;
图5为图1中所示修正模块模拟电路的逻辑原理图。
说明书附图中附图标记说明如下:
1.传感器,2.A/D数据转换卡,3.预处理模块,4.预测模块,5.修正模块,6.存储器,7.平滑模块,8.显示器。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细地描述。
如图1所示,本发明基于卡尔曼滤波器的信号实时时频谱仪的构成是,包括传感器1、A/D数据转换卡2、处理器、控制面板、存储器6和显示器8。所述处理器包括依次连接的预处理模块3、预测模块4、修正模块5和平滑模块7;所述控制面板的主要功能是通过控制面板输入指令,并设定系统正常运行所必需的一些参数,即:通过控制面板输入指令,对系统进行初始化后,将所述传感器1获取的非平稳时变信号经过A/D数据转换卡2进行采样和转换,依次经过预处理模块3、预测模块4、修正模块5和平滑模块7对信号采样点处的参数进行预测、修正和平滑处理,将处理后的参数带入到时频谱函数,通过运算后便可得到高分辨率的时频谱,最终在显示器8上显示该时频谱。
所述A/D数据转换卡2集成有A/D数据采样、转换以及计数功能。
所述预处理模块3的作用包括:(1)运用最大似然法对信号的白噪声方差σ2和参数的高斯随机噪声方差Qr进行估计;(2)根据具体信号来确定最佳阶数,在保证运算速度的前提下又要使拟合残差最小。
所述预处理模块3、预测模块4、修正模块5、平滑模块7分别是由运算放大器组成的模拟电路,所述运算放大器包括加法器、减法器、乘法器、矩阵转置单元和矩阵求逆单元。其中,预测模块4的功能是根据上一采样点处得到的参数来预测下一个采样点处的参数;修正模块5的功能是根据卡尔曼增益对上一环节得到的预测参数进行修正,从而提高参数的准确度。由于卡尔曼滤波中方差不断更新,可能会造成时变参数运算的中断,从而影响到时频谱的精度和分辨率,因此有必要对参数进行固定区间平滑。因此当经由预测模块4、修正模块5得到所有对应N个采样点的参数后,还需要经过平滑模块7进行进一步的平滑处理。
本发明中所述修正模块5的模拟电路逻辑原理图如图5所示,本发明中预处理模块3、预测模块4、平滑模块7的模拟电路逻辑原理图与图5类似。
所述存储器6用来存储运算中的中间数据。
下面用一具体实施例来描述利用本发明基于卡尔曼滤波器的信号实时时频谱仪获得高分频率时频谱的过程,如图4所示。
步骤一、初始化参数,包括:
在传感器进行信号采集前,为了使基于卡尔曼滤波器的信号实时时频谱仪开始正常工作,要通过控制面板事先设定好参数。
需要设定的参数有:
(1)初始时刻的时变参数θ0与其对应的方差C0;因为,随着卡尔曼滤波器的工作,参数θ会逐渐收敛,因此,上述两个参数θ0和C0可任取;本实施例中,设:θ0包含n个元素,全部设为0;C0包含n个元素,全部设为1。
(2)设置信号的白噪声方差σ2和参数的高斯随机噪声方差Qr的初值,σ2设为0.5,Qr含n个元素,全部设为0.1。
(3)设置A/D转换卡采样频率,采样频率一般为信号最高频率的3-4倍,本例中最高频率为350Hz,设定采样频率为1024Hz。
(4)设置采样点N:采样时间设置为1秒,采样频率1024Hz,故N=1024。
当参数设置完毕后,系统即可开始工作。
在本例中,传感器1采集到的原始信号图如图2所示,信号源产生的信号记为x(t),x(t)由三个非平稳信号x1(t),x2(t),x3(t)叠加组成,即x(t)=x1(t)+x2(t)+x3(t)。
其中x1(t)=cos(2πf1t);f1=350-300t;
x1(t)=cos(2πf2t);f1=50+300t;
x1(t)=cos(2πf3t);f1=300。
步骤二、采集数据,并A/D转换
传感器1在t时刻采集到的信号记为xt,xt可以有时变参数表示为xt=-θ1,txt-1-θ2,txt-2…-θn,txt-n+et。
其中,θj,t(j=1,2,…n)是时变参数,et为白噪声序列,其平均值为0,方差为σ2。
引入矩阵Ht和θt,模型可以写成以下矩阵的形式
xt=Htθt+et
其中:Ht=(-xt-1,-xt-1,…-xt-n)
θt=(θ1,t,θ2,t,…θn,t,)=θt-1+ωt
ωt为高斯随机噪声,其方差记为Qr。
将传感器1采集到的信号传递给A/D数据转换卡2,A/D数据转换卡2开始对信号进行采样和转换。
步骤三、利用处理器对信号进行分析处理,包括:
(1)预处理:当采集的样本数量达到N时,系统将长度为N的时间序列x1,x2,…,xN依次输入到预处理模块3,预处理模块3对信号进行分析,自动确定阶数m;同时根据信号的白噪声方差σ2和参数的高斯随机噪声方差Qr的初值对σ2和Qr进行最大似然估计。
(2)得到参数θt|t-1和Ct|t-1:预测模块4将会根据xt-1对应的参数θt-1来预测xt的参数得到θt|t-1,同时其对应的方差Ct|t-1也同步更新。
具体的运算步骤如下:
(3)将上述的θt|t-1和Ct|t-1参数传递到修正模块5中,修正模块5根据卡尔曼增益Kt和预测误差εt对θt|t-1作进一步的修正,使之更准确;同时,方差Ct更新。
具体的运算步骤如下:
本技术领域内技术人员可以根据上述步骤三(3)中所提及到的公式设计该修正模块的模拟电路逻辑原理图,如图5所示,其中,Ht、xt、σ2为从输入端输入的数据,以上各数据按照图5中所示的各基本运算单元进行加法、减法、乘法以及矩阵的转置、求逆等运算后,从输出端输出本发明中预处理模块3、预测模块4和平滑模块7的模拟电路逻辑原理图均与图5类似,在此不再赘述。
(4)存储数据:Ct和修正后得到的参数θt将会被存放在存储器6中。
(5)当t<N时,则返回上述步骤三中(2),继续对参数进行预测和修正,直至完成对N个采样点所对应的N个参数的预测修正处理。
步骤四、数据的平滑处理:
参数θt经由平滑模块7进行固定区间平滑,平滑处理时,θt按照t=N-1到t=1的顺序进行处理,θt|N的值将在θt+1|N的基础上得到。
具体的运算步骤如下:
Lt的初始值设为0。上述过程从t=N-1一直到t=1使得每一时刻经过平滑的θt都可以得到。
步骤五、频谱的显示输出:
当参数θt|N全部处理完毕后,将参数θt|N代入到信号功率谱密度函数中,得到非平稳时变信号的高分辨率时频谱,通过显示器8显示上述时频谱图,如图3所示。
信号功率谱密度函数为:
由图3可以看出,经过基于卡尔曼滤波器的信号实时时频谱仪处理后得到的时频域图的时间和频率分辨率都较高,不存在交叉项干扰,信号中的不同成分可以清晰的分辨出来。当处理完序列长度为N的信号后,本发明时频谱仪会继续对后面的信号进行采样、处理。即信号以N个采样点为一组经由系统进行实时、动态、连续的处理,从而能够获得实时的时频谱图。
本发明的优点是可快速获得高分辨率的非平稳信号时频谱,对于后续的信号分析与处理提供了必要的技术保障。而且卡尔曼滤波的方差由于不断的更新,因此占用储存空间更小,处理速度更快,可以实现时频谱的实时生成。
尽管上面结合图对本发明进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨的情况下,还可以作出很多变形,这些均属于本发明的保护之内。
Claims (2)
1.一种基于卡尔曼滤波器的信号实时时频谱仪,包括传感器(1)、A/D数据转换卡(2)、处理器、控制面板、存储器(6)和显示器(8);其特征在于:
所述处理器包括依次连接的预处理模块(3)、预测模块(4)、修正模块(5)和平滑模块(7);
通过控制面板输入指令,将所述传感器(1)获取的非平稳时变信号经过A/D数据转换卡(2)进行采样和转换后,依次经过预处理模块(3)、预测模块(4)、修正模块(5)和平滑模块(7)对信号采样点处的参数进行预测、修正、平滑处理,通过运算后即可得到高分辨率的时频谱,最终在显示器(8)上显示该时频谱图;
所述预处理模块(3)、预测模块(4)、修正模块(5)和平滑模块(7)分别是由运算放大器构成的模拟电路,所述运算放大器包括加法器、减法器、乘法器、矩阵转置单元和矩阵求逆单元。
2.根据权利要求1所述基于卡尔曼滤波器的信号实时时频谱仪,其特征在于:所述预处理模块(3)运用最大似然法对信号的白噪声方差σ2和参数的高斯随机噪声方差Qr进行估计,并根据具体信号来确定最佳阶数,在保证运算速度的前提下使拟合残差最小;所述预测模块(4)用于根据上一采样点处得到的参数来预测下一个采样点处的参数;所述修正模块(5)用于根据卡尔曼增益对上一环节得到的预测参数进行修正以提高参数的准确度;所述平滑模块(7)用于对参数进行固定区间平滑。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101354382A CN101806834B (zh) | 2010-03-30 | 2010-03-30 | 基于卡尔曼滤波器的信号实时时频谱仪 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101354382A CN101806834B (zh) | 2010-03-30 | 2010-03-30 | 基于卡尔曼滤波器的信号实时时频谱仪 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101806834A CN101806834A (zh) | 2010-08-18 |
CN101806834B true CN101806834B (zh) | 2012-07-04 |
Family
ID=42608718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101354382A Expired - Fee Related CN101806834B (zh) | 2010-03-30 | 2010-03-30 | 基于卡尔曼滤波器的信号实时时频谱仪 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101806834B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102650658B (zh) * | 2012-03-31 | 2014-03-19 | 机械工业第三设计研究院 | 一种时变非平稳信号时频分析方法 |
CN103419083B (zh) * | 2012-11-05 | 2015-10-28 | 上海理工大学 | 数控机床进给系统受力监测方法 |
CN103196450B (zh) * | 2013-04-02 | 2014-05-21 | 武汉大学 | 基于模拟电路形式的卡尔曼滤波方法及模拟电路 |
CN104714088B (zh) * | 2015-02-09 | 2017-12-08 | 苏州蓝色力量信息科技有限公司 | 一种频谱仪及频谱测量方法 |
CN108680787A (zh) * | 2018-05-23 | 2018-10-19 | 成都玖锦科技有限公司 | 基于fpga的实时频谱分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5420501A (en) * | 1993-02-01 | 1995-05-30 | Nsk Ltd. | Frequency spectrum analyzer |
CN101090379A (zh) * | 2007-07-20 | 2007-12-19 | 重庆重邮信科(集团)股份有限公司 | 采用Kalman滤波器进行联合检测的方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3375919B2 (ja) * | 1999-11-11 | 2003-02-10 | アンリツ株式会社 | 信号分析装置 |
-
2010
- 2010-03-30 CN CN2010101354382A patent/CN101806834B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5420501A (en) * | 1993-02-01 | 1995-05-30 | Nsk Ltd. | Frequency spectrum analyzer |
CN101090379A (zh) * | 2007-07-20 | 2007-12-19 | 重庆重邮信科(集团)股份有限公司 | 采用Kalman滤波器进行联合检测的方法及装置 |
Non-Patent Citations (3)
Title |
---|
JP特开2001-141764A 2001.05.25 |
章旭景等.基于子带卡尔曼滤波的语音增强方法.《信号处理》.2009,第25卷(第9期),第1475,1478页及图3. * |
邓自立等.基于Kalman滤波的固定区间观测白噪声平滑器及其在信号估计中的应用.《2002中国控制与决策学术年会论文集》.2002,第201页. * |
Also Published As
Publication number | Publication date |
---|---|
CN101806834A (zh) | 2010-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101806834B (zh) | 基于卡尔曼滤波器的信号实时时频谱仪 | |
Kawahara et al. | Tandem-STRAIGHT: A temporally stable power spectral representation for periodic signals and applications to interference-free spectrum, F0, and aperiodicity estimation | |
CN105391083A (zh) | 基于变分模态分解和相关向量机的风功率区间短期预测方法 | |
CN107167306B (zh) | 基于阶次提取的旋转机械转子运行状态模态分析方法 | |
CN103268519A (zh) | 基于改进Lyapunov指数的电力系统短期负荷预测方法及装置 | |
CN103236825A (zh) | 一种用于高精度数据采集系统的数据校正方法 | |
CN105158200A (zh) | 一种提高近红外光谱定性分析准确度的建模方法 | |
CN116306798A (zh) | 一种超短时风速预测方法及系统 | |
CN104849545A (zh) | 一种电力系统频率测量方法及测量装置 | |
CN106407659A (zh) | 一种空气质量指数预测方法和装置 | |
CN105488335A (zh) | 一种基于李雅普诺夫指数的电力系统负荷预测方法及装置 | |
CN109300484B (zh) | 音频对齐方法、装置、计算机设备以及可读存储介质 | |
CN103279030B (zh) | 基于贝叶斯框架的动态软测量建模方法及装置 | |
Jing et al. | Runoff forecasting model based on variational mode decomposition and artificial neural networks | |
CN102830282A (zh) | 一种2560点分组基2快速傅里叶变换方法 | |
CN109617051B (zh) | 一种新能源电力系统低频振荡参数辨识方法 | |
CN103312337B (zh) | 一种振动信号的稀疏矩阵的自适应获取方法 | |
CN103487653A (zh) | 一种用于实时频谱轨迹处理的装置及方法 | |
CN105137183A (zh) | 一种电力系统谐波分析方法及系统 | |
CN113780356B (zh) | 基于集成学习模型的水质预测方法及系统 | |
CN111551785B (zh) | 基于无迹卡尔曼滤波的频率与谐波检测方法 | |
CN112531710B (zh) | 一种谐波源接入预测评估的方法 | |
CN106096494A (zh) | 一种基于信息聚合构架的行波故障和干扰录波数据辨识方法 | |
CN105425037A (zh) | 一种频谱分析方法及装置 | |
CN104268436B (zh) | 基于趋势拟合的风电日内波动连续时段识别方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120704 |