CN104779934A - 一种用于血压和心率测试的rc滤波方法及rc滤波器 - Google Patents

一种用于血压和心率测试的rc滤波方法及rc滤波器 Download PDF

Info

Publication number
CN104779934A
CN104779934A CN201510179905.4A CN201510179905A CN104779934A CN 104779934 A CN104779934 A CN 104779934A CN 201510179905 A CN201510179905 A CN 201510179905A CN 104779934 A CN104779934 A CN 104779934A
Authority
CN
China
Prior art keywords
preliminary treatment
treatment data
factorial
summation operation
digital 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
Application number
CN201510179905.4A
Other languages
English (en)
Other versions
CN104779934B (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.)
Hefei Sijing Electronic Co ltd
Original Assignee
HEROCHIP MICROELECTRONICS CO Ltd
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 HEROCHIP MICROELECTRONICS CO Ltd filed Critical HEROCHIP MICROELECTRONICS CO Ltd
Priority to CN201510179905.4A priority Critical patent/CN104779934B/zh
Publication of CN104779934A publication Critical patent/CN104779934A/zh
Application granted granted Critical
Publication of CN104779934B publication Critical patent/CN104779934B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种用于血压和心率测试的RC滤波方法及RC滤波器,其中,该方法包括:步骤A:采集待处理的模拟信号,将所述模拟信号转换为数字信号,得到数字信号的原始数据;步骤B:对所述数字信号的原始数据进行滤除高频处理,得到滤除高频后的初步处理数据;步骤C:存储前N个所述初步处理数据,其中N为大于等于20的正整数;步骤D:根据所述N个初步处理数据对应的滤波系数,对所述N个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;步骤E:将N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。本发明RC滤波方法及RC滤波器,剔除了外界因素影响,提高了精度与灵活性,增加了处理带宽。

Description

