CN110196098B - 基于心率变化的车辆声品质评价方法 - Google Patents

基于心率变化的车辆声品质评价方法 Download PDF

Info

Publication number
CN110196098B
CN110196098B CN201910435656.9A CN201910435656A CN110196098B CN 110196098 B CN110196098 B CN 110196098B CN 201910435656 A CN201910435656 A CN 201910435656A CN 110196098 B CN110196098 B CN 110196098B
Authority
CN
China
Prior art keywords
heart rate
vehicle
sound quality
data
evaluator
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
Application number
CN201910435656.9A
Other languages
English (en)
Other versions
CN110196098A (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.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
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 Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201910435656.9A priority Critical patent/CN110196098B/zh
Publication of CN110196098A publication Critical patent/CN110196098A/zh
Application granted granted Critical
Publication of CN110196098B publication Critical patent/CN110196098B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M17/00Testing of vehicles
    • G01M17/007Wheeled or endless-tracked vehicles
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • G10L25/51Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • G10L25/51Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination
    • G10L25/66Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination for extracting parameters related to health condition

Abstract

本发明属于汽车试验技术领域,特别涉及一种基于心率变化的车辆声品质评价方法,包括:一、选择车辆声品质评价员;二、采集待评价工况条件下的待评价车辆的车内音频数据
Figure DDA0002070431710000011
作为待评价样本;三、采集回放车内音频数据A0时评价员心率数据L1,记录平均周期T1、心率R峰均值H1;四、采集不回放车内音频数据A0时评价员的心率数据L2;记录平均周期T2、心率R峰均值H2;五、计算平均周期的权重ω1和心率R峰均值的权重ω2,进而获得车辆声品质量化指标以判断车辆声品质。本发明操作简单,使用仪器价格较便宜,克服了主观评价中评价员因带有个人的情感而造成数据不准确、不合理的缺点,可接受性和可实践性较强。

Description

基于心率变化的车辆声品质评价方法
技术领域
本发明属于汽车试验技术领域,特别涉及一种基于心率变化的车辆声品质评价方法。
背景技术
目前,对车辆声品质的评价方法可分为主观评价和客观评价两种。主观评价是以人为主体,通过问卷调查或评审团主观评价试验的形式进行评价。客观评价是以声音的物理参数以及心理学客观参数为基础,通过测试分析响度、尖锐度、粗糙度、波动度等声音的客观参数来评价声音的优劣。
但是,主观评价的方法主要根据问卷调查方式、成对比较方式、排序方式对车辆声品质进行评价,受评价方式的限制,主观评价方法仅能进行多样本比较或等级划分,而无法给出具体的比例尺度。车辆声品质客观评价方法是以声音固有的指标为研究对象,并着力寻找对心理影响较大的声音指标从而提高车辆声品质客观评价有效性和准确性,在时频域内声音固有的物理指标众多,不同指标排列组合得到的综合指标类别、数量巨大,因而对车辆声品质客观评价指标的选取存在大量分歧。
专利号为ZL201611006119.5的中国专利公开了一种基于瞳孔直径的评价方法,是为避免照度变化对评价数据的影响,需建立照度-瞳孔直径模型,照度-瞳孔直径模型的建立使得模型更加复杂、操作更加繁琐、增大了数据库数据。
为此,人们急需寻求一种能够客观的、有效的、准确的方法评价车辆声品质的方法。一些研究发现,当人们受到外界刺激变化时,情绪发生变化时心率会发生变化。公告号为CN105982678A的专利公开了一种根据心率与呼吸判断情绪的方法,克服了以往测定情绪而出现的浮动状况,过程简单,结果直观。噪音干扰时,人的情绪发生变化,进而人的心率也发生变化。
发明内容
为解决上述技术问题,本发明的目的在于:提供一种基于心率变化的车辆声品质评价方法,操作简单,使用仪器价格较便宜,克服了主观评价中评价员因带有个人的情感而造成数据不准确、不合理的缺点,可接受性和可实践性较强。
本发明为解决其技术问题所采用的技术方案为:
所述基于心率变化的车辆声品质评价方法,包括以下步骤:
步骤一、选择车辆声品质评价员;
步骤二、将音频样本采集装置固定于待评价工况条件下的待评价位置处,采集音频数据
Figure BDA0002070431690000022
作为待评价样本;
步骤三、评价员坐于车辆座椅处,佩戴高保真耳机,回放所采集的车内音频数据A0,使用心率采集装置采集评价员在测试时间t内的心率数据L1,录平均周期T1,心率R峰均值H1
步骤四、评价员坐于车辆座椅处,佩戴高保真耳机,不播放音频,使用心率采集装置采集评价员在测试时间t内的心率数据L2,记录平均周期T2,心率R峰均值H2
步骤五、通过数据处理装置利用心率数据L1与心率数据L2计算平均周期的权重ω1和心率R峰均值的权重ω2,进而获得车辆声品质量化指标以判断车辆声品质。
其中,优选方案为:
所述步骤三中平均周期T1是指播放车内音频数据A0的测试时间t与该时间段内完整最小周期个数n1的比值,T1=t/n1;步骤四中,平均周期T2是指不播放音频数据时测试时间t与该时间段内完整最小周期个数n2的比值,T2=t/n2
所述步骤三中,心率R峰均值H1是指播放音频数据A0时在测试时间t内所有R峰的平均值;所述步骤四中,心率R峰均值H2是指不播放音频时在测试时间t内所有R峰的平均值。
所述步骤三和步骤四中的心率采集装置采用连接在PC端的快速心电检测仪。
所述步骤五中权重ω1和ω2的确定使用熵值法,具体步骤如下:
a.选取心率数据L1记录的平均周期T1和心率R峰均值H1,选取心率数据L2记录的平均周期T2和心率R峰均值H2,定义xij为第i个心率的第j个指标的数值(i=1,2;j=1,2),其中x11=T1、x12=H1、x21=T2、x22=H2
b.指标的归一化处理,各项指标的计量单位并不统一,通过下式对它们先进行标准化处理,即把指标的绝对值转化为相对值:
Figure BDA0002070431690000021
x′ij为第i个心率的第j个指标的相对值(i=1,2;j=1,2);
c.计算第j项指标下第i个心率占该指标的比重:
Figure BDA0002070431690000031
d.计算第j项指标的熵值:
Figure BDA0002070431690000032
其中k=1/ln(2)>0,满足ej≥0;
e.计算信息熵冗余度:
dj=1-ej
f.计算各项指标的权值:
平均周期的权重:
Figure BDA0002070431690000033
心率R峰均值的权重:
Figure BDA0002070431690000034
所述步骤五中,车辆声品质量化指标为p,p=ω1T2/T12H1/H2,当p小于设定阈值时,说明车辆声品质较好;当p大于设定阈值时,说明车辆声品质较差。
与现有技术相比,本发明具有以下有益效果:
1、本发明首次提出基于心率变化的车辆声品质评价方法,该方法操作简单,使用仪器价格较便宜,可接受性、可实践性较强。
2、本发明首次提出基于心率变化的车辆声品质评价方法,该方法克服了主观评价中评价员因带有个人的情感而造成数据不准确、不合理的缺点。
3、本发明所使用仪器在消除误差干扰后均可多次使用,降低车辆声品质评价的成本,且可回收,高效环保。
附图说明
图1是本发明流程图。
图2是本发明实施例1的心率数据L1和心率数据L2监测图。
具体实施方式
下面结合附图对本发明实施例做进一步描述:
实施例1:
如图1-2所示,本发明所述基于心率变化的车辆声品质评价方法,包括以下步骤:
步骤一、选择车辆声品质评价员;
步骤二、将HEAD acoustics公司生产的HMS III仿真人工头固定于待评价机动车待评价工况条件下的驾驶人位置处,使用SQLabII多通道测试系统采集音频数据A0作为待评价样本;
步骤三、评价员坐于车辆座椅处,佩戴开放式高阻抗的(100Ω-300Ω)的HA III高保真耳机,回放所采集的车内音频数据A0,使用力康快速心电检测仪采集评价员在测试时间t内的心率数据L1,记录平均周期T1,心率R峰均值H1,力康快速心电检测仪连接到PC端;
步骤四、评价员坐于车辆座椅处,佩戴开放式高阻抗的(100Ω-300Ω)的HA III高保真耳机,不播放音频,使用力康快速心电检测仪采集评价员在测试时间t内的心率数据L2;记录平均周期T2,心率R峰均值H2,力康快速心电检测仪连接到PC端;
步骤五、通过数据处理装置利用心率数据L1与心率数据L2计算平均周期的权重ω1和心率R峰均值的权重ω2,进而获得车辆声品质量化指标以判断车辆声品质,数据处理装置一般采用MATLAB软件进行仿真计算。
其中,步骤三中平均周期T1是指播放车内音频数据A0的测试时间t与该时间段内完整最小周期个数n1的比值,T1=t/n1;步骤四中,平均周期T2是指不播放音频数据时测试时间t与该时间段内完整最小周期个数n2的比值,T2=t/n2;步骤三中,心率R峰均值H1是指播放音频数据A0时在测试时间t内所有R峰的平均值;所述步骤四中,心率R峰均值H2是指不播放音频时在测试时间t内所有R峰的平均值。
所述步骤五中权重ω1和ω2的确定使用熵值法,具体步骤如下:
a.选取心率数据L1记录的平均周期T1和心率R峰均值H1,选取心率数据L2记录的平均周期T2和心率R峰均值H2,定义xij为第i个心率的第j个指标的数值(i=1,2;j=1,2),其中x11=T1、x12=H1、x21=T2、x22=H2
b.指标的归一化处理,各项指标的计量单位并不统一,通过下式对它们先进行标准化处理,即把指标的绝对值转化为相对值:
Figure BDA0002070431690000051
x′ij为第i个心率的第j个指标的相对值(i=1,2;j=1,2);
c.计算第j项指标下第i个心率占该指标的比重:
Figure BDA0002070431690000052
d.计算第j项指标的熵值:
Figure BDA0002070431690000053
其中,k=1/ln(2)>0,满足ej≥0;
e.计算信息熵冗余度:
dj=1-ej
f.计算各项指标的权值:
平均周期的权重:
Figure BDA0002070431690000054
心率R峰均值的权重:
Figure BDA0002070431690000055
所述步骤五中,车辆声品质量化指标为p,p=ω1T2/T12H1/H2,当p小于设定阈值时,说明车辆声品质较好;当p大于设定阈值时,说明车辆声品质较差,车辆声品质量化指标p的阈值一般取1.5。

Claims (4)

1.一种基于心率变化的车辆声品质评价方法,其特征在于,包括以下步骤:
步骤一、选择车辆声品质评价员;
步骤二、将音频样本采集装置固定于待评价工况条件下的待评价位置处,采集音频数据作为待评价样本;
步骤三、评价员坐于车辆座椅处,佩戴高保真耳机,回放所采集的车内音频数据A0,使用心率采集装置采集评价员在测试时间t内的心率数据L1,记录平均周期T1,心率R峰均值H1
步骤四、评价员坐于车辆座椅处,佩戴高保真耳机,不播放音频,使用心率采集装置采集评价员在测试时间t内的心率数据L2,记录平均周期T2,心率R峰均值H2
步骤五、通过数据处理装置利用心率数据L1与心率数据L2计算平均周期的权重ω1和心率R峰均值的权重ω2,进而获得车辆声品质量化指标以判断车辆声品质;
所述步骤五中权重ω1和ω2的确定使用熵值法,具体步骤如下:
a.选取心率数据L1记录的平均周期T1和心率R峰均值H1,选取心率数据L2记录的平均周期T2和心率R峰均值H2,定义xij为第i个心率的第j个指标的数值,其中i=1,2;j=1,2;x11=T1、x12=H1、x21=T2、x22=H2
b.指标的归一化处理,各项指标的计量单位并不统一,通过下式对它们先进行标准化处理,即把指标的绝对值转化为相对值:
Figure FDA0003241650610000011
x′ij为第i个心率的第j个指标的相对值,其中i=1,2;j=1,2;
c.计算第j项指标下第i个心率占该指标的比重:
Figure FDA0003241650610000012
d.计算第j项指标的熵值:
Figure FDA0003241650610000013
其中k=1/ln(2)>0,满足ej≥0;
e.计算信息熵冗余度:
dj=1-ej
f.计算各项指标的权值:
平均周期的权重:
Figure FDA0003241650610000014
心率R峰均值的权重:
Figure FDA0003241650610000021
所述步骤五中,车辆声品质量化指标为p,p=ω1T2/T12H1/H2,当p小于设定阈值时,说明车辆声品质较好;当p大于设定阈值时,说明车辆声品质较差。
2.根据权利要求1所述的基于心率变化的车辆声品质评价方法,其特征在于,所述步骤三中平均周期T1是指播放车内音频数据A0的测试时间t与该时间段内完整最小周期个数n1的比值,T1=t/n1;步骤四中,平均周期T2是指不播放音频数据时测试时间t与该时间段内完整最小周期个数n2的比值,T2=t/n2
3.根据权利要求1所述的基于心率变化的车辆声品质评价方法,其特征在于,所述步骤三中,心率R峰均值H1是指播放音频数据A0时在测试时间t内所有R峰的平均值;所述步骤四中,心率R峰均值H2是指不播放音频时在测试时间t内所有R峰的平均值。
4.根据权利要求1所述的基于心率变化的车辆声品质评价方法,其特征在于,所述步骤三和步骤四中的心率采集装置采用连接在PC端的快速心电检测仪。
CN201910435656.9A 2019-05-23 2019-05-23 基于心率变化的车辆声品质评价方法 Active CN110196098B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910435656.9A CN110196098B (zh) 2019-05-23 2019-05-23 基于心率变化的车辆声品质评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910435656.9A CN110196098B (zh) 2019-05-23 2019-05-23 基于心率变化的车辆声品质评价方法

Publications (2)

Publication Number Publication Date
CN110196098A CN110196098A (zh) 2019-09-03
CN110196098B true CN110196098B (zh) 2021-10-29

Family

ID=67753034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910435656.9A Active CN110196098B (zh) 2019-05-23 2019-05-23 基于心率变化的车辆声品质评价方法

Country Status (1)

Country Link
CN (1) CN110196098B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111174997B (zh) * 2020-01-14 2021-10-22 合肥工业大学 一种基于心率变化的楼板振动舒适度的初步测试方法
CN112345063B (zh) * 2020-10-14 2022-08-02 合肥工业大学 一种楼板振动舒适度的测试方法
CN115444381A (zh) * 2022-09-19 2022-12-09 吉林大学 一种基于生理信号的车内声品质评价方法和系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928071A (zh) * 2012-10-25 2013-02-13 北京市市政工程研究院 基于心电指标的道路交通噪声检测系统和方法
CN106725289A (zh) * 2016-11-16 2017-05-31 山东理工大学 一种基于瞳孔直径的车辆声品质评价方法
CN107913075A (zh) * 2017-11-15 2018-04-17 重庆邮电大学 一种基于多参数的精神压力评估装置及其评估方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4965162B2 (ja) * 2006-05-10 2012-07-04 トヨタ自動車株式会社 車両用不整脈監視装置
JP5919722B2 (ja) * 2011-10-18 2016-05-18 株式会社豊田中央研究所 生体信号推定装置及びプログラム
CN103091480A (zh) * 2013-01-07 2013-05-08 河北工业大学 基于熵权的地下道路沥青路面使用性能评价方法
CN103471709A (zh) * 2013-09-17 2013-12-25 吉林大学 乘用车车内噪声声品质预测方法
CN107235045A (zh) * 2017-06-29 2017-10-10 吉林大学 考虑生理与操控信息的驾驶员路怒状态车载识别交互系统
CN108630225A (zh) * 2018-03-29 2018-10-09 太原理工大学 基于模糊综合评价的听障儿童元音评估方法
CN108765876A (zh) * 2018-05-31 2018-11-06 东北大学 基于多模信号的驾驶疲劳深度分析预警系统及方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928071A (zh) * 2012-10-25 2013-02-13 北京市市政工程研究院 基于心电指标的道路交通噪声检测系统和方法
CN106725289A (zh) * 2016-11-16 2017-05-31 山东理工大学 一种基于瞳孔直径的车辆声品质评价方法
CN107913075A (zh) * 2017-11-15 2018-04-17 重庆邮电大学 一种基于多参数的精神压力评估装置及其评估方法

Also Published As

Publication number Publication date
CN110196098A (zh) 2019-09-03

Similar Documents

Publication Publication Date Title
CN110196098B (zh) 基于心率变化的车辆声品质评价方法
Matos et al. Detection of cough signals in continuous audio recordings using hidden Markov models
Deliyski et al. Adverse effects of environmental noise on acoustic voice quality measurements
US20130166291A1 (en) Emotional and/or psychiatric state detection
CN112353391A (zh) 一种基于脑电信号的汽车车内声音品质识别方法及装置
CN109285551B (zh) 基于wmfcc和dnn的帕金森患者声纹识别方法
WO2019216320A1 (ja) 機械学習装置、解析装置、機械学習方法および解析方法
CN110234279B (zh) 表征睡眠呼吸障碍的方法
CN111951824A (zh) 一种基于声音判别抑郁症的检测方法
CN103093759A (zh) 一种基于移动终端的嗓音检测评估装置及方法
CN110567575B (zh) 一种汽车门锁闭合声品质评价方法
Noffs et al. Speech metrics, general disability, brain imaging and quality of life in multiple sclerosis
JP2024020321A (ja) 精神・神経系疾患を推定する装置
CN113674767A (zh) 一种基于多模态融合的抑郁状态识别方法
KR102218214B1 (ko) 추정 장치, 추정 프로그램이 기록된 기록 매체, 추정 장치의 작동 방법 및 추정 시스템
CN110068466B (zh) 基于脑电波的车辆声品质评价方法
CN107910019B (zh) 一种人体声音信号处理及分析方法
CN106725289B (zh) 一种基于瞳孔直径的车辆声品质评价方法
Banks et al. Sports related concussion impacts speech rate and muscle physiology
JP5414160B2 (ja) 感性評価装置及び方法
CN113571050A (zh) 一种基于Attention与Bi-LSTM的语音抑郁状态识别方法
CN111798109B (zh) 一种材料摩擦异响匹配方法
Kala et al. Reference free auscultation quality metric and its trends
CN110530653B (zh) 一种汽车声品质的主观评价方法
Bayani et al. Analysis of sound characteristics to design an annoyance metric for rattle sounds in the automotive industry

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