clear?all;
close?all;
clc;
num=100000;
Rb=1000;
tnum=8;?%每個碼元用8個點來表示
Fs=tnum*Rb;
r=0.5;%滾降系數
a=randint（1num）*2-1;
aa=upsample（atnum）;
[NUMDEN]=rcosine（RbFs‘sqrt‘r）;%NUM對應分子系數，DEN是a（1）=1？升余弦濾波器，8倍內插，滾降系數=1；
sender=conv（NUMaa）;%？NUM=delay*（Fs/Rb）*2+1=3*8*2+1=49;卷積長度8e5+49-1
%sender=filter（NUM1aa）;
Eb=1;
s=1;
pe_syx=zeros（113）;%基帶誤碼率矩陣
for?snr_db=-5:10
????snr=10^（snr_db/10）;
????N0=Eb/snr;
????sigma=sqrt（N0/2）;
????noise=randn（1（num+6）*tnum）;%標準正態分布，與sender長度相同
???