matlab仿真高斯信道下的QPSK误码率，有详细的解释，很不错奥

T=100;        %原始信息序列长度
samp=100;     %采样点数
N=samp*T;
m=randint（1T）;     %生成序列
data1=2*m-1;                  %单/双极性变换

for i=1:T
    data1_orig（samp*（i-1）+1:samp*i）=data1（i）;  %信息序列扩展
end

% 信源序列图及其频谱图
figure（1）;
subplot（212）;
pwelch（data1_orig[][][]10）;%data1_orig功率谱密度
title（‘信源序列功率谱密度‘）;
subplot（211）;
x1=0+0.01:0.01:100;
plot（x1data1_orig）;axis（[0 100 -1.2 1.2]）;
title（‘信源序列‘）;

% 将信源分成两路
data11=zeros（1N/2）;            
for i1=1:T/2  
    data11（samp*（i1-1）+1:samp*i1）=data1（2*i1-1）;
end
data12=zeros（1N/2）;              
for i2=1:T/2  
    data12（samp*（i2-1）+1:samp*i2）=data1（2*i2）;
end

%  对载波抽样
for j1=1:N/2                      
    a1（j1）=cos（（j1-1）*pi/20）;
end
for j2=1:N/2                      
    a2（j2）=sin（（j2-1）*pi/20）;
end

% 调制
data21=data11.*a1;  %用余弦调制                         
data22=data12.*a2;  %用正弦调制

data3_modul=data21+data22;                 %经过理想信道的QPSK信号
data3_gauss=awgn（data3_modul10‘measured‘）;  %经过AWGN信道后的QPSK信号


%调制后的信号波形和频谱图
figure（2）;      
subplot（211）;
x2=0+0.01:0.01:50;

plot（x2data3_modul）;   title（‘调制后无衰落QPSK信号波形图‘）;

subplot（212）;
pwelch（data3_modul[][][]10）;  title（‘调制后无衰落QPSK信号频谱图‘）;                                    


%只经过高斯信道后的信号波形和频谱图

figure（3）; 

subplot（211）;
x2=0+0.01:0.01:50;

plot（x2data3_gauss）;title（‘经AWGN信道后的信号波形图（信噪比10dB）‘）; 
                               
subplot（212）;
pwelch（data3_gauss[]