ГОСТ Р ИСО/МЭК 10373-6—2015
/»♦♦«*♦♦«*♦♦♦ Calculate dpOdi
2
ation window and phase factors *»»*/
tor( i-0;i<N_data;i
m
>
i
/* Bartlett winaow */
if ((li_data 4 1) — 0)
i
/* N_data is even */
if <1 < (int) <N_data /2))
I
Wb-2.0F*i/(double)(N_data - l>;
i
else
F
Wb-2.0F*(1C data-i-1)/(double)(N data - 1);
i
1
else
■
/*M_deta is odd */
1Г (1 <- Unt) (O.OOlf(N_data-1)
/2))
Г
i
l
Wb-2.OFM/(double) (K data - 1);
)
else
I
Wb-2.OF-2.OFM/(double)(N dota - 1);
1
i
Yuc|il-(double)cos (wuM*dt)*Wb;
Yus|i]-(double)sin(wuM’dt)’Wb;
Ylc1
1
i-(aouble)cos(wl’i’at)*Wb;
Yls[i]-(double)sin(wlM*dt)’Wb;
sum_Wb i- Wb;
I
c1-N_dalbfsurn_Wb;
у******************* DFT ♦♦***♦**♦♦♦**♦♦**♦♦**♦***♦♦**♦♦*♦♦**♦♦♦**♦♦**♦♦*♦♦
for <J*0;’<count-N_data;$-’*N over)
F
Ol_real-0;
/*
real part O! the up. sideband fourier coefficient
*/
cl_imag-0;
/•
imag part. ol the up. sideband fourier coefficient
*/
c2_real-0; /* real part of the lo. sideband fourier coefficient */
c2_imag-0; /* lir.agpart ol the lo. sideband fourier coefficient *7
i
for( 1-0;KNdata;i♦♦)
I
cl_real-cl_real-*vd[i♦J
1*Yue|iJ
;
cl_imaq-cl_imag*va[i-‘Ji
]*Yus |i)
;
c2_real-c2_real-*vd[i♦J)’Ylc[ i;
c2_imaq-c2_imaq-*vdli-*j|’Yls(ij;
cl_real-2.0F*cf*cl_real/ «double) N_dat.s;
cl_l(rag-2.0F*cr*cl_lmag/(double) N_data;
c2_real-2.0F*cl»c2_real/(double) N_data;
c2_lir>ag-2.0F*ef*c2_irtag/(double) M_daU;
/•••*»*>•*••««»•»*•• j»px scale •*»*•*»*•
/
198