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