%%%產生隨機信源信號
m=randint（110000）;
x=reshape（mlength（m）/44）;
%%%分組編碼
G=[1?0?1?1?0?0?0;1?1?1?0?1?0?0;1?1?0?0?0?1?0;0?1?1?0?0?0?1];%設置編碼矩陣
X=mod（x*G2）;%進行編碼，生成編碼后的碼組
X=X‘;
n=length（m）/4*7;
X=reshape（X1n）;
X=bi2de（（reshape（X2n/2））‘‘left-msb‘）;%二進制數組轉十進制
%%%產生QPSK信號
h=modem.pskmod（4）;%?產生調制句柄?相位偏移默認為0
y=modulate（hX）;%對信號進行調制
P=[];
%%%%加入高斯白噪聲
for?SNR=0:10
z=aw