一种用于血压和心率测试的RC滤波方法及RC滤波器
技术领域
本发明涉及心跳检测的高频滤除技术领域,具体地,涉及一种用于血压和心率测试的RC滤波方法及RC滤波器。
背景技术
在滤波器设计中,由于受到各种因素的影响,滤波器的滤频特性(滤除高频或低频的能力)无法达到理论值,因此在使用硬件电路来设计滤波器就增加了一定的难度。
首先是RC滤波器的驱动能力问题,由于驱动能力较弱,所以必须是有源RC滤波器,这样就增加了电路的功耗。其次,由于电路之间的干扰、电磁场的干扰以及外部时钟的干扰等等,都会使RC滤波器的滤频特性发生改变。最后是温度问题,随着工作时间的推移,温度的升高,RC滤波器的滤频特性也会发生改变。
现有技术中还将RC滤波技术应用到血压和心率测试技术领域,因为心跳声是一种低频信号,所以也可以通过采集心跳信号,滤除高频的方法获取滤频后的心跳信号。但在实际的调试过程,需要根据具体情况更换电阻电容,而且由于外界因素的影响,滤波特性并不是理论值,需要不断更改RC值来进行匹配,这样就增加了一定的麻烦。
发明内容
为了解决现有技术中在进行血压和心率测试时RC滤波器滤波特性不高的问题,本发明提出了一种用于血压和心率测试的RC滤波方法及RC滤波器。其中,该用于血压和心率测试的RC滤波方法包括:
步骤A:采集待处理的模拟信号,将所述模拟信号转换为数字信号,得到数字信号的原始数据;
步骤B:对所述数字信号的原始数据进行滤除高频处理,得到滤除高频后的初步处理数据;
步骤C:存储前N个所述初步处理数据,其中N为大于等于20的正整数;
步骤D:根据所述N个初步处理数据对应的滤波系数,对所述N个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
步骤E:将步骤D中的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
本发明的用于血压和心率测试的RC滤波方法,剔除了外界环境因素的影响,提高了滤波的精度与灵活性,增加了滤波的处理带宽,降低了滤波的成本。不需要再考虑电路功耗问题、外界干扰对滤波特性的影响、温度对滤波特性的影响等因素。而且该RC滤波器设计灵活,发现RC值不满足滤波要求时,只需简单更改滤波系数即可完成RC更换,方便简单易受控制,可以达到滤波特性的理论值的要求。
在一个实施例中,还包括:
步骤F:存储第N+1个初步处理数据,并删除所述前N个初步处理数据中的第一个;
步骤G:对第2个至第N+1个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
步骤H:将步骤G中的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
本实施例还考虑到了在存储的初步处理数据增加时,每新进一个初步处理数据时,将原N个数据中的第一个删除,将剩余的N-1个数据与该新进的数据作为后续N阶乘积求和运算的基础,充分考虑到了根据数字信号的原始数据的变化进行及时更新。
在一个实施例中,所述步骤B具体包括:
根据所述数字信号的原始数据的频率,确定一个基数,所述基数为所述频率的三分之二,并随机设置一个放大系数,所述放大系数为正数;
将所述数字信号的原始数据的频率减去所述基数再乘以所述放大系数,并进行滤除高频处理,根据数字信号的滤除结果加大或减小所述基数和放大系数,直至滤出的数字信号能够显示低频信号,将此时的数据作为所述初步处理数据。
在本实施例中,首选根据原始数据的频率,确定一个基数,通常的,该基数为频率的三分之二,并随机设置一个放大系数。其中基数和系数的选择是在实际运用中调试完成的,基数先选择原始数据的三分之二根据结果再逼近理想基数,系数先随意设定根据结果再逼近理想系数,如果滤出的信号不够明显(即无法显示出低频信号),则可以适当加大基数加大系数或者减小基数减小系数,慢慢调试直至最佳,能够清晰的看出低频信号时为止。
在一个实施例中,所述步骤D具体包括:
计算所述N个初步处理数据分别对应的滤波系数;
将所述N个初步处理数据分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到N阶乘积求和运算的结果;
将所述N个初步处理数据中的前20个分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到20阶乘积求和运算的结果。
本实施例对乘积求和的运算过程进行了限定,能够快速准确的能到对初步处理数据进行N阶乘积求和运算和20阶求和运算的结果。
在一个实施例中,每个时间点对应一个初步处理数据,利用RC放电公式计算,RC常数对应的截止频率即为通过的最高频率,公式中每个时间点通过电容的电压值,即为每个时间点对应初步处理数据的滤波系数。
本发明的用于血压和心率测试的RC滤波器包括:
采集转换模块,用于采集待处理的模拟信号,将所述模拟信号转换为数字信号,得到数字信号的原始数据;
初步处理模块,用于对所述数字信号的原始数据进行滤除高频处理,得到滤除高频后的初步处理数据;
第一存储模块,用于存储前N个所述初步处理数据,其中N为大于等于20的正整数;
第一运算模块,用于根据所述N个初步处理数据对应的滤波系数,对所述N个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
第一滤频模块,用于将第一运算模块计算的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
在一个实施例中,还包括:
第二存储模块,用于存储第N+1个初步处理数据,并删除所述前N个初步处理数据中的第一个;
第二运算模块,用于对第2个至第N+1个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
第二滤频模块,用于将第二运算模块计算的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
在一个实施例中,所述初步处理模块包括:
参数设定子模块,用于根据所述数字信号的原始数据的频率,确定一个基数,所述基数为所述频率的三分之二,并随机设置一个放大系数,所述放大系数为正数;
信号滤除子模块,用于将所述数字信号的原始数据的频率减去所述基数再乘以所述放大系数,并进行滤除高频处理,根据数字信号的滤除结果加大或减小所述基数和放大系数,直至滤出的数字信号能够显示低频信号,将此时的数据作为所述初步处理数据。
在一个实施例中,所述第一运算模块包括:
系数计算子模块,用于计算所述N个初步处理数据分别对应的滤波系数;
N阶运算子模块,用于将所述N个初步处理数据分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到N阶乘积求和运算的结果;
20阶运算子模块,用于将所述N个初步处理数据中的前20个分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到20阶乘积求和运算的结果。
在一个实施例中,每个时间点对应一个初步处理数据,利用RC放电公式计算,RC常数对应的截止频率即为通过的最高频率,公式中每个时间点通过电容的电压值,即为每个时间点对应初步处理数据的滤波系数。
本发明的用于血压和心率测试的RC滤波器,剔除了外界环境因素的影响,提高了滤波的精度与灵活性,增加了滤波的处理带宽,降低了滤波的成本。不需要再考虑电路功耗问题、外界干扰对滤波特性的影响、温度对滤波特性的影响等因素。而且该RC滤波器设计灵活,发现RC值不满足滤波要求时,只需简单更改滤波系数即可完成RC更换,方便简单易受控制,可以达到滤波特性的理论值的要求。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1为本发明实施例的一种用于血压和心率测试的RC滤波方法的流程图;
图2为本发明实施例的另一种用于血压和心率测试的RC滤波方法的流程图;
图3为本发明实施例的一种用于血压和心率测试的RC滤波方法中步骤S102的流程图;
图4为本发明实施例的一种用于血压和心率测试的RC滤波方法中步骤S104的流程图;
图5为本发明实施例的一种用于血压和心率测试的RC滤波器的结构示意图;
图6为本发明实施例的一种用于血压和心率测试的RC滤波器的初步处理模块的结构示意图;
图7为本发明实施例的一种用于血压和心率测试的RC滤波器的第一运算模块的结构示意图
图8为本发明实施例的60阶乘积求和运算和20阶乘积求和运算的示意图。
具体实施方式
下面结合附图,对本发明的具体实施方式进行详细描述,但应当理解本发明的保护范围并不受具体实施方式的限制。
为了解决现有技术中在进行血压和心率测试时RC滤波器滤波特性不高的问题,本发明提出了一种用于血压和心率测试的RC滤波方法及RC滤波器。
实施例一
如图1所示,本实施例的用于血压和心率测试的RC滤波方法包括:
步骤S101:采集待处理的模拟信号,将所述模拟信号转换为数字信号,得到数字信号的原始数据;
具体的,可以通过SPI协议对需要处理的心跳模拟信号进行采样,并转换成数字信号。
步骤S102:对所述数字信号的原始数据进行滤除高频处理,得到滤除高频后的初步处理数据;
具体的,就是滤除数字信号的中的高频信号,以至能够清晰地显示低频信号。
步骤S103:存储前N个所述初步处理数据,其中N为大于等于20的正整数;
具体的,数据应当存储在单片机资源中。在实际调试中,20是最低标准,小于20滤波效果就无法体现出来,因此,N≥20。
步骤S104:根据所述N个初步处理数据对应的滤波系数,对所述N个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
具体的,本实施例中的乘积求和运算类似于卷积运算。卷积是分析数学中一种重要的运算。设:f(x),g(x)是R1上的两个可积函数,作积分:
∫ - ∞ ∞ f ( τ ) g ( x - τ ) dτ
可以证明,关于几乎所有的实数x,上述积分是存在的。这样,随着x的不同取值,这个积分就定义了一个新函数h(x),称为函数f与g的卷积,记为h(x)=(f*g)(x)。容易验证,(f*g)(x)=(g*f)(x),并且(f*g)(x)仍为可积函数。这就是说,把卷积代替乘法,L1(R1)空间是一个代数,甚至是巴拿赫代数。
卷积与傅里叶变换有着密切的关系。利用一点性质,即两函数的傅里叶变换的乘积等于它们卷积后的傅里叶变换,能使傅里叶分析中许多问题的处理得到简化。由卷积得到的函数f*g一般要比f和g都光滑。特别当g为具有紧致集的光滑函数,f为局部可积时,它们的卷积f*g也是光滑函数。利用这一性质,对于任意的可积函数f,都可以简单地构造出一列逼近于f的光滑函数列fs,这种方法称为函数的光滑化或正则化。
步骤S105:将步骤S104中计算得出的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
步骤S104和步骤S105对乘积求和运算的过程及得到最终滤频后的信号的过程进行了具体限定,假设N取60,则前N个初步处理数据分别为数据1,数据2,……,数据60,其对应的滤波系数分别为α1,α2,……,α60。则步骤S104中的N阶乘积求和运算过程为数据1×α1,数据2×α2,……,数据60×α60,并将上述各乘积相加。20阶乘积求和运算过程为数据1×α1,数据2×α2,……,数据20×α20,并将上述各乘积相加。步骤S105的过程即为将上述60阶乘积求和运算的结果减去20阶乘积求和运算的结果,如图8所示为乘积求和运算的示意图。
本发明实施例的用于血压和心率测试的RC滤波方法,剔除了外界环境因素的影响,提高了滤波的精度与灵活性,增加了滤波的处理带宽,降低了滤波的成本。不需要再考虑电路功耗问题、外界干扰对滤波特性的影响、温度对滤波特性的影响等因素。而且该RC滤波器设计灵活,发现RC值不满足滤波要求时,只需简单更改滤波系数即可完成RC更换,方便简单易受控制,可以达到滤波特性的理论值的要求。
实施例二
如图2所示,当存储的初步处理数据增加时,该用于血压和心率测试的RC滤波方法还包括:
步骤S106:存储第N+1个初步处理数据,并删除所述前N个初步处理数据中的第一个;
具体的,利用堆栈的思想,每采集到一个新的原始数据时,便删除第一个进来并存储的数据,接着放大处理当前原始数据后存储起来,此处需要使用到单片机的资源,单片机的存储空间与处理速度完全能满足需求,数据应当存储在单片机资源中。
步骤S107:对第2个至第N+1个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
具体的,当每增加一个初步处理时间时,例如,原N个数据分别为1,2,3……,60,当第61个初步处理数据存储进单片机资源中时,删除第1个数据,对第2,3,4……,61总共60个数据分别进行60阶乘积求运算和20阶乘积求和运算。
步骤S108:将步骤S107中计算得出的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
实施例三
本实施例是对步骤S102的具体描述,如图3所示,包括:
步骤S301:根据所述数字信号的原始数据的频率,确定一个基数,所述基数为所述频率的三分之二,并随机设置一个放大系数,所述放大系数为正数;
步骤S302:将所述数字信号的原始数据的频率减去所述基数再乘以所述放大系数,并进行滤除高频处理,根据数字信号的滤除结果加大或减小所述基数和放大系数,直至滤出的数字信号能够显示低频信号,将此时的数据作为所述初步处理数据。
实施例四
本实施例是对步骤S104的具体描述,如图4所示,包括:
步骤S401:计算所述N个初步处理数据分别对应的滤波系数;
其中,每个时间点(时刻)对应一个初步处理数据,利用RC放电公式计算,RC时间常数对应的截止频率即为通过的最高频率,公式中每个时间点通过电容的电压值,即为每个时间点对应初步处理数据的滤波系数。如下式中Uc(0)设置一个初始值即电容两端电压(此值与选择的阶数N有关,阶数越高初始电压值越高,假设N为60阶,那就要保证在60个时刻后电容电压基本全部衰减,一个时刻为一秒,以此来设定各个时刻的初始电压值),经实践在心跳声检测的高频领域,这样的设置计算方法尤为有效,RC放电公式如下:
Uc(t)=[Uc(0)]e^-t/RC,
也即Uc(t)=[Uc(0)]*exp(-t/RC),其中,Uc(t)为每个时间点通过电容的电压值,也即本发明实施例的滤波系数,Uc(0)为电容两端的初始电压,t为放电时间。
步骤S402:将所述N个初步处理数据分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到N阶乘积求和运算的结果;
步骤S403:将所述N个初步处理数据中的前20个分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到20阶乘积求和运算的结果。
实施例五
本实施例提供了一种用于血压和心率测试的RC滤波器,如图5所示,包括:
采集转换模块51,用于采集待处理的模拟信号,将所述模拟信号转换为数字信号,得到数字信号的原始数据;
初步处理模块52,用于对所述数字信号的原始数据进行滤除高频处理,得到滤除高频后的初步处理数据;
第一存储模块53,用于存储前N个所述初步处理数据,其中N为大于等于20的正整数;
第一运算模块54,用于根据所述N个初步处理数据对应的滤波系数,对所述N个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
第一滤频模块55,用于将第一运算模块54计算的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
第二存储模块56,用于存储第N+1个初步处理数据,并删除所述前N个初步处理数据中的第一个;
第二运算模块57,用于对第2个至第N+1个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
第二滤频模块58,用于将第二运算模块57计算的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
在一个实施例中,如图6所示,所述初步处理模块52包括:
参数设定子模块521,用于根据所述数字信号的原始数据的频率,确定一个基数,所述基数为所述频率的三分之二,并随机设置一个放大系数,所述放大系数为正数;
信号滤除子模块522,用于将所述数字信号的原始数据的频率减去所述基数再乘以所述放大系数,并进行滤除高频处理,根据数字信号的滤除结果加大或减小所述基数和放大系数,直至滤出的数字信号能够显示低频信号,将此时的数据作为所述初步处理数据。
在一个实施例中,如图7所示,所述第一运算模块54包括:
系数计算子模块541,用于计算所述N个初步处理数据分别对应的滤波系数;
N阶运算子模块542,用于将所述N个初步处理数据分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到N阶乘积求和运算的结果;
20阶运算子模块543,用于将所述N个初步处理数据中的前20个分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到20阶乘积求和运算的结果。
在一个实施例中,每个时间点(时刻)对应一个初步处理数据,利用RC放电公式计算,RC常数对应的截止频率即为通过的最高频率,公式中每个时间点通过电容的电压值,即为每个时间点对应初步处理数据的滤波系数。如下式中Uc(0)设置一个初始值即电容两端电压(此值与选择的阶数N有关,阶数越高初始电压值越高,假设N为60阶,那就要保证在60个时刻后电容电压基本全部衰减,一个时刻为一秒,以此来设定各个时刻的初始电压值),经实践在心跳声检测的高频领域,这样的设置计算方法尤为有效,RC放电公式如下:
Uc(t)=[Uc(0)]e^-t/RC,
也即Uc(t)=[Uc(0)]*exp(-t/RC),其中,Uc(t)为每个时间点通过电容的电压值,也即本发明实施例的滤波系数,Uc(0)为电容两端的初始电压,t为放电时间。
本发明实施例的用于血压和心率测试的RC滤波器,剔除了外界环境因素的影响,提高了滤波的精度与灵活性,增加了滤波的处理带宽,降低了滤波的成本。不需要再考虑电路功耗问题、外界干扰对滤波特性的影响、温度对滤波特性的影响等因素。而且该RC滤波器设计灵活,发现RC值不满足滤波要求时,只需简单更改滤波系数即可完成RC更换,方便简单易受控制,可以达到滤波特性的理论值的要求。
本发明能有多种不同形式的具体实施方式,上面以图1-图8为例结合附图对本发明的技术方案作举例说明,这并不意味着本发明所应用的具体实例只能局限在特定的流程或实施例结构中,本领域的普通技术人员应当了解,上文所提供的具体实施方案只是多种优选用法中的一些示例,任何体现本发明权利要求的实施方式均应在本发明技术方案所要求保护的范围之内。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种用于血压和心率测试的RC滤波方法,其特征在于,包括:
步骤A:采集待处理的模拟信号,将所述模拟信号转换为数字信号,得到数字信号的原始数据;
步骤B:对所述数字信号的原始数据进行滤除高频处理,得到滤除高频后的初步处理数据;
步骤C:存储前N个所述初步处理数据,其中N为大于等于20的正整数;
步骤D:根据所述N个初步处理数据对应的滤波系数,对所述N个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
步骤E:将步骤D中的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
2.根据权利要求1所述的用于血压和心率测试的RC滤波方法,其特征在于,还包括:
步骤F:存储第N+1个初步处理数据,并删除所述前N个初步处理数据中的第一个;
步骤G:对第2个至第N+1个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
步骤H:将步骤G中的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
3.根据权利要求1所述的用于血压和心率测试的RC滤波方法,其特征在于,所述步骤B具体包括:
根据所述数字信号的原始数据的频率,确定一个基数,所述基数为所述频率的三分之二,并随机设置一个放大系数,所述放大系数为正数;
将所述数字信号的原始数据的频率减去所述基数再乘以所述放大系数,并进行滤除高频处理,根据数字信号的滤除结果加大或减小所述基数和放大系数,直至滤出的数字信号能够显示低频信号,将此时的数据作为所述初步处理数据。
4.根据权利要求1所述的用于血压和心率测试的RC滤波方法,其特征在于,所述步骤D具体包括:
计算所述N个初步处理数据分别对应的滤波系数;
将所述N个初步处理数据分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到N阶乘积求和运算的结果;
将所述N个初步处理数据中的前20个分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到20阶乘积求和运算的结果。
5.根据权利要求4所述的用于血压和心率测试的RC滤波方法,其特征在于,每个时间点对应一个初步处理数据,利用RC放电公式计算,RC常数对应的截止频率即为通过的最高频率,公式中每个时间点通过电容的电压值,即为每个时间点对应初步处理数据的滤波系数,RC放电公式如下:
Uc(t)=[Uc(0)]e^-t/RC,
其中,Uc(t)为每个时间点通过电容的电压值,也即滤波系数,Uc(0)为电容两端的初始电压,t为放电时间。
6.一种用于血压和心率测试的RC滤波器,其特征在于,包括:
采集转换模块,用于采集待处理的模拟信号,将所述模拟信号转换为数字信号,得到数字信号的原始数据;
初步处理模块,用于对所述数字信号的原始数据进行滤除高频处理,得到滤除高频后的初步处理数据;
第一存储模块,用于存储前N个所述初步处理数据,其中N为大于等于20的正整数;
第一运算模块,用于根据所述N个初步处理数据对应的滤波系数,对所述N个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
第一滤频模块,用于将第一运算模块计算的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
7.根据权利要求6所述的用于血压和心率测试的RC滤波器,其特征在于,还包括:
第二存储模块,用于存储第N+1个初步处理数据,并删除所述前N个初步处理数据中的第一个;
第二运算模块,用于对第2个至第N+1个初步处理数据分别进行N阶乘积求和运算和20阶乘积求和运算;
第二滤频模块,用于将第二运算模块计算的N阶乘积求和运算的结果与20阶乘积求和运算的结果相减,得到滤频后的数字信号。
8.根据权利要求6所述的用于血压和心率测试的RC滤波器,其特征在于,所述初步处理模块包括:
参数设定子模块,用于根据所述数字信号的原始数据的频率,确定一个基数,所述基数为所述频率的三分之二,并随机设置一个放大系数,所述放大系数为正数;
信号滤除子模块,用于将所述数字信号的原始数据的频率减去所述基数再乘以所述放大系数,并进行滤除高频处理,根据数字信号的滤除结果加大或减小所述基数和放大系数,直至滤出的数字信号能够显示低频信号,将此时的数据作为所述初步处理数据。
9.根据权利要求6所述的用于血压和心率测试的RC滤波器,其特征在于,所述第一运算模块包括:
系数计算子模块,用于计算所述N个初步处理数据分别对应的滤波系数;
N阶运算子模块,用于将所述N个初步处理数据分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到N阶乘积求和运算的结果;
20阶运算子模块,用于将所述N个初步处理数据中的前20个分别与对应的滤波系数进行乘积运算,并将各乘积运算结果相加得到20阶乘积求和运算的结果。
10.根据权利要求6所述的用于血压和心率测试的RC滤波器,其特征在于,每个时间点对应一个初步处理数据,利用RC放电公式计算,RC常数对应的截止频率即为通过的最高频率,公式中每个时间点通过电容的电压值,即为每个时间点对应初步处理数据的滤波系数,RC放电公式如下:
Uc(t)=[Uc(0)]e^-t/RC,
其中,Uc(t)为每个时间点通过电容的电压值,也即滤波系数,Uc(0)为电容两端的初始电压,t为放电时间。
CN201510179905.4A 2015-04-15 2015-04-15 一种用于血压和心率测试的 rc 滤波方法及rc 滤波器 Active CN104779934B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510179905.4A CN104779934B (zh) 2015-04-15 2015-04-15 一种用于血压和心率测试的 rc 滤波方法及rc 滤波器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510179905.4A CN104779934B (zh) 2015-04-15 2015-04-15 一种用于血压和心率测试的 rc 滤波方法及rc 滤波器

