CN104506163B - 电压信号的滤波方法及其滤波装置 - Google Patents
电压信号的滤波方法及其滤波装置 Download PDFInfo
- Publication number
- CN104506163B CN104506163B CN201410825332.3A CN201410825332A CN104506163B CN 104506163 B CN104506163 B CN 104506163B CN 201410825332 A CN201410825332 A CN 201410825332A CN 104506163 B CN104506163 B CN 104506163B
- Authority
- CN
- China
- Prior art keywords
- voltage
- voltage signal
- kalman
- filtered
- digitized
- 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
Links
Abstract
本发明公开了电压信号的滤波方法,包括:模数转换器对待滤波的电压信号进行采样,并测量出噪声的峰‑峰值Vpp,然后发给微处理器;微处理器根据Vpp计算电压变化阈值Vt,将数字化的电压信号的电压值与Vt进行比较,如果小于Vt,则用卡尔曼数字滤波器对当前时刻的数字化的电压信号进行滤波,并将卡尔曼数字滤波器在当前时刻的X(k|k)作为下一时刻的X(k‑1|k‑1);如果大于Vt,则不进行滤波,将当前采集到的数字化电压信号的电压值作为卡尔曼数字滤波器在下一时刻的X(k‑1|k‑1),卡尔曼数字滤波器的Q远小于1,R≥ 1。本发明还公开了电压信号的滤波装置。本发明能获得动态响应快、收敛稳定性好的滤波效果。
Description
技术领域
本发明涉及一种电压信号的滤波方法及其滤波装置。
背景技术
自适应前大灯系统(Automatic Front-Lamp system,简称AFS)能够根据路况做动态调节以提高行车安全。目前,AFS控制器中电源管理模块的滤波普遍采用滑动平均法,为具有动态响应快同时不可避免的使收敛稳定性变坏。
卡尔曼数字滤波器是一个最优化自回归数据处理算法,标准卡尔曼数字滤波器的5个基本公式为:
X(k|k-1)=A X(k-1|k-1) ……… (1)
P(k|k-1)=A P(k-1|k-1) A’+Q ……… (2)
X(k|k)= X(k|k-1)+Kg(k) (Z(k)-H X(k|k-1)) ……… (3)
Kg(k)= P(k|k-1) H’ / (H P(k|k-1) H’ + R) ……… (4)
P(k|k)=(I-Kg(k) H)P(k|k-1) ……… (5)
公式(1)中,X(k|k-1)是利用上一状态预测的结果,X(k-1|k-1)是上一状态最优的结果。
公式(2)中,P(k|k-1)是X(k|k-1)对应的协方差(covariance),P(k-1|k-1)是X(k-1|k-1)对应的协方差,A’表示A的转置矩阵,Q是系统过程的协方差,Q增大,动态响应变快,收敛稳定性变坏;Q减小,动态响应变慢,收敛稳定性变好。
公式(3)、(4)中,Kg为卡尔曼增益(Kalman Gain)。
公式(4)中,R是测量噪声的协方差,R增大,动态响应变慢,收敛稳定性变好;R减小,动态响应变快,收敛稳定性变坏
公式(5)中, I 为1的矩阵,对于单模型单测量,I=1。当系统进入k+1状态时,P(k|k)就是式子(2)的P(k-1|k-1)。
卡尔曼数字滤波器的特征为动态响应慢、收敛稳定性好或动态响应快、收敛稳定性差。对于需要动态响应快、收敛稳定性好的滤波效果,使用卡尔曼滤波便无法实现。因此,现行的滤波方法均不能同时保证动态响应时间快和收敛稳定性好。
发明内容
本发明所要解决的技术问题在于提供一种电压信号的滤波方法及其滤波装置,其能获得动态响应快、收敛稳定性好的滤波效果。
为解决上述技术问题,本发明所采用的技术方案是:
一种电压信号的滤波方法, 包括:
模数转换器对待滤波的电压信号进行采样,并测量出该电压信号的噪声的峰-峰值Vpp;
模数转换器将数字化的电压信号及测量到的噪声峰-峰值Vpp发送给微处理器;
微处理器根据测量到的噪声峰-峰值Vpp计算电压变化阈值Vt,Vt=N*Vpp,1.5≤N≤2;
微处理器将数字化的电压信号的电压值与该电压变化阈值Vt进行比较,如果数字化的电压信号的电压值小于电压变化阈值Vt,则用卡尔曼数字滤波器对当前时刻的数字化的电压信号进行滤波,并将卡尔曼数字滤波器在当前时刻的X(k|k)作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1);如果数字化的电压信号的电压值大于电压变化阈值Vt,则不对当前时刻的数字化的电压信号进行滤波,将当前采集到的数字化电压信号的电压值作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),其中,所述卡尔曼数字滤波器的Q远小于1,R大于等于1。
本发明还提供了一种电压信号的滤波装置,其特点在于,包括模数转换器和微处理器,模数转换器的输出端与微处理器的输入端电连接;其中:
模数转换器用于对待滤波的电压信号进行采样,测量出该电压信号的噪声的峰-峰值Vpp,并将数字化的电压信号及测量到的噪声峰-峰值Vpp发送给微处理器;
微处理器用于根据测量到的噪声峰-峰值Vpp计算电压变化阈值Vt,并将数字化的电压信号的电压值与该电压变化阈值Vt进行比较,如果数字化的电压信号的电压值小于电压变化阈值Vt,则用卡尔曼数字滤波器对当前时刻的数字化的电压信号进行滤波,将卡尔曼数字滤波器在当前时刻的X(k|k)作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),如果数字化的电压信号的电压值大于电压变化阈值Vt,则不对当前时刻的数字化的电压信号进行滤波,将当前采集到的数字化电压信号的电压值作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),其中,Vt=N*Vpp,1.5≤N≤2,所述卡尔曼数字滤波器的Q远小于1,R大于等于1。
本发明通过将采集到的数字电压信号与设置的电压变化阀值进行比较,可以分辨出信号与噪声:若采集的电压信号的电压值大于电压变化阈值,则判断采集到的信号更接近于实际的信号,因此不做滤波;若采集的电压信号的电压值小于电压变化阈值,则判断滤波后的信号更接近于实际的信号,因此要使用卡尔曼数字滤波器进行滤波,如此缩短了动态响应时间,提高了响应速率;同时,本发明通过调整Q、R的值使收敛稳定性变好。
附图说明
图1是根据本发明一实施例的电压信号滤波装置的原理框图。
图2示出了采用图1的电压信号滤波装置进行滤波的效果示意图。
图3示出了单纯采用卡尔曼数字滤波器进行滤波的效果示意图。
具体实施方式
根据本发明一实施例的一种电压信号的滤波方法, 包括以下步骤:
模数转换器对待滤波的电压信号进行采样,并测量出该电压信号的噪声的峰-峰值Vpp;
模数转换器将数字化的电压信号及测量到的噪声峰-峰值Vpp发送给微处理器;
微处理器根据测量到的噪声峰-峰值Vpp计算电压变化阈值Vt,Vt=N*Vpp,1.5≤N≤2;
微处理器将数字化的电压信号的电压值与该电压变化阈值Vt进行比较,如果数字化的电压信号的电压值小于电压变化阈值Vt,则用卡尔曼数字滤波器对当前时刻的数字化的电压信号进行滤波,并将卡尔曼数字滤波器在当前时刻的X(k|k)作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1);如果数字化的电压信号的电压值大于电压变化阈值Vt,则不对当前时刻的数字化的电压信号进行滤波,将当前采集到的数字化电压信号的电压值作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),其中,所述卡尔曼数字滤波器的Q远小于1,R大于等于1。上述的当前时刻及下一时刻与模数转换器的采样频率有关。
如果待滤波的电压信号的电压值超过模数转换器的最大允许采样电压,则在用模数转换器对待滤波的电压信号进行采样之前,先采用电压衰减器对待滤波的电压信号进行电压衰减,使待滤波的电压信号的电压值小于模数转换器的最大允许采样电压,然后再用模数转换器对待滤波的电压信号进行采样。
在一具体的实施例中,卡尔曼数字滤波器通过以下步骤进行滤波:
步骤a、给定初始值X(0);X(0)为0或任意大于0的整数;
步骤b、X(k|k-1)= X(k-1|k-1);
步骤c、P(k|k-1)= P(k-1|k-1) +Q;
步骤d、Kg(k)= P(k|k-1) / (P(k|k-1) + R)
步骤e、X(k|k)= X(k|k-1)+Kg(k) (Z(k)- X(k|k-1));Z(k)为k时刻采集到的数字化的电压值;
步骤f、P(k|k)=(1-Kg(k))P(k|k-1)。
其中,优选卡尔曼数字滤波器中的Q小于等于0.00001,R大于等于1。更优选地,Q等于0.00001,R等于1。
本发明还提供了一种实现上述电压信号滤波方法的电压信号滤波装置。请参考图1,以下以本发明的电压信号滤波装置在自适应前大灯系统的电源管理中的应用来说明本发明的结构及工作原理。
该电压信号的滤波装置1包括电压衰减器11、模数转换器13和微处理器15。模数转换器13的输出端与微处理器15的输入端电连接。
电压衰减器11的输出端与模数转换器13的输入端电连接。由于直流电压源2的输出信号的电压大于模数转换器的最大允许采样电压,需要电压衰减器11对直流电压源2输出的电压信号进行电压衰减,使之小于模数转换器13的最大允许采样电压。
模数转换器13用于对待滤波的电压信号进行采样,测量出该电压信号的噪声的峰-峰值Vpp,并将数字化的电压信号及测量到的噪声峰-峰值Vpp发送给微处理器15。
微处理器15用于根据测量到的噪声峰-峰值Vpp计算电压变化阈值Vt,并将数字化的电压信号的电压值与该电压变化阈值Vt进行比较,如果数字化的电压信号的电压值小于电压变化阈值Vt,则用卡尔曼数字滤波器对当前时刻的数字化的电压信号进行滤波,将卡尔曼数字滤波器在当前时刻的X(k|k)作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),如果数字化的电压信号的电压值大于电压变化阈值Vt,则不对当前时刻的数字化的电压信号进行滤波,将当前采集到的数字化电压信号的电压值作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),其中,Vt=N*Vpp,1.5≤N≤2,卡尔曼数字滤波器的Q远小于1,R大于等于1。优选地,Q小于等于0.00001,R大于等于1。更优选地,Q等于0.00001,R等于1。
本发明通过将采集到的数字电压信号与设置的电压变化阀值进行比较,可以分辨出信号与噪声:若采集的电压信号的电压值大于电压变化阈值,则判断采集到的信号更接近于实际的信号,因此不做滤波;若采集的电压信号的电压值小于电压变化阈值,则判断滤波后的信号更接近于实际的信号,因此要使用卡尔曼数字滤波器进行滤波,如此缩短了动态响应时间,提高了响应速率。
图2示出了采用图1的电压信号滤波装置进行滤波的效果示意图。其中,横坐标为时间(ms),纵坐标为电压(V),从图2可以看出,动态响应时间大约为35ms左右,而如果没有本发明的阈值比较步骤,单纯采用卡尔曼数字滤波器进行滤波,如图3所示,动态响应时间为2s左右。图3的横坐标为时间(s),纵坐标为电压(V)。
Claims (6)
1.一种电压信号的滤波方法,其特征在于, 包括:
模数转换器对待滤波的电压信号进行采样,并测量出该电压信号的噪声的峰-峰值Vpp;
模数转换器将数字化的电压信号及测量到的噪声峰-峰值Vpp发送给微处理器;
微处理器根据测量到的噪声峰-峰值Vpp计算电压变化阈值Vt,Vt=N*Vpp,1.5≤N≤2;
微处理器将数字化的电压信号的电压值与该电压变化阈值Vt进行比较,如果数字化的电压信号的电压值小于电压变化阈值Vt,则用卡尔曼数字滤波器对当前时刻的数字化的电压信号进行滤波,并将卡尔曼数字滤波器在当前时刻的X(k|k)作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1);如果数字化的电压信号的电压值大于电压变化阈值Vt,则不对当前时刻的数字化的电压信号进行滤波,将当前采集到的数字化电压信号的电压值作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),其中,所述卡尔曼数字滤波器的Q远小于1,Q小于等于0.00001,R大于等于1;X(k|k)为当前预测结果,X(k-1|k-1) 是上一状态最优的结果。
2.如权利要求1所述的电压信号的滤波方法,其特征在于,如果待滤波的电压信号的电压值超过模数转换器的最大允许采样电压,则在用模数转换器对待滤波的电压信号进行采样之前,先采用电压衰减器对待滤波的电压信号进行电压衰减,使待滤波的电压信号的电压值小于模数转换器的最大允许采样电压。
3.如权利要求1所述的电压信号的滤波方法,其特征在于,所述的卡尔曼数字滤波器通过以下步骤进行滤波:
步骤a、给定初始值X(0);X(0)为0或任意大于0的整数;
步骤b、X(k|k-1)= X(k-1|k-1);X(k|k-1)是利用上一状态预测的结果;
步骤c、P(k|k-1)= P(k-1|k-1) +Q;P(k|k-1)是X(k|k-1)对应的协方差,P(k-1|k-1)是X(k-1|k-1)对应的协方差,Q是系统过程的协方差;
步骤d、Kg(k)= P(k|k-1) / (P(k|k-1) + R);Kg(k)为当前卡尔曼增益,R是测量噪声的协方差;
步骤e、X(k|k)= X(k|k-1)+Kg(k) (Z(k)- X(k|k-1));Z(k)为k时刻的数字化的电压值;
步骤f、P(k|k)=(1-Kg(k))P(k|k-1);P(k|k) 是X(k|k) 对应的协方差。
4.一种电压信号的滤波装置,其特征在于,包括模数转换器和微处理器,所述的模数转换器的输出端与所述微处理器的输入端电连接;其中:
所述的模数转换器用于对待滤波的电压信号进行采样,测量出该电压信号的噪声的峰-峰值Vpp,并将数字化的电压信号及测量到的噪声峰-峰值Vpp发送给微处理器;
所述的微处理器用于根据测量到的噪声峰-峰值Vpp计算电压变化阈值Vt,并将数字化的电压信号的电压值与该电压变化阈值Vt进行比较,如果数字化的电压信号的电压值小于电压变化阈值Vt,则用卡尔曼数字滤波器对当前时刻的数字化的电压信号进行滤波,将卡尔曼数字滤波器在当前时刻的X(k|k)作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),如果数字化的电压信号的电压值大于电压变化阈值Vt,则不对当前时刻的数字化的电压信号进行滤波,将当前采集到的数字化电压信号的电压值作为卡尔曼数字滤波器在下一时刻的X(k-1|k-1),其中,Vt=N*Vpp,1.5≤N≤2,所述卡尔曼数字滤波器的Q远小于1,Q小于等于0.00001,R大于等于1;X(k|k)为当前预测结果,X(k-1|k-1) 是上一状态最优的结果。
5.如权利要求4所述的电压信号的滤波装置,其特征在于,还包括一电压衰减器,所述电压衰减器的输出端与所述模数转换器的输入端电连接,该电压衰减器用于对待滤波的电压信号进行电压衰减,使待滤波的电压信号的电压值小于模数转换器的最大允许采样电压。
6.如权利要求4所述的电压信号的滤波装置,其特征在于,所述的卡尔曼数字滤波器通过以下步骤进行滤波:
步骤a、给定初始值X(0);X(0)为0或任意大于0的整数;
步骤b、X(k|k-1)= X(k-1|k-1);X(k|k-1)是利用上一状态预测的结果;
步骤c、P(k|k-1)= P(k-1|k-1) +Q;P(k|k-1)是X(k|k-1)对应的协方差,P(k-1|k-1)是X(k-1|k-1)对应的协方差,Q是系统过程的协方差;
步骤d、Kg(k)= P(k|k-1) / (P(k|k-1) + R);Kg(k)为当前卡尔曼增益,R是测量噪声的协方差;
步骤e、X(k|k)= X(k|k-1)+Kg(k) (Z(k)- X(k|k-1));Z(k)为k时刻的数字化的电压值;
步骤f、P(k|k)=(1-Kg(k))P(k|k-1);P(k|k) 是X(k|k) 对应的协方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410825332.3A CN104506163B (zh) | 2014-12-27 | 2014-12-27 | 电压信号的滤波方法及其滤波装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410825332.3A CN104506163B (zh) | 2014-12-27 | 2014-12-27 | 电压信号的滤波方法及其滤波装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104506163A CN104506163A (zh) | 2015-04-08 |
CN104506163B true CN104506163B (zh) | 2017-05-17 |
Family
ID=52947887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410825332.3A Active CN104506163B (zh) | 2014-12-27 | 2014-12-27 | 电压信号的滤波方法及其滤波装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104506163B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104868937B (zh) * | 2015-05-26 | 2018-07-17 | 贵州电网有限责任公司电力科学研究院 | 一种远程高速电力线信号采集还原装置及采集和还原方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101800505A (zh) * | 2010-03-12 | 2010-08-11 | 北京航空航天大学 | 一种磁悬浮飞轮转速控制方法 |
CN103941196A (zh) * | 2014-05-07 | 2014-07-23 | 吉林大学 | 锂离子电池荷电状态估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4478584B2 (ja) * | 2005-01-17 | 2010-06-09 | 株式会社ミツトヨ | 位置制御装置、測定装置および加工装置 |
-
2014
- 2014-12-27 CN CN201410825332.3A patent/CN104506163B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101800505A (zh) * | 2010-03-12 | 2010-08-11 | 北京航空航天大学 | 一种磁悬浮飞轮转速控制方法 |
CN103941196A (zh) * | 2014-05-07 | 2014-07-23 | 吉林大学 | 锂离子电池荷电状态估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104506163A (zh) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104614756B (zh) | 一种核素能谱数字化采集系统及其采集方法 | |
CN107015127A (zh) | 一种新能源汽车及其绝缘检测方法、系统 | |
CN103185837B (zh) | 一种电力系统频率测量的方法 | |
CN101030760A (zh) | 一种栅压自适应补偿装置及其栅压自适应补偿方法 | |
CN105301627B (zh) | 一种能谱分析方法、能谱分析系统及伽马射线探测系统 | |
KR20230096130A (ko) | 배터리 셀의 자기 방전 전류 검출 방법, 장치, 기기 및 컴퓨터 저장 매체 | |
WO2006031334A2 (en) | A self-powered vehicle speed sensor | |
CN104506163B (zh) | 电压信号的滤波方法及其滤波装置 | |
CN104679375A (zh) | 一种优化信噪比参数的方法及装置 | |
CN102722627A (zh) | 一种电力一次系统数字物理混合仿真通用接口系统 | |
CN106612084B (zh) | 电流采样处理电路及采样处理方法 | |
CN103634003A (zh) | Osc频率自动校准电路及自动校准方法 | |
CN104168067A (zh) | 一种判断光接收电路中的光功率信号强度的方法及其电路 | |
CN104518755A (zh) | 数字电路噪声滤波器及数字滤波方法 | |
CN107515329B (zh) | 一种基于数字滤波的igbt综合过电流保护方法及系统 | |
CN102288835A (zh) | 信号检测电路及方法 | |
CN112904431B (zh) | 直线电机牵引系统车辆无感应板区的检测方法、装置、系统及介质 | |
CN109450755B (zh) | 增强存储模块总线接口的稳定性余量的方法及系统 | |
CN114034973A (zh) | 一种配电线路接地故障的故障区域识别方法及其装置和系统 | |
CN104820135B (zh) | 电抗器测试方法及测试装置 | |
CN110412356B (zh) | 一种新能源汽车控制器主电路杂散电感测量方法及装置 | |
CN109660256A (zh) | 固定参数的a/d误差修正电路和修正方法 | |
CN103105513A (zh) | 消除被检设备输出随机大信号对标准设备损害的保护装置 | |
CN219498934U (zh) | 一种剩余电流动作保护器及电子设备 | |
CN114141602B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 201203 China Zuchongzhi Road (Shanghai) Free Trade Zone No. 2388 building 1-2 Co-patentee after: Zhejiang Keboda Industrial Co., Ltd. Patentee after: Koboda technology, Limited by Share Ltd Address before: 201203 Shanghai city Zuchongzhi road in Zhangjiang High Tech Park of Pudong New Area No. 2388 building 1-2 Co-patentee before: Zhejiang Keboda Industrial Co., Ltd. Patentee before: Keboda Technology Corp. |