MATLAB实现BPSK误码率曲线

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长度相同