bpsk_qpsk.m

clear
clc
t1=cputime;

Es=2;%功率约束值
N=1000000;%点数

src=zeros（1N）;
src1=zeros（1N/2）;
%sink1=zeros（1N）;
%sink2=zeros（1N）;
%sink22=zeros（1N）;%初始化

clear ij;
SNR0=0;
SNR1=14;
for j=SNR0:SNR1
    %error1=0;
    %error2=0;
    %error3=0;
    snr=j;
    snr1=10^（snr/10）;%将信噪比的值由dB转化为数值
    src=round（rand（1N））;%产生原始发送0、1比特
    %--------将原始0、1比特映向射为0、1、2、3----
    k=1;
    for i=1:2:N-1
        src1（k）=2*src（i）+src（i+1）;
        k=k+1;
    end 
    %-------产生BPSK调制信号---------------------
    M1=modem.pskmod（2）;
    s1=modulate（M1src）;
    pow1=norm（s1‘fro‘）^2/N;
    s1=s1*sqrt（Es/pow1）;
    %-------产生QPSK调制信号---------------------
    M2=modem.pskmod（4）; 
    s2=modulate（M2src1）;
    pow2=norm（s2‘fro‘）^2/N*2;
    s2=s2*sqrt（Es/pow2）;
    %-------产生噪声----------------------------
    noise1=randn（1N）*（sqrt（Es/snr1））;
    noise2=（（randn（1N/2）+sqrt（-1）*randn（1N/2）））*（sqrt（（Es/snr1/2）））;
    %-------产生接收信号-------------------------
    y1=s1+noise1;
    y2=s2+noise2;
    %-------BPSK、QPSK解调----------------------
    M11=modem.pskdemod（2）;
    M22=modem.pskdemod（4）;
    sink