标准规范下载简介

GY／T 3*9-2021 感知音频质量的客观测量方法.pdf

GY／T 3*9-2021 感知音频质量的客观测量方法.pdf简介：

"GY/T 3*9-2021 感知音频质量的客观测量方法.pdf" 是一份由中国广播电视行业标准制定的文件，全称为《广播电视音频质量客观评价方法》，发布于2021年。这份指南主要针对音频质量的客观评估，提供了一套系统和科学的方法，用于量化和测量音频信号的主观感知质量。

该标准可能包括音频的参数分析（如信噪比、失真度、频率响应等）、主观评价技术（如ABX测试、MOS/MOS-LQO等）以及音频编码和传输质量的测量。它的目的是为了帮助研究者、工程师和标准制定者在音频技术的发展和应用中，更准确地评价和优化音频信号的质量，确保广播、电视和在线音频内容的主观听感质量。

这份文件对于音频处理、通信工程、广播与电视技术、音频内容制作等领域具有重要的参考价值。

GY／T 3*9-2021 感知音频质量的客观测量方法.pdf部分内容预览：

GY/T 3*92021

10.2.2.10添加内部噪声

《企业质量文化建设指南 GB/T32230-2015》10.2.2.10添加内部噪

EThres[k] = 100.*x0.3**(00)

IThres[k]l = 100*x0.3**(10)

该处理阶段的模式E2[k,n]见公式（*0），将用于调制模式的计算，并被称为“未被抹除的激励

10.2.2.11时域分布前向掩蔽

E2[k,n] = E,[k,n] + EThres[k, n

为了模拟前向掩蔽，每个滤波器声道中的能量随时间用一阶低通滤波器进行处理抹除。时间常数 每个滤波器的中心频率（见表*）并按照公式（*1）计算。

式中 min 取值为0.00*s。 一阶低通滤波器按照公式（*2）进行计算

式中： T100 取值为0.020s； Tmin 取值为0.00*s。 一阶低通滤波器按照公

该处理阶段的模式E[k，n]被称为激励模式。

10.3激励模式的预处理

10.3激励模式的预处理

10.3.1激励模式的预处理概述

本条中的大多数算法均适用 。由于这两种耳朵的二次米样 因子和频带数量不同，与该因子相关的常数可采用与耳朵相关的变量步长（StepSize）和z进行推 述。对于FFT耳朵，步长值取102*，z的值可以取55（高级模式）或者109（基础模式）。对于滤波 器组耳朵，步长值取192，z的值取*0。若无其他说明，所有变量和递归滤波器的初始值均为0。

10.3.2电平和模式调整

10. 3. 2. 1电平调整和模式调整概述

为了补偿测试信号和参考信号之间的电平差异和线性失真，测试信号和参考信号的电平均调整到两 个信号的平均电平。 第一步，每个滤波器声道中的能量采用一阶低通滤波器进行平滑处理。时间常数取决于滤波器组的 中心频率，计算见公式（**）

式中： T100——取值为0.050s; min 一取值为0.008s。 一阶低通滤波器的计算见公式（*5）和公式（**）

式中： Tmin 取值为0.008s。

GY/T 3*92021

GY/I3*9202

ETest和ERef是相互调整后的激励模式，而时间常数a可根据公式（*7）进行计算得到

10. 3. 2. 2电平调整

如果校正系数大于1，参考信号应除以校正系数，否则测试信号应乘以该校正系数，见公式 公式（50）。

10.3.2.3模式调整

EL,Ref[k,n] = Eref[k,n]/LevCorr[n] |LevCorr[n] > 1 · EL.Test[k, n] = Erest[k,n] . LevCorr[n] LevCorr[n] ≤ ....

RTest[k,n] = RRer[k,n] = 1R[k,n] ≥ 1 R[k,n] RTesr[k, n] = 1, Rrer [k, n] = R[k, n] I R[k,n] < 1 +

1 .***.*. (52)

如果公式（*8）的分母为0（则Rk，n]将不存在）且分子天于0，设置RTestLk，n)为0并设置RRef[k，n』 为1。当公式（*8）中的分子为0，RTest［k，n]和RRef［k，n]的比率可从后续频带计算得到。如果后续频 带不存在（如k=0），则RTest［k，n]和RRef［k，n]比率设置为l。 校正系数在M个滤波器声道上进行平均，并通过公式（**）和公式（*7）中给出的相同时间常数随 时间进行平滑处理，见公式（53）。对滤波器组耳朵，频率窗口M的宽度是3，对FFT耳朵，频 率窗口M的宽度是*（高级模式）或8（基础模式）

GY/T 3*92021

GY/T 3*92021

在频率标度的边界，即频率窗口将超过滤波器组范围的地方，频率窗口的宽度相应降低，见公式（5*）

为了获得频谱自适应模式，电平自适应输入模式采用对应的校正系数进行加权，见公式（ 式（5*）。

