CN110680349B - 一种基于线性调频的脉搏测谎方法与装置 - Google Patents
一种基于线性调频的脉搏测谎方法与装置 Download PDFInfo
- Publication number
- CN110680349B CN110680349B CN201911036322.0A CN201911036322A CN110680349B CN 110680349 B CN110680349 B CN 110680349B CN 201911036322 A CN201911036322 A CN 201911036322A CN 110680349 B CN110680349 B CN 110680349B
- Authority
- CN
- China
- Prior art keywords
- signal
- pulse
- unit
- amplitude
- layer
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 56
- 238000012545 processing Methods 0.000 claims abstract description 42
- 230000006835 compression Effects 0.000 claims abstract description 24
- 238000007906 compression Methods 0.000 claims abstract description 24
- 239000000203 mixture Substances 0.000 claims abstract description 7
- 239000000523 sample Substances 0.000 claims description 60
- 230000006870 function Effects 0.000 claims description 36
- 238000013500 data storage Methods 0.000 claims description 19
- 230000003321 amplification Effects 0.000 claims description 18
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 18
- 238000007781 pre-processing Methods 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 238000012805 post-processing Methods 0.000 claims description 14
- 230000005540 biological transmission Effects 0.000 claims description 11
- 238000002604 ultrasonography Methods 0.000 claims description 11
- 230000004044 response Effects 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 7
- 230000035485 pulse pressure Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000008520 organization Effects 0.000 claims description 3
- 230000026676 system process Effects 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims 1
- 238000013139 quantization Methods 0.000 claims 1
- 210000001519 tissue Anatomy 0.000 description 14
- 238000000034 method Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 238000011835 investigation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 208000004756 Respiratory Insufficiency Diseases 0.000 description 2
- 206010038678 Respiratory depression Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000036772 blood pressure Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035479 physiological effects, processes and functions Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 230000028327 secretion Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000035900 sweating Effects 0.000 description 2
- 206010006322 Breath holding Diseases 0.000 description 1
- 206010049816 Muscle tightness Diseases 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 201000007201 aphasia Diseases 0.000 description 1
- 238000013475 authorization Methods 0.000 description 1
- 210000003403 autonomic nervous system Anatomy 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 230000001079 digestive effect Effects 0.000 description 1
- 206010013781 dry mouth Diseases 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 235000011389 fruit/vegetable juice Nutrition 0.000 description 1
- 230000002496 gastric effect Effects 0.000 description 1
- 210000004247 hand Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
- 238000010223 real-time analysis Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 238000006748 scratching Methods 0.000 description 1
- 230000002393 scratching effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000007920 subcutaneous administration Methods 0.000 description 1
- 210000000106 sweat gland Anatomy 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000003857 wrist joint Anatomy 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
- A61B5/164—Lie detection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/02—Measuring pulse or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Cardiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Child & Adolescent Psychology (AREA)
- Developmental Disabilities (AREA)
- Educational Technology (AREA)
- Hospice & Palliative Care (AREA)
- Psychiatry (AREA)
- Psychology (AREA)
- Social Psychology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于线性调频的脉搏测谎方法与装置，传感器采用非阵列超声探头发射线性调频的连续信号，线性调频信号具有精确、抗多径、去伪影的优点。结合预失真处理、脉冲压缩、分层反演回波算法，通过各反射回波之间的时延，得到各反射面的位置(深度)；利用超声波在体内的衰减函数，通过迭代的方式反演出发射信号与回波信号在每个反射层的幅值，并计算出每一层反射面的反射系数；反射系数可确定各反射层的组织构成，得到脉搏的深度、频率、幅度等多维度的脉搏信息，从而可为测谎结果的判断提供更多维度的参考。
Description
技术领域
本发明涉及声学探测、信号处理等技术领域，具体涉及一种基于线性调频的脉搏测谎方法与装置。
背景技术
测谎技术是一项犯罪心理测试技术，在侦查工作前期，侦查人员根据初步侦查和调查掌握的情况，发现了嫌疑对象，但尚无确凿证据认定该对象是否犯罪分子，通常采用传唤或拘传的方法，进行正面审查，以排除或肯定犯罪嫌疑。但作案者和无辜者都会极力辩解，否认犯罪，其供述有的一时难以查证或虽能查证但需花费大量时间和经费。这时，若使用测谎技术，就能迅速地排除大量的无辜嫌疑人，筛选出重点嫌疑对象，在条件比较好的情况下还可以直接认定犯罪分子，然后围绕重点嫌疑对象开展审讯和调查，可以事半功倍，大大提高破案效率。
现代科学证实，人在说谎时生理上的确发生着一些变化，有一些肉眼可以观察到，如出现抓耳挠腮、腿脚抖动等一系列不自然的人体动作。还有一些生理变化是不易察觉的，如：呼吸速率和血容量异常，出现呼吸抑制和屏息；脉搏加快，血压升高，血输出量增加及成分变化，导致面部、颈部皮肤明显苍白或发红；皮下汗腺分泌增加，导致皮肤出汗，双眼之间或上嘴唇首先出汗，手指和手掌出汗尤其明显；眼睛瞳孔放大；胃收缩，消化液分泌异常，导致嘴、舌、唇干燥；肌肉紧张、颤抖，导致说话结巴。这些生理参量由于受植物神经***支配，所以一般不受人的意识控制，而是自主的运动，在外界刺激下会出现一系列条件反射现象。这一切都逃不过测谎仪的“眼睛”。而测谎一般是从三个方面测定一个人的生理变化，即脉搏、呼吸和皮肤电阻，其中脉搏为首位；因此若能得到多维度的脉搏信息，将会提高对测谎结果判断的准确度。
目前的脉搏检测装置主要采用压力、光学传感器感知脉搏，并配合不同类型的机械装置，例如：微型马达、电磁阀、气泵、液压泵装置等，模拟人手获取脉搏的信息。多个专利公开了基于压力、温度传感装置的脉搏检测装置，例如：中国专利，授权号：ZL86107657.3、授权号：CN102411615、申请号：301310629076.6等。其核心均是通过固定在手上的压力传感器感知患者的脉搏信息，采用温度传感器校准压力传感器的误差，然而其仍具有触感类型传感器的一些不足，例如：精度不高；传感器与手腕贴合层度对结果影响较大，需要用高精密度、高成本的气泵和液压装置减小误差。
而超声因其安全、无损等特点被长期用来采集脉搏数据，例如：中国专利，申请号：200610083339.8，借助医用阵列超声三维成像设备，通过获取脉搏的三维图形得到脉搏各种参数；但该方法采用脉冲信号作为探测信号，脉冲信号一般用于成像，并需结合阵列探头使用、且运算过程复杂。
发明内容
本发明的目的是为了实现脉搏信息的准确获取，提供一种基于线性调频的脉搏测谎方法与装置。
本发明的第一个目的可以如下技术方案达到：
一种基于线性调频的脉搏测谎方法，包括如下步骤：
S1、信号预处理：采用非阵列探头发射线性调频的连续信号，线性调频信号具有精确、抗多径、去伪影等特点；并通过幅度加权探头的脉冲响应，使得回波信号的带宽不局限于探头，从而提高轴向分辨力；最后使用Lanczos窗函数作为幅度削减函数降低脉压后的旁瓣；
S2、脉冲压缩：对接收到的反射回波做脉冲压缩处理，以获得峰值功率较高的窄脉冲；并通过加Hanmming窗，降低脉冲压缩后的旁瓣；
S3、分层反演：通过分析脉冲压缩后各反射回波的时延信息，反演出各反射层的反射系数，并对比已有组织的反射系数，定位到脉搏(血管)的位置信息，从而获得脉搏的深度、频率、幅度和组织状态，为测谎提供更多维度的参考。
进一步地，所述的步骤S1中信号预处理具体包括以下步骤：
S1.1、确定信号形式：由于线性信号具有很好的自相关特性，原始信号为线性调频信号,由用户的输入得
其中f0是信号的中心频率，K是信号带宽，T是信号的时间宽度，A是信号的幅度；
S1.2、幅度加权，对线性调频作幅度加权的预处理：
xpre(t)＝Mprob(t)Mdecr(t)x(t),0≤t≤T
其中，Mprob(t)是基于探头的幅度加权函数；
由于超声探头的脉冲响应可近似为高斯包络的余弦波：
其中，Δfu是探头的带宽，f0是探头的中心频率，μ是一个常量；因此幅度加权函数可表示为高斯函数的倒数：
其中，D是加权系数，表示幅度加权的程度。D值越小，对信号幅度加权的程度越大；D的取值与探头的脉冲响应、信号的带宽和时宽有关系；
而Mdecr(t)是Lanczos窗函数，目的是减小脉压旁瓣幅度：
其中，
i＝1,…,m+1,m为幅度削减长度。
进一步地，步骤S2中脉冲压缩具体包括以下步骤：
S2.1、得匹配滤波函数：使得输出信噪比最大化的准则，可得到匹配滤波器的表达式：
hmatch(t)＝kx*(t0-t)
其中，k为任意常数，x*为线性调频信号的共轭，t0为预期偏移量；
S2.2、在匹配滤波器的基础上加Hanmming窗函数，得失配滤波器的表达式如下：
hmismatch(t)＝w(t)hmatch(t)
其中，Hanmming窗w(n)的定义为：
其中β＝0.46，N为Hanmming窗的总长度；
S2.3、脉冲压缩：接收探头收到的反射回波经过失配滤波器后，可以获得峰值功率较高的窄脉冲：
y(t)＝conv(r(t),F-1(X*(f)Hmismatch(f)))
其中，r(t)为反射回波，F-1为逆傅里叶变换，X*(f)为线性调频信号经过傅里叶变换后的共轭，Hmismatch(f)为失配滤波器的***函数，conv(*)为卷积运算。
进一步地，步骤S3中分层反演具体包括以下步骤：
S3.1、计算各反射层的相对距离：将各窄脉冲峰值对应的时延结果ni根据采样频率与声速转变为距离di：
其中，fs为***采样频率，C为超声波在人体组织内的平均声速；di即表示各反射面相对于皮肤的距离；
S3.2、衰减函数：超声波信号在人体内衰减函数为：
其中，U0为信号的初始幅值，f为信号的频率，αN为体内第N层组织的衰减系数，d为信号在体内移动的距离；
S3.3、原始发射信号的幅值X0经过了已知初始反射层的透射，再经过人体组织衰减函数U，得第一层反射层处的幅值X1：
其中，R0为已知初始反射层的反射系数，α1为初始反射层到第一层反射层之间的衰减系数，m1＝d1-d0为第一层反射层与初始反射层之间的距离；
S3.4、而反射回波Y0则经过人体组织衰减，再经过已知初始反射层的透射，得第一层反射层处的幅值Y1：
S3.5、因此第一层反射层的反射系数即为：
其中，X1原始发射信号X0在第一层反射层处的幅值，Y1为反射回波Y0在第一层反射层处的幅值；
S3.6、以此类推，通过迭代得到发射信号在第i层反射面处的幅值Xi与反射回波信号在第i层反射面处的幅值Yi：
其中，n为反射面的总数，mi＝di-di-1为第i层与第i-1层之间的距离，Ri-1为第i-1层的反射系数，αi为第i-1层反射层到第i层反射层之间的衰减系数；
S3.7、因此第i层反射面的反射系数为：
S3.8、得到了各层的反射系数以后，通过对比已有组织的反射系数，可以确定各反射层的组织构成，从而得到脉搏深度、频率、幅度信息等多维度信息；从而为测谎结果的判断，提供多维度的脉搏参考。
本发明的另一个目的可以如下技术方案达到：
一种基于线性调频的脉搏测谎装置，采用数据处理端与数据采集端分离的方式，具体原理如下：
该基于线性调频的脉搏测谎装置由数据采集端101、数据处理端102两部分组成，两个部分之间是通过无线电磁波相关联。其中数据采集端101负责超声波的收发，模拟信号与数字信号之间的转换，无线发射回波信号；数据处理端102负责通过对发射信号的预处理，提高轴向分辨力和降低旁瓣；同时对反射回波信号进行后处理，得到脉搏深度、频率、幅度等多维度的信息，以为测谎结果的判断提供更全面的参考。
其中，所述的基于线性调频的脉搏测谎装置的数据采集端101负责对数据的采集、数据的转发、***流程的时序控制。其中包括收发逻辑单元201、无线发射单元202、数据存储单元203、AD/DA转换单元204、功率放大单元205、非阵列超声探头单元206、低噪声放大单元207。其组成连接关系如下：
收发逻辑单元201与无线发射单元202、数据存储单元203、AD/DA转换单元204相连接，负责收发流程的时序控制，如发送信号的间隔控制、外部设备(存储器、AD/DA芯片、无线模块)的时序控制等；无线发射单元202与AD/DA转换单元204相连接，负责将数字量化后的超声波反射回波以无线电磁波的形式发送至数据处理端102用于分析处理；数据存储单元203与AD/DA转换单元204、收发逻辑单元201相连接，负责存储数据处理端102中预处理后的信号数据，收发逻辑单元201可通过一定的流程逻辑去读取使用其中的数据；AD/DA转换单元204与功率放大单元205、低噪声放大单元207、无线发射单元202、数据存储单元203相连接，负责将超声探头转换而来的模拟信号与处理芯片能识别的数字信号进行相互转换；功率放大单元205与非阵列超声探头单元206相连接，负责将低电压信号转换为可以驱动非阵列超声探头单元206的驱动信号；低噪声放大单元207与非阵列超声探头单元206相连接，负责通过低噪声放大器将由非阵列超声探头单元206转换过来的电信号进行放大；非阵列超声探头单元206负责实现电信号与非阵列超声信号之间的转换，即实现了电能与声能之间的转换。
其中，所述的基于线性调频的脉搏测谎装置的数据处理端102负责数据的分析处理处理，其包括信号预处理单元301、无线接收单元302、信号后处理单元303。其组成连接关系如下：
若数据存储单元203中无数据存储，或者需要更新发射信号时，信号预处理单元301与数据存储单元203相连接，负责对信号做预失真处理，从而提高轴向分辨力并减低旁瓣；无线接收单元302与信号后处理单元303相连接，负责接收来自数据采集端的超声回波数据；信号后处理单元303负责对无线接收单元302的数据进行后处理，得到脉搏深度、频率、幅度等多维度的信息，以为测谎结果的判断提供更全面的参考。
本发明相对于现有技术具有如下的优点及效果：
(1)本发明使用线性调频的连续信号作为发射信号，通过处理反射回波得到脉搏的多种信息，相比用于成像的脉冲信号，具有所需探头数量少(非阵列探头)、精准度高、抗多径的优点；而相比于压力传感器，则无需考虑装置对探测点的接触力度；
(2)本发明采用预失真处理、脉冲压缩与分层反演回波算法相结合，可在提高轴向分辨力与降低旁瓣的矛盾中取得一个较好的平衡，获得各反射回波的时延信息，从而推得各反射层的反射系数，通过对比已有组织的反射系数，可确定各反射层的组织构成，从而得到脉搏深度、频率、幅度等多维度的信息，为测谎结果的判断提供更多的参考；
(3)本发明采用无线数据回传的方式，因此所佩戴的数据采集端仅进行简单的数据采集，复杂的算法实现放在数据处理端，从而减小了数据采集端的体积、降低了功耗；
综上，本发明可测得脉搏深度、频率、幅度等多维度信息，测谎结果的判断提供更多的参考，同时具有精准度高、装置简单、使用便捷等优点；本发明对脉搏信息进行量化，并且经过实时分析和存储后，可作为大数据的样本用于聚类分析，为未来测谎技术的大数据分析奠定了基础。
附图说明
图1是本发明公开的基于线性调频的脉搏测谎装置的结构图；
图2是本发明公开的数据处理端的结构图；
图3是本发明公开的数据采集端的结构图；
图4是本发明公开的分层反演回波的示例图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例
测谎技术可以通过人在说谎时生理上产生的一系列变化，迅速地排除大量的无辜嫌疑人，筛选出重点嫌疑对象，在条件比较好的情况下还可以直接认定犯罪分子，然后围绕重点嫌疑对象开展审讯和调查，可以事半功倍，大大提高破案效率。
而说谎时的生理变化体现在，呼吸速率和血容量异常，出现呼吸抑制和屏息；脉搏加快，血压升高，血输出量增加及成分变化，导致面部、颈部皮肤明显苍白或发红等。因此对脉搏的多维度检测可以为测谎结果的判断提供更全面的参考。
有鉴于此，本实施例提出了一种基于线性调频的脉搏测谎方法。基于超声波的物理特性，采用非阵列超声探头发射线性调频的连续信号，并在信号处理方法上使用预失真处理、脉冲压缩、分层反演回波算法，其中包括：1)信号发射通过预失真技术抵消发射探头的脉冲响应对发射信号的影响，使得回波信号的带宽不局限于探头，从而提高轴向分辨力；增加幅度削减函数降低旁瓣；2)信号接收采用脉冲压缩，恢复轴向分辨力，从而得到各反射回波之间的时延；3)由时延可得到各反射面的位置(深度)；结合超声波在体内的衰减函数，可迭代地计算出每一层反射面的反射系数；通过反射系数，可确定各反射层的组织构成，从而得到脉搏深度、频率、幅度等多维度的信息，从而可为测谎结果的判断提供更多维度的参考。
本实施例公开的一种基于线性调频的脉搏测谎方法，，包括以下步骤：
S1、信号预处理：采用非阵列探头发射线性调频的连续信号，线性调频信号具有精确、抗多径、去伪影等特点；并通过幅度加权探头的脉冲响应，使得回波信号的带宽不局限于探头，从而提高轴向分辨力；最后使用Lanczos窗函数作为幅度削减函数降低脉压后的旁瓣；
该步骤S1中信号预处理具体包括以下步骤：
S1.1、确定信号形式：由于线性信号具有很好的自相关特性，原始信号为线性调频信号,由用户的输入得：
其中f0是信号的中心频率，K是信号带宽，T是信号的时间宽度；
S1.2、幅度加权：为了使得回波信号的带宽不局限于探头，并减少脉压后的旁瓣，需对线性调频作幅度加权的预处理：
xpre(t)＝Mprob(t)Mdecr(t)x(t)，tC[0，T]
其中，Mprob(t)是基于探头的幅度加权函数，为抵消探头对线性调频信号的影响；由于超声探头的脉冲响应可近似为高斯包络的余弦波：
其中，Δfu是探头的-6dB带宽，f0是探头的中心频率，α是一个常量1.2ln10；因此幅度加权函数可表示为高斯函数的倒数：
其中，D是加权系数，表示幅度加权的程度。D值越小，对信号幅度加权的程度越大；D的取值与探头的脉冲响应、信号的带宽和时宽有关系；
而Mdecr(t)是Lanczos窗函数，目的是减小脉压旁瓣幅度：
其中，
i＝1，...，m+1，m为幅度削减长度。
S2、脉冲压缩：对接收到的反射回波做脉冲压缩处理，以获得峰值功率较高的窄脉冲；并通过加Hanmming窗，降低脉冲压缩后的旁瓣；
该步骤S2中脉冲压缩具体包括以下步骤：
S2.1、得匹配滤波函数：使得输出信噪比最大化的准则，可得到匹配滤波器的表达式：
hmatch(t)＝kx*(t0-t)
其中，k为任意常数，x*为线性调频信号的共轭，t0为预期偏移量；
S2.2、在匹配滤波器的基础上加Hanmming窗函数，得失配滤波器的表达式如下：
hmismatch(t)＝w(t)hmatch(t)
其中，Hanmming窗w(n)的定义为：
其中β＝0.46，N为Hanmming窗的总长度；
S2.3、脉冲压缩：接收探头收到的反射回波经过失配滤波器后，可以获得峰值功率较高的窄脉冲：
y(t)＝conv(r(t)，F-1(X*(f)Hmismatch(f)))
其中，r(t)为反射回波，F-1为逆傅里叶变换，X*(f)为线性调频信号经过傅里叶变换后的共轭，Hmismatch(f)为失配滤波器的***函数，conv(*)为卷积运算。
S3、分层反演：通过分析脉冲压缩后各反射回波的时延信息，反演出各反射层的反射系数，并对比已有组织的反射系数，定位到脉搏(血管)的位置信息，从而获得脉搏的深度、频率、幅度和组织状态，为测谎提供更多维度的参考；
该步骤S3中分层反演具体包括以下步骤：
S3.1、计算各反射层的相对距离：将各窄脉冲峰值对应的时延结果ni根据采样频率与声速转变为距离di：
其中，fs为***采样频率，C为超声波在人体组织内的平均声速；di即表示各反射面相对于皮肤的距离；
S3.2、衰减函数：超声波信号在人体内衰减函数为：
其中，U0为信号的初始幅值，f为信号的频率，αN为体内第N层组织的衰减系数，d为信号在体内移动的距离；
S3.3、原始发射信号的幅值X0经过了已知初始反射层的透射，再经过人体组织衰减函数U，得第一层反射层处的幅值X1：
其中，R0为已知初始反射层的反射系数，α1为初始反射层到第一层反射层之间的衰减系数，m1＝d1-d0为第一层反射层与初始反射层之间的距离；
S3.4、而反射回波Y0则经过人体组织衰减，再经过已知初始反射层的透射，得第一层反射层处的幅值Y1：
S3.5、因此第一层反射层的反射系数即为：
其中，X1原始发射信号X0在第一层反射层处的幅值，Y1为反射回波Y0在第一层反射层处的幅值；
S3.6、以此类推，通过迭代得到发射信号在第i层反射面处的幅值Xi与反射回波信号在第i层反射面处的幅值Yi：
其中，n为反射面的总数，mi＝di-di-1为第i层与第i-1层之间的距离，Ri-1为第i-1层的反射系数，αi为第i-1层反射层到第i层反射层之间的衰减系数；
S3.7、因此第i层反射面的反射系数为：
S3.8、得到了各层的反射系数以后，通过对比已有组织的反射系数，可以确定各反射层的组织构成，从而得到脉搏深度、频率、幅度信息等多维度信息；从而为测谎结果的判断，提供多维度的脉搏参考。
实施例二
本实施例提出一种基于线性调频的脉搏测谎装置，该装置将声学探测、电力电子、信号处理等技术有机地结合在一起，具有装置体积小、成本低的特点。并且根据超声波的物理特性，采用非阵列超声探头发射线性调频的连续信号，并在信号处理方法上使用预失真处理、脉冲压缩、分层反演回波算法相结合，最终得到脉搏的深度、频率、幅度和组织状态。下面结合图1-图4对本装置作进一步详细的描述说明。
本实施例中基于线性调频的脉搏测谎装置结构如图2所示，由数据采集端101、数据处理端102两部分组成，两个部分之间是通过无线电磁波相关联。其中数据采集端101负责超声波的收发，模拟信号与数字信号之间的转换，无线发射回波信号；数据处理端102负责通过对发射信号的预处理，提高轴向分辨力和降低旁瓣；同时对反射回波信号进行后处理，得到脉搏深度、频率、幅度等多维度的信息，以为测谎结果的判断提供更全面的参考。
其中，基于线性调频的脉搏测谎装置的数据采集端101的结构如图3所示，负责对数据的采集、数据的转发、***流程的时序控制。其中包括收发逻辑单元201、无线发射单元202、数据存储单元203、AD/DA转换单元204、功率放大单元205、非阵列超声探头单元206、低噪声放大单元207。其组成连接关系如下：
收发逻辑单元201与无线发射单元202、数据存储单元203、AD/DA转换单元204相连接，负责收发流程的时序控制，如发送信号的间隔控制、外部设备(存储器、AD/DA芯片、无线模块)的时序控制等；无线发射单元202与AD/DA转换单元204相连接，负责将数字量化后的超声波反射回波以无线电磁波的形式发送至数据处理端102用于分析处理；数据存储单元203与AD/DA转换单元204、收发逻辑单元201相连接，负责存储数据处理端102中预处理后的信号数据，收发逻辑单元201可通过一定的流程逻辑去读取使用其中的数据；AD/DA转换单元204与功率放大单元205、低噪声放大单元207、无线发射单元202、数据存储单元203相连接，负责将超声探头转换而来的模拟信号与处理芯片能识别的数字信号进行相互转换；功率放大单元205与非阵列超声探头单元206相连接，负责将低电压信号转换为可以驱动非阵列超声探头单元206的驱动信号；低噪声放大单元207与非阵列超声探头单元206相连接，负责通过低噪声放大器将由非阵列超声探头单元206转换过来的电信号进行放大；非阵列超声探头单元206负责实现电信号与非阵列超声信号之间的转换，即实现了电能与声能之间的转换。
其中，基于线性调频的脉搏测谎装置的数据处理端102的结构如图4所示，负责数据的分析处理处理，其包括信号预处理单元301、无线接收单元302、信号后处理单元303。其组成连接关系如下：
若数据存储单元203中无数据存储，或者需要更新发射信号时，信号预处理单元301与数据存储单元203相连接，负责对信号做预失真处理，从而提高轴向分辨力并减低旁瓣；无线接收单元302与信号后处理单元303相连接，负责接收来自数据采集端的超声回波数据；信号后处理单元303负责对无线接收单元302的数据进行后处理，得到脉搏深度、频率、幅度等多维度的信息，以为测谎结果的判断提供更全面的参考。
综上所述，上述实施例公开了一种基于线性调频的脉搏测谎方法与装置，传感器采用非阵列超声探头发射线性调频的连续信号，线性调频信号具有精确、抗多径、去伪影的优点。结合预失真处理、脉冲压缩、分层反演回波算法，通过各反射回波之间的时延，得到各反射面的位置(深度)；利用超声波在体内的衰减函数，通过迭代的方式反演出发射信号与回波信号在每个反射层的幅值，并计算出每一层反射面的反射系数；反射系数可确定各反射层的组织构成，得到脉搏的深度、频率、幅度等多维度的脉搏信息，从而可为测谎结果的判断提供更多维度的参考。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。
Claims (8)
1.一种基于线性调频的脉搏测谎方法，其特征在于，所述的脉搏测谎方法包括以下步骤：
S1、信号预处理：采用非阵列探头发射线性调频的连续信号，并通过幅度加权探头的脉冲响应，最后使用Lanczos窗函数作为幅度削减函数降低脉压后的旁瓣；
S2、脉冲压缩：对接收到的反射回波做脉冲压缩处理，以获得峰值功率较高的窄脉冲；并通过加Hanmming窗，降低脉冲压缩后的旁瓣；
S3、分层反演：通过分析脉冲压缩后各反射回波的时延信息，反演出各反射层的反射系数，并对比已有组织的反射系数，定位到脉搏的位置信息，获得脉搏的深度、频率、幅度和组织状态，进而得到多维度的测谎参考信息。
2.根据权利要求1所述的一种基于线性调频的脉搏测谎方法，其特征在于，所述的步骤S1中信号预处理的过程如下：
S1.1、确定信号形式，根据原始信号为线性调频信号，由用户的输入得到输入信号如下：
其中f0是信号的中心频率，K是信号带宽，T是信号的时间宽度，A是信号的幅度；
S1.2、幅度加权，对线性调频作幅度加权的预处理：
xpre(t)＝Mprob(t)Mdecr(t)x(t),0≤t≤T
其中，Mprob(t)是基于探头的幅度加权函数；
由于超声探头的脉冲响应可近似为高斯包络的余弦波：
其中，Δfu是探头的带宽，f′0是探头的中心频率，μ是一个常量；因此幅度加权函数可表示为高斯函数的倒数：
其中，D是加权系数，表示幅度加权的程度；
Mdecr(t)是Lanczos窗函数，表达式如下：
其中，
m为幅度削减长度。
3.根据权利要求1所述的一种基于线性调频的脉搏测谎方法，其特征在于，所述的步骤S2中脉冲压缩的过程如下：
S2.1、根据输出信噪比最大化的准则，得到匹配滤波器的表达式如下：
hmatch(t)＝kx*(t0-t)
其中，k为任意常数，x*为线性调频信号的共轭，t0为预期偏移量；
S2.2、在匹配滤波器的基础上加Hanmming窗，得失配滤波器的表达式如下：
hmismatch(t)＝w(t)hmatch(t)
其中，Hanmming窗w(n)的定义为：
其中β＝0.46，N为Hanmming窗的总长度；
S2.3、接收探头收到的反射回波经过失配滤波器进行脉冲压缩后，获得峰值功率较高的窄脉冲：
y(t)＝conv(r(t),F-1(X*(f)Hmismatch(f)))
其中，r(t)为反射回波，F-1为逆傅里叶变换，X*(f)为线性调频信号经过傅里叶变换后的共轭，Hmismatch(f)为失配滤波器的***函数，conv(*)为卷积运算。
4.根据权利要求1所述的一种基于线性调频的脉搏测谎方法，其特征在于，所述的步骤S3中分层反演的过程如下：
S3.1、计算各反射层的相对距离，将各窄脉冲峰值对应的时延结果ni根据采样频率与声速转变为距离di：
其中，fs为***采样频率，C为超声波在人体组织内的平均声速，di即表示各反射面相对于皮肤的距离；
S3.2、通过下式得到超声波信号在人体内衰减函数为：
其中，U0为信号的初始幅值，f为信号的频率，αN为体内第N层组织的衰减系数，d为信号在体内移动的距离；
S3.3、原始发射信号的幅值X0经过已知初始反射层的透射，再经过人体组织衰减函数U，得第一层反射层处的幅值X1：
其中，R0为已知初始反射层的反射系数，α1为初始反射层到第一层反射层之间的衰减系数，m1＝d1-d0为第一层反射层与初始反射层之间的距离；
S3.4、反射回波Y0经过人体组织衰减，再经过已知初始反射层的透射，得第一层反射层处的幅值Y1：
S3.5、计算得到第一层反射层的反射系数如下：
其中，X1原始发射信号X0在第一层反射层处的幅值，Y1为反射回波Y0在第一层反射层处的幅值；
S3.6、以此类推，通过迭代得到发射信号在第i层反射面处的幅值Xi与反射回波信号在第i层反射面处的幅值Yi：
其中，n为反射面的总数，mi＝di-di-1为第i层与第i-1层之间的距离，Ri-1为第i-1层的反射系数，αi为第i-1层反射层到第i层反射层之间的衰减系数；
S3.7、计算得到第i层反射面的反射系数如下：
S3.8、得到各层的反射系数以后，通过对比已有组织的反射系数，确定各反射层的组织构成，得到包括脉搏深度、频率、幅度信息在内的多维度信息，从而为测谎结果的判断提供多维度的脉搏参考。
5.根据权利要求2所述的一种基于线性调频的脉搏测谎方法，其特征在于，D的取值与探头的脉冲响应、信号的带宽和时宽相关。
6.一种基于线性调频的脉搏测谎装置，其特征在于，所述的脉搏测谎装置包括数据采集端和数据处理端，二者之间通过无线电磁波进行信息交互，其中，所述的数据采集端负责超声波的收发，模拟信号与数字信号之间的转换，无线发射回波信号；所述的数据处理端负责通过对发射信号的预处理，提高轴向分辨力和降低旁瓣；同时对反射回波信号进行后处理，得到包括脉搏深度、频率、幅度在内的多维度信息，以为测谎结果的判断提供参考。
7.根据权利要求6所述的一种基于线性调频的脉搏测谎装置，其特征在于，所述的数据采集端实现对数据的采集、数据的转发、***流程的时序控制，包括收发逻辑单元、无线发射单元、数据存储单元、AD/DA转换单元、功率放大单元、非阵列超声探头单元和低噪声放大单元，其组成连接关系如下：
所述的收发逻辑单元分别与无线发射单元、数据存储单元、AD/DA转换单元相连接，负责发送信号的间隔控制、外部设备的时序控制；所述的无线发射单元与AD/DA转换单元相连接，负责将数字量化后的超声波反射回波以无线电磁波的形式发送至数据处理端用于分析处理；所述的数据存储单元分别与AD/DA转换单元、收发逻辑单元相连接，负责存储数据处理端中预处理后的信号数据，所述的收发逻辑单元通过指定的流程逻辑去读取使用其中的数据；所述的AD/DA转换单元分别与功率放大单元、低噪声放大单元、无线发射单元、数据存储单元相连接，负责将超声探头转换而来的模拟信号与处理芯片能识别的数字信号进行相互转换；所述的功率放大单元与非阵列超声探头单元相连接，负责将低电压信号转换为驱动非阵列超声探头单元的驱动信号；所述的低噪声放大单元与非阵列超声探头单元相连接，负责通过低噪声放大器将由非阵列超声探头单元转换过来的电信号进行放大；所述的非阵列超声探头单元负责实现电信号与非阵列超声信号之间的转换，即实现电能与声能之间的转换。
8.根据权利要求6所述的一种基于线性调频的脉搏测谎装置，其特征在于，所述的数据处理端负责数据的分析处理，包括信号预处理单元、无线接收单元、信号后处理单元，其组成连接关系如下：
若数据存储单元中无数据存储，或者需要更新发射信号时，信号预处理单元与数据存储单元相连接，负责对信号做预失真处理，提高轴向分辨力并减低旁瓣；所述的无线接收单元与信号后处理单元相连接，负责接收来自数据采集端的超声回波数据；所述的信号后处理单元负责对无线接收单元的数据进行后处理，得到包括脉搏深度、频率、幅度在内的多维度信息，以为测谎结果的判断提供参考。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911036322.0A CN110680349B (zh) | 2019-10-29 | 2019-10-29 | 一种基于线性调频的脉搏测谎方法与装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911036322.0A CN110680349B (zh) | 2019-10-29 | 2019-10-29 | 一种基于线性调频的脉搏测谎方法与装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110680349A CN110680349A (zh) | 2020-01-14 |
| CN110680349B true CN110680349B (zh) | 2021-07-20 |
Family
ID=69114569
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911036322.0A Active CN110680349B (zh) | 2019-10-29 | 2019-10-29 | 一种基于线性调频的脉搏测谎方法与装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110680349B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114305225B (zh) * | 2021-12-31 | 2023-02-28 | 福建汉特云智能科技有限公司 | 一种扫地机器人垃圾箱的垃圾容量检测方法及其*** |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0240573A1 (en) * | 1985-09-27 | 1987-10-14 | Yokogawa Medical Systems, Ltd | Ultrasonic wave pulse doppler apparatus |
| CN101474079A (zh) * | 2009-01-06 | 2009-07-08 | 华南理工大学 | 基于超声稳态反射波的生物体声阻抗的测量方法与装置 |
| CN204336935U (zh) * | 2014-12-24 | 2015-05-20 | 重庆大学 | 一种基于脉搏、呼吸和皮肤电阻的多道生理心理测试仪 |
| CN105105739A (zh) * | 2015-04-01 | 2015-12-02 | 上海兆观信息科技有限公司 | 短距离无线心率及心率变异性检测方法 |
| CN105338908A (zh) * | 2013-05-09 | 2016-02-17 | 爱飞纽医疗机械贸易有限公司 | 超声波优化方法和用于该超声波优化方法的超声波医疗装置 |
| CN106908768A (zh) * | 2017-02-24 | 2017-06-30 | 四川九洲空管科技有限责任公司 | 一种超低旁瓣的脉冲压缩方法 |
| CN109106383A (zh) * | 2017-06-22 | 2019-01-01 | 罗杰谊 | 情绪感测***及方法 |
-
2019
- 2019-10-29 CN CN201911036322.0A patent/CN110680349B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0240573A1 (en) * | 1985-09-27 | 1987-10-14 | Yokogawa Medical Systems, Ltd | Ultrasonic wave pulse doppler apparatus |
| CN101474079A (zh) * | 2009-01-06 | 2009-07-08 | 华南理工大学 | 基于超声稳态反射波的生物体声阻抗的测量方法与装置 |
| CN105338908A (zh) * | 2013-05-09 | 2016-02-17 | 爱飞纽医疗机械贸易有限公司 | 超声波优化方法和用于该超声波优化方法的超声波医疗装置 |
| CN204336935U (zh) * | 2014-12-24 | 2015-05-20 | 重庆大学 | 一种基于脉搏、呼吸和皮肤电阻的多道生理心理测试仪 |
| CN105105739A (zh) * | 2015-04-01 | 2015-12-02 | 上海兆观信息科技有限公司 | 短距离无线心率及心率变异性检测方法 |
| CN106908768A (zh) * | 2017-02-24 | 2017-06-30 | 四川九洲空管科技有限责任公司 | 一种超低旁瓣的脉冲压缩方法 |
| CN109106383A (zh) * | 2017-06-22 | 2019-01-01 | 罗杰谊 | 情绪感测***及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110680349A (zh) | 2020-01-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cespedes et al. | Theoretical bounds on strain estimation in elastography | |
| O'Donnell et al. | Internal displacement and strain imaging using ultrasonic speckle tracking | |
| Kuc | Clinical application of an ultrasound attenuation coefficient estimation technique for liver pathology characterization | |
| Varghese et al. | Noise reduction in elastograms using temporal stretching with multicompression averaging | |
| US7535797B2 (en) | High-resolution ultrasound displacement measurement apparatus and method | |
| CN102458260B (zh) | 测量心脏收缩功能的方法和装置 | |
| Varghese et al. | Performance optimization in elastography: Multicompression with temporal stretching | |
| KR101048554B1 (ko) | 초음파를 이용한 생체 조직 탄성도 및 경화 측정 시스템 | |
| US20060074320A1 (en) | Home ultrasound system | |
| US4881549A (en) | Apparatus and method for obtaining ultrasonic backcatter measurement from tissue | |
| CN107928674B (zh) | 一种基于声波测距的非接触式呼吸检测方法 | |
| Ranta et al. | Digestive activity evaluation by multichannel abdominal sounds analysis | |
| EP0154869A1 (en) | Ultrasonic measurement apparatus | |
| JP2004520094A (ja) | 超音波トモグラフ | |
| CN101884526A (zh) | 基于超声血流信息的动脉血压测量装置 | |
| KR970701017A (ko) | 골조송증 진단장치 및 방법 | |
| CN111012319B (zh) | 皮肤血流和血管的监测成像方法、***及存储介质 | |
| Chernyakova et al. | Fourier-domain beamforming and structure-based reconstruction for plane-wave imaging | |
| US20140187951A1 (en) | Method of remove blackground noise in short pulse excitation of arfi | |
| Shi et al. | Neural network based real-time heart sound monitor using a wireless wearable wrist sensor | |
| CN110680349B (zh) | 一种基于线性调频的脉搏测谎方法与装置 | |
| CA3105427A1 (en) | Apparatus and method for characterization of a ductile membrane, surface and sub-surface properties | |
| Steinberg et al. | Continuous artery monitoring using a flexible and wearable single-element ultrasonic sensor | |
| Fang et al. | Wrist pulse recording with a wearable piezoresistor-piezoelectret compound sensing system and its applications in health monitoring | |
| CN111685760B (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 |