wfcreate u 1 10000
genr x=nrnd
smpl 2 10000
x=0.5*x(-1)+nrnd
smpl @all
genr y=0.5*x(-1)+nrnd
group group01
group01.add x y
group01.cspec(all)
smpl @all
genr w=@acos(-1)*@trend/10000
genr fxw=(1/(1+0.25-@cos(w)))*(1/(2*@acos(-1)))
genr co=0.5*@cos(w)*fxw
genr qu=0.5*@sin(w)*fxw
genr am=0.5*fxw
genr phase=@atan(-qu/co)
genr phase_shift=phase/w
genr gain=am/fxw
genr fe=1/(2*@acos(-1))
genr coherence=((1+fe/(0.25*fxw)))^(-1)