%用于生成QPSK序列
N=32; %信息序列個數
Fs=32; %采樣頻率
Fc=2; %載波頻率
SNR=-2; %信噪比

%生成隨機信息序列
source=randsrc(1,N,[0,1;.5,.5]);

%映射矩陣
A=[0 1;3,2];


%生成Qpsk序列
 for i=1:N,
    for j=1:Fs,
        modsrc((i-1)*Fs+j)=cos(2*pi*Fc/Fs*j+source(i)*pi);
    end
end

%加入高斯白噪聲
modsrcnoise=awgn(modsrc,SNR,'measured','db');