Моделирование особенностей речи диктатора

Для создания системы верификации диктора в статье была предложена методика верификации на основе шипящих звуков, которая использует авторские методы, осуществляющие обобщенную классификацию звуков посредством дискретного и непрерывного вейвлет-преобразования. Для шипящих звуков, выделенн...

Full description

Saved in:
Bibliographic Details
Date:2008
Main Authors: Криводубский, О.А., Федоров, Е.Е.
Format: Article
Language:Russian
Published: Інститут проблем математичних машин і систем НАН України 2008
Subjects:
Online Access:https://nasplib.isofts.kiev.ua/handle/123456789/755
Tags: Add Tag
No Tags, Be the first to tag this record!
Journal Title:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Cite this:Моделирование особенностей речи диктатора / О.А. Криводубский, Е.Е. Федоров // Математические машины и системы. — 2008. — № 1. — С. 163-170.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
_version_ 1860198706554339328
author Криводубский, О.А.
Федоров, Е.Е.
author_facet Криводубский, О.А.
Федоров, Е.Е.
citation_txt Моделирование особенностей речи диктатора / О.А. Криводубский, Е.Е. Федоров // Математические машины и системы. — 2008. — № 1. — С. 163-170.
collection DSpace DC
description Для создания системы верификации диктора в статье была предложена методика верификации на основе шипящих звуков, которая использует авторские методы, осуществляющие обобщенную классификацию звуков посредством дискретного и непрерывного вейвлет-преобразования. Для шипящих звуков, выделенных этими методами, проведен количественный анализ систем признаков, основанных на линейном предсказании, нормированном количестве импульсов равной длины и мел-частотных кепстральных коэффициентах (MFCC). Полученные признаки используются в методе верификации, основанном на алгоритме DTW. Табл.: 1. Библиогр.: 10 назв. Для створення системи верифікації диктора у статті була запропонована методика верифікації на підставі шиплячих звуків, що використовує авторські методи, які здійснюють узагальнену класифікацію звуків за допомогою дискретного та безперервного вейвлет-перетворення. Для шиплячих звуків, виділених цими методами, проведений кількісний аналіз систем ознак, заснованих на лінійному пророкуванні, нормованої кількості імпульсів рівної довжини та мел-частотних кепстральних коефіцієнтах (MFCC). Отримані ознаки використовуються в методі верифікації, заснованому на алгоритмі DTW. Табл.: 1. Бібліогр.: 10 назв. For creation of a system of verification of the speaker in the article the technique of verification on the basis of unvoiced fricative sounds was offered which uses author's methods executing a generalized classification of sounds by means of discrete and continuous wavelet-transformation. For unvoiced fricative sounds chosen by these methods, the quantitative analysis of systems of indications based on a linear prediction, normalized quantity of pulses of equal length and mel-frequency cepstral coefficients (MFCC) is conducted. The obtained indications are used in a method of verification based on the algorithm DTW. Tabl.: 1. Refs.: 10 titles.
first_indexed 2025-12-07T18:09:35Z
format Article
fulltext ISSN 1028-9763. Математичні машини і системи, 2008, № 1 163 УДК 004.934.1’1 О.А. КРИВОДУБСКИЙ, Е.Е. ФЕДОРОВ МОДЕЛИРОВАНИЕ ОСОБЕННОСТЕЙ РЕЧИ ДИКТОРА Abstract: For creation of a system of verification of the speaker in the article the technique of verification on the basis of unvoiced fricative sounds was offered which uses author's methods executing a generalized classification of sounds by means of discrete and continuous wavelet-transformation. For unvoiced fricative sounds chosen by these methods, the quantitative analysis of systems of indications based on a linear prediction, normalized quantity of pulses of equal length and mel-frequency cepstral coefficients (MFCC) is conducted. The obtained indications are used in a method of verification based on the algorithm DTW. Key words: verification of the speaker, generalized classification of sounds, wavelet-transformation, linear prediction, normalized quantity of pulses of equal length, mel-frequency cepstral coefficients, algorithm DTW. Анотація: Для створення системи верифікації диктора у статті була запропонована методика верифікації на підставі шиплячих звуків, що використовує авторські методи, які здійснюють узагальнену класифікацію звуків за допомогою дискретного та безперервного вейвлет-перетворення. Для шиплячих звуків, виділених цими методами, проведений кількісний аналіз систем ознак, заснованих на лінійному пророкуванні, нормованої кількості імпульсів рівної довжини та мел-частотних кепстральних коефіцієнтах (MFCC). Отримані ознаки використовуються в методі верифікації, заснованому на алгоритмі DTW. Ключові слова: верифікація диктора, узагальнена класифікація звуків, вейвлет-перетворення, лінійне пророкування, нормована кількість імпульсів рівної довжини, мел-частотні кепстральні коефіцієнти, алгоритм DTW. Аннотация: Для создания системы верификации диктора в статье была предложена методика верификации на основе шипящих звуков, которая использует авторские методы, осуществляющие обобщенную классификацию звуков посредством дискретного и непрерывного вейвлет-преобразования. Для шипящих звуков, выделенных этими методами, проведен количественный анализ систем признаков, основанных на линейном предсказании, нормированном количестве импульсов равной длины и мел- частотных кепстральных коэффициентах (MFCC). Полученные признаки используются в методе верификации, основанном на алгоритме DTW. Ключевые слова: верификация диктора, обобщенная классификация звуков, вейвлет-преобразование, линейное предсказание, нормированное количество импульсов равной длины, мел-частотные кепстральные коэффициенты, алгоритм DTW. 1. Введение Постановка проблемы. В настоящее время актуальной является разработка систем, предназначенных для верификации диктора. Эти системы имеют широкую область применения: криминалистика, криптография, охранные системы и др. При разработке подобных систем важную роль играет выбор системы признаков и методов верификации, использующих данные признаки. Анализ исследований. В работах [1–2] приведены системы верификации, дающие в большинстве случаев вероятность распознавания ниже 90%. Верификация обычно проводится на целых словах, получение которых не всегда возможно. Цель и задачи исследования. Для повышения надежности верификации необходимо разработать методику верификации диктора на основе шипящих звуков. 2. Методика верификации диктора и методы классификации звуков В статье рассматриваются: – методика верификации дикторов; – авторские методы ФЕОК-ДВП и ФЕОК-НВП, осуществляющие обобщенную классификацию звуков посредством дискретного и непрерывного вейвлет-преобразования; – проведенный количественный анализ для выбора системы признаков шипящих звуков. Методика верификации диктора на основе шипящих звуков включает три этапа. ISSN 1028-9763. Математичні машини і системи, 2008, № 1 164 На первом этапе с целью выделения шипящих звуков осуществляется обобщенная классификация звуков посредством авторских методов ФЕОК-ДВП и ФЕОК-НВП. Авторский метод ФЕОК-ДВП Метод ФЕОК-ДВП осуществляет обобщенную классификацию звуков посредством дискретного вейвлет-преобразования и заключается в следующем. Производится декомпозиция сигнала )(nx на P уровней с вычислением высоко- ( imd ) и низкочастотных ( imc ) составляющих (субполос) посредством свертки на текущем i -том уровне ( Pi ,1∈ ) сигнала с полосовыми фильтрами с коэффициентами ng , nh [3–5]: ∑ − = +− − = 12/ 0 2,1 2/1 1 2 iN n mnniim gcd , ∑ − = +− − = 12/ 0 2,1 2/1 1 2 iN n mnniim hcc , где )(0 nsc n = , 12/,0 1 −∈ −iNm . Затем сигнал разбивается на фреймы длиной N∆ . Для s -го фрейма на i -м уровне разложения вычисляется энергия ∑ ∆ ∆−= = Ns Nsm imsi dDE )1( 2 , . )1( 2∑ ∆ ∆−= = Ns Nsm PmsP cCE Для s -го фрейма на i -м уровне разложения вычисляется мера контрастности . 1 ∑ = = i j sjsisi DEDEContr Тип каждого s -го фрейма определяется следующим образом:      ≤ ≥ << = ,, , , 2 1 21 тональныйContr шипящийContr шумContr n s s s γγ γγ γγγ α α αα где γα1 , γα2 – пороги, которые автоматически вычисляются в подсистеме адаптации и представляют собой минимальное и максимальное значения контрастности шума на уровне γ . Если 2=n , то определяются левая 1l и правая 2l границы шипящих звуков. Достоинством метода ФЕОК-ДВП является малая вычислительная сложность. Недостатком – сложность проведения классификации в зашумленном сигнале. Авторский метод ФЕОК-НВП Метод ФЕОК-НВП осуществляет обобщенную классификацию звуков посредством аппроксимированного непрерывного вейвлет-преобразования и заключается в следующем. Производится декомпозиция сигнала )(nx на P уровней с вычислением вейвлет- коэффициентов ( imd ) [4–5]. ∑ − = ∆= 1 0 )()( N n ilil tnnxd ψ , 1,0 −∈ Nl , ,,1 Pi ∈ ISSN 1028-9763. Математичні машини і системи, 2008, № 1 165 где t∆ – величина, обратная частоте дискретизации; ( )lbnaan ii il 00 2/ 0)( −= −− ψψ , – вейвлет, 10 >a , 0≠b . Затем сигнал разбивается на фреймы длиной N∆ . Для s -го фрейма на i -м уровне разложения вычисляется энергия . )1( 2∑ ∆ ∆−= = Ns Nsm imsi dDE Для s -го фрейма по всем уровням разложения вычисляется мера контрастности ∑ = = i j sjsisi DEDEContr 1 , .,1 Pi ∈ Тип каждого k -го фрейма определяется следующим образом:       = ∑ Θ∈γ γδ tk t n maxarg , },...,1{ P⊂Θ ,    << = иначе Contrk k ,0 ,1 21 1 γγγ γ ααδ ,    > = иначе Contrk k ,0 ,1 1 2 γγ γ αδ ,    > = , ,0 ,1 2 3 иначе Contrk k γγ γ αδ где γα1 , γα2 – пороги, которые автоматически вычисляются в подсистеме адаптации и представляют собой минимальное и максимальное значения контрастности шума на уровне γ , Θ – подмножество уровней, получаемое в результате численного исследования. . ,3 ,2 ,1      = тональный шипящий шум n Если 2=n , то определяются левая 1l и правая 2l границы шипящих звуков. Достоинством метода ФЕОК-НВП является возможность проведения классификации в зашумленном сигнале. Недостатком – большая вычислительная сложность. На втором этапе вычисляются значения признаков шипящих звуков. Сигнал )(ms предварительно разбивался на фреймы (участки равной длины). В качестве признаков использовались: 1. Нормированная автокорреляция [6–7], вычисляемая на n -м фрейме сигнала )(msn , )()()( mwmsms nn =) , где N n nw ∆ += π2 cos46,054,0)( – окно Хемминга, ISSN 1028-9763. Математичні машини і системи, 2008, № 1 166 ∑ −− = += iN m nnn imsmsiR 1 0 )()()( )) , pn ,1∈ , где p – порядок линейного предсказателя, )0( )( )( n n n R iR iR = , ))(,...,)1((1 pRRQ nnin = , )(,1 1Qη∈i , Ln ,1∈ , где )( 1Qη – количество эталонов шипящих фонем, L – количество фреймов. 2. Коэффициенты линейного предсказания, вычисляемые с помощью алгоритма Дарбина [6–7] на n -м фрейме сигнала: )0(:)0( nn RE = , piEjiRiRk i n i j n i njnni ≤≤      −−= − − = −∑ 1,)()(: )1( 1 1 )1(α , ni i ni k=:)(α , 11,: )1( , )1()( −≤≤−= − − − ijk i jinni i nj i nj ααα , )1(2)( )1(: −−= i nni i n EkE , pja p njnj ≤≤= 1,: )(α , где )(iRn – автокорреляционная функция, )(i njα – j -й коэффициент линейного предсказателя порядка i , nik – i -й коэффициент отражения, )(i nE - среднеквадратичная погрешность предсказания для линейного предсказателя порядка i . ),...,( 1 2 npnin aaQ = , )(,1 2Qη∈i , Ln ,1∈ . 3. Коэффициенты отражения КЛП (или PARCOR-коэффициенты) [6–7], определяемые по алгоритму Дарбина на n -м фрейме сигнала: ),...,( 0 3 npnin kkQ = , )(,1 3Qη∈i , Ln ,1∈ . 4. Нормированная автокорреляция КЛП [6–7], получаемая по коэффициентам линейного предсказания: ∑ = = p j njn ar 0 2)0( , ∑ − = += kp j kjnnjn aakr 0 ,2)( , )0( )( )( n n n r ir ir = , ))(,...,)1((4 prrQ nnin = , )(,1 4Qη∈i , Ln ,1∈ . 5. Кепстр импульсной характеристики системы линейного предсказания [6–7], вычисляемый по коэффициентам линейного предсказания: ISSN 1028-9763. Математичні машини і системи, 2008, № 1 167 1)0(ˆ 0 == nn ah , ∑ − = ⋅+= 1 1 ,)(ˆ)(ˆ j k j-knnnjn akh j k ajh , pn ,1∈ , ))(ˆ),...,1(ˆ(5 phhQ nnin = , )(,1 5Qη∈i , Ln ,1∈ . 6. Площади поперечных сечений кусочно-постоянной акустической трубы [6–7], содержащей )1( +p цилиндрическую секцию фиксированной длины, вычисляемые с помощью коэффициентов отражения: 11 =nA , ni ni ni in A k k A + −=+ 1 1 1, , 1,2 +∈ pi , ),...,( 1,2 6 += pnnin AAQ , )(,1 6Qη∈i , Ln ,1∈ . 7. Нормированный энергетический спектр КЛП [1]: ∑ ∑ = =      − − = p s nn p k nnkn n sk N srr kRR kW 1 1 2 cos)()0( )()0( )( π α , _________ 12/,0 −∈ Nk , ( ) ( )∑ − = = 12/ 0 )( N i n n n iW kW kW , 12/0 −≤≤ Nk , ))12/(,...,)0((7 −= NWWQ nnin , )(,1 7Qη∈i , Ln ,1∈ . 8. Нормированный энергетический спектр, вычисленный на основе энергетического спектра ДПФ [1]: )()()( mwmsms nn =) , где N n nw ∆ += π2 cos46,054,0)( – окно Хемминга. ( ) ( )∑ −∆ = ∆ − = 1 0 2 ˆ N n N nk j enskS π , 12/0 −∆≤≤ Nk , ( ) ( )∑ − = = 12/ 0 )( N i n n n iS kS kS , 12/0 −≤≤ Nk , ))12/(,...,)0((8 −= NSSQ nnin , )(,1 8Qη∈i , Ln ,1∈ . 9. Нормированное количество импульсов равной длины. Для каждого n -го фрейма вычисляется nzd – количество импульсов длины z [8], lenz ,1∈ , которое далее нормируется ∑ = = len s ns nz nz d d d 1 |||| , ISSN 1028-9763. Математичні машини і системи, 2008, № 1 168 ),...,( ,1 9 lennnin ddQ = , )(,1 9Qη∈i , Ln ,1∈ . 10. Мел-частотные кепстральные коэффициенты (MFCC) [9], вычисляемые с использованием обратного дискретного косинусного преобразования на n -м фрейме сигнала: ( )∑ = −= P l nlnk PlkEMFCC 1 /)5.0(cos π , ,,1 Pk ∈ где ( )         ∆+−= ∑ = i i k kk llnnl KkkwkSE 2 1 2 ))2/1(()(lg ) – логарифмированная энергия l -той мел- частотной полосы. Для ее вычисления использовалась треугольная оконная функция Бартлета:          ∆> ∆≤≤ ∆ − ≤≤∆− ∆ + ∆−< = 2/,0 2/0, 2 1 02/, 2 1 2/,0 )( l l l l l l Kk Kk K k kK K k Kk kw , где lll kkK 12 −=∆ – четное число, ll kk 2,1 – границы частотных диапазонов l -той полосы. ),...,( 1 10 nPnin MFCCMFCCQ = , )(,1 10Qη∈i , Ln ,1∈ . На основе каждой s -й системы признаков формируются эталоны шипящих звуков. )),...,,...,(),...,,...,,...,(),...,,...,,...,(( 111111 s iLK s iLk s iL s ijK s ijk s ij s Ki s ki s i s i QQQQQQQQQQ = , 10,1∈s , )(,1 si Qη∈ , где )( sQη – количество эталонов шипящих фонем для s -й системы признаков; L – количество участков равной длины (фреймов), на которые разбивается сигнал; K – количество признаков, описывающих один фрейм. На третьем этапе верификации диктора, используя полученные векторы признаков шипящих звуков, производится собственно верификация. При этом используется алгоритм DTW [10]. В основе алгоритма DTW лежит рекуррентная формула s i s i DC 1111 = , ),,min( 1,1,1,,,1, −−−−+= jmijmijmi s imj s imj CCCDC , LjLm ,1,,1 ∈∈ , где s imjC – расстояние между левыми частями фонемы (фреймы от 1 до m ) и i -м эталоном (фреймы от 1 до j ); s imjD – расстояние между m -м фреймом фонемы и j -м фреймом i -го эталона. В качестве s imjD выбрана евклидова метрика ∑ = −= K k s ijk s mk s imj QQD 1 2)( ) , ISSN 1028-9763. Математичні машини і системи, 2008, № 1 169 где s mkQ ) – k -й признак m -го фрейма фонемы; s ijkQ – k -й признак j -го фрейма i -го эталона; s – тип системы признаков. Количественная оценка сопоставления шипящего звука верифицируемого диктора sQ ) с эталонами s iQ для s -й системы признаков вычислялась в соответствии с функционалом ошибки minmin →=ℑ s iLL i C . Результат верификации диктора определялся в соответствии с правилом    ∈ ∈ = 2 1 ,"" ,"" Nnчужой Nnсвой result , s iLL i Cn minarg= . Для выбора системы признаков было проведено численное исследование, в котором участвовало 100 дикторов (50 «своих», 50 «чужих»): мужчины и женщины с разными голосовыми данными. В табл.1 приведены результаты верификации диктора по шипящему звуку («ш»). Численное исследование позволяет сделать вывод, что наиболее перспективными являются признаки MFCC и коэффициенты отражения КЛП. Таблица 1. Результаты численного исследования систем признаков Система признаков Вероятность верификации Нормированная автокорреляция 0,44 Коэффициенты КЛП 0,78 Коэффициенты отражения КЛП 0,96 Нормированная автокорреляция КЛП 0,78 Кепстр КЛП 0,68 Площади поперечных сечений акустической трубы КЛП 0,42 Нормированный энергетический спектр КЛП 0,4 Нормированный энергетический спектр ДПФ 0,32 Нормированное количество импульсов равной длины 0,56 MFCC 0,96 3. Выводы Новизна. В статье предложены методика верификации диктора на основе шипящих звуков, авторские методы ФЕОК-ДВП и ФЕОК-НВП, осуществляющие обобщенную классификацию звуков в дискретном сигнале с целью выделения шипящих. Было проведено численное исследование систем признаков шипящих звуков, при этом в качестве метода распознавания был выбран алгоритм DTW. Преимуществом предлагаемого подхода является то, что для верификации диктора вместо всего слова достаточно использовать входящий в него шипящий звук, что расширяет область применения системы верификации. Практическое значение. Основные положения работы были использованы при разработке системы верификации диктора, которая может использоваться в криминалистике и охранных системах. ISSN 1028-9763. Математичні машини і системи, 2008, № 1 170 СПИСОК ЛИТЕРАТУРЫ 1. Рабинер Л.Р., Шафер Р.В. Цифровая обработка речевых сигналов. – М.: Радио и связь,1981. – 496 с. 2. Атал Б.С. Автоматическое опознавание дикторов по голосам // ТИИЭР. – 1976. – Т. 64, № 4. – С. 48–66. 3. Чуи К. Введение в вэйвлеты. – М.: Мир, 2001. – 412 с. 4. Малла С. Вэйвлеты в обработке сигналов. – М.: Мир, 2005. – 671 с. 5. Добеши И. Десять лекций по вейвлетам. – М.: РХД, 2004. – 464 с. 6. Rabiner L.R., Jang B.H. Fundamentals of speech recognition // New Jersey: Prentice Hall PTR, Englewood Cliffs, 1993. – P. 507. 7. Маркел Д.Д., Грэй А.Х. Линейное предсказание речи. – М.: Связь,1980. – 308 с. 8. Молдокулова Н.В., Трунин-Донской В.Н. Лингво-акустические проблемы создания системы распознавания слитной речи на ЭВМ / Под ред. Ю.Й. Журавлева; АН КиргССР, Вычислительный центр АН СССР. – Фрунзе: Илим, 1989. – 136 с. 9. Davis S.B., Mermelstein P. Comparison of Parametric Representations for Monosyllabic Word Recognition in Continuously Spoken Sentences // IEEE Trans. on Acoustic, Speech and Signal Proc. – 1980. – Vol. 28, N 4. – Р. 357–366. 10. Винцюк Т.К. Анализ, распознавание и интерпретация речевых сигналов. – К.: Наукова думка, 1987. – 261 с. Стаття надійшла до редакції 30.08.2007
id nasplib_isofts_kiev_ua-123456789-755
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
issn 1028-9763
language Russian
last_indexed 2025-12-07T18:09:35Z
publishDate 2008
publisher Інститут проблем математичних машин і систем НАН України
record_format dspace
spelling Криводубский, О.А.
Федоров, Е.Е.
2008-06-25T10:55:24Z
2008-06-25T10:55:24Z
2008
Моделирование особенностей речи диктатора / О.А. Криводубский, Е.Е. Федоров // Математические машины и системы. — 2008. — № 1. — С. 163-170.
1028-9763
https://nasplib.isofts.kiev.ua/handle/123456789/755
004.934.1’1
Для создания системы верификации диктора в статье была предложена методика верификации на основе шипящих звуков, которая использует авторские методы, осуществляющие обобщенную классификацию звуков посредством дискретного и непрерывного вейвлет-преобразования. Для шипящих звуков, выделенных этими методами, проведен количественный анализ систем признаков, основанных на линейном предсказании, нормированном количестве импульсов равной длины и мел-частотных кепстральных коэффициентах (MFCC). Полученные признаки используются в методе верификации, основанном на алгоритме DTW. Табл.: 1. Библиогр.: 10 назв.
Для створення системи верифікації диктора у статті була запропонована методика верифікації на підставі шиплячих звуків, що використовує авторські методи, які здійснюють узагальнену класифікацію звуків за допомогою дискретного та безперервного вейвлет-перетворення. Для шиплячих звуків, виділених цими методами, проведений кількісний аналіз систем ознак, заснованих на лінійному пророкуванні, нормованої кількості імпульсів рівної довжини та мел-частотних кепстральних коефіцієнтах (MFCC). Отримані ознаки використовуються в методі верифікації, заснованому на алгоритмі DTW. Табл.: 1. Бібліогр.: 10 назв.
For creation of a system of verification of the speaker in the article the technique of verification on the basis of unvoiced fricative sounds was offered which uses author's methods executing a generalized classification of sounds by means of discrete and continuous wavelet-transformation. For unvoiced fricative sounds chosen by these methods, the quantitative analysis of systems of indications based on a linear prediction, normalized quantity of pulses of equal length and mel-frequency cepstral coefficients (MFCC) is conducted. The obtained indications are used in a method of verification based on the algorithm DTW. Tabl.: 1. Refs.: 10 titles.
ru
Інститут проблем математичних машин і систем НАН України
Моделювання і управління великими системами
Моделирование особенностей речи диктатора
Моделювання особливостей мови диктатора
The method of verification of the speaker on the basis of unvoiced fricative sounds
Article
published earlier
spellingShingle Моделирование особенностей речи диктатора
Криводубский, О.А.
Федоров, Е.Е.
Моделювання і управління великими системами
title Моделирование особенностей речи диктатора
title_alt Моделювання особливостей мови диктатора
The method of verification of the speaker on the basis of unvoiced fricative sounds
title_full Моделирование особенностей речи диктатора
title_fullStr Моделирование особенностей речи диктатора
title_full_unstemmed Моделирование особенностей речи диктатора
title_short Моделирование особенностей речи диктатора
title_sort моделирование особенностей речи диктатора
topic Моделювання і управління великими системами
topic_facet Моделювання і управління великими системами
url https://nasplib.isofts.kiev.ua/handle/123456789/755
work_keys_str_mv AT krivodubskiioa modelirovanieosobennosteirečidiktatora
AT fedorovee modelirovanieosobennosteirečidiktatora
AT krivodubskiioa modelûvannâosoblivosteimovidiktatora
AT fedorovee modelûvannâosoblivosteimovidiktatora
AT krivodubskiioa themethodofverificationofthespeakeronthebasisofunvoicedfricativesounds
AT fedorovee themethodofverificationofthespeakeronthebasisofunvoicedfricativesounds