通信系统设计—数字带通传输系统仿真

clear all;
clc;
nsymbol=100;              %每种信噪比下的发送符号数
fs=1;                     %每个符号的采样点数
ts=1/fs;                    %采样时间间隔
M=2;                        %BPSK
graycode=[0 1];             %Gray编码规则   
EbN0=0:15;
msg=randint（1nsymbolM）;               %消息数据
figure（1）;
subplot（121）;stem（msg）;title（‘调制前的消息序列‘）;


msg1=graycode（msg+1）;                   %Gray映射
msgmod=pskmod（msg1M）.‘;                %基带BPSK调制
subplot（122）;stem（msgmod）;title（‘调制后的消息序列‘）;
scatterplot（msgmod）;title（‘星座图‘）


msgmod2=zeros（19*nsymbol）;      %8倍零插值
for i=0:nsymbol-1
    msgmod2（9*i+1）=msgmod（i+1）;  %插值实现
end
figure（3）;
subplot（121）;stem（real（msgmod2））;title（‘8倍零插值后实部图‘）;
subplot（122）;stem（imag（msgmod2））;title（‘8倍零插值后虚部图‘）;


fsl=10;    %滤波器的采样频率
rolloff=0.25;   %滤波器的滚降系数
delay=5;    %滤波器时延
rrcfilter=rcosine（fsfsl‘sqrt‘rolloffdelay）;%设计根升余弦滤波器
s=rcosflt（msgmod2fsfsl‘filter‘rrcfilter）;%通过根升余弦滤波器进行脉冲成形
figure（4）;
subplot（121）;plot（real（s））;title（‘经过发送滤波器后的实部信号‘）;
subplot（122）;plot（imag（s））;title（‘经过发送滤波器后的虚部信号‘）;

jrg=awgn（s1‘measured‘）;
eyediagram（jrg50）;title（‘SNR=1dB时的正交信号眼图‘）;
jrg1=awgn（s15‘measured‘）;
eyediagram（jrg1