function Md_m=Modulate(Md_f,phcode,fd,i)
phcode = phcode;
fd = fd;
Loop_time = i;
BOC_Code = BOCcode(phcode,fd,Loop_time);

 BOC_Mapping=BOC_Code;
if Md_f == 1
    Md_m = BOC_Mapping;
else
    Md_m = -BOC_Mapping;
end