Publications (2)

Publication Number Publication Date
CN104779934A true CN104779934A (zh) 2015-07-15
CN104779934B CN104779934B (zh) 2018-04-24

Family

ID=53621216

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510179905.4A Active CN104779934B (zh) 2015-04-15 2015-04-15 一种用于血压和心率测试的 rc 滤波方法及rc 滤波器

Country Status (1)

Country Link
CN (1) CN104779934B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105832323A (zh) * 2016-06-12 2016-08-10 电子科技大学 一种基于rc电路的光电血管容积图信号处理方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005040223A (ja) * 2003-07-24 2005-02-17 Olympus Corp 医療用画像表示装置
CN101032396A (zh) * 2007-04-05 2007-09-12 上海交通大学 基于心率变异信号的脑缺血缺氧监护系统
CN103239222A (zh) * 2013-04-17 2013-08-14 深圳市科曼医疗设备有限公司 心电信号采集装置与采集方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005040223A (ja) * 2003-07-24 2005-02-17 Olympus Corp 医療用画像表示装置
CN101032396A (zh) * 2007-04-05 2007-09-12 上海交通大学 基于心率变异信号的脑缺血缺氧监护系统
CN103239222A (zh) * 2013-04-17 2013-08-14 深圳市科曼医疗设备有限公司 心电信号采集装置与采集方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
于金: "《机电一体化系统设计及实践》", 29 February 2008 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105832323A (zh) * 2016-06-12 2016-08-10 电子科技大学 一种基于rc电路的光电血管容积图信号处理方法

