本发明涉及听障儿童的汉语语音研究领域。
背景技术：
对于听障儿童的发音问题一直是病理语音方向重要的研究内容，由于受到听障患者语音语料难采集的限制、还有一些语音分析仪器价格昂贵难普及现状，以及目前听障患者的发音评估主要依靠主观专业医师听辨，而专业医师对于数目庞大的听障儿童来说远远不能满足等情况。
作者何维佳等人在文献《成人痉挛型构音障碍普通话元音声学分析》中，观察痉挛型构音障碍患者普通话元音/a/、/i/、/u/的声学特征，通过多维语音分析系统mdvp提取出病例组与对照组的元音声学参数，包括元音基频、基频微扰值、振幅微扰值、相对平均微扰，音高周期扰动商，平滑音高周期扰动商、振幅微扰商、平滑振幅微扰商、谐噪比等九项声学参数，之后利用spss统计学软件对两组元音之间的声学特征参数进行方差分析，结果讨论了病例组与病例组声学参数差异的原因，并得出汉语元音/a/可以成为病例的检测声样结论。该实验的研究对象是痉挛型构音障碍患者，且只考虑了基本的声学特征参数，没有发音时运动学方面的参数，评价指标方面不够全面综合，且在结论中，仅将/a/作为病例的检测声样，没有继续评估其它元音的发音质量。
作者万勤等人在文献《3～6岁听障儿童与键听儿童嗓音声学特征比较》中，随机选取3～6岁听障儿童与键听儿童各100名进行对照实验，分别提取两者发音
技术实现要素：
本发明所要解决的技术问题是：如何解决背景技术中不能够全面评价听障儿童的发音的问题。
本发明所采用的技术方案是：基于模糊综合评价的听障儿童元音评估方法，按照如下步骤进行
步骤一、使用三维电磁发音记录仪获得多个正常人和听障儿童阅读元音语料时的基本声学参数u1、舌位的声学参数u2、唇舌运动参数u3建立模糊综合评价数据库，基本声学参数u1包括特征参数发音时长u11、基频u12、短时振幅u13，舌位的声学参数u2包括特征参数u21即共振峰f1、u22即共振峰f2，唇舌运动参数u3包括特征参数上唇z方向上的位移u31、下唇z方向上的位移u32、舌尖z方向上的位移u33、舌中z方向上的位移u34，以竖直方向为z方向并且从下向上递增；
步骤二、建立二级评价指标u1＝{u11、u12、u13}、u2＝{u21、u22}、u3＝{u31、u32、u33、u34}，将二级评价指标的每个特征参数根据发音最异样到发音标准建立隶属函数，给隶属函数设置六个得分等级，根据每个元音语料的六个得分等级分数建立二级评价指标判断矩阵即基本声学参数判断矩阵r1、舌位的声学参数判断矩阵r2、唇舌运动参数判断矩阵r3，二级评价指标得分为，基本声学参数评价指标得分b1＝r1*a1，舌位的声学参数评价指标得分b2＝r2*a2，唇舌运动参数评价指标得分b3＝r3*a3，其中，a1为基本声学参数权重，a2为舌位的声学参数权重，a3为唇舌运动参数权重；
步骤三、建立一级评价指标为u＝{u1、u2、u3}，一级评价指标判断矩阵r为{基本声学参数评价指标得分b1、舌位的声学参数评价指标得分b2、舌位的声学参数评价指标得分b2}，一级评价指标得分b＝r*a，a为一级评价指标权重，用一级评价指标得分评价听障儿童阅读元音语料时的发音情况。
作为一种优选方式：元音语料为单元音/a/，/i/，/u/，/o/，/e/，/ü/。
作为一种优选方式：基本声学参数权重a1＝{0.1876、0.1812、0.631}，舌位的声学参数权重a2＝{0.605、0.395}、唇舌运动参数权重a3＝{0.3906、0.1518、0.2258、0.2317}，一级评价指标权重a＝{0.33825、0.34771、0.34767}。
作为一种优选方式：隶属函数设置的六个得分等级{0.9、0.8、0.6、0.5、0.4、0.2}。
本发明的有益效果是：将基本声学参数、舌位声学参数、唇舌运动学参数共同作为评价普通话单元音发音指标因素，更加全面地考虑了发音的评价因素，且分别通过熵权法，标准差法客观地计算参数指标所占权重，并结合模糊综合评价法计算出得出了听障儿童汉语单元音的得分情况，为评估听障儿童汉语单元音发音状况提供了一种客观的评估方法，通过听障儿童的发音数据，提供了一种客观的评价听障儿童汉语单元音方法。
具体实施方式
元音语料为普通话单元音/a/，/i/，/u/，/o/，/e/，/ü/，实验被试人员为9名听障儿童，年龄为14～20岁，除听力受损外无其他障碍，可以自主发音。正常键听学生9名，均为普通话标准的在校大学生，让两组被试人员在自然放松的状态下，依次阅读测试语料，选取最好的一次发音作为实验数据，本实验采集数据的仪器型号为ag501，以250帧每秒的采样率录制发音动作，在测试者产生语音的同时并记录各个器官的运动数据，同时记录下音频数据，实验中采集的运动数据所涉及的坐标为三维立体坐标，本发明选取z轴方向上的唇舌运动数据，三维坐标以每个被试人员上下方向为z轴并且方向是从下向上递增(以竖直方向为z方向并且从下向上递增)。
首先对语音信号进行预处理，包括端点检测、加窗、分帧，再提取听障儿童和正常人语音各项特征参数，包括发音时长、基频、短时振幅平均值、共振峰f1和f2、上唇z方向上的位移、下唇z方向上的位移、舌尖z方向上的位移、以及舌中z方向上的位移。发音时长、共振峰f1和f2通过语音分析软件praat提取出，基频和短时振幅平均值通过matlabr2010b提取出，其中基频采用平均幅度差法提取，取其平均值，上唇z方向上的位移、下唇z方向上的位移、舌尖z方向上的位移、以及舌中z方向上的位移通过运动学数据可视化软件visartico提取，所有唇舌运动数据为与发音时长相对应的运动数据，取其平均值。通过以上方法可以获得听障儿童与正常人的各项特征参数，将正常人的特征参数数据作为标准，正常人数据如下表一所示：
表一各项参数参考标准值
建立二级模糊综合评价模型，本实施例中，听障儿童普通话单元音/a/，/i/，/u/，/o/，/e/，/ü/的一级评价指标因素有三个，u＝{u1、u2、u3},u1基本声学参数、u2舌位声学参数、u3唇舌运动参数，二级评价指标因素有九个，u1＝{u11、u12、u13}、u2＝{u21、u22}、u3＝{u31、u32、u33、u34}、u11为发音时长，u12为基频、u13为短时振幅平均值，u21为共振峰f1、u22为共振峰f2,u31为上唇z方向上的位移、u32为下唇z方向上的位移、u33为舌尖z方向上的位移、u34为舌中z方向上的位移。在matlabr2010b环境中，首先对所有的特征参数进行归一化处理，去除指标因素的量纲影响，再通过熵权法计算出各个二级评价指标因素所占的权重，经计算后：u1＝{u11、u12、u13}所占的权重为a1＝{0.1876、0.1812、0.631}、u2＝{u21、u22}所占的权重为a2＝{0.605、0.395}、u3＝{u31、u32、u33、u34}所占的权重为a3＝{0.3906、0.1518、0.2258、0.2317}，建立评判矩阵，对每一个指标因素进行评价，本实施例将听障儿童的特征参数与正常人标准的特征参数之间距离作为评判标准，建立隶属度函数，函数根据距离的大小设置为六个分数{0.90.80.60.50.40.2}，所有的因素指标均依照距离建立隶属度函数，通过计算听障儿童的特征参数与正常人标准的特征参数距离后，所有二级评价指标因素的分数依次为r1＝{0.60.40.80.50.90.2}、r2＝{0.50.20.80.60.90.4}、r3＝{0.60.20.40.50.80.9}、r4＝{0.80.40.60.50.20.9}、r5＝{0.90.20.50.60.80.4}、r6＝{0.50.40.80.90.20.6}、r7＝{0.20.50.60.40.80.9}、r8＝{0.20.80.90.60.50.4}、r9＝{0.60.20.90.80.50.4}，故对u1＝{u11u12u13}建立判断矩阵为
，将二级模糊综合评价的b1、b2和b3作为一级评价指标u＝{u1、u2、u3}的评价矩阵，记为：
最后通过标准差法，将听障儿童的普通话单元音特征参数到正常人的标准元音特征参数的标准差作为一级评价指标u＝{u1、u2、u3}的权重，所有数据需先进行归一化处理，经计算后得到得到u＝{u1、u2、u3}的权重a＝{0.338250.347710.34767}，为兼顾到所有因素指标的作用，一级模糊综合评判，
在将每个元音的得分除以元音总得分，即求出了每个元音的得分百分比，故最终对听障儿童的汉语单元音/a/，/i/，/u/，/o/，/e/，/ü/的评价结果为b(vowel)＝[0.179710.100020.189030.175970.165940.18933]，从结果可以看出，本发明评价听障儿童的汉语单元音评估结果为听障儿童的发音得分为/ü/高于/u/高于/a/高于/o/高于/e/高于/i/，得分越高，表明此发音越好，在对听障儿童的汉语单元音进行评估时，可以通过模糊综合评价法，判断听障儿童的汉语元音发音情况。