dpsk_pe1.m

一个基于matlab的dpsk仿真

snr=[0:10 20:10:50];
for(q=1:15)
m=0;
for(l=1:1000)
    a=rand(1,11);
    a=2*a-1;
    a=ceil(a);
    da(1)=1;
    for(i=1:11)
        da(i+1)=xor(a(i),da(i));
    end;
    da=2*da-1;
    n=[0:0.001:12];
    for(i=1:12)
        d(1000*(i-1)+1:1000*i+1)=da(i)*cos(4*pi*n(1000*(i-1)+1:1000*i+1));
    end;
    y=awgn(d,snr(q)-27,'measured');
    z1=y.*cos(4*pi*n);
    for(i=1:12)
        for(p=1:1001)
            z11(1000*(i-1)+p)=sum(z1(1000*(i-1)+1:1000*(i-1)+p))*0.001;
        end;
    end;
    z11=hilbert(z11);
    for(i=2:11)
        z(i-1)=z11(1000*i+1)*z11(1000*(i-1)+1);
    end;
    for(i=1:10)
        if(z(i)<0)
            o(i)=0;
        else
            o(i)=1;
        end;
        if(o(i)~=a(i))
            m=m+1;
        end;
    end;
end;
pe(q)=m/10000;
end;
snr=[0:10 20:10:50];
plot(snr,pe,'V-')
hold on;
plot(snr,0.5*exp(-10.^(snr/10)),'r.-')