t=linspace(0,5,j); %产生基带信号

　　a=round(rand(1,i));

　　st1=t;

　　for n=1:10

　　if a(n)<1;

　　for m=j/i*(n-1)+1:j/i*n

　　st1(m)=0;

　　end

　　else

　　for m=j/i*(n-1)+1:j/i*n

　　st1(m)=1;

　　end

　　end

　　end

　　subplot(321);

　　plot(t,st1);

　　title('基带信号');

　　axis([0,5,-1,2]);

　　st2=t; %基带信号求反

　　for k=1:j;

　　if st1(k)>=1;

　　st2(k)=0;

　　else

　　st2(k)=1;

　　end

　　end;

　　subplot(322);

　　plot(t,st2);

　　title('基带信号反码');

　　axis([0,5,-1,2]); %载波信号

　　s1=sin(2*pi*fc*t);

　　s2=sin(2*pi*fc*t+pi);

　　subplot(323);

　　plot(s1);

　　title('载波信号'); %调制

　　p1=st1.*s1;

　　p2=st2.*s2;

　　subplot(324);

　　plot(t,p1);

　　title('st1*s1');

　　subplot(325);

　　plot(t,p2);

　　title('st2*s2');

　　e_psk=p1+p2;

　　subplot(326);

　　plot(t,e_psk);

　　title('e_2psk');

　　noise=rand(1,j);

　　psk=e_psk+noise;%加入噪声

　　figure(2);

　　subplot(411);

　　plot(t,psk);

　　title('加噪后波形'); %相干解调

　　psk=psk.*s1;%与载波相乘

　　subplot(412);

　　plot(t,psk);

　　title('与载波相乘后波');

　　[f,af] = T2F(t,psk);%通过低通滤波器

　　[t,psk] = lpf(f,af,B);

　　subplot(413);

　　plot(t,psk);

　　title('通过低通滤波器波形'); %抽样判决

　　for m=0:i-1;

　　if psk(1,m*500+250)<0;

　　for j=m*500+1:(m+1)*500;

　　psk(1,j)=0;

　　end

　　else

　　for j=m*500+1:(m+1)*500;

　　psk(1,j)=1;

　　end

　　end

　　end

　　subplot(414);

　　plot(t,psk);

　　axis([0,5,-1,2]);

　　title('抽样判决后波形')

　　通过以上的仿真分析，我们可以知道，在误码率相同的情况下，所需要的信噪比2ASK比2FSK高3dB，2fsk比2psk高3db，2ask比2psk高6db。也就是说，若信噪比r一定，2psk系统的误码率比2fsk的小，2fsk系统的误码率比2ask的小，由此看来，在抗加性高斯白噪声方面，想干2psk性能最好，2fsk次之，2ask最差。当然，后面还有DPSK,等等改进型的。(如果你对本文的代码和仿真图像有任何疑问可联系本文作者QQ:1529620855)