ГОСТ Р 53556.4—2013
pit [ismp] = (pot
0
[idiv] *pit_cvO [index
0
[idiv][icv]] *
po
/1
[idiv] *prt_cv
1
[index
1
[idiv]][icvj) /
2
;
;
}
else {
for (icv =
0
: /cv <lengthp pdiv]-
1
; lev ++J {
ismp = ((icv*idiv) %N_DIV_P) * icv * N_DIV_P:
ismp = ismp/2 * (ismp%
2
) *
20
;
pit [ismp] = (polO [idiv] *pit_cvO [indoxO [idiv] [icv]] *
po
/1
[idiv] ’pit_cv
1
[index
1
[idiv]][icv]) / 2;
}
icv = lengthp [idiv]A;
ismp = idiv * icv* N_DIV_P;
ismp = ismp/2 * (ismp%2) *
20
;
pit [ismp] = (polO [idiv] *pit_cvO [indexO [idiv] [tcvj] *
po
/1
[idiv] *prt_cv
1
[indexl fidivj] [icv]) / 2;
}
}■
6.10.3.3.4 Проектирование периодических пиковых компонентов в линейной шкале
Сначала параметры вычисляются следующим образом:
femin = к>д2 ((N_FR / (double) ISAMPF) -0.2);
femax = Юд2 ((N_FR / (double) ISAMPF) *2.4);
if (ISAMPF ==
8
) bandwidth = 1.5;
в/se. if (ISAMPF> = 11J bandwidth = 2.0;
e/se. if (ISAMPF> = 22) bandwidth = 4.0;
if(bandwidth <1 ./UPPER_BOUNDARY [0] [0]) bandwidth = 1JUPPER_BOUNDARY[0][0].
где femin является минимальной частотой при квантовании, /стах является максимальной частотой
и ISAMPFявляется целочисленной частотой дискретизации, усеченной из значений стандартной частоты,
перечисленных в правой графе таблицы 82.
Затем базовая частота периодических пиковых компонентовдекодируется согласно следующей про
цедуре;
for (i_ch = 0; i_ch <N_CH; i_ch **) {
p o w j = (int) (pow (1.009792f. (float) index_pit[i_sup]) *4096. * 0.5):
bi_i = (int) ((float) bk>ck_size_sampies / (float) isampf *0,2*1024+0.5);
p itch j =p o w j *b/_i-256.;
bfreq [l_ch] = p/fc/jj/16384.;
}■
Прежде, чем спроектировать периодические пиковые компоненты в линейную шкалу, все элементы
целевого массива pit_seq [][] обнуляются;
for (i_ch = 0; i_ch <N_CH; i_ch **) {
for (ismp = 0; ismp <N_FR; ismp **) {
pit_seq p_ch] [ismp] =
0
.;
}
}■
Затем восстановленные периодические пиковые компоненты проектируются в линейную шкалу сле
дующим образом:
for (i_ch = 0; i_ch <N_CH; i_ch **) {
if(bandwidth * uppedimitj <16384) {
tmpnpOJ =p itc h j’16384. / (upperlimitj):
tm p n p ij =tmpnpOJ *N_FR_P:
tmpnpOJ = tmpnp1J / N_FR:
npcount =tmpnp0_i/16384.;
} else {
tmpnpOJ =p itch j *bandwidth*2:
tmpnp1J =tmpnpOJ *N_FR_P;
158