Teza doctorat (pdf) - Universitatea Tehnică

Estimarea perioadei actuale se face conform relaţiilor:
113
Cap. 4. Analiza semnalului vocal
D
1 d
( PIV , M S ( i))
D
2 d
( PIV , M D ( j))
(4.67)
PER
(
D1
D
2 ) / 2
unde:
- MS(i) este punctul de maxim local din stânga pivotului de amplitudine comparabilă cu
pivotul: ( ( i),
PIV ) S
M S
- MD(j) este punctul de maxim local din dreapta pivotului de amplitudine comparabilă cu
acesta: ( ( j),
PIV ) S
M D
Algoritmul de estimare a perioadei este explicat în cele ce urmează :
Algoritmul Estimare_Perioadă :
(1) Se iniţializează valorile de perioadă din stânga şi dreapta pivotului: Per[Stânga] =
Per[Dreapta] = 0 , Stânga =0 , Dreapta = 1 , ind_per = 0;
(2) Pentru direcţia Stânga, apoi Dreapta se execută paşii (3) – (9) ;
(3) Pentru raportul de amplitudine R_AMP = 70%, scăzând cu 10% la fiecare iteraţie, se
execută paşii (4) – (9) ;
(4) Pentru fiecare punct P(i) de tip maxim local (MAX) ce urmează după pivotul PIV în
direcţia curentă se execută paşii (5) - (9).
(5) Se calculează distanţa – în număr de eşantioane - dintre pivot şi punctul curent : DIST
= D( P(i) , PIV )
(6) Dacă DIST < D_MIN (se ia D_MIN corespunzător valorii de 2,5 ms), atunci P(i) nu este
un maxim de perioadă şi se reia pasul (4) cu următoarea iteraţie ;
(7) Dacă DIST > D_MAX (se ia D_MAX corespunzător valorii de 12,5 ms), atunci P(i) nu
este un maxim de perioadă şi se reia pasul (4) cu următoarea iteraţie;
(8) Se calculează COND_GASIT = P(i) / PIV > R_AMP;
(9) Dacă COND_GASIT e adevărată, înseamnă că s-a găsit următorul maxim de perioadă:
se reţine valoarea perioadei corespunzătoare: Per[ind_per] = | i – i_pivot | ; se reia pasul (2) cu
următoarea iteraţie, căutându-se următorul maxim de perioadă în cealaltă direcţie;
(10) În finalul iteraţiilor, dacă Per[Stânga] = 0, se ia pentru valoarea perioadei estimate
PER = Per[Dreapta] ; dacă Per[Dreapta] = 0, se ia pentru valoarea perioadei estimate PER =
Per[Stânga] ; se returnează True. Acest caz se întâlneşte când pivotul a fost detectat în
extremitatea stânga sau dreapta a segmentului considerat;
(11) Dacă valorile celor două perioade diferă cu o marjă mică de eroare: |Δ(Per[Stânga] -
Per[Dreapta])| < Δ_ER, atunci se ia pentru valoarea perioadei estimate: PER = (Per[Stânga] +
Per[Dreapta]) / 2; se returnează True;
(12) Dacă |Δ(Per[Stânga] - Per[Dreapta])| > Δ_ER, atunci se ia pentru valoarea perioadei
estimate perioada cea mai apropiată de valoarea medie Per_Med a perioadelor vocalice
calculată până în acel moment: PER = Per[k], unde k minimizează expresia |Δ(Per[k] –
Per_Med)| ; se returnează True.