output.m

对学好matlab很有用

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,100);%随机序列
for i=1:100 
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:220);
output2=output1/max(output1);
for i=1:100;
    if output2(i)<0.5
        output2(i)=0;
    else 
        output2(i)=1;
    end;
end;
figure(6)
stem(output2);