R=0.5;
N_T=15;
RATE=4;
T=5;
B = RCOSFIR(R, N_T, RATE, T,'sqrt');
stem(B);
figure(1)
stem(B);%發送濾波器時域波形
x=rand(1,400);%隨機序列
i=1;
while(i<=400)
    x(i)=x(i)
    i=i+1
    x(i)=0
    i=i+1
    x(i)=0
    i=i+1
    x(i)=0
    i=i+1
end
for i=1:400 
if x(i)<0.5 
x(i)=0;
else 
x(i)=1;
end; 
end;
figure(2)
stem(x);%0、1序列
y=conv(x,B);%發送濾波器與0、1序列卷積
figure(3)
stem(y);
z=awgn(y,10);%加SNR=10的高斯噪聲
figure(4)
stem(z);
xt2=rcosfir(R,N_T,RATE,T,'sqrt');%接收濾波器
output=conv(xt2,z);
figure(5)
stem(output);%接收濾波器卷積
output1=output(121:520);
output2=output1/max(output1);
i=1;
while(i<=400)
    output2(i)=output2(i)
    i=i+1
    output2(i)=0
    i=i+1
    output2(i)=0
    i=i+1
    output2(i)=0
    i=i+1
end
for i=1:400;
    if output2(i)<0.5
        output2(i)=0;
    else 
        output2(i)=1;
    end;
end;
figure(6)
stem(output2);