MATLAB环境下QPSK信号AWGN信道误码率仿真

信号源：数字基带信号、根升余弦脉冲成型（上采样8倍，即每个符号8个采样点，滚降：0.2）； 传输：AWGN信道（信噪比范围可调） 接收器：匹配滤波，相关解调，判决 画出接收信号的眼图、星座图；计算误码率并与理论值比较。

clear all
nsymbol=1000;         %每种信噪比下的发送符号数

T=1;                        %符号周期
fs=100;                     %每个符号的采样点数
ts=1/fs;                    %采样时间间隔
t=0:ts:T-ts;                %时间向量
fc=10;  
Fd=1;
Fc=0;
Fs=Fd*8;%载波频率



c=sqrt（2/T）*exp（j*2*pi*fc*t）;   %载波信号
c1=sqrt（2/T）*cos（2*pi*fc*t）;    %同相载波
c2=-sqrt（2/T）*sin（2*pi*fc*t）;   %正交载波     

M=4;  %8-PSK
%graycode=[0 1 2 3 6 7 4 5];             %Gray编码规则                      
EsN0=0:10;                              %信噪比，Es/N0
snr1=10.^（EsN0/10）;                     %信噪比转换为线性值


msg=randint（1nsymbolM）;               %消息数据
%msg1=graycode（msg+1）;                   %Gray映射
msgmod1=pskmod（msgM）.‘;                %基带Q-PSK调制


% 经过滤波器   成型滤波器
rolloff=0.2;
delay=3;
msgmod=rcosflt（msgmod1FdFs‘fir/sqrt‘rolloffdelay）;

msgmod=msgmod（（delay*Fs/Fd+1:（end - delay*Fs/Fd））1）;

scatterplot（msgmod）

 %载波调制
tx=real（msgmod*c）;                     
tx1=reshape（tx.‘1length（msgmod）*length（c））;   

spow=norm（tx1）.^2/nsymbol;              %求每个符号的平均功率

for indx=1:length（EsN0）
    sigma=sqrt（spow/（2*snr1（indx）））;                %根据符号功率求噪声功率
    rx=tx1+sigma*randn（1length（tx1））;              %加入高斯白噪声
    rx1=reshape（rxlength（c）length（msgmod））;     
    
    r1=（c1*rx1）/length（c1）;                         %相关运算
    r2=（c2*rx1）/length（c2）;
    r=r1+j*r2;
    eyediagram（r2）;
    scatterplot（r）;
    
    
     %接受端  信号经过低通滤波器  匹配滤波器
     rolloff=0.2;     
     %delay=3;
     k=rcosflt（rFdFs‘fir/sqrt/Fs‘rolloff）;
     % g= downsample（k8）;
     k=k（（delay*Fs/Fd+1:（end - delay*Fs/Fd））1）;
%     h = resample（k 1000 8048）;

    y=pskdemod（kM）;        %PSK解调;
    %h = resample（y 1000 8000）;
     g= downsample（y8）;
    
    u=g‘;
 
    %[errber（indx）]=biterr（msgylog2（M））;     %误比特率
    [errser（indx）]=symerr（msgu）;             %误符号率
end
figure（）;
ser1=2*qfunc（sqrt（2*snr1）*sin（pi/M））;               %理论误符号率
%ber1=1/log2（M）*ser1;                                %理论误比特率
semilogy（EsN0ser‘-k*‘EsN0ser1‘-k.‘）;
title（‘Q-PSK载波调制信号在AWGN信道下的性能‘）
xlabel（‘Es/N0‘）;ylabel（‘误比特率‘）
legend（‘误比特率‘‘理论误比特率‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       43475  2014-04-11 00:09  QQ截图20140411000906.png
     文件        2341  2014-04-10 23:45  qpsk过高斯信道.m
     文件       14579  2014-04-11 00:07  QQ截图20140411000753.png