车内音响音质的客观评价方法、系统、设备及存储介质与流程

1.本发明涉及汽车音响技术领域，尤其涉及一种车内音响音质的客观评价方法、一种车内音响音质的客观评价系统、一种计算机设备及一种计算机可读存储介质。


背景技术：

2.随着汽车行业的高速发展，顾客对汽车的舒适性要求也随之提高，声品质已成为产品竞争力的重要因素，因此，声品质的研究与应用成为汽车声学系统的研究热点。
3.整车音效作为音响最终的表现形式，受到消费者最为直观地关注，其好坏已成为定位一款车型品质的重要指标，因此，如何快速客观地评价车内音响声品质成为了急需解决的问题。
4.然而汽车内声场环境较为复杂，影响了标准客观指标的有效性，使其缺少鲁棒性、可重复性和感知相关性。同时，汽车车内声场与一般房间环境不同，车内声场环境属于小空间，具有高反射面，几何形状复杂，次优布置的扬声器，声学路径不对称等特点，不能直接沿用厅堂声学的客观指标。
5.因此，建立一套适用于车内音响音质的评价方法，具有重要的研究意义和经济价值。


技术实现要素：

6.本发明所要解决的技术问题在于，提供一种车内音响音质的客观评价方法、系统、计算机设备及计算机可读存储介质，可在车内测量并计算音质感客观参数，能够快速有效地评价整车音响音质的好坏。
7.为了解决上述技术问题，本发明提供了一种车内音响音质的客观评价方法，包括：获取车内录制信号及播放源信号；根据所述车内录制信号及播放源信号计算音质感客观参数，所述音质感客观参数包括音质失真参数和频率响应参数；根据预设权重系数对所述音质失真参数和所述频率响应参数进行加权处理，生成并输出评价结果。
8.作为上述方案的改进，所述根据车内录制信号及播放源信号计算音质感客观参数的步骤包括：将所述车内录制信号与播放源信号进行对齐处理；根据对齐后的车内录制信号与播放源信号提取阶跃扫频信号；将所述阶跃扫频信号的每一段单频信号，使用平顶窗进行加窗并进行fft变换，计算频率及谐波的幅值；根据所述幅值及预设参数定义公式计算所述音质失真参数，所述音质失真参数包括低频谐波失真参数、中频谐波失真参数及高频谐波失真参数；根据对齐后的车内录制信号与播放源信号提取白噪声信号；从所述白噪声信号中截取目标信号，使用汉宁窗进行加窗并进行短时fft变换，计算频率响应；根据所述频率响应及预设参数定义公式计算所述频率响应参数，所述频率响应参数包括绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心。
9.作为上述方案的改进，所述根据幅值及预设参数定义公式，计算低频谐波失真参数、中频谐波失真参数及高频谐波失真参数的步骤包括：根据下述公式，计算低频谐波失真
参数、中频谐波失真参数及高频谐波失真参数；
[0010][0011][0012][0013]
其中，thd_l为低频谐波失真参数，thd_m为中频谐波失真参数，thd_h为高频谐波失真参数，n为相应频段测量的频率点的总数，fn为第n个频率点对应的频率值，thd(fn)为fn下的谐波失真。
[0014]
作为上述方案的改进，所述根据频率响应及预设参数定义公式计算绝对平均偏差的步骤包括：根据下述公式，计算绝对平均偏差；
[0015][0016]
其中，aad为绝对平均偏差，n为300hz-16khz之间的1/20倍频程带的总数，y
band_n
为第n个倍频程带对应的平均幅度，y
ref@1-3khz
为1khz-3khz的参考频带的平均幅度。
[0017]
作为上述方案的改进，所述根据频率响应及预设参数定义公式计算低频延展参数的步骤包括：根据下述公式，计算低频延展参数；
[0018][0019]
其中，lfx为计算低频延展参数，y
ref
为参考频带内的平均幅度再加上10db，x-6db
为在300hz以下并低于y
ref
值6db的第一个频率。
[0020]
作为上述方案的改进，所述根据频率响应及预设参数定义公式计算高频延展参数的步骤包括：根据下述公式，计算高频延展参数；
[0021][0022]
其中，hfx为高频延展参数，y
ref
为参考频带内的平均幅度，x-6db
为在5khz以上并低于y
ref
值6db的第一个频率。
[0023]
作为上述方案的改进，所述根据频率响应及预设参数定义公式计算频谱平坦度的步骤包括：根据下述公式，计算频谱平坦度；
[0024][0025]
其中，spf为频谱平坦度，sk为信号幅度谱在频率点k上的值，k为频带计算范围的
上限频率。
[0026]
作为上述方案的改进，所述根据频率响应及预设参数定义公式计算频谱重心的步骤包括：根据下述公式，计算频谱重心；
[0027][0028]
spc＝log
10
(spectral_centroid)，
[0029]
其中，spc为频谱重心，sk(tm)为在时间tm的信号加窗后的幅度谱在频率点k上的值，fk为k对应的频率值，k为频带计算范围的上限频率点，spectral_centroid为每个时间tm计算结果的平均。
[0030]
相应地，本发明还提供了一种车内音响音质的客观评价系统，包括：采集模块，用于获取车内录制信号及播放源信号；计算模块，用于根据所述车内录制信号及播放源信号计算音质感客观参数，所述音质感客观参数包括音质失真参数和频率响应参数；评价结果输出模块，用于根据预设权重系数对所述音质失真参数和所述频率响应参数进行加权处理，生成并输出评价结果。
[0031]
作为上述方案的改进，所述计算模块包括：对齐单元，用于将所述车内录制信号与播放源信号进行对齐处理；第一提取单元，用于根据对齐后的车内录制信号与播放源信号提取阶跃扫频信号；第一变换单元，用于将所述阶跃扫频信号的每一段单频信号，使用平顶窗进行加窗并进行fft变换，计算频率及谐波的幅值；第一计算单元，用于根据所述幅值及预设参数定义公式计算所述音质失真参数，所述音质失真参数包括低频谐波失真参数、中频谐波失真参数及高频谐波失真参数；第二提取单元，用于根据对齐后的车内录制信号与播放源信号提取白噪声信号；第二变换单元，用于从所述白噪声信号中截取目标信号，使用汉宁窗进行加窗并进行短时fft变换，计算频率响应；第二计算单元，用于根据所述频率响应及预设参数定义公式计算所述频率响应参数，所述频率响应参数包括绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心。
[0032]
相应地，本发明还提供了一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现上述车内音响音质的客观评价方法的步骤。
[0033]
相应地，本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述车内音响音质的客观评价方法的步骤。
[0034]
实施本发明，具有如下有益效果：
[0035]
本发明在车内音响客观评价的音质感上，选择8个音质感客观参数进行计算，可在不同的汽车音响系统之间体现出显著差异，适用于车内小空间声场环境。
[0036]
进一步，本发明结合一般车内音响系统分频方式，将整车频率响应曲线分为3个频段，分别为低频(50-300hz)，中频(300-3000hz)，高频(3000-20000hz)，从而在非线性的谐波失真中提取出低频谐波失真参数、中频谐波失真参数及高频谐波失真参数。同时，本发明参考上述频段划分方式，形成5个由频率响应二次计算获取的音质感客观参数(绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心)。通过在车内测量并计算音质感
客观参数，从能够通过客观测量进一步预估主观评价结果，可快速有效地评价整车音响音质感的好坏。
附图说明
[0037]
图1是本发明车内音响音质的客观评价方法的第一实施例流程图；
[0038]
图2是本发明车内音响音质的客观评价方法的第二实施例流程图；
[0039]
图3是本发明车内音响音质的客观评价系统的结构示意图；
[0040]
图4是本发明车内音响音质的客观评价系统中计算单元的结构示意图。
具体实施方式
[0041]
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步地详细描述。
[0042]
参见图1，图1显示了本发明车内音响音质的客观评价方法的第一实施例流程图，其包括：
[0043]
s101，获取车内录制信号及播放源信号。
[0044]
工作时，在车内播放事先准备好的“播放源信号”，同时通过录制设备录制车内的声音信号以形成“车内录制信号”。也就是说，播放源信号是指输入车辆音响设置中的音频文件；车内录制信号是指通过车内音响播放后，再录制的音频文件。
[0045]
s102，根据所述车内录制信号及播放源信号计算音质感客观参数。
[0046]
通过将车内录制信号与播放源信号进行对比分析，从而根据预设参数定义公式，计算出音质感客观参数。
[0047]
需要说明的是，车内音响系统与一般环境不同，车内反射物较多，导致频响的峰谷值较多，同时车内音响系统一般会通过均衡的方式提升低频辐射性能，因此，选择的音质感客观参数定义需要考虑车内声场情况，以进一步提升音质感客观参数的有效性。
[0048]
现有技术中，一般沿用厅堂声学相关参数在车内的声场进行评价。与现有技术不同的是，本发明在车内音响客观评价的音质感上，选择8个音质感客观参数进行计算，可在不同的汽车音响系统之间体现出显著差异，适用于车内小空间声场环境。
[0049]
具体地，所述音质感客观参数包括低频谐波失真参数、中频谐波失真参数、高频谐波失真参数、绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心。其中，低频谐波失真参数是指在低频频段下的平均谐波失真；中频谐波失真参数是指在中频频段下的平均谐波失真；高频谐波失真参数是指在高频频段下的平均谐波失真；绝对平均偏差用于表示各倍频程频带与参考频带的绝对幅度差，绝对平均偏差的值越高，则表示幅度与参考频带的偏差越大；频谱平坦度即信号功率谱的几何平均值与算术平均值的比值，该值越接近1，表明平坦度越高。
[0050]
s103，根据预设权重系数对所述音质感客观参数进行加权处理，生成并输出评价结果。
[0051]
通过预设权重系数与上述音质感客观参数构建模型，即可得到评价结果。需要说明的是，所述权重系数可根据实际情况进行预设，灵活性强。
[0052]
例如，基础权重系数为a，绝对平均偏差对应的权重系数为b，低频延展参数对应的
权重系数为c，高频延展参数对应的权重系数为d，频谱平坦度对应的权重系数为e，频谱重心对应的权重系数为f，低频谐波失真参数对应的权重系数为g，中频谐波失真参数对应的权重系数为h、高频谐波失真参数对应的权重系数为i，则评价结果s＝a b
×
绝对平均偏差 c
×
低频延展参数 d
×
高频延展参数 e
×
频谱平坦度 f
×
频谱重心 g
×
低频谐波失真参数 h
×
中频谐波失真参数 i
×
高频谐波失真参数。
[0053]
因此，本发明创新地将音质感从基础测量中分离进行评价，在车内测量并计算音质感客观参数，从而能够通过客观测量进一步预估主观评价结果，可快速有效地评价整车音响音质感的好坏。
[0054]
参见图2，图2显示了本发明车内音响音质的客观评价方法的第二实施例流程图，其包括：
[0055]
s201，获取车内录制信号及播放源信号。
[0056]
s202，将所述车内录制信号与播放源信号进行对齐处理。
[0057]
将车内录制信号与播放源信号进行对比分析，并根据播放源信号进行信号对齐。
[0058]
s203，根据对齐后的车内录制信号与播放源信号提取阶跃扫频信号。
[0059]
s204，将所述阶跃扫频信号的每一段单频信号，使用平顶窗进行加窗并进行fft变换，计算频率及谐波的幅值。
[0060]
将对齐后的阶跃扫频信号的每一段单频信号，使用时长1s(离散采样点48000)的平顶窗进行加窗，计算fft变换，计算该频率下的幅值，以及各个谐波的幅值。其中，fft的点数与平顶窗的窗长相等。
[0061]
s205，根据所述幅值及预设参数定义公式，计算低频谐波失真参数、中频谐波失真参数及高频谐波失真参数。
[0062]
谐波失真的相关参数包括低频谐波失真参数、中频谐波失真参数及高频谐波失真参数，其中，低频谐波失真参数是指在低频频段下的平均谐波失真，中频谐波失真参数是指在中频频段下的平均谐波失真，高频谐波失真参数是指在高频频段下的平均谐波失真。具体地，可根据下述公式，计算低频谐波失真参数thd_l、中频谐波失真参数thd_m及高频谐波失真参数thd_h；
[0063][0064][0065][0066]
其中，n为相应频段测量的频率点的总数，fn为第n个频率点对应的频率值，thd(fn)为fn下的谐波失真。
[0067]
因此，本发明结合一般车内音响系统分频方式，将整车频率响应曲线分为3个频段，分别为低频(50-300hz)，中频(300-3000hz)，高频(3000-20000hz)，从而在非线性的谐
波失真中提取出低频谐波失真参数、中频谐波失真参数及高频谐波失真参数。
[0068]
s206，根据对齐后的车内录制信号与播放源信号提取白噪声信号。
[0069]
s207，从所述白噪声信号中截取目标信号，使用汉宁窗进行加窗并进行短时fft变换，计算频率响应。
[0070]
将对齐后的白噪声信号，截取一段目标信号(10s)进行短时fft变换。使用时长1s(离散采样点48000)的汉宁窗，其中，窗移动重叠率为66.7％，fft的点数与窗长相等；移动窗25次，进行25次fft变换，对得到的25次结果进行平均，再根据测量量程调整幅度为真实的幅度值，即得到其频率响应。
[0071]
s208，根据所述频率响应及预设参数定义公式，计算绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心。
[0072]
本发明参考上述步骤s205的频段划分方式，形成5个由频率响应二次计算获取的音质感客观参数(绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心)。
[0073]
下面分别结合具体的计算方式分别对绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心进行说明。
[0074]
(1)绝对平均偏差aad
[0075]
绝对平均偏差用于表示各倍频程频带与参考频带的绝对幅度差；其中，绝对平均偏差的值越高，则表示幅度与参考频带的偏差越大。具体地，可根据下述公式，计算绝对平均偏差aad，
[0076][0077]
其中，n为300hz-16khz之间的1/20倍频程带的总数，y
band_
n为第n个倍频程带对应的平均幅度，y
ref@1-3khz
为1khz-3khz的参考频带的平均幅度。
[0078]
(2)低频延展参数lfx
[0079]
根据下述公式，计算低频延展参数lfx，
[0080][0081]
其中，y
ref
为参考频带内的平均幅度再加上10db(即，y
ref
＝参考频带内的平均幅度 10db)，x-6db
为在300hz以下并低于y
ref
值6db的第一个频率。
[0082]
需要说明的是，由于车内音响均衡器对低频的提升，y
ref
需在参考频带内的平均幅度的基础上 10db。其中，参考频带一般取为300hz-10khz，而对于车内音响系统，本发明采用的参考频带取为300hz
–
3khz。
[0083]
(3)高频延展参数hfx
[0084]
根据下述公式，计算高频延展参数hfx；
[0085][0086]
其中，y
ref
为参考频带内的平均幅度，x-6db
为在5khz以上并低于y
ref
值6db的第一个频率，本发明采用的参考频带取为300hz
–
3khz。
[0087]
(4)频谱平坦度spf
[0088]
频谱平坦度即信号功率谱的几何平均值与算术平均值的比值；该比值在0到1之
间，对于一个白噪声其计算结果为1，对于一个纯音信号，计算结果为0；该值越接近1，表明平坦度越高；可通过取对数扩大取值范围，取对数后结果越大(越接近于0)，平坦度越高。
[0089]
具体地，可根据下述公式，计算频谱平坦度spf；
[0090][0091]
其中，sk为信号幅度谱在频率点k上的值，k为频带计算范围的上限频率，本发明采用的频带计算范围选择300hz
–
20khz的中高频。
[0092]
(5)频谱重心spc
[0093]
根据下述公式，计算频谱重心spc；
[0094][0095]
spc＝log
10
(spectral_centroid)，
[0096]
其中，sk(tm)为在时间tm的信号加窗后的幅度谱在频率点k上的值，fk为k对应的频率值，k为频带计算范围的上限频率点，spectral_centroid为每个时间tm计算结果的平均。
[0097]
s209，根据预设权重系数对所述音质感客观参数进行加权处理，生成并输出评价结果。
[0098]
具体地，评价结果用s表示，a、b、c、d、e、f、g、h、i分别为权重系数，则s＝a b
×
aad c
×
lfx d
×
hfx e
×
spf f
×
spc g
×
thd_l h
×
thd_m i
×
thd_h。
[0099]
因此，本发明采用8个与音质感相关的客观参数，包含5个由频率响应二次计算获取的音质感客观参数(绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心)，3个从非线性谐波失真中提取的音质感客观参数(低频谐波失真参数、中频谐波失真参数及高频谐波失真参数)。通过在车内测量并计算上述音质感客观参数，从而能够通过客观测量进一步预估主观评价结果。
[0100]
下面以4辆某型号车辆的车内音响系统为例，对本发明做进一步的详细说明。
[0101]
步骤一、获取车内录制信号及播放源信号；
[0102]
步骤二、将所述车内录制信号与播放源信号进行对齐处理。
[0103]
步骤三、根据对齐后的车内录制信号与播放源信号提取阶跃扫频信号。
[0104]
步骤四、将所述阶跃扫频信号的每一段单频信号，使用平顶窗进行加窗并进行fft变换，计算频率及谐波的幅值。
[0105]
步骤五、根据所述幅值及预设参数定义公式，计算低频谐波失真参数、中频谐波失真参数及高频谐波失真参数。
[0106]
步骤六、根据对齐后的车内录制信号与播放源信号提取白噪声信号；
[0107]
步骤七、从所述白噪声信号中截取目标信号，使用汉宁窗进行加窗并进行短时fft变换，计算频率响应；
[0108]
步骤八、根据所述频率响应及预设参数定义公式，计算绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心。
[0109]
具体的音质感客观参数如表1所示：
[0110]
表1
[0111]
车型车型1车型2车型3车型4绝对平均偏差5.846.597.125.3低频延展参数1.041.231.461.58高频延展参数4.214.213.953.78频谱平坦度-2.2-0.63-1.73-1.32频谱重心3.9933.9863.9833.987低频谐波失真参数0.751.040.533.49中频谐波失真参数0.110.30.390.43高频谐波失真参数0.030.020.090.05
[0112]
步骤九、根据预设权重系数对所述音质感客观参数进行加权处理，生成并输出评价结果，具体的权重系数如表2所示：
[0113]
表2
[0114]
权重系数取值a3b0.5c-3d1e1f1g-1h-1i-1
[0115]
计算出评价结果如表3所示：
[0116]
表3
[0117]
车型车型1车型2车型3车型4s7.9138.8116.3733.387
[0118]
因此，本发明中音质感客观参数，在不同车上均体现出了明显差异，从而可以得到评价结果为：车型2》车型1》车型3》车型4。
[0119]
参见图3，图3显示了本发明车内音响音质的客观评价系统100的具体结构，其包括：
[0120]
采集模块1，用于获取车内录制信号及播放源信号。其中，播放源信号是指输入车辆音响设置中的音频文件；车内录制信号是指通过车内音响播放后，再录制的音频文件。
[0121]
计算模块2，用于根据所述车内录制信号及播放源信号计算音质感客观参数。所述音质感客观参数包括低频谐波失真参数、中频谐波失真参数、高频谐波失真参数、绝对平均
偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心。其中，低频谐波失真参数是指在低频频段下的平均谐波失真；中频谐波失真参数是指在中频频段下的平均谐波失真；高频谐波失真参数是指在高频频段下的平均谐波失真；绝对平均偏差用于表示各倍频程频带与参考频带的绝对幅度差，绝对平均偏差的值越高，则表示幅度与参考频带的偏差越大；频谱平坦度即信号功率谱的几何平均值与算术平均值的比值，该值越接近1，表明平坦度越高。
[0122]
评价结果输出模块3，用于根据预设权重系数对所述音质感客观参数进行加权处理，生成并输出评价结果。通过预设权重系数与上述音质感客观参数构建模型，即可得到评价结果。需要说明的是，所述权重系数可根据实际情况进行预设，灵活性强。
[0123]
例如，基础权重系数为a，绝对平均偏差对应的权重系数为b，低频延展参数对应的权重系数为c，高频延展参数对应的权重系数为d，频谱平坦度对应的权重系数为e，频谱重心对应的权重系数为f，低频谐波失真参数对应的权重系数为g，中频谐波失真参数对应的权重系数为h、高频谐波失真参数对应的权重系数为i，则评价结果s＝a b
×
绝对平均偏差 c
×
低频延展参数 d
×
高频延展参数 e
×
频谱平坦度 f
×
频谱重心 g
×
低频谐波失真参数 h
×
中频谐波失真参数 i
×
高频谐波失真参数。
[0124]
因此，本发明创新地将音质感从基础测量中分离进行评价，在车内测量并计算音质感客观参数，从而能够通过客观测量进一步预估主观评价结果，可快速有效地评价整车音响音质感的好坏。
[0125]
如图4所示，所述计算模块2包括：
[0126]
对齐单元21，用于将所述车内录制信号与播放源信号进行对齐处理。具体地，将车内录制信号与播放源信号进行对比分析，并根据播放源信号进行信号对齐。
[0127]
第一提取单元22，用于根据对齐后的车内录制信号与播放源信号提取阶跃扫频信号。
[0128]
第一变换单元23，用于将所述阶跃扫频信号的每一段单频信号，使用平顶窗进行加窗并进行fft变换，计算频率及谐波的幅值。具体地，第一变换单元23将对齐后的阶跃扫频信号的每一段单频信号，使用时长1s(离散采样点48000)的平顶窗进行加窗，计算fft变换，计算该频率下的幅值，以及各个谐波的幅值。其中，fft的点数与平顶窗的窗长相等。
[0129]
第一计算单元24，用于根据所述幅值及预设参数定义公式，计算低频谐波失真参数、中频谐波失真参数及高频谐波失真参数。
[0130]
具体地，第一计算单元24执行以下计算：
[0131]
根据下述公式，计算低频谐波失真参数thd_l、中频谐波失真参数thd_m及高频谐波失真参数thd_h；
[0132][0133]
[0134][0135]
其中，n为相应频段测量的频率点的总数，fn为第n个频率点对应的频率值，thd(fn)为fn下的谐波失真。
[0136]
第二提取单元25，用于根据对齐后的车内录制信号与播放源信号提取白噪声信号。
[0137]
第二变换单元26，用于从所述白噪声信号中截取目标信号，使用汉宁窗进行加窗并进行短时fft变换，计算频率响应。
[0138]
第二计算单元27，用于根据所述频率响应及预设参数定义公式，计算绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心。
[0139]
具体地，第二计算单元27执行以下计算：
[0140]
(1)绝对平均偏差aad
[0141]
根据下述公式，计算绝对平均偏差aad，
[0142][0143]
其中，n为300hz-16khz之间的1/20倍频程带的总数，y
band_n
为第n个倍频程带对应的平均幅度，y
ref@1-3khz
为1khz-3khz的参考频带的平均幅度。
[0144]
(2)低频延展参数lfx
[0145]
根据下述公式，计算低频延展参数lfx，
[0146][0147]
其中，y
ref
为参考频带内的平均幅度再加上10db(即，y
ref
＝参考频带内的平均幅度 10db)，x-6db
为在300hz以下并低于y
ref
值6db的第一个频率，参考频带一般取为300hz-10khz，而对于车内音响系统，本发明采用的参考频带取为300hz
–
3khz。
[0148]
(3)高频延展参数hfx
[0149]
根据下述公式，计算高频延展参数hfx；
[0150][0151]
其中，y
ref
为参考频带内的平均幅度，x-6db
为在5khz以上并低于y
ref
值6db的第一个频率，本发明采用的参考频带取为300hz
–
3khz。
[0152]
(4)频谱平坦度spf
[0153]
根据下述公式，计算频谱平坦度spf；
[0154][0155]
其中，sk为信号幅度谱在频率点k上的值，k为频带计算范围的上限频率，本发明采
用的频带计算范围选择300hz
–
20khz的中高频。
[0156]
(5)频谱重心spc
[0157]
根据下述公式，计算频谱重心spc；
[0158][0159]
spc＝log
10
(spectral_centroid)，
[0160]
其中，sk(tm)为在时间tm的信号加窗后的幅度谱在频率点k上的值，fk为k对应的频率值，k为频带计算范围的上限频率点，spectral_centroid为每个时间tm计算结果的平均。
[0161]
因此，本发明结合一般车内音响系统分频方式，将整车频率响应曲线分为3个频段，分别为低频(50-300hz)，中频(300-3000hz)，高频(3000-20000hz)，从而在非线性的谐波失真中提取出低频谐波失真参数、中频谐波失真参数及高频谐波失真参数。同时，本发明参考上述频段划分方式，形成5个由频率响应二次计算获取的音质感客观参数(绝对平均偏差、低频延展参数、高频延展参数、频谱平坦度及频谱重心)。因此，通过在车内测量并计算上述音质感客观参数，从而能够通过客观测量进一步预估主观评价结果。
[0162]
相应地，本发明还提供了一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现上述车内音响音质的客观评价方法的步骤。同时，本发明还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述车内音响音质的客观评价方法的步骤。
[0163]
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。