bits=70000;
msg=randint（bits12）;
SNR=0:1:12;
L=length（SNR）;
BER1=zeros（1L）;
BER2=zeros（1L）;
modbit1=pskmod（msg2）;
for?k=1:L
????y1=awgn（modbit1SNR（k）‘measured‘）;
????demmsg1=pskdemod（y12）;
????recode=reshape（demmsg1‘1[]）;
????error1=（recode~=msg‘）;
????errorbits=sum（error1）;
????BER1（k）=errorbits/length（msg）;
end
code=encode（msg157‘cyclic‘[1?0?0?0?1?0?1?1?1]）;
modbit2=pskmod（code2）;
f