ГОСТ Р ИСО/МЭК10373-7— 2011
{
cO
_
real=cO
_
real+vd[i]‘(ftoat)cos(wO’vtime[i]):
c0Jmag=cQJmag+vd[i]*(fk>at)sin(v/O*vtime[iJ);
c1
_
real=c1
_
real+vd[i]‘(fk>at)cos(wu,vlime[i]):
с1Jmag=c1Jmag+vd[i]*(fk>at}sin(v/u*vtime[i]);
c2
_
real=c2
_
real+vd[i]‘(fk>at)cos(wrvtime[i]):
c2Jmag=c2Jmag+vd[i]*(fk>at)sin(wl’vtimeli}):
}
/**‘ *
.......
•*••••«••• OFT scale......................•.«»•/
c0
_
real=2.0F’c0
_
real/(floatXN
_
data);
c0Jmag=2.0F’c0Jmag/(float)(NjJata):
c1
_
real=2.0F’c1
_
real/(float){N
_
data);
c1Jmag=2.0F*c1
_
imag/(lloat)(N
_
data):
c2
_
real=2.0F’c2
_
real/(floatXN
_
data);
c2Jmag=2.0F*c2Jmag/(fIoat)(N
_
data):
/****’ ........... abSO|uie fourier coefficient ..........*’V
cO
_
abs=(float)sqrt(cO
_
real*cO
_
real + cOJmag’cOjmag);
c1
_
abs=(float)sqrt(c1
_
rearc1
_
real + c1
_
imag*c1Jmag);
c2
_
abs=(float)sqrt(c2
_
real*c2
_
real+c2
_
imag*c2Jm ag);
/•*•*” ****“ ** Phase of fourier coefficient
............
**/
cO
_
phase=(fIoat)atan2(c0Jmag ,c0
_
real):
c1
_
phase=(float)atan2(c1 Jmag.c1
_
real):
c2
_
phase=(float)atan2(c2Jmag,c2
_
real):
/* "•* ’
.......
.Result Display
.
........................................*/
printfC4n\nResults: Vi*):
pmtffCamer’);
printffAbs: %7.3fmV ”.1000’c0
_
abs);
printffPhase: %3.0fdeg\n".c0
_
phase/pi‘ 18Q):
printffUpper sideband "):
printf("Abs: %7.3fmV ".1000*c1
_
abs);
printf(’Phase: %3.0fdegW.c1
_
phase/pi*18Q);
pnntffLower sideband ");
printf(’Abs: %7.3fmV ".1000*c2
_
abs):
printffPhase: %3.0fdeg\n\n".c2
_
phase/p)*180):
printf(*\n*****•******’■,••****,*•
...........
...................**W):
return;
V......................... End DFT ......................*/
/
.....
* ..............................................................................................*
...................
/***MAIN LOOP" 7
/**..................................
.......
.................................
int main(unsigned short paramCount.char ’paramListQ)
{
char fname[256J;
unsigned int i.sample
_
count;
pi = (float)atan(l )*4; l‘ calculate pi V
printf(‘\n
.......
........................................................................\n")
printfpn"” ISO/IEC 10373-7 VICC Test-Program •***\n‘)
printf(
Л
п**** Version: 1.1 JUL 2000" “ Vi";
printfCVi"**
.......
.............................................
23