CN109793519A - 呼吸速率检测方法、装置、存储介质及终端 - Google Patents
呼吸速率检测方法、装置、存储介质及终端 Download PDFInfo
- Publication number
- CN109793519A CN109793519A CN201910057439.0A CN201910057439A CN109793519A CN 109793519 A CN109793519 A CN 109793519A CN 201910057439 A CN201910057439 A CN 201910057439A CN 109793519 A CN109793519 A CN 109793519A
- Authority
- CN
- China
- Prior art keywords
- signal
- respiratory rate
- segment
- frequency
- rate detection
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 56
- 230000029058 respiratory gaseous exchange Effects 0.000 title claims abstract description 33
- 230000004044 response Effects 0.000 claims abstract description 51
- 230000005236 sound signal Effects 0.000 claims abstract description 47
- 230000036387 respiratory rate Effects 0.000 claims description 55
- 239000011159 matrix material Substances 0.000 claims description 42
- 238000005520 cutting process Methods 0.000 claims description 38
- 230000015654 memory Effects 0.000 claims description 23
- 238000012790 confirmation Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 230000011218 segmentation Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 238000000034 method Methods 0.000 abstract description 21
- 238000005516 engineering process Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明实施例公开了一种呼吸速率检测方法、装置、存储介质及终端,所述方法包括:分割原始音频信号为多个第一预设长度的第一信号片段;分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道;根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;根据多个信道冲激响应计算目标呼吸速率。本发明实施例的呼吸速率检测方案使用更加简便。
Description
技术领域
本发明实施例涉及通信技术领域,尤其涉及一种呼吸速率检测方法、装置、存储介质及终端。
背景技术
呼吸速率(respiratory rate)是在确定身体基本生命力时使用的最基本的生命体征。现有技术可以使用各种方法来测量呼吸速率,例如可以通过使用肺活量计测量进出肺的空气流量来计算呼吸速率,还可以通过测量每次呼吸气体中二氧化碳的浓度或分压来计算呼吸速率,还可以使用佩戴在用户胸部上的加速度传感器来估计呼吸速率。这些现有技术基本都需要专业的测量仪器对用户进行呼吸速率的测量,普通用户需要前往医院使用专门仪器进行测量,不方便在家或办公场所随时进行呼吸速率的测量。
因此,一种更加简便的呼吸速率检测方法成为了现有技术急需解决的技术问题。
发明内容
本发明实施例提供一种更加简便的呼吸速率检测方法、装置、存储介质及终端。
第一方面,本发明实施例提供了一种呼吸速率检测方法,包括:
分割原始音频信号为多个第一预设长度的第一信号片段;
分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道;
根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;
根据多个信道冲激响应计算目标呼吸速率。
第二方面,本发明实施例还提供了一种呼吸速率检测装置,所述呼吸速率检测装置包括:
第一分割模块,用于分割原始音频信号为多个第一预设长度的第一信号片段;
第二分割单元,用于分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道;
冲激响应模块,用于根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;
速率计算模块,用于根据多个信道冲激响应计算目标呼吸速率。
第三方面,本发明实施例提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述的呼吸速率检测方法。
第四方面,本发明实施例提供了一种终端,其特征在于,包括存储器,处理器及存储在存储器上并可在处理器运行的计算机程序,所述处理器执行所述计算机程序时实现上述的呼吸速率检测方法。
相对于现有技术,本发明实施例中提供的呼吸速率检测的方案根据预先定义的多个信道将原始音频信号和音频采集信号分割为多个信号片段,即第一信号片段和对应的第二信号片段;根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;根据多个信道冲激响应计算目标呼吸速率。本实施例的方案使用常见的终端例如手机进行呼吸速率检测,用户随时可以根据需求进行呼吸速率的检测,无需额外的专门检测设备,检测过程也非常方便。另外,在检测的同时用户会像平时听音乐一样感受终端扬声器发出的音乐,在用户倾听音乐的时候就同步实现了呼吸速率的检测,避免了现有技术进行呼吸速率检测需要使用扬声器发送超声波信号,用户无法正常使用扬声器播放音乐的缺陷,增强了用户的使用体验。
附图说明
图1为本发明实施例一提供的呼吸速率检测方法的流程图;
图2为本发明实施例二提供的呼吸速率检测方法的流程图;
图3为本发明实施例三提供的一种呼吸速率检测装置的结构示意图;
图4为本发明实施例四提供的一种终端的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
在更加详细地讨论示例性实施例之前应当提到的是,一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各步骤描述成顺序的处理,但是其中的许多步骤可以被并行地、并发地或者同时实施。此外,各步骤的顺序可以被重新安排。当其操作完成时处理可以被终止,但是还可以具有未包括在附图中的附加步骤。处理可以对应于方法、函数、规程、子例程、子程序等等。
此外,术语“第一”、“第二”等可在本文中用于描述各种方向、动作、步骤或元件等,但这些方向、动作、步骤或元件不受这些术语限制。这些术语仅用于将第一个方向、动作、步骤或元件与另一个方向、动作、步骤或元件区分。举例来说,在不脱离本申请的范围的情况下,可以将第一速度差值为第二速度差值,且类似地,可将第二速度差值称为第一速度差值。第一速度差值和第二速度差值两者都是速度差值,但其不是同一速度差值。术语“第一”、“第二”等而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
实施例一
图1为本发明实施例一提供的呼吸速率检测方法的流程图,用于在用户倾听音乐的同时实现呼吸速率检测,本实施例的方法可以适用于进行呼吸速率检测功能的终端设备,本实施例的方法可以由终端设备的处理器执行,具体包括:
S101、分割原始音频信号为多个第一预设长度的第一信号片段;
S102、分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道;
S103、根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;
S104、根据多个信道冲激响应计算目标呼吸速率。
本实施例终端设备可以是手机、平板、笔记本电脑、台式机等固定终端或移动终端。
针对步骤S101,本实施例中,原始音频信号为驱动终端扬声器的音频信号,本实施例中驱动终端扬声器的音频信号可以是终端针对用户选择的音乐文件进行解码播放过程中产生的音频信号。用户选择的音乐文件可以是存储在终端内的音乐文件也可以是访问网络服务器下载的音乐文件,当用户使用终端对选择的音乐文件进行播放时,可以通过音频解码芯片或处理器获取所述原始音频信号。本实施例中根据预先设置的切割数量、切割步长以及预设切割长度对原始音频信号进行切割,以生成第一预设长度(时间长度)的多个第一信号片段。本实施例的多个第一信号片段的切割长度为第一预设长度,多个第一信号片段之间的间隔为切割步长。本实施例的切割数量可以是100-1000之间的任意数值;切割步长Tstep可以为0.0.5-0.2秒之间的任一时间长度,例如0.1秒;第一预设长度可以是0.3-1.0之间的任意时间长度,例如0.4秒或0.6秒。
针对步骤S102,本实施例中,所述音频采集信号为终端麦克风采集的音频信号。本实施例中,当用户执行步骤S101使用终端播放音乐文件时,终端处理器同时控制麦克风采集扬声器播放的音乐声音以及其他环境噪音,以生成音频采集信号。本实施例中根据预先设置的切割数量、切割步长以及预设切割长度对音频采集信号进行切割,以生成第二预设长度(时间长度)的多个第二信号片段。本实施例的多个第二信号片段的切割长度为第二预设长度,多个第二信号片段之间的间隔为切割步长。本实施例的切割数量可以是100-1000之间的任意数值;切割步长Tstep可以为0.0.5-0.2秒之间的任一时间长度,例如0.1秒;第二预设长度可以是0.3-1.0之间的任意时间长度,例如0.4秒或0.6秒。本实施例中,切割过程中,每个信道的第一信号片段和对应的所述第二信号片段的切割的起点相同,且第一预设长度等于所述第二预设长度,本实施例根据预设的切割数量确定切割步长Tstep,即两次切割之间的时间间隔,切割步长Tstep和第一预设长度/第二预设长度之间的比例可以为1:4。
针对步骤S103,一实施例中,针对相同的信道,通过对第一信号片段进行快速傅里叶变换(Fast Fourier Transform,FFT),获得第一频域信号X=FFT(x),对第二信号片段进行快速傅里叶变换,获得第二频域信号Y=FFT(y),然后对每个信道对应的第一频域信号X和第二频域信号Y进行相除运算以获得多个信道的频率响应H=Y/X。随后,对每个信道的频率响应H进行快速傅里叶逆变换(Inverse Fast Fourier Transform,IFFT),以得到当前信道的冲激响应h=IFFT(H)。本实施例中对多对第一信号片段与第二信号片段计算得出多个信道的冲激响应,例如对n多个信道对应的n对第一信号片段与第二信号片段计算得出n帧信道冲激响应h1,h2,h3,…,hn。
针对步骤S104,本实施例中,将步骤S103获取的多个信道冲激响应处理为预设矩阵;对预设矩阵的每行数据进行快速傅里叶变换,以获得距离频率矩阵;根据距离频率矩阵内最大值元素确认所述目标呼吸速率。本实施例中,预设矩阵为CIR(Channel ImpulseResponse,CIR)矩阵,例如将步骤S103中获取的n帧信号h1,h2,h3,…,hn排列为矩阵MCIR后进行傅里叶变换处理为距离频率矩阵F,矩阵MCIR满足下面公式:
一实施例中,根据距离频率矩阵F内最大值元素的确认所述目标呼吸速率具体包括以下步骤:
S1041、确认距离频率矩阵内最大值元素的位置;
S1042、根据所述最大值元素的位置确认所述最大值元素的列数;
S1044、根据所述列数和预设的呼吸系数的乘积确认所述目标呼吸速率。
一实施例中,确认距离频率矩阵F中最大值元素的方式可以为:即通过两两比较遍历矩阵中所有元素后确认最大值元素。
设通过两两比较确认距离频率矩阵F最大值元素为F(f,d)后,可以根据F(f,d)所在位置确定最大值元素为hfk的列数为距离频率矩阵F的第f列,本实施例预设的呼吸系数的计算满足公式其中Tstep为第一预设长度,K为预设矩阵的总列数,则目标呼吸速率的计算方式为
相对于现有技术,本发明实施例中提供的呼吸速率检测的方法根据预先定义的步长、第一预设长度/第二预设长度以及切割数量将原始音频信号和音频采集信号分割为多个信号片段,即第一信号片段和对应的第二信号片段;根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;根据多个信道冲激响应计算目标呼吸速率。本实施例的方法使用常见的终端例如手机进行呼吸速率检测,用户随时可以根据需求进行呼吸速率的检测,无需额外的专门检测设备,检测过程也非常方便。另外,在检测的同时用户会像平时听音乐一样感受终端扬声器发出的音乐,在用户倾听音乐的时候就同步实现了呼吸速率的检测,避免了现有技术进行呼吸速率检测需要使用扬声器发送超声波信号,用户无法正常使用扬声器播放音乐的缺陷,增强了用户的使用体验。
实施例二
图2为本发明实施例二提供的呼吸速率检测方法的流程图,本实施例进一步提供了获取并同步原始音频信号和音频采集信号的方法,该方法包括:
S201、获取原始音频信号和音频采集信号;
S202、计算所述原始音频信号和所述音频采集信号的同步系数;
S203、根据所述同步系数将所述原始音频信号和音频采集信号同步对齐;
S204、分割原始音频信号为多个第一预设长度的第一信号片段;
S205、对第一信号片段进行快速傅里叶变换,获得第一频域信号X;
S206、分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道;
S207、对第二信号片段进行快速傅里叶变换,获得第二频域信号Y;
S208、对每个信道对应的第一频域信号X和第二频域信号Y进行相除运算以获得多个信道的频率响应;
S209、对所多个信道的频率响应进行快速傅里叶逆变换,以获得当前多个信道的冲激响应;
S210、将所述多个信道冲激响应处理为预设矩阵;
S211、对预设矩阵的每行数据进行快速傅里叶变换,以获得距离频率矩阵;
S212、确认距离频率矩阵内最大值元素的位置;
S213、根据所述最大值元素的位置确认所述最大值元素的列数;
S214、根据所述列数和预设的呼吸系数的乘积确认所述目标呼吸速率。
针对步骤S201-S203,本实施例中,原始音频信号为驱动终端扬声器的音频信号,本实施例中驱动终端扬声器的音频信号可以是终端针对用户选择的音乐文件进行解码播放过程中产生的音频信号。用户选择的音乐文件可以是存储在终端内的音乐文件也可以是访问网络服务器下载的音乐文件,当用户使用终端对选择的音乐文件进行播放时,可以通过音频解码芯片或处理器获取所述原始音频信号。所述音频采集信号为终端麦克风采集的音频信号。本实施例中,当用户使用终端播放音乐文件时,终端处理器同时控制麦克风采集扬声器播放的音乐声音以及其他环境噪音,以生成音频采集信号。
将原始音频信号和音频采集信号做互相关,以相关系数最高的点开始截取音频采集信号,使得接收到的音频采集信号与发出的原始音频信号对齐。本实施例对原始音频信号和音频采集信号做互相关计算满足以下公式:
其中,f为原始音频信号和音频采集信号之一,g为原始音频信号和音频采集信号之另一,本实施例原始音频信号和音频采集信号为离散的两个信号,星号表示f信号的复共轭。
针对步骤S204-S205,本实施例中根据预先设置的切割数量、切割步长以及预设切割长度对原始音频信号进行切割,以生成第一预设长度(时间长度)的多个第一信号片段。对第一信号片段进行快速傅里叶变换(Fast Fourier Transform,FFT),获得第一频域信号X=FFT(x)。本实施例的多个第一信号片段的切割长度为第一预设长度,多个第一信号片段之间的间隔为切割步长。本实施例的切割数量可以是100-1000之间的任意数值;切割步长Tstep可以为0.0.5-0.2秒之间的任一时间长度,例如0.1秒;第一预设长度可以是0.3-1.0之间的任意时间长度,例如0.4秒或0.6秒。
针对步骤S206-S207,本实施例中根据预先设置的切割数量、切割步长以及预设切割长度对音频采集信号进行切割,以生成第二预设长度(时间长度)的多个第二信号片段,针对相同信道对应的第一信号片段的第二信号片段,对第二信号片段进行快速傅里叶变换,获得第二频域信号Y=FFT(y)。本实施例的多个第二信号片段的切割长度为第二预设长度,多个第二信号片段之间的间隔为切割步长。本实施例的切割数量可以是100-1000之间的任意数值;切割步长Tstep可以为0.0.5-0.2秒之间的任一时间长度,例如0.1秒;第二预设长度可以是0.3-1.0之间的任意时间长度,例如0.4秒或0.6秒。本实施例中,切割过程中,每个信道对应的第一信号片段和所述第二信号片段的切割的起点相同,且第一预设长度等于所述第二预设长度,本实施例根据预设的切割数量确定切割步长Tstep,即两次切割之间的时间间隔,切割步长Tstep和第一预设长度/第二预设长度之间的比例可以为1:4。
步骤S208中,对每个信道对应的第一频域信号X和第二频域信号Y进行相除运算以获得多个信道的频率响应H=Y/X。本实施例中同时对多个信道中的每个信道进行频率响应计算,通过分别获得每个信道的频率响应H来确定多个信道的频率响应。
步骤S209中,对所有信道的频率响应H进行快速傅里叶逆变换(Inverse FastFourier Transform,IFFT),以得到多个当前信道的冲激响应h=IFFT(H),例如n帧信号h1,h2,h3,…,hn。
步骤S210中,将所述多个信道冲激响应处理为预设矩阵,例如CIR(ChannelImpulse Response,CIR)矩阵:
本实施例的矩阵MCIR的列数等于分割数量或信道数量,一个信道的采集数据对应为一帧数据,矩阵MCIR的行数为每帧数据中对应的采集频点的数量。
步骤S211中,对预设矩阵的每行数据进行快速傅里叶变换,以获得距离频率矩阵可以为:
步骤S212中,确认距离频率矩阵F中最大值元素的方式可以为:即通过两两比较遍历距离频率矩阵F中所有元素后确认最大值元素。
步骤S213中,设通过两两比较确认距离频率矩阵F的最大值元素为F(f,d)后,可以根据F(f,d)所在位置确定最大值元素为hfk的列数为距离频率矩阵F的第f列。
步骤S214中,本实施例预设的呼吸系数的计算满足公式其中Tstep为切割步长,K为预设矩阵的总列数,则目标呼吸速率的计算方式为
相对于现有技术,本实施例的呼吸速率检测方法进一步提供了获取并同步原始音频信号和音频采集信号的方法,用于在原始音频信号和音频采集信号同步不理想的时候进行信号处理,增加了呼吸速率检测的准确率。
实施例三
图3为本发明实施例三提供的一种呼吸速率检测装置300的结构示意图,该装置可由软件和/或硬件实现,一般集成在终端中,可通过执行呼吸速率检测方法实现对属性参数的自动存取。
本实施例的呼吸速率检测装置300包括:
第一分割模块301,用于分割原始音频信号为多个第一预设长度的第一信号片段;所述原始音频信号为驱动终端扬声器的音频信号,所述音频采集信号为终端麦克风采集的音频信号。
第二分割单元302,用于分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道,所述第一预设长度等于所述第二预设长度;
冲激响应模块303,用于根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;
速率计算模块304,用于根据多个信道冲激响应计算目标呼吸速率。
进一步地,呼吸速率检测装置300还包括信号同步模块305,用于获取原始音频信号和音频采集信号;计算所述原始音频信号和所述音频采集信号的同步系数;根据所述同步系数将所述原始音频信号和音频采集信号同步对齐。
进一步地,冲激响应模块303还包括:
第一变换模块,用于对第一信号片段进行快速傅里叶变换,获得第一频域信号X;
第二变换模块,用于对第二信号片段进行快速傅里叶变换,获得第二频域信号Y;
频响计算模块,用于对每个信道对应的第一频域信号X和第二频域信号Y进行相除运算以获得多个信道的频率响应。
进一步地,速率计算模块304还包括:
矩阵转换模块,用于将所述多个信道冲激响应处理为预设矩阵;
距离矩阵模块,用于对预设矩阵的每行数据进行快速傅里叶变换,以获得距离频率矩阵;
速率计算模块,用于根据距离频率矩阵内最大值元素的确认所述目标呼吸速率。
进一步地,速率计算模块还包括:
位置确认模块,用于确认距离频率矩阵内最大值元素的位置;
列数确认模块,用于根据所述最大值元素的位置确认所述最大值元素的列数;
速率计算模块,用于根据所述列数和预设的呼吸系数的乘积确认所述目标呼吸速率。
相对于现有技术,本实施例呼吸速率检测装置用于执行前述本实施例呼吸速率检测方法,因此还可以实现前述实施例的所有优点,在此不再赘述。
实施例四
图4为本发明实施例四提供的一种终端的结构示意图,如图4所示,该种终端包括存储器401、处理器402,终端中处理器402的数量可以是一个或多个,图4中以一个处理器402为例;终端中的存储器401、处理器402可以通过总线或其他方式连接,图4中以通过总线连接为例。
存储器401作为一种计算机可读存储介质,可用于存储软件程序、计算机可执行程序以及模块,如本发明实施例中的呼吸速率检测方法对应的程序指令/模块,处理器402通过运行存储在存储器401中的软件程序、指令以及模块,从而执行设备/终端/终端的各种功能应用以及数据处理,即实现上述的呼吸速率检测方法。
存储器401可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据终端的使用所创建的数据等。此外,存储器401可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中,存储器401可进一步包括相对于处理器402远程设置的存储器,这些远程存储器可以通过网络连接至设备/终端/终端。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
存储介质——任何的各种类型的存储器设备或存储设备。术语“存储介质”旨在包括:安装介质,例如CD-ROM、软盘或磁带装置;计算机系统存储器或随机存取存储器,诸如DRAM、DDRRAM、SRAM、EDORAM,兰巴斯(Rambus)RAM等;非易失性存储器,诸如闪存、磁介质(例如硬盘或光存储);寄存器或其它相似类型的存储器元件等。存储介质可以进一步包括其它类型的存储器或其组合。另外,存储介质可以位于程序在其中被执行的第一计算机系统中,或者可以位于不同的第二计算机系统中,第二计算机系统通过网络(诸如因特网)连接到第一计算机系统。第二计算机系统可以提供程序指令给第一计算机用于执行。术语“存储介质”可以包括可以驻留在不同位置中(例如在通过网络连接的不同计算机系统中)的两个或更多存储介质。存储介质可以存储可由一个或多个处理器执行的程序指令(例如具体实现为计算机程序)。
一实施例中,本发明实施例提供的一种包含计算机可执行指令的存储介质计算机可执行指令在由计算机处理器执行时用于执行一种呼吸速率检测方法,该方法包括:
分割原始音频信号为多个第一预设长度的第一信号片段;
分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道;
根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;
根据多个信道冲激响应计算目标呼吸速率。
当然,本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上的方法操作,还可以执行本发明任意实施例所提供的呼吸速率检测方法中的相关操作。
通过以上关于实施例的描述,所属领域的技术人员可以清楚地了解到,本发明可借助软件及必需的通用硬件来实现,当然也可以通过硬件实现,但很多情况下前者是更佳的实施例。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如计算机的软盘、只读存储器(Read-OnlyMemory,ROM)、随机存取存储器(RandomAccessMemory,RAM)、闪存(FLASH)、硬盘或光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,终端,或者网络设备等)执行本发明各个实施例的方法。
值得注意的是,上述呼吸速率检测装置的实施例中,所包括的各个单元和模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。
Claims (10)
1.一种呼吸速率检测方法,其特征在于,包括:
分割原始音频信号为多个第一预设长度的第一信号片段;
分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道;
根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;
根据多个信道冲激响应计算目标呼吸速率。
2.根据权利要求1所述的呼吸速率检测系统,其特征在于:所述分割原始音频信号为多个第一预设长度的第一信号片段,之前包括:
获取原始音频信号和音频采集信号;
计算所述原始音频信号和所述音频采集信号的同步系数;
根据所述同步系数将所述原始音频信号和音频采集信号同步对齐。
3.根据权利要求1所述的呼吸速率检测方法,其特征在于:所述原始音频信号为驱动终端扬声器的音频信号,所述音频采集信号为终端麦克风采集的音频信号。
4.根据权利要求1所述的呼吸速率检测方法统,其特征在于:所述第一预设长度等于所述第二预设长度。
5.根据权利要求1所述的呼吸速率检测方法,其特征在于:所述根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应,包括:
对第一信号片段进行快速傅里叶变换,获得第一频域信号X;
对第二信号片段进行快速傅里叶变换,获得第二频域信号Y;
对每个信道对应的第一频域信号X和第二频域信号Y进行相除运算以获得多个信道的频率响应;
对多个信道的频率响应进行快速傅里叶逆变换,以获得当前多个信道的冲激响应。
6.根据权利要求1所述的呼吸速率检测系统,其特征在于:所述根据多个信道冲激响应计算目标呼吸速率,包括:
将所述多个信道冲激响应处理为预设矩阵;
对预设矩阵的每行数据进行快速傅里叶变换,以获得距离频率矩阵;
根据距离频率矩阵内最大值元素的确认所述目标呼吸速率。
7.根据权利要求6所述的呼吸速率检测系统,其特征在于:所述根据距离频率矩阵内最大值元素的确认所述目标呼吸速率,包括:
确认距离频率矩阵内最大值元素的位置;
根据所述最大值元素的位置确认所述最大值元素的列数;
根据所述列数和预设的呼吸系数的乘积确认所述目标呼吸速率。
8.一种呼吸速率检测装置,其特征在于,包括:
第一分割模块,用于分割原始音频信号为多个第一预设长度的第一信号片段;
第二分割单元,用于分割音频采集信号为多个第二预设长度的第二信号片段,所述多个第一信号片段和所述多个第二信号片段一一对应定义多个信道;
冲激响应模块,用于根据多个第一信号片段和对应的第二信号片段确认多个信道的冲激响应;
速率计算模块,用于根据多个信道冲激响应计算目标呼吸速率。
9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-7中任一所述的呼吸速率检测方法。
10.一种呼吸速率检测终端,其特征在于,包括存储器,处理器及存储在存储器上并可在处理器运行的计算机程序,所述处理器执行所述计算机程序时实现如权利要求1-7任一所述的呼吸速率检测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910057439.0A CN109793519B (zh) | 2019-01-22 | 2019-01-22 | 呼吸速率检测方法、装置、存储介质及终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910057439.0A CN109793519B (zh) | 2019-01-22 | 2019-01-22 | 呼吸速率检测方法、装置、存储介质及终端 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109793519A true CN109793519A (zh) | 2019-05-24 |
CN109793519B CN109793519B (zh) | 2021-12-24 |
Family
ID=66559943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910057439.0A Active CN109793519B (zh) | 2019-01-22 | 2019-01-22 | 呼吸速率检测方法、装置、存储介质及终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109793519B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007264567A (ja) * | 2006-03-30 | 2007-10-11 | Railway Technical Res Inst | 音声中の非発話音声の判別処理方法 |
CN102215740A (zh) * | 2008-11-17 | 2011-10-12 | 索尼爱立信移动通讯有限公司 | 用于从生理声信号检测生理测量结果的设备、方法和计算机程序 |
US8652040B2 (en) * | 2006-12-19 | 2014-02-18 | Valencell, Inc. | Telemetric apparatus for health and environmental monitoring |
CN106308801A (zh) * | 2015-06-29 | 2017-01-11 | 无锡智感星际科技有限公司 | 一种利用智能手机检测人体呼吸频率的方法 |
CN107928673A (zh) * | 2017-11-06 | 2018-04-20 | 腾讯科技(深圳)有限公司 | 音频信号处理方法、装置、存储介质和计算机设备 |
CN108417225A (zh) * | 2018-02-11 | 2018-08-17 | 广州市碳码科技有限责任公司 | 一种呼吸音监测方法、装置、终端及计算机可读存储介质 |
WO2018205013A1 (en) * | 2017-05-10 | 2018-11-15 | Ecole De Technologie Superieure | System and method for determining cardiac rhythm and/or respiratory rate |
-
2019
- 2019-01-22 CN CN201910057439.0A patent/CN109793519B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007264567A (ja) * | 2006-03-30 | 2007-10-11 | Railway Technical Res Inst | 音声中の非発話音声の判別処理方法 |
US8652040B2 (en) * | 2006-12-19 | 2014-02-18 | Valencell, Inc. | Telemetric apparatus for health and environmental monitoring |
CN102215740A (zh) * | 2008-11-17 | 2011-10-12 | 索尼爱立信移动通讯有限公司 | 用于从生理声信号检测生理测量结果的设备、方法和计算机程序 |
CN106308801A (zh) * | 2015-06-29 | 2017-01-11 | 无锡智感星际科技有限公司 | 一种利用智能手机检测人体呼吸频率的方法 |
WO2018205013A1 (en) * | 2017-05-10 | 2018-11-15 | Ecole De Technologie Superieure | System and method for determining cardiac rhythm and/or respiratory rate |
CN107928673A (zh) * | 2017-11-06 | 2018-04-20 | 腾讯科技(深圳)有限公司 | 音频信号处理方法、装置、存储介质和计算机设备 |
CN108417225A (zh) * | 2018-02-11 | 2018-08-17 | 广州市碳码科技有限责任公司 | 一种呼吸音监测方法、装置、终端及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109793519B (zh) | 2021-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9357306B2 (en) | Multichannel audio calibration method and apparatus | |
Park et al. | Sound-field analysis by plane-wave decomposition using spherical microphone array | |
EP3320692B1 (en) | Spatial audio processing apparatus | |
KR101591220B1 (ko) | 공간적 전력 밀도에 기초하여 마이크 위치 결정을 위한 장치 및 방법 | |
US7590248B1 (en) | Head related transfer function filter generation | |
US10334357B2 (en) | Machine learning based sound field analysis | |
TWI530201B (zh) | 經由自抵達方向估值提取幾何資訊之聲音擷取技術 | |
CN109313909B (zh) | 评估麦克风阵列一致性的方法、设备、装置和系统 | |
JP5957446B2 (ja) | 音響処理システム及び方法 | |
US9351070B2 (en) | Positional disambiguation in spatial audio | |
US8144896B2 (en) | Speech separation with microphone arrays | |
JP2006519406A5 (zh) | ||
KR101882423B1 (ko) | 적어도 제1 쌍의 룸 임펄스 응답에 기초하여, 믹싱 시간 전체를 추정하는 장치 및 방법, 대응하는 컴퓨터 프로그램 | |
CN104937955A (zh) | 自动的扬声器极性检测 | |
JP5000297B2 (ja) | 音場の表現を求めるためのシステムおよび方法 | |
Traa et al. | Blind multi-channel source separation by circular-linear statistical modeling of phase differences | |
CN113687305A (zh) | 声源方位的定位方法、装置、设备及计算机可读存储介质 | |
CN109793519A (zh) | 呼吸速率检测方法、装置、存储介质及终端 | |
WO2014083380A1 (en) | A shared audio scene apparatus | |
EP3203760A1 (en) | Method and apparatus for determining the position of a number of loudspeakers in a setup of a surround sound system | |
Di Carlo et al. | dEchorate: a calibrated room impulse response database for echo-aware signal processing | |
Defrance et al. | Finding the onset of a room impulse response: Straightforward? | |
Qiao et al. | Isolation performance metrics for personal sound zone reproduction systems | |
CN111899749A (zh) | 一种监听混凝土泵送管道作业声音的降噪方法 | |
CN114303392A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |