2dpsk非相干解调

function y=fxg_2DPSK
data=rand_data（20000）;
d（1:）=0.14:0.01:0.45;
n=length（d）;
pe=zeros（1n）; 
for j=1:n
for i=1:20000
     if（data（i）==1）
         r=0.5*（（d（1j）+wgn（11-20））*（-d（1j）+wgn（11-20））+wgn（11-20）*wgn（11-20））;
     else
          r=0.5*（（d（1j）+wgn（11-20））*（d（1j）+wgn（11-20））+wgn（11-20）*wgn（11-20））;
     end
     if（r>0）
         data2（i）=0;
     e