眼图.m

对学好matlab很有用

R=0.5;
N_T=15;
RATE=4;
T=5;
B = RCOSFIR(R, N_T, RATE, T);
stem(B);
figure(1)
stem(B);%发送滤波器时域波形
i=1;
while(i<=249)
    x(i)=1;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    i=i+1;
    x(i)=0;
    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:256 
if x(i)<0.5 
x(i)=0;
else 
x(i)=1;
end; 
end;
x=x(1:256);
figure(2)
stem(x);%0、1序列
y=conv(x,B);%发送滤波器与0、1序列卷积
figure(3)
stem(y);
for i=1:24
    xf(:,i)=5;
    for i=25:121
        xf(:,i)=0;
    end
end
xt2=ifft(xf);%接收滤波器
output=conv(xt2,y);
figure(4)
stem(output);%接收滤波器卷积
output1=output(121:376);
output2=output1/max(output1);
i=1;
while(i<=256)
    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:256;
    if output2(i)<0.5
        output2(i)=0;
    else 
        output2(i)=1;
    end;
end;
figure(5)
stem(output2);
eyediagram(output2,256);