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)