Also Published As

Publication number Publication date
CN104779934B (zh) 2018-04-24

Similar Documents

Publication Publication Date Title
Usui et al. Digital low-pass differentiation for biological signal processing
Chang et al. On-line empirical mode decomposition biomedical microprocessor for Hilbert Huang transform
Chand et al. FPGA implementation of fast FIR low pass filter for EMG removal from ECG signal
CN113568032B (zh) 一种基于z变换的负指数核脉冲信号处理方法及系统
CN104779934A (zh) 一种用于血压和心率测试的rc滤波方法及rc滤波器
CN104539264A (zh) 应用于eps电源系统的滤波方法及滤波电路
Vainio et al. A class of predictive analog filters for sensor signal processing and control instrumentation
Turulin et al. Analysis of controlled digital recursive high-pass filters structures with infinite non-negative impulse response
CN104579239B (zh) 一种滤波系统的过滤方法
Burleson Polynomial evaluation in VLSI using distributed arithmetic
Mishra et al. Optimized first-order s-to-z mapping function for IIR filter designing
CN107659290A (zh) 带宽扩展滤波器及其设计方法
Reddy et al. Comparison of FIR and IIR Filters using ECG Signal with Different Sampling Frequencies
CN102710237B (zh) 一种一阶数字低通滤波方法、滤波器及电子设备
CN112737546A (zh) 信号跟踪方法、滤波方法、装置及医疗设备
CN106896251B (zh) 一种示波器带内平坦度的自动修正方法
US10114509B2 (en) Multi-touch device having impedances positioned at electrode cross points and operation method therefor for detecting touch inputs
CN102868385B (zh) 一种可配置系数的数字滤波器和实现方法
KR101898534B1 (ko) 고주파 신호 필터링을 위한 디지털 체인 필터
Vityazeva et al. Optimal design of heart rate variability analyzing filter
CN103378821A (zh) 滤波器系统
Wang Design of digital signal filtering algorithm based on Bluetooth ECG monitoring system
Lei et al. Efficient 2D linear-phase IIR filter design and application in image filtering
CN116865720A (zh) 一种数字周期信号低通滤波方法
CN116388750A (zh) 高精度数字锁相放大器以及微弱信号检测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 214111 Jiangsu city of Wuxi province Wuxi hi tech District Branch IC design building block B Room 408

Applicant after: WUXI HEYOO TECHNOLOGY CO.,LTD.

Address before: 214111 Jiangsu city of Wuxi province Wuxi hi tech District Branch IC design building block B Room 408

Applicant before: WUXI HEROCHIP MICROELECTRONICS CO.,LTD.

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20180321

Address after: 214000 block No. 1, block A, block No. 34, block 34, block No. 400-2, Changjiang Road, Xin Wu District, Jiangsu Province, 400-2

Applicant after: Wuxi Si Rong Microelectronics Co.,Ltd.

Address before: 214111 Jiangsu city of Wuxi province Wuxi hi tech District Branch IC design building block B Room 408

Applicant before: WUXI HEYOO TECHNOLOGY CO.,LTD.

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200107

Address after: Room 516, building A3, innovation industrial park, No. 800, Wangjiang West Road, hi tech Zone, Hefei City, Anhui Province 230000

Patentee after: Hefei Sijing Electronic Co.,Ltd.

Address before: 214000 block No. 1, block A, block No. 34, block 34, block No. 400-2, Changjiang Road, Xin Wu District, Jiangsu Province, 400-2

Patentee before: Wuxi Si Rong Microelectronics Co.,Ltd.