// описываем sinc
function sinc(t)
{
 if(t==0.0) return 1.0;
 return Math.sin(t)/t;
}

// задаём параметры кодированного сигнала
var t0 = 1.0;       // длительность такта
var T0 = t0/2.0;   // середина первого такта
var N = 7;          // количество кодовых элементов

// задаём кодовую последовательность
var a = new Array(N+10);
a[0]=1.0; a[1]=1.0; a[2]=1.0; 
a[3]=-1.0; a[4]=-1.0; a[5]=1.0; a[6]=-1.0;

// рассчитываем действительную и мнимую часть спектра
var Re = 0.0; 
var Im = 0.0;
for(var n=0; n<N; n++)
{
 Re = Re + a[n]*Math.cos(n*x*t0);
 Im = Im - a[n]*Math.sin(n*x*t0);
}
Re = t0*Re*sinc(x*t0/2.0);
Im = t0*Im*sinc(x*t0/2.0);

// рассчитываем амплитудный спектр (при расчёте фазового спектра это убрать)
return Math.sqrt(Re*Re+Im*Im);

// рассчитываем фазовый спектр
return Math.atan2(Im,Re);