Ep,Ref[k,n] = EL,Ref[k,n] · PattCorrRef[k,n] Ep,Test[k,n] = EL.Test[k,n] : PattCorTest[k,n

通过对未抹除的激励模式E2[k，n]取0.3次幂，并通过公式（57）和公式（58），计算得到简化 该值和其时域得到的绝对值随时间进行抹除

根据Ede和E按照公式（*0）计算每个滤波器输出包络的调制测量值

值E在下文中还将用于调制差异的计算！

被测信号和参考信号的指定响度模式通过公式（*1）进行计算

100 Hz T100 0.050 S = To + To 0.008s

Ntotai[n] = · = max(N[k,n],0)

为计算1kHz*0dBsPL正弦波1宋时的整体响度，FFT耳朵的比例常数const为1.07***，滤波器组 耳朵的比例常数const为1.2*539。阈值s和阅值点的激励EThres分别根据公式（*3）和公式（**） 进行计算。

[kl = 100.3**(1)

Erha[kl = 1003**(100)

GY/T 3*92021

(22.05atn(0)0.75atn()

注：由于存在不同的耳朵周边，此处所指的响度与IS0532（声学－计算响度电平的方法1975)中定

10.3.5误差信号的计算

误差信号只在FFT耳朵模式中计算，在频域里计算经外耳和中耳过滤后的参考信号和测试信号 频谱之间的差异得到该值，见公式（*5），其中参考信号和测试信号的外耳加权FFT输出见10.2

采用10.2.1.5描述的算法，将Fnoise映射到音高域。 该算法的输出Pnoise[n,kl，被称为噪声模式。

10.*输出变量的计算

用于预测基本音频质量的MOV值应与表5相符合

表5用于预测基本音频质量的MOV值

10. *. 2调制差异

10. *. 2调制差异

10. *. 2. 1概述

GY/T 3*92021

被测信号和参考信号的时域包络调制间的差异，通过计算每个滤波器声道的本地调制差异后再进行

ModDiff[n] = 10 Z= ModDiff[k,n]

静音阈值在计算时考虑了电平相关的加权系数，见公式（*8）。该电平相关的加权系数是根据公式 （55）给出的参考信号的修正激励模式以及内部噪声函数计算得到，滤波器组耳朵的内部噪声定义 见公式（39），FFT耳朵的内部噪声定义见公式（15）。

采用加权系数TempWt[n]的瞬时调制差异ModDiff[n]的时域平均值将在10.5.2中描述。常数negWt offset和levWt的值应与表*相符合。

表*估计总体调制差异的输出变量

10.*.2.2 RmsModDiff

10.*.2.2RmsModDiff.

输出变量RmsModDiff.调制差异的方均值。该差异从滤波器组耳朵中计算得到。瞬时平均 值参考10.5.2，常数应与表*相符合。 10.*.2.3WinModDiff1B 输出变量WinModDiff1s调制差异的窗口平均值。调制差异根据FFT耳朵进行计算得到。时 域平均法参考10.5.2.*，常数应与表*相符合。公式（*5）给出时间加权系数不能应用于该MOV

输出变量RmsModDiff.调制差异的方均值。该差异从滤波器组耳朵中计算得到。 值参考10.5.2，常数应与表*相符合

10.*. 2. 3WinModDiff1

出变量WinModDiff1s调制差异的窗口平均值。调制差异根据FFT耳朵进行计算得到。时 考10.5.2.*，常数应与表*相符合。公式（*5）给出时间加权系数不能应用于该MOV

10.*.2.*AvgModDiff1.和AvgModDiff2

莫型输出变量AvgModDiff1:和AvgModDiff2:是依据FFT耳朵计算得到的调制差异的线性平均值 dDiff1:和AvgModDiff2:之间的差异是因为采用了不同常数。时间平均法见10.5.2.2，常数应与表

10. *. 3噪声响度

10.*. 3. 1概述

GY/T 3*92021

噪声响度的计算公式（见公式（*9））用于计算噪声的指定响度；如果与掩蔽信号相比噪声非常小的情 况下，噪声响度的计算公式用于计算噪声掩蔽比

其中Eo恒等于1，Erres指内部噪音函数ETres［k]，定义见公式（3*）土木工程施工-第17讲-结构安装工程，s可通过公式（70）进行 s = ThresFaco × Mod[k,n] + So

瞬时噪声响度值需测试信号和参考信号的左右任一声道的噪声响度超过NThres=0.1sone，50ms后开 始计算，见10.5.2.5.2。 在谱平均中，瞬时值由每个临界频带滤波器组的数量进行归一化，而不是由滤波器组的总数进行归 化，也就是，谱平均的结果乘以值为2*的系数。 如果瞬时噪声响度低于阈值NLain，设置其为0。估算整体噪声响度的MOV值应与表7相符合

表7估算整体噪声响度的MOV

10. *. 3.2 RmsNoiseLoud.

0. *. 3.2 RmsNoiseLou

NoiseLoud.指滤波器组耳朵计算得到的噪声响度方均值。时间平均法见10.5.2.3JT∕T 975-2015 沥青比重瓶，常数见

10.*. 3.3RmsMissingComponents.

10.*.3.3RmsMissingComponents