clc
close all
clear
SNR=0:1:20;                 %信噪比变化范围
SNR1=0.5*(10.^(SNR/10));    %将信噪比转化成直角坐标
N=1000000;                  %仿真点数
M1=2;                       %BPSK
M2=4;                       %QPSK
x1=randi([0,1],1,N);         %产生随机信号
x2=randi([0,3],1,N);
R=raylrnd(0.5,1,N);         %产生瑞利信号
h1=pskmod(x1,M1);            %BPSK调制
h2=pskmod(x2,M2);            %QPSK调制
H1=h1.*R;                   %BPSK经瑞利信道
H2=h2.*R;                   %QPSK经瑞利信道