CN107340405A - 一种气固两相流交流电信号递归最小二乘自适应滤波方法 - Google Patents
一种气固两相流交流电信号递归最小二乘自适应滤波方法 Download PDFInfo
- Publication number
- CN107340405A CN107340405A CN201710439594.XA CN201710439594A CN107340405A CN 107340405 A CN107340405 A CN 107340405A CN 201710439594 A CN201710439594 A CN 201710439594A CN 107340405 A CN107340405 A CN 107340405A
- Authority
- CN
- China
- Prior art keywords
- signal
- dual
- phrase distribution
- gas olid
- upstream
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/08—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring variation of an electric variable directly affected by the flow, e.g. by using dynamo-electric effect
- G01P5/083—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring variation of an electric variable directly affected by the flow, e.g. by using dynamo-electric effect by using electronic circuits for measuring the dynamoelectric effect
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H21/0043—Adaptive algorithms
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H2017/0072—Theoretical filter design
- H03H2017/0081—Theoretical filter design of FIR filters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H21/0043—Adaptive algorithms
- H03H2021/0049—Recursive least squares algorithm
- H03H2021/005—Recursive least squares algorithm with forgetting factor
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H2021/0085—Applications
- H03H2021/0089—System identification, i.e. modeling
- H03H2021/009—System identification, i.e. modeling with recursive filters
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明公开了气固两相流交流电信号递归最小二乘自适应滤波方法,步骤1:采集上、下游气固两相流交流电信号X1(k)、X2(k);步骤2:运用递归最小二乘自适应滤波方法对步骤1采集的上、下游气固两相流交流电信号X1(k)、X2(k)进行气固两相流交流电信号递归最小二乘自适应滤波处理,得到滤波后的上、下游信号Y1(k)、Y2(k);步骤3:对步骤2滤波后的上、下游信号Y1(k)、Y2(k)进行互相关测速运算,得出气固两相流固相颗粒流速的测量值。本发明应用到气固两相流交流电信号处理上,运用该方法处理后得到的上、下游气固两相流交流电信号进行互相关测速算法运算分析,使得所测气固两相流参数更准确、更稳定。
Description
技术领域
本发明属于信号检测技术领域,具体涉及一种气固两相流交流电信号递归最小二乘自适应滤波方法。
背景技术
气固两相流流动参数主要包括气固两相流固体颗粒的流速、浓度、质量流量等,气固两相流参数的准确测量对提高气力输送相关行业以及部门的生产安全性、生产效率,以及降低能耗、节约能源具有十分积极的意义。
静电传感器所检测到的气固两相流信号是微伏级别的微弱信号,复杂、恶劣的工业现场环境中,微弱信号极易受噪声干扰,导致气固两相流交流电信号信噪比偏低,运用互相关测速算法分析计算出来的气固两相流固体颗粒流速精度低,运用传统的滤波方法对气固两相流交流电信号进行处理,效果不太理想,并且静电传感器检测到的气固两相流信号先验知识未知,噪声复杂、多变,针对这类信号特征,寻求一种适用于气固两相流交流电信号的滤波方法,对于气固两相流参数的准确、稳定测量有着重要的意义。
发明内容
本发明的目的是提供一种气固两相流交流电信号递归最小二乘自适应滤波方法,以解决目前气固两相流参数测量精度低的问题。
本发明采用以下技术方案:一种气固两相流交流电信号递归最小
二乘自适应滤波方法,包括以下步骤:
步骤1:采集上、下游气固两相流交流电信号X1(k)、X2(k);
步骤2:运用递归最小二乘自适应滤波方法对步骤1采集的上、下游气固两相流交流电信号X1(k)、X2(k)进行气固两相流交流电信号递归最小二乘自适应滤波处理,得到滤波后的上、下游信号Y1(k)、Y2(k);
步骤3:对步骤2滤波后的上、下游信号Y1(k)、Y2(k)进行互相关测速运算,得出气固两相流固相颗粒流速的测量值。
作为本发明进一步的方案,所述步骤1利用上、下游静电传感器
所检测到的荷电信号经过信号调理模块、信号采集模块,得到上、
下游气固两相流交流电信号X1(k)、X2(k)。
作为本发明进一步的方案,下游气固两相流交流电信号X2(k)的处理与上游气固两相流交流电信号X1(k)处理方法相同,以上游气固两相流交流电信号X1(k)作为滤波信号进行说明,所述步骤2具体方法为:
步骤2.1:设置滤波器阶数N、数据长度n、遗忘因子λ、初始权值序列w0(k)、自相关矩阵逆矩阵SD(0);
步骤2.2:计算k时刻有限长单位冲激响应(FIR)数字滤波器输出信号y(k),y(k)=X1(k)TW(k),其中,X1(k)=[x(k) x(k-1) L x(k-N+1)]Τ为k时刻气固两相流交流电信号序列的前N个值,W(k)=[w0(k) w1(k) L wN-1(k)]Τ是k时刻滤波器的权值序列;
步骤2.3:计算k时刻的误差e(k),即k时刻气固两相流交流电信号与FIR滤波器输出信号之差e(k)=x(k)-y(k);
步骤2.4:按照权值序列更新公式w(k+1)=w(k)+e(k)SD(k)X1(k)更新下一时刻权值序列;
步骤2.5:根据自相关矩阵逆矩阵更新公式
更新k+1时刻的自相关矩阵逆矩阵SD(k+1);
步骤2.6:k=k+1,若k≤n,返回步骤2.2,否则,完成滤波,
得到滤波后的信号Y1(k)。
本发明的有益效果:本发明应用到气固两相流交流电信号处理上,与传统的滤波方法相比较,递归最小二乘自适应滤波方法对于气固两相流交流电信号的处理具有较强针对性,运用该方法处理后得到的上、下游气固两相流交流电信号进行互相关测速算法运算分析,使得所测气固两相流参数更准确、更稳定。
附图说明
图1是本发明中气固两相流交流电信号递归最小二乘自适应滤波器基本结构示意图;
图2是本发明中气固两相流交流电信号递归最小二乘自适应滤波方法对气固两相流交流电信号的处理流程图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种气固两相流交流电信号递归最小二乘自适应滤波方法,如图1所示,x(k),x(k-1),x(k-2),…,x(k-N+1)是k时刻气固两相流交流电信号序列X(k)的前N个值,y(k)是x(k),x(k-1),x(k-2),…,x(k-N+1)经过有限长单位冲激响应(FIR)滤波器得到的输出信号,w0(k),w1(k),…,wN-1(k)是k时刻滤波器权值,误差信号e(k)是k时刻气固两相流交流电信号x(k)与FIR滤波器输出信号y(k)之差,根据递归最小二乘的目标函数推导得到下一时刻权值序列更新公式w(k+1)=w(k)+e(k)SD(k)X1(k),通过权值更新公式来调节FIR滤波器的权值序列w0(k),w1(k),…,wN-1(k),以达到最佳滤波,具体按照以下步骤实施:
步骤1:采集上、下游气固两相流交流电信号X1(k)、X2(k):上、下游静电传感器所检测到的荷电信号经过信号调理模块、信号采集模块,得到上、下游气固两相流交流电信号X1(k)、X2(k);
步骤2:运用递归最小二乘自适应滤波方法对步骤1采集的上、下游气固两相流交流电信号X1(k)、X2(k)进行气固两相流交流电信号递归最小二乘自适应滤波处理,得到滤波后的信号Y1(k)、Y2(k),以上游信号X1(k)作为滤波信号进行说明,下游信号X2(k)的处理同上游信号处理方法一样,具体流程图如图2所示:
步骤2.1:设置滤波器阶数N、数据长度n、遗忘因子λ、初始权值序列w0(0)、自相关矩阵逆矩阵SD(0);
步骤2.2:计算k时刻有限长单位冲激响应(FIR)数字滤波器输出信号y(k),y(k)=X1(k)TW(k),其中,X1(k)=[x(k) x(k-1) … x(k-N+1)]Τ为k时刻气固两相流交流电信号序列的前N个值,W(k)=[w0(k) w1(k) … wN-1(k)]Τ是k时刻滤波器的权值序列;
步骤2.3:计算k时刻的误差e(k),即k时刻气固两相流交流电信号与FIR滤波器输出信号之差e(k)=x(k)-y(k);
步骤2.4:按照权值序列更新公式w(k+1)=w(k)+e(k)SD(k)X1(k)更新下一时刻权值序列;
步骤2.5:根据自相关矩阵逆矩阵更新公式
更新k+1时刻的自相关矩阵逆矩阵SD(k+1);
步骤2.6:k=k+1,若k≤n,返回步骤2.2,否则,完成滤波,得到滤波后的信号Y1(k)。
自适应滤波器具有学习和跟踪的能力,气固两相流交流电信号先验知识未知,自适应滤波器通过自身的学习能力,获得最优权值,实现最佳滤波,在噪声随着周围的环境变化而变化时,自适应滤波器通过自身具有的跟踪能力,调节滤波器权值,最终实现最佳滤波。
步骤3:对步骤2滤波后的上、下游信号Y1(k)、Y2(k)进行互相关测速运算,得出气固两相流固相颗粒流速的测量值。
本发明具有快速的收敛速度和跟踪的能力,对气固两相流交流电信号的处理具有很强的针对性,通过对气固两相流交流电信号进行递归最小二乘自适应滤波处理,提高信号的信噪比,使用该方法处理后的信号进行运算分析,能够提高气固两相流参数测量准确性以及稳定性。
以上所述为本发明较佳实施例,对于本领域的普通技术人员而言,根据本发明的教导,在不脱离本发明的原理与精神的情况下,对实施方式所进行的改变、修改、替换和变型仍落入本发明的保护范围之内。
Claims (3)
1.一种气固两相流交流电信号递归最小二乘自适应滤波方法,其特征在于,包括以下步骤:
步骤1:采集上、下游气固两相流交流电信号X1(k)、X2(k);
步骤2:运用递归最小二乘自适应滤波方法对步骤1采集的上、下游气固两相流交流电信号X1(k)、X2(k)进行气固两相流交流电信号递归最小二乘自适应滤波处理,得到滤波后的上、下游信号Y1(k)、Y2(k);
步骤3:对步骤2滤波后的上、下游信号Y1(k)、Y2(k)进行互相关测速运算,得出气固两相流固相颗粒流速的测量值。
2.根据权利要求1所述的一种气固两相流交流电信号递归最小二乘自适应滤波方法,其特征在于,所述步骤1利用上、下游静电传感器所检测到的荷电信号经过信号调理模块、信号采集模块,得到上、下游气固两相流交流电信号X1(k)、X2(k)。
3.根据权利要求1所述的一种气固两相流交流电信号递归最小二乘自适应滤波方法,其特征在于,下游气固两相流交流电信号X2(k)的处理与上游气固两相流交流电信号X1(k)处理方法相同,以上游气固两相流交流电信号X1(k)作为滤波信号进行说明,所述步骤2具体方法为:
步骤2.1:设置滤波器阶数N、数据长度n、遗忘因子λ、初始权值序列w0(k)、自相关矩阵逆矩阵SD(0);
步骤2.2:计算k时刻有限长单位冲激响应(FIR)数字滤波器输出信号y(k),y(k)=X1(k)TW(k),其中,X1(k)=[x(k) x(k-1) L x(k-N+1)]Τ为k时刻气固两相流交流电信号序列的前N个值,W(k)=[w0(k) w1(k) L wN-1(k)]Τ是k时刻滤波器的权值序列;
步骤2.3:计算k时刻的误差e(k),即k时刻气固两相流交流电信号与FIR滤波器输出信号之差e(k)=x(k)-y(k);
步骤2.4:按照权值序列更新公式w(k+1)=w(k)+e(k)SD(k)X1(k)更新下一时刻权值序列;
步骤2.5:根据自相关矩阵逆矩阵更新公式
更新k+1时刻的自相关矩阵逆矩阵SD(k+1);
步骤2.6:k=k+1,若k≤n,返回步骤2.2,否则,完成滤波,得到滤波后的信号Y1(k)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710439594.XA CN107340405A (zh) | 2017-06-12 | 2017-06-12 | 一种气固两相流交流电信号递归最小二乘自适应滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710439594.XA CN107340405A (zh) | 2017-06-12 | 2017-06-12 | 一种气固两相流交流电信号递归最小二乘自适应滤波方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107340405A true CN107340405A (zh) | 2017-11-10 |
Family
ID=60219965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710439594.XA Pending CN107340405A (zh) | 2017-06-12 | 2017-06-12 | 一种气固两相流交流电信号递归最小二乘自适应滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107340405A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828261A (zh) * | 2018-07-04 | 2018-11-16 | 天津大学 | 用于气固两相流静电相关测速的调理电路及信号复原方法 |
CN110379440A (zh) * | 2019-07-19 | 2019-10-25 | 宁波奥克斯电气股份有限公司 | 语音降噪方法、装置、语音空调及计算机可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975799A (zh) * | 2010-10-29 | 2011-02-16 | 西安科瑞自动化有限责任公司 | 交流电荷感应式煤粉浓度测量装置 |
CN103078640A (zh) * | 2013-01-16 | 2013-05-01 | 华南理工大学 | 一种用于adc的rls自适应滤波校准算法 |
CN103776743A (zh) * | 2013-12-24 | 2014-05-07 | 中电投河南电力有限公司技术信息中心 | 交流电荷感应式煤粉细度在线监测方法及装置 |
CN105929238A (zh) * | 2016-04-13 | 2016-09-07 | 西安理工大学 | 一种气固两相流交流电信号最小均方误差自适应滤波方法 |
-
2017
- 2017-06-12 CN CN201710439594.XA patent/CN107340405A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975799A (zh) * | 2010-10-29 | 2011-02-16 | 西安科瑞自动化有限责任公司 | 交流电荷感应式煤粉浓度测量装置 |
CN103078640A (zh) * | 2013-01-16 | 2013-05-01 | 华南理工大学 | 一种用于adc的rls自适应滤波校准算法 |
CN103776743A (zh) * | 2013-12-24 | 2014-05-07 | 中电投河南电力有限公司技术信息中心 | 交流电荷感应式煤粉细度在线监测方法及装置 |
CN105929238A (zh) * | 2016-04-13 | 2016-09-07 | 西安理工大学 | 一种气固两相流交流电信号最小均方误差自适应滤波方法 |
Non-Patent Citations (1)
Title |
---|
许国威 等: "自适应滤波RLS算法研究与DSP实现", 《计算机技术与发展》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828261A (zh) * | 2018-07-04 | 2018-11-16 | 天津大学 | 用于气固两相流静电相关测速的调理电路及信号复原方法 |
CN110379440A (zh) * | 2019-07-19 | 2019-10-25 | 宁波奥克斯电气股份有限公司 | 语音降噪方法、装置、语音空调及计算机可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105067061B (zh) | 一种基于静电与电容传感器阵列的粉体质量流量测量装置及方法 | |
CN106679741B (zh) | 基于涡街流量计抗干扰信号的处理方法及系统 | |
CN107271768A (zh) | 一种最小二乘拟合动态频率测量方法 | |
Gao et al. | Improved Oustaloup approximation of fractional-order operators using adaptive chaotic particle swarm optimization | |
CN103808380B (zh) | 一种用于超声波流量计量仪表的流量快速跟踪方法 | |
CN104535890B (zh) | 一种电网电压跌落的快速检测方法 | |
CN107340405A (zh) | 一种气固两相流交流电信号递归最小二乘自适应滤波方法 | |
CN103344825A (zh) | 基于交流采样的电能计量系统 | |
Doroshin et al. | Open research issues on nonlinear dynamics, dynamical systems and processes | |
CN103675629A (zh) | 基于电场特征量的电极起晕电压预测方法 | |
CN107102189A (zh) | 基于s函数的变步长lms谐波电流检测方法 | |
Kellman et al. | Barker-Coded node-pore resistive pulse sensing with built-in coincidence correction | |
CN104539264B (zh) | 应用于eps电源系统的滤波方法及滤波电路 | |
CN105894014B (zh) | 基于多因素不一致度量的异常行为序贯检测方法 | |
CN106767362A (zh) | 一种电容式微位移传感器检测电路 | |
CN107356783A (zh) | 一种基于小波变换自适应滤波的气固两相流交流电信号处理方法 | |
CN105929238A (zh) | 一种气固两相流交流电信号最小均方误差自适应滤波方法 | |
CN107228886A (zh) | 水的电导率电阻率检测的测量装置和方法 | |
CN106094972B (zh) | 一种基于函数模型的光伏发电系统最大功率点跟踪方法 | |
CN105606893B (zh) | 基于空间平滑修正music的电力间谐波检测方法 | |
CN208805238U (zh) | 一种抗振型双涡街传感器流体控制装置 | |
CN106323589A (zh) | 一种密相气力输送两相流流型的归类方法 | |
CN105676143A (zh) | 蓄电池出厂参数在线检测装置 | |
CN106767952B (zh) | 一种电感式位移传感器的干扰消除方法 | |
CN109839412A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171110 |
|
RJ01 | Rejection of invention patent application after publication |